Article Text

Lung volume measurement using chest CT in COVID-19 patients: a cohort study in Japan
  1. Shiro Otake1,
  2. Yusuke Shiraishi2,
  3. Shotaro Chubachi1,
  4. Naoya Tanabe2,
  5. Tomoki Maetani2,
  6. Takanori Asakura3,
  7. Ho Namkoong3,
  8. Takashi Shimada1,
  9. Shuhei Azekawa1,
  10. Kensuke Nakagawara1,
  11. Hiromu Tanaka1,
  12. Takahiro Fukushima1,
  13. Mayuko Watase4,
  14. Hideki Terai1,
  15. Mamoru Sasaki5,
  16. Soichiro Ueda5,
  17. Yukari Kato6,
  18. Norihiro Harada6,
  19. Shoji Suzuki7,
  20. Shuichi Yoshida7,
  21. Hiroki Tateno7,
  22. Yoshitake Yamada8,
  23. Masahiro Jinzaki8,
  24. Toyohiro Hirai9,
  25. Yukinori Okada10,11,
  26. Ryuji Koike12,
  27. Makoto Ishii13,
  28. Naoki Hasegawa14,
  29. Akinori Kimura15,
  30. Seiya Imoto16,
  31. Satoru Miyano17,
  32. Seishi Ogawa18,19,
  33. Takanori Kanai20 and
  34. Koichi Fukunaga1
  1. 1ivision of Pulmonary Medicine, Department of Medicine, Keio University School of Medicine, Tokyo, Japan
  2. 2Department of Respiratory Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
  3. 3Division of Pulmonary Medicine, Department of Medicine, Keio University School of Medicine, Tokyo, Japan
  4. 4Department of Respiratory Medicine, National Hospital Organization Tokyo Medical Centre, Tokyo, Japan
  5. 5Department of Internal Medicine, Saitama Medical Center, Tokyo, Japan
  6. 6Division of Respiratory Medicine, Juntendo University School of Medicine Graduate School of Medicine, Bunkyo-ku, Japan
  7. 7Department of Pulmonary Medicine, Saitama City Hospital, Saitama, Japan
  8. 8Keio University Department of Radiology, Shinjuku-ku, Tokyo, Japan
  9. 9Respiratory Medicine, Graduate School of Medicine, Kyoto University, Kyoto, Japan
  10. 10Department of Statistical Genetics, Osaka University Graduate School of Medicine, Suita, Japan
  11. 11Department of Genome Informatics, The University of Tokyo Graduate School of Medicine Faculty of Medicine, Bunkyo-ku, Japan
  12. 12Department of Pharmacovigilance, Tokyo Medical and Dental University, Tokyo, Japan
  13. 13Faculty of Medicine, Nagoya University Graduate School of Medicine, Nagoya, Japan
  14. 14Center for Infectious Diseases and Infection Control, Keio University, School of Medicine, Tokyo, Japan
  15. 15Medical Research Institute, Tokyo Medical and Dental University, Tokyo, Japan
  16. 16The University of Tokyo, Bunkyo-ku, Japan
  17. 17Tokyo Medical and Dental University, Bunkyo-ku, Japan
  18. 18Department of Pathology and Tumor Biology, Kyoto University Graduate School of Medicine Faculty of Medicine, Kyoto, Japan
  19. 19Department of Medicine, Regenerative Medicine Karolinska Institute, Stockholm, Sweden
  20. 20Division of Gastroenterology and Hepatology Department of Internal Medicine, Keio University School of Medicine, Shinjuku-ku, Japan
  1. Correspondence to Dr Shotaro Chubachi; bachibachi472000{at}z6.keio.jp

Abstract

Objective This study aimed to investigate the utility of CT quantification of lung volume for predicting critical outcomes in COVID-19 patients.

Methods This retrospective cohort study included 1200 hospitalised patients with COVID-19 from 4 hospitals. Lung fields were extracted using artificial intelligence-based segmentation, and the percentage of the predicted (%pred) total lung volume (TLC (%pred)) was calculated. The incidence of critical outcomes and posthospitalisation complications was compared between patients with low and high CT lung volumes classified based on the median percentage of predicted TLCct (n=600 for each). Prognostic factors for residual lung volume loss were investigated in 208 patients with COVID-19 via a follow-up CT after 3 months.

