Skip to main content

Temporal trend and climate factors of hemorrhagic fever with renal syndrome epidemic in Shenyang City, China

Abstract

Background

Hemorrhagic fever with renal syndrome (HFRS) is an important infectious disease caused by different species of hantaviruses. As a rodent-borne disease with a seasonal distribution, external environmental factors including climate factors may play a significant role in its transmission. The city of Shenyang is one of the most seriously endemic areas for HFRS. Here, we characterized the dynamic temporal trend of HFRS, and identified climate-related risk factors and their roles in HFRS transmission in Shenyang, China.

Methods

The annual and monthly cumulative numbers of HFRS cases from 2004 to 2009 were calculated and plotted to show the annual and seasonal fluctuation in Shenyang. Cross-correlation and autocorrelation analyses were performed to detect the lagged effect of climate factors on HFRS transmission and the autocorrelation of monthly HFRS cases. Principal component analysis was constructed by using climate data from 2004 to 2009 to extract principal components of climate factors to reduce co-linearity. The extracted principal components and autocorrelation terms of monthly HFRS cases were added into a multiple regression model called principal components regression model (PCR) to quantify the relationship between climate factors, autocorrelation terms and transmission of HFRS. The PCR model was compared to a general multiple regression model conducted only with climate factors as independent variables.

Results

A distinctly declining temporal trend of annual HFRS incidence was identified. HFRS cases were reported every month, and the two peak periods occurred in spring (March to May) and winter (November to January), during which, nearly 75% of the HFRS cases were reported. Three principal components were extracted with a cumulative contribution rate of 86.06%. Component 1 represented MinRH0, MT1, RH1, and MWV1; component 2 represented RH2, MaxT3, and MAP3; and component 3 represented MaxT2, MAP2, and MWV2. The PCR model was composed of three principal components and two autocorrelation terms. The association between HFRS epidemics and climate factors was better explained in the PCR model (F = 446.452, P < 0.001, adjusted R 2 = 0.75) than in the general multiple regression model (F = 223.670, P < 0.000, adjusted R 2 = 0.51).

Conclusion

The temporal distribution of HFRS in Shenyang varied in different years with a distinctly declining trend. The monthly trends of HFRS were significantly associated with local temperature, relative humidity, precipitation, air pressure, and wind velocity of the different previous months. The model conducted in this study will make HFRS surveillance simpler and the control of HFRS more targeted in Shenyang.

Peer Review reports

Background

Hemorrhagic fever with renal syndrome (HFRS), with characteristics of fever, hemorrhage, kidney damage and hypotension, is an important infectious disease caused by different species of hantaviruses. In hantavirus-endemic areas, HFRS outbreaks have occurred among farmers and others who have close contact with excreta of infected rodents [1, 2]. HFRS has been recognized as a notable public health problem in China [3]. It is considered that the number of HFRS cases in China accounts for 90% of the total cases worldwide [4, 5]. At present, HFRS is endemic in 28 of 31 provinces, autonomous regions, and metropolitan areas in mainland China [6]. Although some prevention and control measures such as scientific rodent control, vaccination and environmental management have been performed, HFRS remains a serious public health problem with about 20,000-50,000 human cases annually in mainland China [7].

As a rodent-borne disease with a seasonal distribution, external environmental factors including climate factors may play a significant role in its transmission. Studies in different areas of China and other countries have suggested that climate factors, such as temperature, precipitation, and relative humidity, may influence the incidence of HFRS [8–10]. The role of climate factors in such transmission may vary, given the large variations in climate types, ecological characteristics, population immunity, public health intervention measures, and socioeconomic status in the different regions.

Liaoning Province is one of the most seriously affected areas with the highest incidence during the years 2004 and 2005 [11]. Since the first HFRS case in 1958, Shenyang, the capital city of Liaoning Province, has been among the forefront of the province in recent years, with all its 13 districts/counties being susceptible to the disease [12]. Through the use of a 6-year (2004-2009) record of HFRS cases in Shenyang, this study aimed to characterize the seasonal patterns of HFRS distribution and explore associations between climate factors and HFRS transmission.

Methods

Study area

