Article Text

Original research
Development of influenza-associated disease burden pyramid in Shanghai, China, 2010–2017: a Bayesian modelling study
  1. Jing Li1,2,3,
  2. Chunfang Wang4,
  3. Luanqi Ruan5,
  4. Shan Jin4,
  5. Chuchu Ye5,
  6. Huiting Yu4,
  7. Weiping Zhu5,
  8. Xiling Wang1,6
  1. 1School of Public Health, Fudan University, Shanghai, Shanghai, China
  2. 2Renal Division, Peking University Shenzhen Hospital, Shenzhen, Guangdong, China
  3. 3Clinical Research Academy, Peking University Shenzhen Hospital, Shenzhen, Guangdong, China
  4. 4Department of Vital Statistics, Shanghai Municipal Center for Disease Control and Prevention, Shanghai, China
  5. 5Research Base of Key Laboratory of Surveillance and Early Warning on Infectious Disease, Pudong New Area Center for Disease Control and Prevention, Shanghai, China
  6. 6Shanghai Key Laboratory of Meteorology and Health, Shanghai, China
  1. Correspondence to Xiling Wang; erinwang{at}fudan.edu.cn; Weiping Zhu; wpzhu{at}pdcdc.sh.cn

Abstract

Objectives Negative estimates can be produced when statistical modelling techniques are applied to estimate morbidity and mortality attributable to influenza. Based on the prior knowledge that influenza viruses are hazardous pathogens and have adverse health outcomes of respiratory and circulatory disease (R&C), we developed an improved model incorporating Bayes’ theorem to estimate the disease burden of influenza in Shanghai, China, from 2010 to 2017.

Design A modelling study using aggregated data from administrative systems on weekly R&C mortality and hospitalisation, influenza surveillance and meteorological data. We constrained the regression coefficients for influenza activity to be positive by truncating the prior distributions at zero.

Setting Shanghai, China.

Participants People registered with R&C deaths (450 298) and hospitalisations (2621 787, from 1 July 2013), and with influenza-like illness (ILI) outpatient visits (342 149) between 4 January 2010 and 31 December 2017.

Primary outcome measures Influenza-associated disease burden (mortality, hospitalisation and outpatient visit rates) and clinical severity (outpatient–mortality, outpatient–hospitalisation and hospitalisation–mortality risks).

Results Influenza was associated with an annual average of 15.49 (95% credibility interval (CrI) 9.06–22.06) excess R&C deaths, 100.65 (95% CrI 48.79–156.78) excess R&C hospitalisations and 914.95 (95% CrI 798.51–1023.66) excess ILI outpatient visits per 100 000 population in Shanghai. 97.23% and 80.24% excess R&C deaths and hospitalisations occurred in people aged ≥65 years. More than half of excess morbidity and mortality were associated with influenza A(H3N2) virus, and its severities were 1.65-fold to 3.54-fold and 1.47-fold to 2.16-fold higher than that for influenza A(H1N1) and B viruses, respectively.

Conclusions The proposed Bayesian approach with reasonable prior information improved estimates of influenza-associated disease burden. Influenza A(H3N2) virus was generally associated with higher morbidity and mortality, and was relatively more severe compared with influenza A(H1N1) and B viruses. Targeted influenza prevention and control strategies for the elderly in Shanghai may substantially reduce the disease burden.

  • epidemiology
  • statistics & research methods
  • public health

Data availability statement

Data are available upon reasonable request. Data are available upon reasonable request to the corresponding author.

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.

Strengths and limitations of this study

  • We comprehensively assessed the impact of influenza on both morbidity and mortality in subtropical Shanghai by integrating multiple data sources of influenza surveillance, hospital records and death registration.

  • The developed ‘influenza pyramid’, together with the risks between layers, provided key parameters for handy estimation of influenza-associated disease burden in Shanghai as well as in other high-income countries/regions/cities.

  • The proposed Bayesian approach with reasonable prior information improved statistical estimates of influenza-associated disease burden.

  • Our results may be confounded by other cocirculating respiratory viruses which were not included in the national virological surveillance scheme in China.

Introduction