Results The incidence of critical outcomes was higher in the low TLCct (%pred) group than in the high TLCct (%pred) group (14.2% vs 3.3%, p<0.0001). Multivariable analysis of previously reported factors (age, sex, body mass index and comorbidities) demonstrated that CT-derived lung volume was significantly associated with critical outcomes. The low TLCct (%pred) group exhibited a higher incidence of bacterial infection, heart failure, thromboembolism, liver dysfunction and renal dysfunction than the high TLCct (%pred) group. TLCct (%pred) at 3 months was similarly divided into two groups at the median (71.8%). Among patients with follow-up CT scans, lung volumes showed a recovery trend from the time of admission to 3 months but remained lower in critical cases at 3 months.

Conclusion Lower CT lung volume was associated with critical outcomes, posthospitalisation complications and slower improvement of clinical conditions in COVID-19 patients.

  • COVID-19
  • imaging/CT MRI etc

Data availability statement

Data are available on reasonable request. The datasets used and/or analysed in the current study are available from the corresponding author on reasonable request.

http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

WHAT IS ALREADY KNOWN ON THIS TOPIC

  • The measurement of lung volume using chest CT has been reported to be a simple and useful indicator of various lung diseases.

  • However, few reports have examined lung volume in COVID-19 patients.

WHAT THIS STUDY ADDS

  • CT-derived lung volume measurement has significant value as an easily assessable indicator for predicting the severity of COVID-19.

  • CT-derived lung volume can identify populations at risk of COVID-19 sequelae by evaluating residual lung volume reduction.

HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY

  • Longitudinal studies should examine the association between CT-derived lung volume and the symptoms of long COVID.

Introduction

The first case of COVID-19 occurred in October 2019; since then, 764 473 623 cases and 6 915 273 deaths have been reported worldwide.1 Although initially feared as an unidentified pathogen, comprehensive research has been conducted on a global scale, resulting in advances in vaccines and therapeutic agents, thereby improving prognosis.2–4 However, the virus persists with mutations, and the epidemic remains unrelenting in various regions of the world.5

Chest CT is extensively used to diagnose COVID-19 and predict its severity.6 7 The measurement of lung volume using chest CT has been reported to be a simple and useful indicator of various lung diseases, including chronic obstructive disease (COPD) and idiopathic pulmonary fibrosis (IPF), where CT-derived lung volume correlates with pulmonary function test results.8–11 Recent studies have indicated its usefulness in assessing disease activity, especially IPF because it enables the evaluation of longitudinal lung volume reduction resulting from lung fibrosis.12 13 COVID-19 manifests as ground glass opacities and consolidation shadows during the acute phase, which may improve or lead to fibrotic lesions over time.14 15 Fibrotic scarring following acute inflammation may also reduce lung volume in patients with COVID-19 pneumonia. However, few reports have examined lung volume in COVID-19 patients,16–19 and only one report has specifically investigated its association with severity.17

COVID-19 causes long-term lung sequelae.20 Chest CT serves as a valuable tool for assessing structural sequelae of the lungs, and many studies using qualitative evaluation by radiologists have confirmed the presence of residual opacities at the 3-month time point following COVID-19 diagnosis, particularly in severe cases.21–24 Therefore, CT-derived lung volume measurements may be advantageous for evaluating long-term lung fibrosis in COVID-19 patients.16 However, the frequency of patients with a long-term reduction in lung volume and clinical characteristics has not been previously reported.

Based on this background, we hypothesised that automated quantification of lung volume using chest CT would be useful for predicting the severity of COVID-19 and the pulmonary sequelae at 3 months. Therefore, using the world’s largest sample size of CT analysis data in a multicentre, long-term patient cohort, we aimed to (1) compare detailed clinical characteristics, such as the relationship between lung volume and critical outcomes and complications during hospitalisation and (2) examine the 3-month course of lung volume, clinical characteristics and predictive factors in patients with residual lung volume loss.

Methods

Study design and settings

Data for all COVID-19 cases used for the present secondary analysis of a multicentre retrospective study were obtained by the Japan COVID-19 Task Force.4 All patients who participated in the present study provided written or oral informed consent. The Japan COVID-19 Task Force collected clinical information on patients with COVID-19 aged >18 years who were diagnosed by PCR or antigen testing from four hospitals nationwide in Japan. Of the 1410 patients identified, 210 were excluded due to a lack of CT imaging or height data; thus, 1200 patients were included in the baseline chest CT analysis. The characteristics of the excluded cases are described in online supplemental table 1. Among all patients, there were 102 critical cases, of these, 20 cases (19.6%) were critical at the time the CT was taken and 82 cases (80.4%) subsequently became critical. Of these, 208 patients for whom CT data were available at 3 months were included in the 3-month CT analysis (online supplemental figure 1).

