- Research article
- Open Access
- Open Peer Review
Rates and risk factors for drug resistance tuberculosis in Northeastern China
BMC Public Healthvolume 13, Article number: 1171 (2013)
Drug-resistant tuberculosis (TB) has emerged as a major challenge toward TB control and prevention. In Lianyungang city, the extent and trend of drug resistant TB is not well known. The objective of the survey was to assess drug resistance pattern of MTB and risk factors for drug resistant TB, including multidrug resistance tuberculosis (MDR-TB) in this area.
We performed drug susceptibility testing on Mycobacterium tuberculosis (MTB) isolates with first- and second-line anti-tuberculosis drugs of 1012 culture positive TB cases by using the proportion method, who were consecutively enrolled from January 2011 to December 2012 in Lianyungang city, China. The patterns of drug resistance in MTB were investigated and multiple logistic regression analysis was performed to assess the risk factors for drug resistant TB.
Among the 1012 strains tested, 308 (30.4%) strains were resistant to at least one first-line drug; the prevalence of MDR-TB was 88 (8.7%), 5 (0.5%) strains were found to be extensively drug-resistant tuberculosis (XDR-TB). Female gender was a risk factor for MDR-TB (adjusted odds ratio (aOR) 1.763, 95% CI (1.060-2.934). The aged 28–54 years was significantly associated with the risk of MDR-TB with an aOR: 2.224, 95% CI (1.158-4.273) when compared with those 65 years or older. Patients with previous treatment history had a more than 7-fold increased risk of MDR-TB, compared with those never previously treated.
The burden of drug resistant TB cases is sizeable, which highlights an urgent need to reinforce control, detection and treatment strategies for drug resistant TB.
Tuberculosis (TB) remains a major global health problem, especially in high TB burden countries with large numbers of TB patients and poorer sanitation. TB ranks as the second leading cause of death among infectious diseases worldwide, following only the human immunodeficiency virus (HIV) . According to the 2012 WHO global TB report, China ranks as 2nd among the world’s 22 high burden countries with a TB incidence around 1 million, and the prevalence rate of TB was 104 per 100,000 patients (95% confidence interval [95% CI] 91–119) and the incidence rate was 75 per 100,000 patients (95% CI 66–85) in 2011 .
In a nationwide survey in China in 2007, the estimated multidrug-resistant tuberculosis (MDR-TB) rate was 5.7% for new cases and 25.6% for previously treated cases. Approximately 8% of the patients with MDR-TB had extensively drug-resistant tuberculosis (XDR-TB) . However, information on the prevalence of MDR and XDR-TB remains scant in the region. Drug-resistant Mycobacterium tuberculosis(MTB) strains, MDR-TB and XDR-TB strains may be the potential propellers for the spread of TB, and through now we put great efforts in TB control and treatment, the situation has not greatly improved. Despite poor compliance to chemotherapy duration, drug resistance of MTB would be another important factor for treatment failure among new TB cases , which may lead to a fall in successful TB cure rates .
Lianyungang city is located in the northeastern part of China, and covers an area of 7.5 thousand square kilometers with a population around 5 million in 2012. Information on the prevalence of MDR-TB and XDR-TB remains unknown in this region. Thus, we conducted this study to assess drug resistance patterns of MTB on first- and second-line anti-TB drugs in this region, and evaluated related risk factors of drug resistance among new and previously treated smear-positive TB patients.
Study population and isolates
In Lianyungang city, all newly registered patients with sputum smear-positive pulmonary tuberculosis were collected for strain identification and drug susceptibility testing (DST). From January 2011 to December 2012, a total of 1170 clinical isolates were collected from the sputum samples of pulmonary tuberculosis patients. In all surveys, each newly registered TB patient, positive on sputum smear microscopy, was interviewed by the clinician using the medical records to obtain the treatment history. The medical records include personal information, physical examinations, present illness, TB-related complaints, previous medical history, family history and PPD result. The treatment history of included cases was classified by medical staff into new and previously treated cases. New cases were defined as patients with tuberculosis who have never been treated with anti-TB drugs or received them for less than one month. Previously treated cases were defined as patients who have been treated for tuberculosis for at least one month. The definitions of new/previously treated cases referred to the WHO guidelines .
Strain identification and drug susceptibility test
The sputum samples were cultured and isolated on Lowenstein-Jensen (LJ) culture media which were prepared by the Baso biotechnology Ltd (Zhuhai, China). All the MTB isolates from positive culture were identified and subjected to DST. Identification of MTB was done using the p-Nitrobenzoic acid (PNB) method and the growth in LJ media containing PNB indicates that the bacilli do not belong to the MTB complex. Species other than MTB were excluded from the final analysis.
LJ media impregnated one anti-tuberculosis drug was used for DST, and the corresponding drug concentrations were as follows: 0.2 μg/ml for INH, 40 μg/ml for RMP, 4 μg/ml for SM, 2 μg/ml for EMB, 30 μg/ml for KM and 2 μg/ml for OFX. The LJ culture media were incubated at 37°C.They were read twice at the first week to detect contaminations and/or fast growth of atypical mycobacteria and again every week thereafter for slower growing bacteria; if no bacteria grew by 6 weeks, the result was recorded as negative. Contaminated slant cultures were discarded, and additional decontamination and culture were undertaken using a stored portion of the original sample. Resistance was expressed as the percentage of colonies that grew on the drug-containing media compared to those on control media. The growth of colonies in the drug containing plate was compared to the control plate as a proportion. If the bacterial growth on the media with the specific drug was ≥1% compared to the control, the strain was declared resistant to the specific drug; or it was defined as sensitive when the growth rate was < 1% compared to the control. Strains isolation, identification and DST were performed at the fourth people’s hospital of Lianyungang city. Procedures for external quality assurance for smear and culture were based on WHO guidelines . For internal quality assurance of DST, a standard H37Rv strain was included with each new batch of LJ media. External quality control for culture and DST was conducted by the provincial TB reference laboratory, which participates in the annual proficiency review of DST organized by the Hong Kong Supranational Tuberculosis Reference Laboratory and has passed each review since 2010.
MDR and XDR
The following WHO-recommended definitions were used for resistant cases:
MDR-TB was defined as TB with resistance to at least INH and RMP. XDR-TB defined as MDR-TB plus resistance to a fluoroquinolone and at least one second-line injectable agent: amikacin, kanamycin and/or capreomycin . In our study the XDR-TB was defined as TB with resistance to at least INH, RMP, KM and OFX.
For comparison of categorical variables, significance testing was done by χ2 test with by 2-sided Fisher exact test as appropriate. Associations between selected factors were estimated by computing odds ratios (ORs) and their 95% confidence intervals (CIs) from an unconditional logistic regression model. The criterion for significance was set at P < 0.05 based on a two-sided test. Analyses were conducted with SPSS version 13.0 (SPSS Inc., Chicago, IL, USA).
Informed written consent was obtained from all patients and the study was approved by the Institutional Review Board of Jiangsu Provincial Center for Disease Control and Prevention. Ethics were respected throughout the study period.
Bacterial strains and patient characteristics
During the study period, 1170 sputum smear positive TB cases were enrolled and 1093 clinical isolates reported a culture positive MTB. Meanwhile, 66 cultured negative cases and 11 contaminated cases were excluded. Of the 1093 culture positive samples, 59 cultures were distinguished as Non-MTB (NTM). 22 samples failed for DST assay were excluded. Thus, 1012 samples were included in the final analysis (Figure 1). The age of the 1012 patients ranged from 14 to 90 years (mean ± SD), 49.8 ± 20.6, media 54), 748 (73.9%) of them were male and 264 (26.1%) were female. Of the 1012 case-patients included in the study, 816 (80.6%) were new incidence and 196 (19.4%) were previously treated. 9 (0.9%) cases were migrants.859 (84.9%) of the 1012 cases came from rural areas (Table 1).
Prevalence of resistance to anti-tuberculosis drugs
This study observed most resistance of isoniazid 184 (18.2%) followed by streptomycin 141 (13.9%), rifampicin 114 (11.3%), ofloxacin 78 (7.7%), ethambutol 72 (7.1%), and kanamycin 42 (4.2%). Among the new cases, drug resistance was reported in 191(23.4%); INH resistance was reported in 91 (11.2%), RMP resistance in 48 (5.9%), EMB resistance in 30 (3.7%), SM resistance in 90 (11.0%), KM resistance in 28 (3.4%), OFX resistance in 42 (5.1%), Mono-resistance was found in 113 patients (13.8%), MDR in 34 (4.2%) and XDR in 0 (0%). Among the 196 previously treated cases, drug resistance was reported in 117 (59.7%); INH resistance was reported in 93 (47.4%), RMP resistance in 66 (33.7%), EMB resistance in 42 (21.4%), SM resistance in 51 (26.0%), KM resistance in 14 (7.1%), OFX resistance in 36 (18.4%), Mono-resistance was present in 22 (11.2%), MDR in 54 (27.6%) and XDR in 5 (2.6%) patients, respectively (Table 2).
Factors associated with MDR-TB
The univariate analysis revealed a higher proportion of female gender (cOR 1.704, 95% CI 1.077–2.697), individuals aged 28–54 years (cOR 2.914, 95% CI 1.502–5.652), and patients with previous TB treatment (cOR 8.746, 95% CI 5.495–13.922) amongst MDR-TB patients. The final multivariable logistic regression model identified female gender (aOR 1.763, 95% CI 1.060–2.934), aged 28–54 years (aOR 2.224, 95% CI 1.158–4.273) and previous history of TB treatment (aOR 8.910, 95% CI 5.478–14.492), to be associated with MDR-TB (Table 3).
Information on anti-tuberculosis drug resistance levels is an essential management tool for evaluating the performance of national TB control programmers (NTPs). Resistance in previously treated cases is an indicator of current treatment practices in the community. Drug resistance in new cases reflects transmission of disease with resistant bacilli .
The prevalence of drug resistance in previously treated cases was higher than new cases for each drug alone as well as for all six drugs, as reported previously [2, 7, 8]. This suggests that retreatment is deficient and poses a threat to continued transmission, which has not yet manifested itself among new patients. The prevalence of MDR-TB was 4.2% in new cases and 27.6% in previously treated cases. As reported in a national survey of China, the rates for MDR-TB in new and previously treated cases were 5.7% and 25.6%, respectively . The MDR-TB prevalence among previously treated cases was modestly higher than the national average, and among new cases the rate was slightly lower. MDR-TB was observed 2.8% -14.7% of new cases and 9.7% -34.3% of previously treated cases in China [2, 9–12]. In Korea, resistance to at least one first-line drug was identified in 11.7% of new cases and 41.6% of previously treated cases. MDR-TB was detected in 3.9% of new cases and 27.2% of previously treated cases. The proportion of XDR-TB among MDR-TB patients was 16.7% (9/54) . In Japan, The prevalence of MDR in new and previously treated cases was 0.7% and 9.8% respectively . The overall MDR-TB rate was 4.3%; 2.5% in new cases and 13.9% in previously treated cases in Pakistan . MDR-TB was 0.2% of new cases and 3.4% of previously treated cases in Madagascar . In Cambodia, no single MDR-TB was found among the new cases and 3.1% in previously treated cases . In Yemen, The prevalence of MDR in new and previously treated cases was 3% and 9.4% respectively .
Compared MDR-TB patients, patients with XDR-TB are more likely to die or have treatment failure [19, 20]. The prevalence of XDR-TB among combined cases was 0.5%; it was 0% and 2.6% (5 cases) in new and previously treated cases, respectively. XDR-TB was distributed widely, albeit sporadically. This might imply that China’s treatment strategy has improved, but that community TB control is still insufficient . The prevalence of XDR-TB in new cases was slightly lower than the national average (0.5%). However, XDR-TB prevalence among previously treated cases was a little higher than the national average (2.1%) . These findings points to the need for interventions that will increase continuity of treatment and reduce the rate of treatment default, especially among patients treated within the hospital system .
As reported, 11% of new cases of tuberculosis and 16% of previously treated cases were resistant to either INH or RMP (but not both) before they received standard first-line short-course treatment . Resistance to these first-line drugs were usually not detected because culturing and drug susceptibility testing are not routinely performed at the local tuberculosis clinics. The use of standard first-line drugs in the treatment of these patients may increase the risk of relapse, treatment failure and acquired drug resistance [23–25].
Male or female TB patients could have different levels of risk for drug resistance due to differences in access to health-care services or exposure to other risk factors . In the present study, the resistance rate to any drug was similar in female and male, but females were more likely to have MDR-TB than males, with an OR of 1.763 (95% CI: 1.060-2.934), also reported in previous studies [8, 26–29]. The reasons for the association between female gender and MDR-TB are not well known. We hypothesize that this association could be related to the fact that women spend a long period of time caring for men and others with MDR-TB both in households and in healthcare settings in China where the majority of health care workers are female [26, 27]. TB diagnosis and DST may be delayed among female patients, making treatment more difficult and inappropriate chemotherapeutic regimens more likely [30, 31]. It is also possible that referring patients to the hospital may be different for male and female patients and the providers may feel reluctant to refer female patients with less complicated infection to a tertiary care center . Discovering gender disparities associated with the risks of MDR-TB could provide insight into the development of targeted measures and improve access to health care and reduce the risk of acquiring drug-resistance.
We also found that age was a significant factor in the development of drug resistance. Resistance to any drug and MDR was highest in those aged 28–54 years. The association between age and the risk of MDR-TB is not well established in the literature as different studies use different cut-off points for age groups. However, it was also reported that MDR-TB patients were more likely to be younger than 65 years [33, 34]. We assumed that age-related difference in treatment adherence might be a possible explanation, as patients at 28–54 years old were often occupied by study, work or other activities on a daily basis, in contrast with the more sedentary lifestyle of elderly patients .
Many other risk factors for drug resistance tuberculosis and MDR-TB have been identified in recent studies: such as irregular treatment , urban residence , non-permanent residents , urban migration , lack of a sewage system in the home , frequent travelers , alcoholism plus smoking, and lung cavities . Other important risk factors are associated with the supply or quality of the drugs, the possible inadequate drug intake by patients, and others, such as the lack of full supervision during the intensive phase of treatment , poor NTPs [14, 17], and the possible high prevalence of highly virulent MDR-TB strains of MTB [17, 18, 37].
There are several limitations of this survey. First, these is potential misclassification of the new and previously treated cases when some cases registered as new but may actually have had TB treatment in the past. Second, the burden of XDR-TB was underestimated because only resistance to ofloxacin and kanamycin not resistance to capreomycin or other aminoglycosides. Third, classification was based on patient history of prior treatment for TB and reviewing of medical records (which were not available for all patients enrolled). Fourth, this survey did not collect information on HIV infection status because patients with tuberculosis in China are not routinely tested for HIV.
This is the first prevalence survey on resistance to the six major anti-tuberculosis drugs in Lianyungang city. Our study showed that the high prevalence of drug resistance is a major challenge for TB control. Patients with previous treatment history and female gender are risk factors for MDR-TB. Prevention and control of drug-resistant TB should be emphasized by the revised DOTS program through prompt case detection as well as routine and quality-assured DST.
WHO: World Health Organization: Global tuberculosis control. 2012, WHO report, http://wwww.hoint/tb/publications/global_report/gtbr12_mainpdf 2012, website WHO report
Zhao Y, Xu S, Wang L, Chin DP, Wang S, Jiang G, Xia H, Zhou Y, Li Q, Ou X: National survey of drug-resistant tuberculosis in China. N Engl J Med. 2012, 366 (23): 2161-2170.
Iwamoto T, Yoshida S, Suzuki K, Tomita M, Fujiyama R, Tanaka N, Kawakami Y, Ito M: Hypervariable loci that enhance the discriminatory ability of newly proposed 15‒loci and 24‒loci variable‒number tandem repeat typing method on Mycobacterium tuberculosis strains predominated by the Beijing family. FEMS Microbiol Lett. 2007, 270 (1): 67-74.
Glynn J, Vyonycky E, Fine P: Influence of sampling on estimates of clustering and recent transmission of Mycobacterium tuberculosis derived from DNA fingerprinting techniques. Am J Epidemiol. 1999, 149 (4): 366-371.
Organization WH: Multidrug and extensively drug-resistant TB (M/XDR-TB): 2010 global report on surveillance and response. 2010, Geneva: World Health Organization
Chonde T, Basra D, Mfinanga S, Range N, Lwilla F, Shirima R, Van Deun A, Zignol M, Cobelens F, Egwaga S: National anti-tuberculosis drug resistance study in Tanzania. Int J Tuberc Lung Dis. 2010, 14 (8): 967-972.
Yu C-C, Chang C-Y, Liu C-E, Shih L-F, Hsiao J-H, Chen C-H: Drug resistance pattern of <i> mycobacterium tuberculosis</i> complex at a Medical Center in Central Taiwan, 2003–2007. J Microbiol Immunol Infect. 2010, 43 (4): 285-290.
Shao Y, Yang D, Xu W, Lu W, Song H, Dai Y, Shen H, Wang J: Epidemiology of anti-tuberculosis drug resistance in a chinese population: current situation and challenges ahead. BMC Public Health. 2011, 11: 110-
Liu Q, Yang D, Xu W, Wang J, Lv B, Shao Y, Song H, Li G, Dong H, Wan K, et al: Molecular typing of mycobacterium tuberculosis isolates circulating in Jiangsu Province, China. BMC Infect Dis. 2011, 11: 288-
Li X, Wang H, Jing H, Wang Y, Yu C, Wang J, Liu Z, Graviss E, Ma X: Population-based surveillance of extensively drug-resistant tuberculosis in Shandong Province, China [Short communication]. Int J Tuberc Lung Dis. 2012, 16 (5): 612-614.
Hu Y, Mathema B, Wang W, Hoffner S, Kreiswirth B, Xu B: Prevalence of multidrug-resistant pulmonary tuberculosis in counties with different duration of DOTS implementation in rural China. Microb Drug Resist. 2008, 14 (3): 227-232.
Shen X, DeRiemer K, Yuan Z-A, Shen M, Xia Z, Gui X, Wang L, Gao Q, Mei J: Drug-resistant tuberculosis in Shanghai, China, 2000–2006: prevalence, trends and risk factors. Int J Tuberc Lung Dis. 2009, 13 (2): 253-259.
Choi JC, Lim SY, Suh GY, Chung MP, Kim H, Kwon OJ, Lee NY, Park YK, Bai GH, Koh W-J: Drug resistance rates of Mycobacterium tuberculosis at a private referral center in Korea. J Korean Med Sci. 2007, 22 (4): 677-681.
Chacón L, Laínez M, Rosales E, Mercado M, Caminero J: Evolution in the resistance of Mycobacterium tuberculosis to anti-tuberculosis drugs in Nicaragua. Int J Tuberc Lung Dis. 2009, 13 (1): 62-67.
Yoshikawa T: Tuberculosis in aging adults. J Am Geriatr Soc. 1992, 40 (2): 178-
Wright A, Zignol M, Van Deun A, Falzon D, Gerdes SR, Feldman K, Hoffner S, Drobniewski F, Barrera L, van Soolingen D: Epidemiology of antituberculosis drug resistance 2002–07: an updated analysis of the Global Project on Anti-Tuberculosis Drug Resistance Surveillance. Lancet. 2009, 373 (9678): 1861-1873.
Caminero J: Likelihood of generating MDR-TB and XDR-TB under adequate National Tuberculosis Control Programme implementation [Review article]. Int J Tuberc Lung Dis. 2008, 12 (8): 869-877.
Dye C: Doomsday postponed? Preventing and reversing epidemics of drug-resistant tuberculosis. Nat Rev Microbiol. 2009, 7 (1): 81-87.
Kim H-R, Hwang SS, Kim HJ, Lee SM, Yoo C-G, Kim YW, Han SK, Shim Y-S, Yim J-J: Impact of extensive drug resistance on treatment outcomes in non-HIV-infected patients with multidrug-resistant tuberculosis. Clin Infect Dis. 2007, 45 (10): 1290-1295.
Banerjee R, Allen J, Westenhouse J, Oh P, Elms W, Desmond E, Nitta A, Royce S, Flood J: Extensively drug-resistant tuberculosis in California, 1993–2006. Clin Infect Dis. 2008, 47 (4): 450-457.
Jou R, Chiang CY, Huang WL: Distribution of the Beijing family genotypes of Mycobacterium tuberculosis in Taiwan. J Clin Microbiol. 2005, 43 (1): 95-100.
Organization WH: Engaging all health care providers in TB control: guidance on implementing public-private mix approaches. 2006, Geneva: World Health Organization
Yoshiyama T, Yanai H, Rhiengtong D, Palittapongarnpim P, Nampaisan O, Supawitkul S, Uthaivorawit W, Mori T: Development of acquired drug resistance in recurrent tuberculosis patients with various previous treatment outcomes. Int J Tuberc Lung Dis. 2004, 8 (1): 31-38.
Lew W, Pai M, Oxlade O: Initial drug resistance and tuberculosis treatment outcomes. Ann Intern Med. 2008, 149: 123-134.
Mak A, Thomas A, del Granado M, Zaleskis R, Mouzafarova N, Menzies D: Influence of multidrug resistance on tuberculosis treatment outcomes with standardized regimens. Am J Respir Crit Care Med. 2008, 178 (3): 306-312.
Vashakidze L, Salakaia A, Shubladze N, Cynamon M, Barbakadze K, Kikvidze M, Papitashvili L, Nonikashvili M, Solomonia N, Bejanishvili N: Prevalence and risk factors for drug resistance among hospitalized TB patients in Georgia. Int J Tuberc Lung Dis. 2009, 13 (9): 1148-
Mdivani N, Zangaladze E, Volkova N, Kourbatova E, Jibuti T, Shubladze N, Kutateladze T, Khechinashvili G, del Rio C, Salakaia A: High prevalence of multidrug-resistant tuberculosis in Georgia. Int J Infect Dis. 2008, 12 (6): 635-
Toungoussova O, Caugant D, Sandven P, Mariandyshev A, Bjune G: Drug resistance of Mycobacterium tuberculosis strains isolated from patients with pulmonary tuberculosis in Archangels, Russia. Int J Tuberc Lung Dis. 2002, 6 (5): 406-414.
Lockman S, Kruuner A, Binkin NJ, Levina K, Wang YC, Danilovitsh M, Hoffner SE, Tappero JW: Clinical outcomes of Estonian patients with primary multidrug-resistant versus drug-susceptible tuberculosis. Clin Infect Dis. 2001, 32 (3): 373-380.
Chiang C, Chang C, Chang R, Li C, Huang R: Patient and health system delays in the diagnosis and treatment of tuberculosis in Southern Taiwan. Int J Tuberc Lung Dis. 2005, 9 (9): 1006-1012.
Long NH, Johansson E, Lonnroth K, Eriksson B, Winkvist A, Diwan V: Longer delays in tuberculosis diagnosis among women in Vietnam. Int J Tuberc Lung Dis. 1999, 3 (5): 388-393.
Jeon CY, Hwang SH, Min JH, Prevots DR, Goldfeder LC, Lee H, Eum SY, Jeon DS, Kang HS, Kim JH: Extensively drug‒resistant tuberculosis in South Korea: risk factors and treatment outcomes among patients at a tertiary referral hospital. Clin Infect Dis. 2008, 46 (1): 42-49.
Faustini A, Hall AJ, Perucci CA: Risk factors for multidrug resistant tuberculosis in Europe: a systematic review. Thorax. 2006, 61 (2): 158-163.
Law W, Yew W, Chiu Leung C, Kam K, Tam C, Chan C, Leung C: Risk factors for multidrug-resistant tuberculosis in Hong Kong. Int J Tuberc Lung Dis. 2008, 12 (9): 1065-1070.
Barroso EC, Mota RMS, Santos RO, Sousa ALO, Barroso JB, Rodrigues JLN: Risk factors for acquired multidrug-resistant tuberculosis. J Pneumol. 2003, 29 (2): 89-97.
Gandhi NR, Moll A, Sturm AW, Pawinski R, Govender T, Lalloo U, Zeller K, Andrews J, Friedland G: Extensively drug-resistant tuberculosis as a cause of death in patients co-infected with tuberculosis and HIV in a rural area of South Africa. Lancet. 2006, 368 (9547): 1575-1580.
Amor YB, Nemser B, Singh A, Sankin A, Schluger N: Underreported threat of multidrug-resistant tuberculosis in Africa. Emerg Infect Dis. 2008, 14 (9): 1345-
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/13/1171/prepub
This study was supported by the National Nature Science Foundation of China (81302480), Key personnel of Science and Education Industrial Engineering (JKRC2011005), Provincial Department of Public Health of Jiangsu (Y201030), Key Medical Disciplines of Jiangsu Province (XK201150) and Ministry of Health (W201208). The sponsors have no role in the study design or data analysis.
We also thank Christopher Whalen, M.D., M.S. and Leonardo Martinez, MPH of Department of Epidemiology and Biostatistics, Miller Hall College of Public Health University of Georgia for their careful reviewing of this manuscript.
The authors declare that they have no competing interests.
QL, LZ, WL conceived the study, analyzed the data and drafted the manuscript; YS, HS, GL, JS participated in the study design, implemented the field investigation and performed DST tests; YZ, CZ and CC participated in the study design and helped draft the manuscript. All authors contributed to the study and have read and approved the final manuscript.
Qiao Liu, Limei Zhu contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.