The study area covers Shenyang, the capital city of Liaoning Province, located between latitude 41°48' and 41°50' north and longitude 123°23' and 123°40' east (Figure 1). The city comprises 9 districts and four counties with a total land area of 12,948 km2 and a population of about 7,860,000 in 2009. Shenyang has a temperate climate with dry winters and wet, hot summers. The annual temperature ranges from -29 to +36°C and the annual mean temperature is about 7.5°. The annual rainfall is typically between 600 and 800 mm. HFRS is a rodent-borne disease. Therefore, rodent control is important for disease control. The major rodent species in Shenyang are Apodemus agrarius and Rattus norvegicus. The widespread distribution of rodents makes their effective control difficult.

Figure 1
figure 1

The black dot indicates the location of Shenyang in China.

Data collection and management

Surveillance data

In this study, records on HFRS cases from 2004 to 2009 were obtained from the National Notifiable Disease Surveillance System. The disease surveillance data used in this study was already anonymous. We state that the China CDC provides the permission to use the data. All HFRS cases were first diagnosed by clinical symptoms. Patient blood samples were collected in hospitals and sent to the laboratory of Liaoning Provincial Center for Disease Control and Prevention (CDC) for serological confirmation. Finally, the data were collected by case number according to the sampling results. There might be admission rate bias in the disease report, but this has been reduced as much as possible. According to the legislation, physicians in hospitals must report every case of HFRS to the local health authority within 12 hours. The health authority reports these cases to the next level of the organization every month. Therefore, it is believed that the degree of compliance in disease notification over the study period was consistent. Demographic data were obtained from the statistical yearbooks of Liaoning Province, which were compiled by the Liaoning Provincial Bureau of Statistics.

Meteorological data

From 2004 to 2009, monthly climate data from approximately 700 surveillance stations in mainland China were collected from the China Meteorological Data Sharing Service System http://cdc.cma.gov.cn/index.jsp. The climate factors included monthly mean temperature (MT), monthly mean maximum temperature (MaxT), monthly mean minimum temperature (MinT), monthly mean relative humidity (RH), monthly mean minimum relative humidity (MinRH), monthly accumulative precipitation (AP), monthly mean air pressure (MAP), monthly mean wind velocity (MWV), and monthly sunshine duration (SD).

Temporal distribution analysis

The annual HFRS incidence from 2004 to 2009 was calculated and plotted to observe annual fluctuations in Shenyang. Cumulative HFRS cases for each month from 2004 to 2009 were also calculated to observe seasonal fluctuations.

Cross-correlation and autocorrelation analysis

A cross-correlation analysis was conducted to detect the effects of climate factors on HFRS transmission with a lag time of 6 months. The cross-correlation could be observed if the absolute value of cross-correlation coefficient (CCF) was two times larger than the standard error (SE). An autocorrelation analysis was performed to explore whether the monthly HFRS cases were affected by the cases in previous months, by using the Ljung-Box Q test, autocorrelation coefficient (AC) and partial autocorrelation coefficient (PAC).

Multivariate time series analysis

Principal component analysis (PCA) was constructed by using climate data from 2004 to 2009 to extract principal components. The extracted principal components and autocorrelation terms of monthly HFRS cases were added into a multiple regression model called principal components regression model (PCR) to quantify the relationship between climate factors, autocorrelation terms and transmission of HFRS. The PCR model was compared to a general multiple regression model conducted with only climate factors as independent variables.

Ethical review

The present study was reviewed by the research institutional review board of Shandong University and China CDC, and it is found that utilization of disease surveillance and meteorological data did not require oversight by an ethics committee.

Results

Descriptive analysis

There were 1,931 cases in total in Shenyang over the study period. HFRS cases were reported in all districts/counties of Shenyang during this period, with annual average incidences ranging from 0.15 to 34.65 per 100,000.

Temporal distribution of HFRS in Shenyang

The annual HFRS incidence in Shenyang varied from 2004 to 2009. The highest incidence of 9.11 cases per 100,000 persons occurred in 2004, and the lowest incidence of 1.75 cases per 100,000 occurred in 2009. A distinctly declining temporal trend in annual HFRS incidence was identified (Figure 2). HFRS cases were reported every month, and the two peak periods occurred in spring (March to May) and winter (November to January), during which nearly 75% of the HFRS cases were reported (Figure 3).