Patient and public involvement

Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

Data collection

The following information was extracted from the electronic case record form: age, sex, height, weight, clinical signs and symptoms, laboratory findings on admission, comorbidities and treatment details. Critical outcomes were defined as the need for oxygen supplementation via high-flow oxygen therapy, mechanical ventilation, extracorporeal membrane oxygenation or death.2 The follow-up period for the outcomes was throughout the hospitalisation period. All laboratory tests were performed according to the clinical care needs of the patients. Symptoms and signs were not only considered at the time of referral and admission but also throughout the hospitalisation period. Laboratory results were collected within 48 hours of the initial visit or admission. The attending physicians at each facility reviewed and evaluated the chest CT images for qualitative pneumonia.25 A team of respiratory clinicians reviewed the data. If core data were missing, the clinician was contacted for data collection. Missing or absent data on patient background were noted as unknown.

CT acquisition

All CT images were obtained after full inspiration. Images of the entire lung, with a slice thickness of 1–5 mm, were reconstructed using standard kernels. The CT scanners used were the SOMATOM series (Siemens Healthineers), Aquilion series (Canon Medical Systems), Revolution series (GE Healthcare), Discovery series (GE Healthcare) and BrightSpeed (GE Healthcare). All CT images were taken within 48 hours of admission.

AI-based image analysis

Segmentation of pneumonia and the total lung capacity on CT (TLCct) were calculated on full-inspiratory CT using a SYNAPSE VINCENT volume analyzer (Fujifilm Medical, Tokyo, Japan),25 and the predicted TLCct% was calculated based on the predicted values.26 Representative images of the TLCct are shown in online supplemental figure 2.

Statistical analysis

Data are presented as the mean±SE. Data were compared among groups using a t-test and χ2 test. Receiver operating characteristic (ROC) curve analysis was performed to determine the appropriate TLCct% predicted cut-off values for various outcomes: critical, oxygen requirement, invasive mechanical ventilation (IMV) requirement and mortality. To evaluate the relationship between TLCct and each outcome, we performed a logistic regression analysis as a multivariate analysis adjusted for known severity risk factors: age, sex, body mass index (BMI), hypertension, diabetes mellitus, cardiovascular disease, hyperuricaemia, COPD and chronic kidney disease.27–29.30 Statistical significance was set at p<0.05. Data were analysed using JMP V.16 software (SAS Institute).

Results

Comparison of clinical characteristics between high and low TLCct groups

The distribution of the TLCct (%pred) on admission is shown in online supplemental figure 3. Patients were divided into two groups according to the median value, with the top 1/2 defined as high TLCct (%pred) (≧67.1%, n=600) and the bottom 1/2 as low TLCct (%pred) (<67.1%, n=600). A comparison of the clinical characteristics between the two groups is shown in online supplemental table 2. The low TLCct (%pred) group exhibited a higher proportion of men, a higher BMI, a higher incidence of hypertension and chronic kidney disease, and a lower incidence of asthma than the high TLCct (%pred) group. The low TLCct (%pred) group exhibited a higher frequency of lower respiratory symptoms, a lower frequency of upper respiratory symptoms, tachycardia, hyperthermia, tachypnoea and lower SpO2 than the high TLCct (%pred) group (online supplemental table 3). In the low TLCct (%pred) group, leucocyte, aspartate aminotransferase, alanine transaminase, blood urea nitrogen, lactate dehydrogenase (LDH), ferritin, KL-6, HbA1c, D-dimer, procalcitonin and C reactive protein (CRP) levels were higher and lymphocyte and albumin levels were lower than those in the high TLCct (%pred) group (online supplemental table 4).