Annual influenza epidemics result in substantial morbidity and mortality. Influenza virus infections can cause a wide spectrum of diseases, from mild to severe illness requiring hospitalisation, and at times even death.1–4 The influenza burden can be presented using the disease pyramid with multiple tiers, including asymptomatic infections, non-medically attended illnesses, medically attended illnesses, hospitalisations, intensive care unit admissions and deaths.5 The ‘influenza pyramid’ provides a full perspective of the disease burden associated with influenza and allows identification of relationships between each layer of the outcome. However, the construction of the pyramid may require merging data gathered from different surveillance systems. Considering the availability of time series data, very few studies have attempted to assess the overall burden of influenza and the estimation tended to focus on only one or two particular levels of the pyramid.6–9

The disease burden of influenza is difficult to quantify directly because influenza infections are rarely laboratory-confirmed, or specifically coded as influenza-related in hospital discharge records or death certificates. Statistical modelling techniques have been developed to estimate morbidity and mortality attributable to influenza.3 10 11 Pneumonia and influenza (P&I), respiratory and circulatory disease (R&C) and all-cause are often chosen as outcome variables in modelling methods. P&I is deemed to be the most specific outcome measure, but this category excludes other respiratory and circulatory diseases exacerbated by influenza. All-cause is too broad and sacrifices specificity; thus, R&C optimises the balance between sensitivity and specificity.12–14

Negative estimates can be produced when statistical modelling techniques are applied to estimate morbidity and mortality attributable to influenza.15–17 Based on the prior knowledge that influenza viruses are hazardous pathogens and have adverse health outcomes of R&C,3 4 18 we have developed an improved model incorporating Bayes’ theorem to construct the pyramid of influenza-associated disease burden. We presented the estimates of influenza-associated deaths and hospitalisations due to R&C and influenza-like illness (ILI) outpatient visits in Shanghai, China, from 2010 to 2017.

Methods

Influenza surveillance data