Figure 2
figure 2

Annual incidence of HFRS from 2004 to 2009 in Shenyang, China.

Figure 3
figure 3

Monthly distribution of HFRS cases in Shenyang, China.

Cross-correlation between monthly HFRS incidence and climate factors

As shown in Table 1 (See additional file 1), the monthly mean relative humidity (RH0) and monthly mean minimum relative humidity (MinRH0) were negatively correlated with the monthly number of HFRS cases in Shenyang, while monthly mean wind velocity (MWV0) was positively correlated with the monthly number of HFRS cases in the same month. During the previous month, the monthly number of HFRS cases was negatively correlated with climate factors, including the monthly mean temperature (MT1), monthly mean maximum temperature (MaxT1), monthly mean minimum temperature (MinT1), monthly mean relative humidity (RH1), monthly accumulative precipitation (AP1), as well as positively correlated with monthly mean wind velocity (MWV1). During the previous 2 months, the monthly HFRS cases were negatively correlated with the monthly mean temperature (MT2), monthly mean maximum temperature (MaxT2), monthly mean minimum temperature (MinT2), and monthly mean relative humidity (RH2), but positively correlated with the monthly mean air pressure (MAP2) and monthly mean wind velocity (MWV2). During the previous 3 months, the monthly HFRS cases were negatively correlated with the monthly mean temperature (MT3), and monthly mean maximum temperature (MaxT3), but positively correlated with the monthly mean air pressure(MAP3).

Autocorrelation of monthly number of HFRS cases

The P value of the Ljung-Box Q statistic of each lagged month was < 0.05. The absolute value of the autocorrelation coefficient and partial autocorrelation coefficient during the first two lagged months (Lag1 and Lag2) was greater than that of other lagged months, which indicated that there was a strong autocorrelation of monthly HFRS cases during the first two lagged months (Table 2, See additional file 2).

Model evaluation

The PCA was conducted by using variables RH0, MinRH0, MWS0, MT1, MaxT1, MinT1, RH1, AP1, MWS1, MT2, MaxT2, MinT2, RH2, MAP2, MWS2, MT3, MaxT3, and MAP3. Three principal components were extracted with a cumulative contribution rate of 86.06% (Table 3, See additional file 3). The loadings of the three principal components of each variable were calculated (Table 4 See additional file 4). Component 1 represented MinRH0, MT1, RH1, and MWV1; component 2 represented RH2, MaxT3, and MAP3; and component 3 represented MaxT2, MAP2, and MWV2. Consequently, the three principal components and the two autocorrelation terms of monthly HFRS cases (Lag1, and Lag2) were added into the multiple regression model. The following model gave the best results:

Y = - 0 . 078 c o m p o n e n t 1 - 0 . 168 c o m p o n e n t 2 - 0 . 372 c o m p o n e n t 3 + 0 . 395 L a g 1 + 0 . 378 L a g 2

The association between HFRS epidemics and climate factors was better explained in the PCR model (Figure 4, F = 446.452, P < 0.001, adjusted R 2 = 0.75) than in the general multiple regression model (F = 223.670, P < 0.000, adjusted R 2 = 0.51).

Figure 4
figure 4

Fitted (green) and actual (blue) monthly numbers of HFRS cases in Shenyang, China.

Discussion

The results of this study confirmed that the temporal distribution in Shenyang varied by months and years, which indicated that public health resource allocations should focus on the months with the highest HFRS incidence.

Previous studies have debated the association between HFRS epidemics and temperature in areas that are severely affected by HFRS around the world; especially in many areas of China [8, 9, 13, 14]. These studies have indicated that temperature is an important factor for HFRS epidemics in China. Some studies have found that temperature has a positive effect on HFRS [8, 13], whereas a negative effect has been found in other studies [9, 14, 15]. The results of the present study showed that MT1, MaxT2 and MaxT3 were negatively associated with the monthly HFRS cases in Shenyang. Based on the results of this study, effective surveillance of HFRS epidemics can be conducted by monitoring the fluctuations of MT1, MaxT2 and MaxT3 in the area, and targeted countermeasures can be taken ahead of time.