Critical cases exhibited a lower TLCct(%pred) than non-critical cases. (51.4% vs 67.7%, p<0.0001) (figure 1A). ROC curves showed that the optimal cut-off for TLCct for predicting critical outcomes was 51.1 (area under the curve, AUC 0.77, sensitivity 53.9%, specificity 85.8%) (figure 1B). Multivariable analysis adjusted for known poor prognostic factors showed that TLCctwas independently associated with critical outcomes. Also, unadjusted univariate analysis for each variable is shown in online supplemental figure 4. In addition, TLCct predicted oxygen and IMV requirements (online supplemental figures 5 and 6). Also, lung volume was lower in critical cases, both with and without qualitative pneumonia. Multivariable analysis adjusted for AI-based pneumonia volume or the presence of qualitative pneumonia showed an independent association between TLCct and critical outcomes (online supplemental figure 7). Compared with the CURB-65, an existing severity prediction score, the AUC in the ROC curve was higher in TLCct(%pred) (online supplemental figure 8). Then we performed the following sensitivity analyses. In multivariable analysis including the variables found to differ at baseline characteristics (online supplemental table 2), and it was shown that TLC independently predicted critical outcomes. In addition, we excluded cases that were critical at the time the CT was taken and analysed the association between TLCct (%pred) and critical outcomes. TLCct (%pred) was significantly lower in critical cases, and TLCct (%pred) still independently predicted critical outcomes in multivariable analysis including known severity risk factors (online supplemental figure 9).

Figure 1

Relationship between critical outcomes and TLCct (%pred). (A) Comparison of TLCct (%pred) according to disease severity. (B) ROC curve of TLCct (%pred) in predicting critical outcomes. (C) Forest plot showing multivariate logistic regression analysis to evaluate the relationship between each variable and critical outcome. ****p<0.0001. AUC, area under the curve; BMI, body mass index; COPD, chronic obstructive pulmonary disease; ROC, receiver operating characteristic; TLCct (%pred): total lung capacity on CT (%predicted).

Association between complications during hospitalisation and TLCct(%pred)

Comparisons of the TLCct(%pred) with and without complications during hospitalisation are shown in figure 2. The cases with bacterial infection exhibited lower TLCct (%pred) than cases without bacterial infection (58.81% vs 67.42%, p<0.0001). The same applies to cases with and without heart failure (54.68% vs 66.88%, p=0.0037), thromboembolism (59.62% vs 66.97%, p=0.0174), liver dysfunction (61.25% vs 70.84%, p<0.0001) and renal dysfunction (62.66% vs 67.78%, p<0.0001).

Figure 2

Comparison of the incidence of posthospitalisation complications between high and low TLCct (%pred) groups. Patients were divided into two groups according to the median TLCct value (%pred). *p<0.05, **p<0.01, ****p<0.0001. TLCct (%pred), total lung capacity on CT (%predicted).

Association between the clinical outcomes and TLCct(%pred), stratified epidemic waves

The epidemic waves were categorised around the prevalent virus variants as the 1st–3rd wave (conventional variant), 4th wave (alpha variant), 5th wave (delta variant) and 6th–7th wave (omicron variant).31 The critical cases exhibited a lower TLCct(%pred) than non-critical cases for all epidemic waves (1st–3rd waves, 69.39% vs 54.78%, p<0.0001; 4th wave, 69.43% vs 53.03%, p<0.0001; 5th wave, 60.66% vs 47.49%, p<0.0001; and 6th–7th waves, 58.48% vs 47.92%, p=0.0024) (figure 3).

Figure 3

Relationship between critical outcome and TLCct (%pred) adjusted for epidemic waves in Japan. Patients were divided into two groups according to the median TLCct value (%pred). **p<0.01, ****p<0.0001. TLCct (%pred), total lung capacity on CT (%predicted).

Comparison of clinical characteristics between high and low TLCct groups at 3 months

The distribution of TLCct (%pred) at 3 months is shown in online supplemental figure 10. The patients were divided into two groups at the median (71.8%) and then compared in terms of baseline clinical characteristics (online supplemental table 5). During the 3-month CT analysis, the low TLCct (%pred) group exhibited a higher proportion of men and a higher incidence of hypertension, chronic kidney disease, diabetes mellitus and hyperuricaemia. After admission, the low TLCct (%pred) group in the 3-month CT analysis exhibited a higher incidence of critical outcomes, oxygen requirement, IMV requirement and bacterial infection (table 1). The longitudinal changes in TLCct (%pred) stratified by patients with and without critical outcomes are shown in figure 4. Overall, lung volumes showed a recovery trend from the time of admission to 3 months, but they were still lower in critical cases at 3 months (p<0.0001). When restricting the analysis to cases with a 3-month follow-up, 73.6% of those with low lung volumes on admission maintained this classification after 3 months (online supplemental figure 11). In addition, there was a trend towards more severe cases in the group with a larger absolute value of change than those with a lower value of change (online supplemental table 6).