Weekly ILI surveillance and influenza virological surveillance data from 2010 to 2017 were obtained from the Pudong New Area Center for Disease Control and Prevention. Four hospitals conducted year-round surveillance of influenza, including two national sentinel hospitals (Dongfang Hospital and Zhoupu Hospital) and two municipal sentinel hospitals (Pudong Hospital and Eastern Division of Renji Hospital) in Pudong New Area, Shanghai. As we do not have access to raw influenza surveillance data in Shanghai, we used influenza virus activity in Pudong New Area instead. There was a strong correlation of the weekly positive proportion of laboratory-confirmed influenza (LAB%) between the two data sources from 2010 to 2015 (r=0.80–0.93, all p§amp;lt;0.001) (online supplemental figure S1). Weekly LAB% data in Shanghai were digitised from the time series plot of influenza activity in published literature19 using the WebPlotDigitizer software (https://automeris.io/WebPlotDigitizer). ILI surveillance data included weekly numbers of total outpatient visits and age-specific ILI consultations (0–4, 5–14, 15–24, 25–59 and ≥60 years). As the elderly people with an age cutpoint of ≥65 years were used in our study, we reclassified these five age groups into 0–14, 15–64 and ≥65 years according to the age distribution of the permanent population in Pudong New Area. Influenza virological surveillance data included weekly numbers of total specimens tested and specimens positive for influenza A(H1N1) (referring to the 2009 pandemic strain A(H1N1)pdm09), A(H3N2) and B viruses.

Mortality and hospitalisation data

The death registration system of Shanghai is based on its household registration (Hukou) system.20 We obtained the weekly mortality data for registered residents in Shanghai for the years 2010–2017 from the Shanghai Center for Disease Control and Prevention. Weekly inpatient data from 2013 to 2017 were obtained from the Shanghai Municipal Health and Family Planning Commission. All medical institutions in Shanghai are required to apply the new edition of the inpatient medical record home page where full indicators were collected since 1 January 2013. We excluded the unstable data for the first half of the year (from the 1st week of 2013 to the 26th week of 2013) during the transition period between the new and old editions of inpatient electronic medical record. We only included inpatients from secondary and tertiary hospitals which had higher data quality than primary hospitals, and they accounted for more than 90% of the total inpatients in Shanghai. Shanghai provided medical treatments to many patients seeking medical care from other provinces; we only used the data of Shanghai residents for representativeness. Hospital admissions were greatly reduced around the traditional Chinese Spring Festival and National Day, which cannot reflect patients’ true healthcare needs, and thus we replaced the week containing the holidays with the average value of the previous and next weeks.21

We retrieved underlying causes of death and primary discharge diagnoses coded as respiratory disease (J00–J99), circulatory disease (I00–I99) and R&C (J00–J99 and I00–I99) in the International Classification of Diseases Codes V.10. We assumed that the excess R&C deaths or hospitalisations would approximate the totality of the influenza-associated deaths or hospitalisations. Data for death and hospitalisation were stratified by age groups (<65 and ≥65 years; 0–5, 6–17, 18–64 and ≥65 years, respectively).

Population and meteorological data

Age-specific annual population sizes, including registered population and permanent population, were obtained from the Shanghai Statistical Yearbook22 and National Bureau of Statistics of China, respectively. Annual outpatient visits (medicine and paediatrics) were obtained from China’s Health and Family Planning Statistical Yearbook.23 Daily temperature and dew point temperature data from 2010 to 2017 were obtained from the Shanghai Meteorological Bureau and were averaged into a weekly level. Absolute humidity was derived from temperature and dew point temperature.

Statistical analysis

Negative binomial regression models were fitted separately to weekly deaths, hospitalisations and ILI outpatient visits by age groups and diseases in a Bayesian framework using the Markov Chain Monte Carlo approach. We used LAB% as a proxy variable for influenza activity in the models and assumed a multiplicative association between weekly influenza activity and counts of deaths, hospitalisations or ILI outpatient visits, which has been widely applied in previous literature.3 4 24 Natural cubic splines of the calendar week and absolute humidity were added to adjust for time-varying confounders. Viral surveillance data lagged by 0–3 weeks. We selected the dfs and lag weeks based on the minimum deviance information criterion25 and curve fitting (online supplemental table S1, figure S2–S4). We constrained the regression coefficients for influenza activity to be positive by truncating the normal priors at zero as previous studies have demonstrated positive associations between influenza and R&C deaths, hospitalisations and ILI outpatient visits.3 4 26 More details of the statistical model are available in online supplemental appendix S1.

The influenza-associated excess deaths, hospitalisations and ILI outpatient visits were estimated as the difference between the predicted numbers under the model and the baseline numbers when influenza activity proxies were assumed to be zero. The influenza-associated excess mortality rates, hospitalisation rates and ILI outpatient visit rates were estimated as the numbers of excess deaths, hospitalisations and ILI outpatient visits per 100 000 population. Age-standardised rates were derived using the World (WHO 2000–2025) Standard Population as the reference.27 With the assumption that the proportion of ILI outpatient visits among medicine and paediatric outpatient visits in the surveillance hospitals of Pudong New Area is representative of Shanghai, we extrapolated the influenza-associated ILI outpatient visits in surveillance hospitals to the general population in Shanghai (see details in online supplemental appendix S2).28

Each model was run in 60 000 Markov Chain Monte Carlo iterations with the first 30 000 iterations discarded as burn-in. We estimated the clinical severity of influenza with different influenza virus types/subtypes and age groups, including outpatient–mortality risks, outpatient–hospitalisation risks and hospitalisation–mortality risks, as the ratios of influenza-associated excess mortality rates to ILI outpatient visit rates, hospitalisation rates to ILI outpatient visit rates and mortality rates to hospitalisation rates. The 95% credibility intervals (CrIs) for excess crude and age-standardised rates were based on the samples drawn from the posterior distributions. Estimates of the outpatient–mortality risks, outpatient–hospitalisation risks, hospitalisation–mortality risks and their 95% CrIs were also based on these posterior samples.

Analyses were performed using the ‘rjags’ package of R software V.3.5.2 (R Foundation for Statistical Computing, Vienna, Austria), an interface between R and JAGS software V.4.3.0 .29

Patient and public involvement

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

Results

Descriptive statistics

During the study period of 2010–2017, Pudong New Area tested 12 283 specimens for influenza viruses, of which 409 were positive for influenza A(H1N1) virus, 1427 for influenza A(H3N2) virus and 812 for influenza B virus (table 1). Influenza B virus was the predominant circulating virus in 2010–2011, and shifted to A(H3N2) virus in subsequent years (figure 1). In general, R&C mortality rates, R&C hospitalisation rates and ILI outpatient visit rates increased from 2010 to 2017 in Shanghai, with annual average rates of 392.32, 2407.98 and 5247.92 per 100 000 population, respectively (table 1).

Table 1

Annual summary of influenza activity, ILI outpatient visit rate, hospitalisation rate, and mortality rate in Shanghai, 2010–2017

Figure 1

Influenza activity by type/subtype in Pudong New Area, Shanghai, 2010–2017. LAB%, weekly positive proportion of laboratory-confirmed influenza.

Influenza-associated disease burden

During the study period, we estimated an annual average of 2223 (95% CrI 1300–3164) excess R&C deaths, 24 353 (95% CrI 11 805–37 934) excess R&C hospitalisations and 218 733 (95% CrI 1 90 897–244 722) excess ILI outpatient visits, attributable to influenza in Shanghai, corresponding to 3.95%, 4.18% and 17.43% of all R&C deaths, R&C hospitalisations and ILI outpatient visits, respectively. Accordingly, the estimated crude and age-standardised rates per 100 000 population were 15.49 (95% CrI 9.06–22.06) and 7.36 (4.76–10.25) excess R&C deaths, 100.65 (95% CrI 48.79–156.78) and 97.90 (70.74–129.02) excess R&C hospitalisations, and 914.95 (95% CrI 798.51–1023.66) and 974.65 (95% CrI 855.22–1100.21) excess ILI outpatient visits, respectively (table 2). The pyramids of influenza-associated disease burden in Shanghai are shown in figure 2.

Table 2

Mean annual influenza-associated excess mortality rates and hospitalisation rates due to R&Cs, and ILI outpatient visit rates (per 100 000 population) by influenza type/subtype and age group in Shanghai, 2010–2017

Figure 2

Pyramid of influenza burden by influenza virus types/subtypes and age groups in Shanghai, 2010–2017. ILI, influenza-like illness; R&C, respiratory and circulatory disease.

Influenza-associated R&C mortality rate for people aged ≥65 years was 81.03 (95% CrI 50.41–116.99) per 100 000 population, which was substantially higher than that for those aged <65 years (0.70, 95% CrI 0.23–1.37 per 100 000 population) (table 2). The age-specific influenza-associated R&C hospitalisation rates per 100 000 population showed a J-shaped pattern: highest among people aged ≥65 years (696.38, 95% CrI 396.31–1026.38), second highest among children aged 0–5 years (123.21, 95% CrI 46.56–212.08), followed by children aged 6–17 years (69.19, 95% CrI 46.65–91.48) and lowest among people aged 18–64 years (20.73, 95% CrI 7.02–39.76). Children aged 0–14 years had the highest rate of influenza-associated ILI outpatient visits (1430.91, 95% CrI 1096.85–1773.40 per 100 000 population), followed by people aged ≥65 (1096.79, 95% CrI 914.25–1261.59 per 100 000 population) and 15–64 years (781.92, 95% CrI 664.68–894.82 per 100 000 population). Among people aged ≥65 years, 97.23% of excess R&C deaths and 80.24% of excess R&C hospitalisations occurred, whereas of excess ILI outpatient visits, only 13.13% occurred in this age group.

Influenza A(H3N2) virus was generally associated with the highest rates of influenza infections, with all-age annual average rates of 10.69 (95% CrI 5.42–15.80) R&C deaths, 60.96 (95% CrI 19.95–100.85) R&C hospitalisations and 424.74 (95% CrI 340.90–510.55) ILI outpatient visits per 100 000 population, respectively (table 2 and figure 2). When broken down by diseases, for all ages and ≥65 years, the influenza-associated mortality rates for respiratory disease (6.90, 95% CrI 5.27–8.39, and 35.70, 95% CrI 27.78–44.16, per 100 000 population) were 1.40-fold to 1.48-fold lower than those for circulatory disease (9.64, 95% CrI 5.05–14.51, and 52.81, 95% CrI 28.66–79.37, per 100 000 population), while the influenza-associated hospitalisation rates for respiratory disease (72.48, 95% CrI 51.34–95.11, and 501.90, 95% CrI 406.51–600.23, per 100 000 population) were 2.39-fold to 2.40-fold higher than those for circulatory disease (30.20, 95% CrI 9.61–60.12, and 210.20, 95% CrI 63.08–404.70) per 100 000 population) (table 2). A similar J-shaped age pattern was also observed in influenza-associated respiratory hospitalisation rates, while for circulatory hospitalisation, the estimated rates increased with age.