Lin et al. have indicated that relative humidity might make some contribution to the high epidemic occurrence of HFRS in Liaoning Province [11]. Liu et al. have applied a case-crossover design and conditional logistic regression to analyze the relationship between incidence rates of HFRS and meteorological variables in one national surveillance area for HFRS in Shandong Province [9]. They have found that the MRH, AP, and WV were positively associated with HFRS incidence. In contrast, some researchers have reported that the associations between HFRS epidemics and precipitation and MRH are not significant [16]. The present study found that relative humidity factors, including MinRH0, RH1 and RH2, affected HFRS negatively, which was consistent with previous studies in China [8, 14, 17]. In addition, MWV1, MWV2, MAP2 and MAP3 had positive effects on the monthly HFRS cases in cross-correlation analysis, but negative effects in multiple regression, which could have been due to the combined effects with other elements of the same principal component.

Excessive precipitation could have a negative impact on rodents by destroying their habitats in Eastern China [18, 19]. In addition, frequent rainfall may decrease the likelihood of rodent-to-rodent contact, rodent-to-human contact, and virus transmission due to decreased rodent activity and reduced human exposure [18]. The impact of precipitation on HFRS in this study was consistent with previous studies. Although AP1 was not included in the model as a predictor of HFRS, it still had a negative effect on HFRS, according to cross-correlation analysis. The absence of AP in the model was due to its relative small loading. Therefore, when forecasting the trend of HFRS in Shenyang, AP1 still should be taken into consideration. In addition, the monthly SD was not related to the monthly number of HFRS cases. This finding was not consistent with some other studies [9, 14, 17], which could have been due to different study areas. This result indicated that it is not necessary to consider SD as a predictor of HFRS in Shenyang.

In the present study, multiple regression was constructed to detect the relationship between HFRS epidemics and climate factors in Shenyang, by incorporating climate factors of previous months based on a PCA. This model was reliable and had a good fit (adjusted R 2 of fitting = 0.75). Previous studies also have used other methods to identify climate-related risk factors of HFRS, such as Bayesian discriminant analysis [17], structure equation model [20], generalized linear model, and generalized additive model [21]. To obtain the expected results, it was necessary to choose the most suitable method according to available data.

Despite the insights gained, the limitations of our study should also be acknowledged. First, this study examined only the effect of climate factors and autoregression on HFRS epidemics, without taking into account rat density, socioeconomic circumstances, and countermeasures. Transmission of hantaviruses to humans is also related to other factors such as human activities, farming patterns, and rodent abatement strategies [22]. Second, the model in this study was constructed by using relatively short-term data for a 6-year period, which could have made the results less reliable. Third, this study didn't consider the impact of the vaccination. Vaccination plays an important role in the control of hemorrhagic fever with renal syndrome. A comprehensive intervention based on vaccination has been adopted in Shenyang since 2004 [23]. Vaccination partly contributed to the obvious decline of the annual HFRS incidence in 2007 compared to that in 2006. Climatic factors also contributed to the lower incidence in 2007 compared with that in 2006, because the amount of precipitation of 2007 is much higher than that in 2006. As the vaccination intervention is done in specific areas, it is difficult to take the impact of the vaccination into the model. HFRS incidence in the past two decades was much higher in the frigid-temperate zone, mostly in northeastern China, followed by the warm-temperate zone, with a lower incidence seen in southeastern China, where there are higher temperatures, higher humidity, and greater precipitation [6, 24]. This was consistent with the results of our analysis of climate factors. These results lay a foundation for HFRS surveillance and control in this area. As a next step, it will be advantageous to explore the method for assigning fluctuant coefficients to climate factors, to improve upon the current methods. Furthermore, it is necessary to conduct further research on how socioeconomic factors affect HFRS transmission in Shenyang, which will assist in amending prevention and control strategies.

Conclusion