Figure 4

Comparison of TLCct (%pred) at 3 months. Longitudinal changes in TLCct (%pred) stratified by patients with and without critical outcomes. ***p<0.001, ****p<0.0001. TLCct (%pred), total lung capacity on CT (%predicted).

Table 1

Comparison between clinical outcomes of high and low TLCCT groups, as per 3-month CT analysis

Discussion

To the best of our knowledge, this is the largest study to demonstrate the utility of lung volume measurement using chest CT in COVID-19 patients, and it is also the first to investigate the clinical characteristics of patients with decreased lung volumes after the acute phase. This study demonstrated that CT-derived lung volume can predict critical outcomes and complications during hospitalisation. Additionally, the longitudinal analysis revealed a higher incidence of decreased lung volume at 3 months in patients with critical outcomes. These findings revealed that CT-derived lung volume measurement holds significant value as an easily assessable indicator for predicting the severity of COVID-19. Furthermore, this tool identifies populations at risk of COVID-19 sequelae by evaluating residual lung volume reduction.

In this study, CT-derived lung volume on admission was associated with multiple outcomes, such as critical outcome, oxygen requirement, IMV requirement and death. It was associated with critical outcomes independently of AI-based pneumonia and pneumonia qualitative assessment, suggesting that diminished lung volume itself may correlate with the severity of the illness. Qualitative evaluation of lung shadows by radiologists using chest CT can predict infection severity in the early stages of an outbreak.32 However, variability in evaluations among radiologists is problematic.33 An advantage of automated lung volume measurement is its robust quantitative evaluation compared with qualitative evaluation by a radiologist. The strength of this study is that we demonstrated an association between CT-derived lung volume and outcomes in a large sample size and a multicentre long-term patient cohort. A previous report17 examined the association between lung volume and disease severity in COVID-19 but found no association with outcomes. This discrepancy may be partly attributed to substantial differences in sample sizes between the two studies. The large number of patients in this study also allowed for multivariate analysis, thereby demonstrating that CT-derived lung volume predicted the outcomes after adjusting for known risk factors. In routine clinical practice, it is difficult to perform CT imaging under the same conditions, particularly during pandemics. Although this study included multiple CT scanner models and imaging conditions, the results were clinically meaningful for predicting disease severity. SARS-CoV-2 has caused multiple epidemics worldwide.31 The clinical characteristics of patients vary across different periods owing to various factors, including advances in therapeutic agents, vaccines and viral variants. This study is the first to demonstrate the utility of CT-derived lung volume measurements for predicting severe disease, regardless of epidemic waves. These findings have significant clinical implications for future pandemics.

Owing to concerns regarding transmission, COVID-19 complicates the assessment of pulmonary function in the acute phase.34 Previous reports showed that CT-derived lung volume strongly correlated with the pulmonary function test results.8 9 Therefore, CT-derived lung volume could serve as a surrogate marker for the physiological assessment of COVID-19 in the acute phase. In this study, patients with decreased lung volume exhibited a higher incidence of multiorgan complications. Severe COVID-19 is associated with various posthospitalisation complications and poor outcomes.35 Our study demonstrated that lung volume measurements can predict these complications and enable early intervention.

The decreased lung volume in COVID-19 may be caused by several factors: (1) virus-induced damage to the type II alveolar epithelium, resulting in a reduction in surfactants that maintain alveolar expansion; (2) increased vascular permeability triggered by inflammatory cytokines, leading to alveolar collapse due to effusion and (3) inflammatory microvascular thrombi inducing capillary vasospasm in the alveolar space, resulting in augmented dead-space ventilation and subsequent alveolar collapse.16 In this study, the low TLCct (%pred) group exhibited elevated levels of markers suggestive of alveolar epithelial damage (LDH and KL-6), other inflammatory markers (CRP and procalcitonin) and coagulation system markers (D-dimer), which aligns with the proposed mechanism. Additionally, we revealed that the low TLCct (%pred) group exhibited clinical characteristics such as male sex, high BMI, hypertension and chronic kidney disease. These results are consistent with previous reports showing that these factors predict COVID-19 severity.36–39