Influenza-associated clinical severity

Influenza-associated outpatient–mortality risk, outpatient–hospitalisation risk and hospitalisation–mortality risk were estimated to be 1.69% (95% CrI 0.97%–2.48%), 11.09% (95% CrI 5.23%–17.39%) and 15.29% (95% CrI 7.68%–34.85%), respectively (online supplemental table S2); that is to say, we estimated that there were 1.69 excess R&C deaths and 11.09 excess R&C hospitalisations for every 100 excess ILI outpatient visits, and 15.29 excess R&C deaths for every 100 excess R&C hospitalisations in Shanghai. The estimated risks for people aged ≥65 years were higher than those for all ages except for hospitalisation–mortality risks (figure 3). When comparing these risks between different types/subtypes of influenza viruses, we found that influenza A(H3N2) virus had the greatest outpatient–mortality risks, outpatient–hospitalisation risks and hospitalisation–mortality risks across different age groups. Its severities were 1.65-fold to 3.54-fold and 1.47-fold to 2.16-fold higher than those for influenza A(H1N1) and B viruses, respectively (online supplemental table S2). However, the difference did not reach statistical significance. Also, we presented separately the influenza-associated hospitalisation–mortality risks for respiratory disease and circulatory disease (online supplemental figure S5). People with the circulatory disease had a 3.22-fold higher hospitalisation–mortality risk than those with respiratory disease (online supplemental table S2).