This study shows that the temporal distribution of HFRS in Shenyang varied in different years(2004-2009), with a distinctly declining trend. The monthly trends for HFRS were significantly associated with local temperature, relative humidity, precipitation, air pressure, and wind velocity of the previous months. The model conducted in this study will make HFRS surveillance simpler and its control more targeted in Shenyang.

References

  1. Vapalahti K, Paunio M, Brummer-Korvenkontio M, Vaheri A, Vapaahti O: Puumala virus infections in Finland:increased occupational risk for farmers. Am J Epidemiol. 1999, 149: 1142-1151.

    Article  CAS  PubMed  Google Scholar 

  2. Glass GE, Childs JE, Korch GW, LeDuc JW: Association of intraspecific wounding with hantaviral infection in wild rats(Rattus norvegicus). Epidemiol Infect. 1988, 101: 459-472. 10.1017/S0950268800054418.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Zhang YZ, Xiao DL, Wang Y, Wang HX, Sun L, Tao XX, Qu YG: The epidemic characteristics and preventive measures of hemorrhagic fever with renal syndrome in China. Chin J Epidemiol. 2004, 25: 466-469.

    Google Scholar 

  4. Ulrich R, Hjelle B, Pitra C, Kuger DH: Emerging viruses: the case 'hantavirus'. Intervirology. 2002, 45: 318-327. 10.1159/000067924.

    Article  PubMed  Google Scholar 

  5. Simmons JH, Riley LK: Hantavirus: an overview. Com Med. 2002, 52: 97-110.

    CAS  Google Scholar 

  6. Yan L, Fang LQ, Huang HG, Zhang LQ, Feng D, Zhao WJ, Zhang WY, Li XW, Cao WC: Landscape elements and Hantaan virus-related hemorrhagic fever with renal syndrome, People's Republic of China. Emerg Infect Dis. 2007, 13: 1301-1306.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Fang LQ, Yan L, Liang S, de Vlas SJ, Feng D, Han X, Zhao W, Xu B, Bian L, Yang H, Gong P, Richardus JH, Cao WC: Spatial analysis of hemorrhagic fever with renal syndrome in China. BMC Infect Dis. 2006, 6: 77-10.1186/1471-2334-6-77.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Fang LQ, Wang XJ, Liang S, Li YL, Song SX, Zhang WY, Qian Q, Li YP, Wei L, Wang ZQ, Yang H, Cao WC: Spatiotemporal trends and climatic factors of hemorrhagic fever with renal syndrome epidemic in Shandong Province, China. PLoS Negl Trop Dis. 2010, 4: e789-10.1371/journal.pntd.0000789.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Liu J, Wang JZ, Xue FZ, Kang DM, Li SW: Association between incidence of hemorrhagic fever with renal syndrome(HFRS) and meteorological factors. Chinese Journal of Health Statistics. 2006, 23: 326-329.

    Google Scholar 

  10. Klempa B: Hantaviruses and climate change. Clin Microbiol Infect. 2009, 15: 518-523. 10.1111/j.1469-0691.2009.02848.x.

    Article  CAS  PubMed  Google Scholar 

  11. Lin HL, Liu QY, Guo JQ, Zhang JB, Wang JF, Chen HX: Analysis of the geographic distribution of HFRS in Liaoning province between 2000 and 2005. BMC Public Health. 2007, 7: 207-10.1186/1471-2458-7-207.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Chen HZ, Xu T, Dai MY, Zhang J: Study on host change trend of hemorrhagic fever with renal syndrome in Shenyang, 1984-2006. Chin J Dis Control & Prev. 2009, 13: 97-99.

    Google Scholar 

  13. Luo L, Wang M, Zhou RH, Liu YF, Qiu JC, Chen XS: An analysis of the relationship among the climatic indices and the HFRS and the mice density surveillance. Modern Preventive Medicine. 2005, 32: 205-206. 222

    Google Scholar 

  14. Wu RJ, Hu XS, Zheng YF, Liu GZ, Li L: The correlation analysis between hemorrhagic fever with renal syndrome(HFRS) and meteorological factors and forecast of HFRS. Chin J Vector Bio & Control. 2005, 16: 118-120.

    Google Scholar 

  15. Luo CW, Liu QY, Hou JL: Correlation analysis and regression model of epidemic factors of hemorrhagic fever with renal syndrome in Heihe City, Heilongjiang province. Disease Surveillance. 2009, 24: 118-120.

    Google Scholar 

  16. Zhang JJ, Lin T, Wu J, Zhang MY, Hu LY: Analysis of epidemic factors of hemorrhagic fever in Binzhou City from 1984 to 2004. Prev Med Trib. 2005, 11: 725-727.

    Google Scholar 

  17. Shen TF, Huang DS, Wu W, Guan P, Zhou BS: Application of Bayesian discriminant analysis in forecasting hemorrhagic fever with renal syndrome. Chin J Vector Bio & Control. 2009, 20: 147-150.

    Google Scholar 

  18. Bi P, Tong SL, Donald K, Parton K, Ni JF: Climatic, reservoir and occupational variables and the transmission of hemorrhagic fever with renal syndrome in China. Int J Epidemiol. 2002, 31: 189-193. 10.1093/ije/31.1.189.

    Article  PubMed  Google Scholar 

  19. Bi P: El Nino and incidence of hemorrhagic fever with renal syndrome in China. JAMA. 2003, 289: 176-177. 10.1001/jama.289.2.176-d.

    Article  PubMed  Google Scholar 

  20. Guan P, Huang DS, He M, Shen TF, Guo JQ, Zhou BS: Investigating the effects of climatic variables and reservoir on the incidence of hemorrhagic fever with renal syndrome in Huludao City, China: a 17-year data analysis based on structure equation model. BMC Infect Dis. 2009, 9: 109-10.1186/1471-2334-9-109.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Hu W, Mengersen K, Bi P, Tong S: Time-series analysis of the risk factors for haemorrhagic fever with renal syndrome: comparison of statistical models. Epidemiol Infect. 2007, 135: 245-252. 10.1017/S0950268806006649.

    Article  CAS  PubMed  Google Scholar 

  22. Schmaljohn C, Hjelle B: Hantaviruses: a global disease problem. Emerg Infect Dis. 1997, 3: 95-104. 10.3201/eid0302.970202.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Bai S, Li P: Evaluation on effect of HFRS vaccination in Shenyang from 2004 to 2006. Chin J Dis Control Prev. 2009, 13: 211-212.

    Google Scholar 

  24. Chen HX, Qiu FX, Zhao XQ, Luo CW, Li XQ: Characteristics of the distribution of epidemic season of hemorrhagic fever with renal syndrome in different regions and different years in China. Chin J Clin Virol. 1994, 8: 197-203.

    Google Scholar 