COVID-19 is widely known to cause sequelae after the acute phase and is termed long COVID.40 Chest CT is a useful tool for evaluating structural sequelae of the lungs, and a qualitative evaluation by a radiologist revealed that 84% of patients exhibited residual lung shadows at 3 months.23 Additionally, these residual lesions correlate with pulmonary sequelae, worse pulmonary function and decreased exercise tolerability.23 41 This study revealed that the lung volumes of almost all patients recovered gradually from 0 to 3 months but were still lower in critical cases at 3 months. These results are consistent with a previous report showing that qualitative lung shadows improve over time21 and that severely ill patients are more likely to have residual qualitative lung shadows.22 Consequently, clinicians should be aware of a residual reduction in lung volume, especially in critically ill patients.

Our study has several limitations. First, the number of CT follow-up cases after 3 months was small owing to dropouts. In this study, we performed CT analyses after a 3-month interval based on the findings of a previous study, which showed that a 3-month period was optimal for evaluating residual shadows.42 Another study found that residual shadows at 3 months were correlated with physiological sequelae and dyspnoea.22–24 Notably, mild cases frequently interrupt hospital visits, resulting in potential bias. Second, the lung volume measurements were affected by the extent of inspiration. Total lung capacity is affected by sex and body size, and to redress these influences, the predicted TLCct% was calculated based on the predicted values. Nevertheless, in instances characterised by deficient inspiration attributable to hypoxia or other physiological symptoms, the total lung capacity may be underestimated, particularly in severe cases. Although the usefulness of respiratory-synchronised CT has been reported,43 we were unable to study it because COVID-19 is an emergency disease and many patients are admitted to the emergency department for CT during emergency admission. Third, lung volume measurements were affected by pre-existing lung diseases. Although reduced lung volume has been reported in pulmonary diseases such as ILD, COPD and asthma, the present study collected patient history only through interviews, the and we did not make a diagnosis based on imaging or pulmonary function test results of potential background diseases. Consequently, the impact of these background conditions could not be entirely ruled out and might have influenced the results.

In conclusion, low CT-derived lung volume at admission is associated with poor clinical outcomes and posthospitalisation complications in COVID-19 patients. CT-derived lung volumes recovered gradually over 3 months but were still lower in critical cases at 3 months. Future studies are warranted to investigate these effects in patients with decreased lung volume over time.

Data availability statement

Data are available on reasonable request. The datasets used and/or analysed in the current study are available from the corresponding author on reasonable request.

Ethics statements

Patient consent for publication

Ethics approval

This study involves human participants and was approved by the ethics committee of the Keio University School of Medicine (20200061) and related research institutions. All the participants provided informed consent. Participants gave informed consent to participate in the study before taking part.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Contributors Guarantor: SC. Conceptualisation: SOtake, YS, SC, NT, TM and YY. Data curation: SOtake, YS, SC, HN, TS, SA, KN, HTanaka, TF, MW and HTerai. Formal analysis: SOtake, YS, SC, NT, TM and TA. Methodology: SOtake, YS, SC, NT, YS, TM, TA, HN and YY. Project administration: SC and NT. Resource: SOtake, SC, HN, TS, SA, KN, HTanaka, TF, MW, HTerai, MS, SU, YK, NH, SS, SY, HTateno and KF. Software: YS, NT, TM and YY. Supervision: SC, NT, HN, YY, MJ, TH, YO, RK, MI, NH, AK, SI, SM, SOgawa, TK and KF. Visualisation: HTateno, TM, SC, NT and TA. Writing–original draft: HTerai, TM, SC, NT, YS and YY. Writing–review and editing: TA, HN, HTanaka, MS, SU, YK, NH, SS, SY, HTateno, MJ, TH, YO, RK, MI, NH, AK, SI, SM, SOgawa, TK and KF. All the authors have read and approved the final manuscript.

  • Funding This study was supported by the Japan Agency for Medical Research and Development (AMED) (grant numbers JP20nk0101612, JP20fk0108415, JP21jk0210034, JP21km0405211, JP21km0405217, JP21wm0325031, JP21fk0108563, JP21fk0108573, JP20fk0108452, JP21fk0108553, JP21fk0108431, JP22fk0108510, JP22fk0108513, JP22wm0325031, JP22fk0108573, JP23tm0524008), Japan Science and Technology Agency (JST) CREST (grant number JPMJCR20H2], Japan Science and Technology Agency (JST) PRESTO [grant number JPMJPR21R7], and the Ministry of Health, Labour and Welfare (grant number 20CA2054).

  • Competing interests None declared.

  • Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.