Figure 3

Clinical severity of influenza with different influenza virus types/subtypes and age groups. (A) Outpatient–mortality risk, (B) outpatient–hospitalisation risk and (C) hospitalisation–mortality risk.

Discussion

Although there were great year-to-year variations in influenza-associated mortality and morbidity from 2010 to 2017, we estimated that influenza contributed to an annual average of 15.49 (95% CrI 9.06–22.06) excess R&C deaths, 100.65 (95% CrI 48.79–156.78) excess R&C hospitalisations and 914.95 (95% CrI 798.51–1023.66) excess ILI outpatient visits per 100 000 population in Shanghai. Our crude or age-standardised estimates of influenza-associated mortality and morbidity are comparable with the corresponding estimates published for the USA (R&C death: 9.9 per 100 000 population,30 R&C hospitalisation: 88.4 per 100 000 population4 and ILI outpatient visit: 755–870 per 100 000 population)31 32; Argentina (R&C death: 21.3 per 100 000 population and R&C hospitalisation: 57 per 100 000 population)14; Hong Kong SAR (R&C death: 7.7 per 100 000 population)33; and five southern Chinese cities (R&C death: 8.8 per 100 000 population).34 In particular, Hong Kong SAR reported an approximately twofold higher hospitalisation burden than our estimates,33 which could be explained by the result of expanding health service capacities and lowering admission criteria due to the ageing population of Hong Kong in recent years.13 24 35 However, caution is needed when comparing these estimates between countries/regions/cities due to variations in the study period, predominant circulating influenza virus, age structure of population, influenza vaccination coverage and statistical model used. Influenza A(H3N2) virus accounted for the highest disease burden than the other two types, which has also been demonstrated by other studies.4 14 36 Influenza A(H3N2) virus is believed to have more frequent antigen drift and higher virulence than seasonal influenza A(H1N1) and B viruses.37 38 Additionally, our study revealed that influenza A(H3N2) virus was relatively more severe as well.

The elderly people were most severely affected by influenza, with the highest influenza-associated mortality and hospitalisation rates, and the second highest influenza-associated ILI outpatient visit rates. The ageing of the Shanghai population is an important contributor to the high rates of influenza-associated mortality and morbidity. For example, the influenza-associated mortality rates after age standardisation decreased by more than half. Between 2010 and 2017, the registered population aged ≥65 years in Shanghai increased by 28.70%, of about 3.18 million by 2017, accounting for 21.83% of the total population. Interestingly, people aged ≥65 years had lower hospitalisation–mortality risks than for all ages. A possible explanation might be that young adults are less likely to be hospitalised compared with elderly people, but some of them do develop serious illnesses and require hospitalisation, and thus have a greater risk of death.