Pre-publication history

Download references

Acknowledgements

This study was supported by the National Basic Research Program of China (973 Program) (Grant No.2012CB955500) and Special Infectious Diseases Program of the Ministry of Science & Technology, China(Grant No.2008ZX10004-010). We would like to thank Miss Fei Zhang and Miss Jingyi Duan for their earnest assistance in revising and editing this manuscript.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Baofa Jiang or Qiyong Liu.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

XL, QL, BJ were involved in the research design, execution and write up of the first draft of the manuscript. XL, QL and WG contributed to data collection and statistical analysis. All authors contributed to the writing of the manuscript and approved the submitted version of the manuscript.

Electronic supplementary material

12879_2010_1664_MOESM1_ESM.DOC

Additional file 1: Table 1. Cross correlation between monthly HFRS cases and climate factors in Shenyang, China. (DOC 60 KB)

12879_2010_1664_MOESM2_ESM.DOC

Additional file 2: Table 2. Autocorrelation and partial correlation of monthly HFRS cases in Shenyang, China. (DOC 34 KB)

Additional file 3: Table 3. Total variance of relative variables explained by PCA. (DOC 49 KB)

Additional file 4: Table 4. PCA, matrix of component loadings. (DOC 38 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Liu, X., Jiang, B., Gu, W. et al. Temporal trend and climate factors of hemorrhagic fever with renal syndrome epidemic in Shenyang City, China. BMC Infect Dis 11, 331 (2011). https://doi.org/10.1186/1471-2334-11-331

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2334-11-331

Keywords