Currently, population-based studies of influenza-associated mortality and morbidity have been increasingly reported in mainland China. Of these, two nationally representative studies estimated province-level influenza-associated respiratory mortality rates and ILI outpatient visit rates. However, studies on influenza-associated hospitalisation burden were limited to Beijing,39 40 Jiangsu41–45 and Hubei.46 Our study comprehensively assessed the impact of influenza in Shanghai, especially the hospitalisation burden, which has never been assessed before. The estimated ratios between each layer of the influenza pyramid provided key parameters for handy estimation of influenza-associated disease burden in Shanghai as well as in other high-income countries/regions/cities.

Performing regression analyses in a Bayesian framework had two advantages in our study. First, positive priors truncated the posterior distributions of regression coefficients for influenza activity to ensure the positivity of the generated samples. Otherwise, it was difficult to explain that influenza was protective to human health. Second, Bayesian posterior samples allowed us to estimate the 95% CrIs for these ratios, which could be challenging in a frequentist framework. Various proxies for influenza activity have been developed, and the most commonly used proxy variable was LAB%.47 Yu et al19 estimated influenza-associated P&I mortality, considering three different proxies for influenza activity: LAB number (positive number of laboratory-confirmed influenza), LAB% and LAB%×ILI% (the product of positive proportion of laboratory-confirmed influenza and ILI consultation rate). The results of this study demonstrated that LAB% proxy produced slightly higher estimates of influenza-associated mortality than LAB number and LAB%×ILI% proxies. Another study by Chan et al48 evaluated the performance of various virus proxy variables in estimating excess hospitalisation and concluded that LAB%×ILI% proxy provided more reliable estimates if age-specific virus data were not available. However, in general, different influenza activity proxies provided comparable estimates for influenza disease burden. We used LAB% proxy in this study which has been adopted to link influenza and health outcomes in previous literature.3 4 49

Nevertheless, our study had several limitations. First, despite a strong correlation of influenza virus activity between Pudong New Area and Shanghai, LAB% in Shanghai was slightly higher than that in Pudong New Area (median difference: 6.12%). Replacement of influenza activity in Pudong New Area may lead to an underestimation of the influenza burden in Shanghai from 2010 to 2017. Second, national ILI surveillance data were recorded by five predefined age groups; thus, influenza burden estimates for ILI outpatient visits and hospitalisations were difficult to compare across different age groups. Additionally, when extrapolating influenza-associated ILI outpatient visits and reclassifying ILI consultations, we were based on representativeness assumptions, but further evaluation may be needed. Third, R&C was chosen to optimise sensitivity-specificity balance, but influenza-associated R&C events may be slightly lower than the totality of the influenza-associated events. Fourth, we excluded hospitalisation data of inpatients from primary hospitals and with non-Shanghai residential addresses, considering the study’s quality and population representativeness. The former may underestimate and the latter may overestimate the influenza-associated hospitalisation burden. Fifth, our estimates of influenza-associated mortality and morbidity burden were not exactly based on the same population (registered population and permanent population respectively). Deaths of permanent population without Shanghai Hukou are usually registered at their place of household registration. Using the permanent population as the denominator may underestimate the influenza-associated mortality burden. Finally, we did not have data on cocirculating respiratory viruses, such as respiratory syncytial virus, which may have confounded the results.

In conclusion, the proposed Bayesian approach with reasonable prior information improved estimates of influenza-associated disease burden. Our study highlighted the substantial morbidity and mortality burden attributed to influenza in Shanghai. Influenza A(H3N2) virus was generally associated with more morbidity and mortality, and was relatively more severe compared with influenza A(H1N1) and B viruses. Targeted influenza prevention and control strategies for the elderly in Shanghai may substantially impact the disease burden.

Data availability statement

Data are available upon reasonable request. Data are available upon reasonable request to the corresponding author.

Ethics statements

Patient consent for publication

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

  • JL and CW contributed equally.

  • Contributors XW and WZ conceived, designed and supervised the study. JL, CW, WZ, LR, SJ, CY and HY participated in data collection. JL, XW and SJ conducted statistical analyses. JL and CW drafted the manuscript. XW and WZ commented on the data and their interpretation, and revised the content critically. All authors read and approved the final manuscript.

  • Funding XW is supported by the National Nature and Science Foundation of China (grant number 81602936). WZ, CY and LR are supported by the National Science and Technology Major Project (grant number 2018ZX10713001008). The funding sources had no role in the study design, data collection, data analysis or writing of the report.

  • Competing interests None declared.

  • 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.