Skip to main content

Primary drug resistance of mycobacterium tuberculosis in Shandong, China, 2004–2018



Primary drug-resistant tuberculosis (DR-TB) has contributed to a significant health and economic burden on a global scale, especially in China. we sought to estimate epidemiological characteristics of primary DR-TB in China from 2004 to 2018.


Eleven thousand four hundred sixty-seven newly diagnosed and 1981 retreated TB cases with drug susceptibility data were included. Chi-Square test for trends, linear regression, a joinpoint regression model and temporal trend in proportions of the different resistance patterns were carried out.


The proportion of primary DR-TB and mono-resistant TB (MR-TB) in China had reduced by more than 12% since 2004, and were 21.38%, 13.35% in 2018 respectively. Among primary DR-TB cases (2173,18.95%), the percentage of multiresistant TB (MDR-TB, from 5.41 to 17.46%), male (from 77.03 to 84.13%), cavity (from 13.51 to 43.92%), rifampicin(RFP)-resistant TB (from 8.11 to 26.98%), streptomycin(SM)-resistant TB (from 50.00 to 71.43%) increased significantly (P < 0.05). On the contrary, the proportion of female, non-cavity, isoniazide(INH)-resistant TB (from 55.41 to 48.15%) and MR-TB (from 82.43 to 62.43%) decreased significant (P < 0.05). The primary drug resistance rate among female, cavity, smoking, drinking, 15 to 44 year-old TB subgroups increased by 0.16, 6.24, 20.95, 158.85, 31.49%, respectively. The percentage of primary DR-TB, RFP-resistant TB dropped significantly during 2004–2007 in Joinpoint regression model.


The total rate of drug resistance among new TB cases showed a downward trend in Shandong, China, from 2004 to 2018. Primary drug resistance patterns were shifting from female, non-cavity, INH-resistant TB, and MR-TB groups to male, cavity, RFP/SM-resistant TB, and MDR-TB groups. Considering the rising drug resistance rate among some special population, future control of primary DR-TB in China may require an increased focus on female, cavity, smoking, drinking, or 15 to 44 year-old TB subgroups.


Drug resistance of mycobacterium tuberculosis (MTB), especially multidrug resistant tuberculosis (MDR-TB, defined as resistance to at least rifampin and isoniazid) has contributed to a significant health and economic burden on a global scale [1,2,3]. Drug-resistant tuberculosis (DR-TB) can be transmitted (primary resistance, refers to the infection with drug-resistant MTB) or develop during the course of treatment (secondary or acquired resistance) [3]. Plenty of researches [4,5,6] had revealed that the majority of DR-TB were primary instead of acquired drug resistance, in other words, the main mechanism of drug resistance in TB was the transmission of drug-resistant MTB strains from the existing TB patients rather than mismanagement of previous treatment episode. For instance, a molecular epidemiological study on DR-TB strains in Shanghai using highly discriminatory whole genome sequencing (WGS) found that more than 73% of MDR-TB cases were caused by primary transmission [5, 6].

According to the World Health Organization (WHO) Global TB report 2018, an estimated 10 million people developed TB in 2017, of whom 3.5% of newly diagnosed and 18% of retreated TB cases were MDR-TB, and China had the world’s second largest number of TB cases and MDR/rifampicin-resistant TB(RR-TB), comprising around 9 and 13% of the world total, respectively, behind only India (TB cases, 27%; MDR/RR-TB, 24%) [7]. Based on a national survey of drug-resistant tuberculosis (DR-TB) in China, there were 869,092 new TB cases and 54,216 previously treated TB cases in 2010, and the rate of primary and acquired MDR-TB were 5.7% (95% confidence interval [CI], 4.5–7.0) and 25.6% (95% CI, 21.5–29.8), respectively [4]. Great progress in TB control and prevention had been made in China since the implement of DOTS strategy in 1990s, and prevalence rate of TB reduced 3.4% per year [8]. However, previous studies reported that the rate of total DR-TB remained a high level, and the overall MDR rate was 6.2% in Shandong during 2007–2014, overall rifampin (RFP) resistance and rifampin monoresistance (RMR) increased at a yearly rate of 0.2 and 0.1%, respectively [9].

Understanding the burden of primary DR-TB and the factors associated with its transmission may help determine the high-risk population for drug-resistant MTB infection and develop measures for preventing transmission. Moreover, investigating the current epidemiological status of primary DR-TB in China also help to assess the effectiveness of existing TB interventions and provide guidance for TB control in future. However, few studies [4, 10] had focused on the epidemic of primary DR-TB in China. To evaluate the trend, high risk factors and epidemiological characteristics of primary DR-TB in China from 2004 to 2018, we collected 11,467 new TB cases with drug susceptibility test (DST) results for statistical analysis by applying a join-point regression model and chi-square regression for trend, and DR-TB groups were stratified by drug-resistant profiles, age, sex, smoking, drinking and cavity for further analysis.

Materials and methods

Ethics statement

The protocols applied in this study were approved by the Ethics Committee of Shandong Provincial Hospital, affiliated with Shandong University (SPH) and and the Ethic Committee of Shandong Provincial Chest Hospital (SPCH), China. Before analysis, patient records were anonymized and deidentified.


This study was carried out in Shandong, a coastal province of the East China region. Shandong consisted of 17 municipalities and 140 counties (districts) with 95.8 million inhabitants at the 2010 Census. It is located at 36°24′N latitude 118°24′E longitude with an area of 157,100 km2. Shandong has emerged as one of the most populous and most affluent provinces in China since the late nineteenth century.

Study population and data collection

A total of 11,467 newly diagnosed and 1981 retreated MTB cases were collected from 36 TB prevention and control institutions of Shandong Province, China, Jan 1, 2004 to Dec 31, 2018. Two province-level hospitals (SPH and SPCH), 13 municipal-level and 21 county-level local health departments were involved in the surveillance of DR-TB from 2004 to 2018. All MTB cases enrolled in our study were consecutive culture-confirmed and finished DST for first line anti-TB drugs. Demographic and clinical characteristics on age, sex, drinking, smoking, cavity, treatment history and extra-pulmonary TB were available. Smoker (or drinker) refers to those who satisfy at least one of the following two conditions: i) smoking (or drinking) for 6 months or above; ii) those who was still smoking(or drinking) or had stopped smoking (or drinking) for less than 6 months before TB diagnosis. Non-smoker and non-drinker were defined as the person who had never smoked or never drunk, respectively.

Bacteriologic examinations and drug susceptibility testing

Each surveillance site collected two sputum samples from all eligible patients, and then sent all sputum samples to the TB Reference Laboratory of SPCH for further examination including bacteriologic culture, DST, and species identification. Isolates were inoculated into tubes of acidified Löwenstein-Jensen (LJ) medium after conventional pretreatment process and cultured at 37 °C [4]. Cultures with growing colonies were sent for further identification and DST. According to previous published protocol [11], standard traditional biochemical testings such as P-nitrobenzoic acid, 2-thiophene carboxylic acid hydrazide testing and 16S rRNA gene sequence analysis (MicroSeq ID Microbial Indentification Software(version 2.0); Applied Biosystems, Foster City, CA, USA) [12] were used to differentiate M. tuberculosis from other Mycobacteria spp.. DST of M. tuberculosis were performed using absolute concentration method on L-J media, and all procedures were carried out in accordance with the guideline of WHO [13]. The concentration of four first-line anti-TB drugs were as follows: 0.2 μg/mL (isoniazid, INH), 40 μg/mL(rifampin, RFP), 10 μg/mL (streptomycin, SM), 2 μg/mL (ethambutol, EMB). DST of ethionamide, fluoroquinolone, kanamycin and pyrazinamide were not routinely performed.

Quality control

Two professionally trained investigators were independently responsible for quality assessment and data extraction; and all laboratories involved in our study regularly accepted external quality assessment of Superior TB National Reference laboratory in SPCH.


Drug-resistant tuberculosis (DR-TB) are classified as having acquired or primary drug resistance on the basis of a history of previous treatment [13].

Mono-resistance (MR) refers to resistance to one first-line anti-TB drug only [14].

Polydrug resistance (PDR) refers to resistance to more than one first-line anti-TB drug, other than both isoniazid and rifampicin [14].

Multidrug resistance (MDR) refers to resistance to at least both isoniazid and rifampicin [14].

Statistical analysis

Categorical variables including age (0–14, 15–44, 45–64, 65+), sex(male or female), drinking (yes/no/unknown), smoking (yes/no/unknown), cavity (yes/no/unknown), patients type (extra-pulmonary TB/pulmonary TB) of new and relapse TB cases were calculated as counts and proportions, respectively. In addition, odds ratios (ORs) and 95% CIs for the comparisons of these characteristics were estimated between newly diagnosed susceptible and DR-TB cases, retreated susceptible and DR-TB cases, new and relapse TB cases, primary and acquired DR-TB cases. Differences in drug susceptibility profiles between primary DR M. tuberculosis and acquired drug-resistant M. tuberculosis isolates were analyzed using Pearson Chi-square test, and p-value < 0.05 was considered statistically significant. Chi-square test for trends and linear regression in line charts and stacked bar charts were used to assess the changes and temporal trend in quantity and proportions of the different resistance patterns among total DR-TB cases from 2004 to 2018.

We calculated total or annual primary drug resistance rate of tuberculosis as the number of annual primary DR-TB cases divided by the population of newly diagnosed TB cases. The overall rate and the annual rate stratified by age groups, sex, drinking (yes/no/unknown), smoking (yes/no/unknown), cavity (yes/no/unknown) and type (extra-pulmonary TB/pulmonary TB) were estimated as well. Indicators including annual percent changes (APCs) and average annual percent changes (AAPCs) were calculated from 2004 to 2018 inclusive on the basis of a joinpoint regression model. Each segment described a short-term trend (APC). Long-term trends over the entire study period are AAPCs and were estimated as the weighted average of the short-term APCs, with the weights equal to the length of the short-term line segment [15]. Two-sided t-tests at p < 0.05 were used to test whether AAPCs and APCs were statistically significantly different from zero; 95% CI for each segment were calculated. If the AAPC was within one segment, the t-distribution was used. Otherwise, the normal (z) distribution was used [15,16,17]. A non-significant (p ≥ 0.05) APC was described as stable while a significant (p < 0.05) positive or negative APC was termed as increase or decrease. All analyses were carried out using SPSS software (version 20.0) and the Joinpoint Regression Software (version


Patients’ characteristics

Table 1 illustrates the baseline characteristics of the study participants. A total of 11,467 new and 1981 relapse TB cases were reported during the period from Jan 1, 2004 to Dec 31, 2018 in Shandong, China, of which 2173 (18.95%) and 505 (23.24%) were drug-resistant TB cases, respectively, P < 0.01. Among these primary DR-TB cases, there were more males (83.06% vs 16.94%), non-smokers (40.73% vs 13.90%), non-drinkers (42.84% vs 11.50%), pulmonary TB (99.68% vs 0.32%). In addition, most primary DR-TB cases (2168, 99.77%) were from three age groups: 15–44 years (930, 42.80%), 45–64 years (728, 33.50%), > 65 years (486, 22.37%). The distribution of age, sex, smoking, et al. among acquired DR-TB cases were similar to primary DR-TB cases (Table 1).

Table 1 Demographic and clinical characteristics of new and relapse TB patients, Shandong, China, 2004–2018a

Compared with newly diagnosed susceptible TB cases, primary DR-TB cases were more likely to be female (OR: 1.166, 95%CI:1.031–1.320) and less likely to be aged more than 65(OR: 0.840, 95%CI:0.745–0.948). Among retreated TB cases, factors associated with acquired resistance were baseline cavitary diseases (OR:1.544, 95%CI:1.24–1.924) and > 65 years (OR:0.684, 95%CI:0.521–0.900). New TB were less likely to be aged 45–64 years old (OR:0.84, 95%CI: 0.750–0.942) or > 65 years old (OR:0.851, 95%CI: 0.753–0.962) when compared with relapse TB cases. Moreover, patients with pulmonary TB had a statistically significant higher risk of primary drug resistance (vs acquired drug resistance) than those with extra-pulmonary TB (OR:4.349, 95%CI:1.519–12.456) (Table 1).

Resistance patterns

Among 11,467 newly diagnosed and 1981 retreated TB clinical isolates, the number and proportion of any resistance to first-line drugs including INH, RFP, EMB, SM were 1232(10.74%) and 333(16.81%), 508(4.43%) and 223(11.26%), 174(1.52%) and 68(3.43%), 1482(12.92%) and 329(16.61%), respectively; new TB cases had a lower rate of DR-TB(18.95% vs 23.24%, P < 0.01), MDR-TB(3.19% vs 8.43%, P < 0.01), PDR-TB(4.12% vs 5.55%, P < 0.01) than relapse TB cases, but almost the same in MR-TB (11.56% vs 11.36%, P > 0.05). Of 1326 primary MR-TB cases, SM was associated with the highest rate of resistance (6.41%,758), followed by INH (3.88%, 445), RFP (0.78%, 90), and EMB (0.17%, 20). There were four main types of primary MDR-TB, for instance, MDR1 (INH + RFP), MDR2 (INH + RFP + EMB), MDR3 (INH + RFP + EMB + SM), and MDR4 (INH + RFP + SM), which accounted for 0.71%(81), 0.13%(15), 0.72%(83), 1.45%(166) respectively. Similarly, primary PDR-TB mainly consisted of PDR1 (INH + EMB, 0.12%, 14), PDR2 (INH + SM, 3.41%, 391), PDR3 (RFP + EMB, 0.06, 7), PDR4 (RFP + SM, 0.33%, 38), PDR5 (INH + EMB + SM, 0.14%, 16) (Table 2).

Table 2 Primary and acquired drug resistance profiles of Mycobacterium tuberculosis, Shandong, China, 2004–2018a

Disparity in the proportion of primary DR-TB subgroups

As shown in Figs. 1a-f and 2, primary DR-TB cases was sub-divided into different subgroups according to the types of drug resistance (MR-TB/MDR-TB/PDR-TB, INH/RFP/SM/EMB-resistance), age (0–14, 15–44, 45–64, 65+), sex (male or female), drinking history (yes/no/unknown), smoking history (yes/no/unknown), cavity (yes/no/unknown), and the amount and proportions of each subgroup among the total primary DR-TB cases varied every year since 2004 to 2018. Among all primary DR-TB cases, the percentage of MDR-TB ([R2 = 0.2034] from 5.41 to 17.46%; χ2 test for trends: χ2 = 5.376, P = 0.020), males ([R2 = 0.4834] from 77.03 to 84.13%; χ2 test for trends: χ2 = 12.570, P < 0.001), cavity ([R2 = 0.7022] from 13.51 to 43.92%; χ2 test for trends: χ2 = 120.53, P < 0.001), RFP-resistance ([R2 = 0.5355] from 8.11 to 26.98%; χ2 test for trends: χ2 = 16.785, P < 0.001), SM-resistance ([R2 = 0.5365] from 50.00 to 71.43%; χ2 test for trends: χ2 = 22.076, P < 0.001) increased significantly from 2004 to 2018. On the contrary, the proportion of females, non-cavity, INH ([R2 = 0.4269] from 55.41 to 48.15%; χ2 test for trends: χ2 = 14.725, P < 0.001) and MR-TB ([R2 = 0.178] from 82.43 to 62.43%; χ2 test for trends: χ2 = 6.287, P = 0.012) decreased significantly (P < 0.05) (Table 5 in Appendix).

Fig. 1

Trends for the quantity and proportions of different subgroups among total primary DR-TB cases, Shandong, China, 2004–2018*. a Trends for MR-TB, MDR-TB, PDR-TB among total primary cases of TB; b Trends for primary DR-TB cases of different age (0–14, 15–44, 45–64, 65+); c Trends for primary DR-TB cases of different sex(male or female); d Trends for primary DR-TB cases with or without cavity; e Trends for primary DR-TB cases with or without smoking history; f Trends for primary DR-TB cases with or without drinking history; The proportions of each subgroups were calculated as follows: (the quantity of each subgroups/ the quantity of total primary DR-TB subgroups in the same year)*100%. The χ2 and linear regression results are shown in Table 5 in Appendix. TB, tuberculosis; DR-TB, drug-resistant tuberculosis; MR-TB, mono-resistant tuberculosis; MDR-TB, multi-resistant tuberculosis; PDR-TB, polydrug resistant tuberculosis. EMB, ethambutol; INH, isoniazid; RFP, rifampin; SM, streptomycin

Fig. 2

Overall first-line drug resistance for INH, RIF, EMB, and SM in primary cases of tuberculosis in Shandong China, 2004–2018. The proportions of INH-, RIF-, EMB-, and SM-resistance were calculated as follows: (the quantity of each subgroups/ the quantity of total primary DR-TB subgroups in the same year)*100%. The χ2 and linear regression results are shown in Table 5 in Appendix. EMB, ethambutol; INH, isoniazid; RFP, rifampin; SM, streptomycin

Total and annul primary drug resistance rate

Table 3 illustrates the overall and annul primary drug resistance rate of various TB subgroups from 2004 to 2018. During the study period, the primary drug resistance rate among female, cavity, smoking, drinking, 15 to 44 year-old TB subgroups increased by 0.16, 6.24, 20.95, 158.85, 31.49%, and varied from 20.24 to 20.27%, from 20.83 to 22.13%, from 21.57 to 26.09%, from 10.26 to 26.55%, 18.59 to 24.44%, respectively; meanwhile the rate of DR-TB among male, non-cavity, non-smoking, non-drinking, 45 to 64 year-old, > 65 year-old TB subgroups reduced by 16.62, 26.34, 15.79, 22.18, 35.01, 36.66% respectively, decreasing from 20.24 to 20.27%, 24.90 to 18.34%, 21.57 to 26.09%, 25.00 to 19.45%, 32.47 to 21.10%, and 28.17 to 17.84%, respectively. In addition, the rate of MDR-TB, PDR-TB, RFP-resistant TB, and SM-resistant TB among new TB cases increased from 1.32 to 3.73%, 2.96 to 4.30%, 1.97 to 5.77%, and 12.17 to 15.27%, with a percentage change of 183.71, 45.20, 192.31, 25.47%, respectively; the rates of MR-TB, INH-resistant TB, EMB-resistant TB, and total DR-TB dropped by 33.48, 23.67, 25.49, 12.17%, respectively, with a decrease from 20.07 to 13.35%, 13.49 to 10.29%, 1.97 to 1.47%, respectively.

Table 3 The temporal change trend of primary drug resistance rate among newly diagnosed TB cases in Shandong, China, 2004–2018

Temporal trends and Joinpoint regression model

The Joinpoint regression analysis (Table 4) revealed that total primary MR-TB rate had decreased during 2004–2007 (APC = -20.6% [95%CI: − 33.7, − 4.9], Z = -3.00, P < 0.01). In contrast, total primary MDR-TB rate had increased during 2004–2006 (APC = 57.9% [95%CI: 0.3, 148.7], Z = 2.40, P < 0.01). The rate of INH resistant TB showed a significant decline in the trend with an APC of − 2.3% during 2004–2018(P < 0.01), while the rate of RFP resistant TB were rising significantly during 2004–2010 (APC = 16.90% [95%CI: 5.8, 29.1], Z = 3.50, P < 0.01).

Table 4 Annual percentage change in primary drug-resistant rate of Mycobacterium tuberculosis in Shandong, China, 2004–2018

The primary DR-TB rate among smoking and drinking subgroups had been on the rise during 2011–2018 (APC = 7.1% [95%CI: 0.8, 13.8], Z = 2.90, P < 0.01) and 2008–2018 (APC = 6.7% [95%CI: 1.6,12.1], Z = 3.00, P < 0.01), respectively. Of 45–64 year-old TB group, the primary DR-TB rate decline during 2004–2007(APC = -20.5% [95%CI: − 33.2,-5.3], Z = -3.10, P < 0.01). DR-TB among elderly TB cases (> 65 years) showed a sharply drop during 2004–2009 (APC = -12.7% [95%CI: − 21.8,-2.4], Z = -2.70, P < 0.01).


Our study had enrolled 11,467 newly diagnosed TB cases with DST results at 36 TB surveillance centers across Shandong, China from 2004 to 2018 to evaluate the epidemiology and high risk factors of primary TB drug resistance, and we found that ongoing transmission of drug resistant MTB remained to be the major mechanism of DR-TB in China, approximately one fifth of new TB cases were primary DR-TB while 1.3–4.4% were MDR-TB, slightly differed from previous researches in Zhejiang, Jiangsu, Heilongjiang, China [10, 18,19,20]. Trends in annual prevalence of primary DR revealed that the rates of MDR-TB (from 1.32 to 3.73%), PDR-TB (from 2.96 to 4.30%) began to increase slowly along with that of MR-TB (from 20.07 to 13.35%) decreasing, However, MR-TB still dominated the drug resistance rate. In a similar way the proportion of these three sub-types of DR-TB changed during our study. The increment in both MDR-TB and PDR-TB alarmed us that it would continue to be a tough challenge to TB elimination strategy in China [4, 21].

According the temporal trend in proportions of different primary DR-TB subgroups each year, we pointed out that primary resistance patterns were shifting from female, non-cavity, INH resistant TB, and MR-TB groups to male, cavity, RFP/SM resistant TB, and MDR-TB groups in Shandong China from 2004 to 2018, and more than two thirds were male or had a resistance to SM, but more than half of the present primary DR-TB cases were still MR-TB. The percentage of males among MDR-TB cases were more than 70%, which was consistent with previous researches in Bangladesh [22], Mozambique [23], et al., but in contrast with Ethiopia [24] and Lianyungang city, China [25] where females were in majority. The major reasons for this phenomenon were that with genuine gender or behavioral differences [25], men had a significant higher risk of contracting and dying from TB than women [26], and it was reported by WHO that TB incidence around China in 2017 for male patients (67.5%) was twice as much as females (32.5%) [7]. In accordance with Li D [27] and Yeom’s [28] results, we found that cavities had turned to be more common for primary DR TB since 2012 in China. Low effectiveness of anti-TB drug sensitivity towards primary DR-TB, limited drug penetration into cavities, specific microorganism virulence and immune status of patients may contribute to this phenomenon [27,28,29]. However, another reason the more important may be the imbalance between compensatory evolution and compensation cost of drug resistant MTB strains [30, 31]. Four out of ten primary DR-TB cases were in the age group of 15–44 years, followed by groups of 45–64 years (33.5%) and > 65 years (22.37%), suggesting that DR-TB prevention and control should focus more on young and middle aged patients. In recent years, TB patients who were alcohol users and smokers were associated with a higher annual primary DR-TB prevalence than those who were non-smokers and non-drinkers. However, when total TB cases across 2004–2018 were involved in bivariate analysis, personal behaviours like smoking and alcohol use were no longer predictors of primary DR-TB on bivariate analysis, in contrast to findings in India [32] and Botswana [33] where alcohol use was a risk factor for MDR-TB, but consistent with studies in Nepal [34], indicating that it may be distinguished by regions and phases.

Directly observed treatment, short-course (DOTS) strategy had been put forward in China since 1990s, and was expanded to cover all smear-positive and smear-negative TB patients since 2005 instead of previous policy that only patients diagnosed with smear-positive or severe smear-negative TB were available to free treatment [35, 36]. As a consequence, TB prevalence has declined sharply from 1990 to 2010 [36, 37], in accordance with which both the percentage of primary DR-TB, RFP resistant TB and the drug resistant rate of 45–64 year-old, > 65 year-old or female group dropped significantly during 2004–2007. Possible reason for this trend was the reduction of acquired DR-TB due to the free treatment policy in China (35–37), in other words, the source of primary DR-TB had been controlled effectively. An obvious increment in primary MDR-TB and RFP resistant TB during 2004–2008 had been observed, and one hypothesis was that MDR/RR-TB were more difficult to treat than drug-susceptible TB and other types of DR-TB, thus it was hard to prevent the newly transmission of previous MDR/RR-TB as well [4, 7, 38]. Furthermore, evidence has shown that substantial costs associated with TB diagnosis and treatment remained a heavy financial burden for TB patients in spite of the “free” TB care policy [35]. Other mechanisms remained to be further explored. Another issue of concern was that people with smoking or drinking habits tended to be more susceptible to infection of drug-resistant MTB [39, 40], consistent to our research since 2008. In comparison with never-smokers, current smokers had an excess risk of pulmonary TB (adjusted HR, 2.87; 95% CI, 2.00–4.11; p < 0.001) [40].

Our study had several advantages. First, our study covered all DST data of Shandong province from 2004 to 2018, which is the second largest province in China with a population of nearly 90 million, thus the findings of our research were more likely to be rolled out nationally. Second, we distinguished primary DR-TB from acquired DR-TB, and the total TB or DR-TB population were stratified by various factors including sex, age, smoking history, drinking history, cavity and so on, other than most previous epidemiological researches in DR-TB that considered new and relapse DR-TB together [1, 2, 9, 10]. Third, we proposed a Joinpoint regression model to describe the temporal trend and turning point of TB control in Shandong, China.

There were also some limitations of our study. First, drug sensitive tests were not routinely carried out among TB cases, thus our data from TB monitoring stations were affected by screening intensity and local medical conditions, and could be overestimated due to selection bias. Second, diversities in technical levels and experimental conditions in different TB monitoring stations may contributed to unavoidable bias. Third, DST for second-line anti-TB drugs were seldom performed by TB monitoring stations in Shandong, China, thus the epidemiology of primary resistance to second-line anti-TB drugs and extensive drug-resistant TB (XDR-TB) could not be estimated.


Our study had described the epidemiological features and temporal trend of primary DR-TB in Shandong province, China from 2004 to 2018 to evaluate the current situation of TB prevention and control. We found that the proportion of primary DR-TB and MR-TB had reduced by more than 12% since 2004, and were 21.38, 13.35% in 2018 respectively. Moreover, primary drug resistance patterns were shifting from female, non-cavity, INH resistant TB, and MR-TB groups to male, cavity, RFP/SM resistant TB, and MDR-TB groups, even more than two thirds were male or had a resistance to SM, but more than half of the present primary DR-TB cases were still MR-TB. This study indicates the DOTS strategy and “free” policy in China had achieved some effects in TB control, but considering the increment of drug resistance rate among some special population, more attention should be focused on female, cavity, smoking, drinking, 15 to 44 year-old TB subgroups. In addition, more project support, high-quality training for medical staff and enhanced public awareness of TB prevention and control are also necessary for TB elimination goal in China.

Availability of data and materials

Data can be available through contact with the corresponding author.



Average annual percent changes


Annual percent changes


Confidence interval


Drug-resistant tuberculosis


Drug susceptibility test






Multi-resistant tuberculosis


Mono-resistant tuberculosis


Mycobacterium tuberculosis


Odds ratios


Polydrug resistant tuberculosis




Rifampin-resistant tuberculosis






Whole genome sequencing


World health organization


Extensive drug resistant TB


  1. 1.

    Gandhi NR, Nunn P, Dheda K, Schaaf HS, Zignol M, van Soolingen D, et al. Multidrug-resistant and extensively drug-resistant tuberculosis: a threat to global control of tuberculosis. Lancet. 2010;375(9728):1830–43.

    Article  Google Scholar 

  2. 2.

    Yang Y, Zhou C, Shi L, Meng H, Yan H. Prevalence and characterization of drug-resistant tuberculosis in a local hospital of Northeast China. Int J Infect Dis. 2014;22:83–6.

    Article  Google Scholar 

  3. 3.

    World Health Organization. Antimicrobial Resistance. Geneva: World Health Organization; 2018.

    Google Scholar 

  4. 4.

    Zhao Y, Zhao B, Xu S, Wang S, Wang L, Wang Y, et al. National survey of drug-resistant tuberculosis in China. N Engl J Med. 2012;366(23):2161–70.

    CAS  Article  Google Scholar 

  5. 5.

    Yang C, Luo T, Shen X, Wu J, Gan M, Xu P, et al. Transmission of multidrug-resistant Mycobacterium tuberculosis in Shanghai, China: a retrospective observational study using whole-genome sequencing and epidemiological investigation. Lancet Infect Dis. 2017;17(3):275–84.

    CAS  Article  Google Scholar 

  6. 6.

    Marais BJ, Sintchenko V. Epidemic spread of multidrug-resistant tuberculosis in China. Lancet Infect Dis. 2017;17(3):238–9.

    Article  Google Scholar 

  7. 7.

    World Health Organization. Global tuberculosis report 2018. Geneva: World Health Organization; 2018.

    Google Scholar 

  8. 8.

    WHO 2019. Tuberculosis in China.

  9. 9.

    He X, Zhang X, Zhao J, Liu Y, Yu C, Yang G, et al. Epidemiological trends of drug-resistant tuberculosis in China from 2007 to 2014 a retrospective study. Medicine. 2016;95(15):e3336.

    Article  Google Scholar 

  10. 10.

    Yang C, Gao Q. Recent transmission of mycobacterium tuberculosis in China: the implication of molecular epidemiology for tuberculosis control. Front Med. 2018;12(1):76–83.

    Article  Google Scholar 

  11. 11.

    El Amin NM, Hanson HS, Pettersson B, Petrini B, Von Stedingk LV. Identification of non-tuberculous mycobacteria: 16S rRNA gene sequenceanalysis vs. conventional methods. Scand J Infect Dis. 2000;32:47–50.

    Article  Google Scholar 

  12. 12.

    Jing H, Wang H, Wang Y, Deng Y, Li X, Liu Z, et al. Prevalence of nontuberculous mycobacteria infection, China, 2004-2009. Emerging Infect Dis. 2012;18(3):527–8.

    Article  Google Scholar 

  13. 13.

    Van Rie A, Warren R, Richardson M, Gie RP, Enarson DA, Beyers N, et al. Classification of drug-resistant tuberculosis in an epidemic area. Lancet. 2000;356(9223):22–5.

    Article  Google Scholar 

  14. 14.

    Falzon D, Schünemann HJ, Harausz E, González-Angulo L, Lienhardt C, Jaramillo E, et al. World Health Organization treatment guidelines for drug-resistant tuberculosis, 2016 update. Eur Respir J. 2017;49(3):1602308.

    Article  Google Scholar 

  15. 15.

    Nistal-Nuño B. Joinpoint regression analysis to evaluate traffic public health policies by national temporal trends from 2000 to 2015. Int J Inj Control Saf Promot. 2018;25(2):128–33.

    Article  Google Scholar 

  16. 16.

    Bandi P, Silver D, Mijanovich T, Macinko J. Temporal trends in motor vehicle fatalities in the United States, 1968 to 2010 - a joinpoint regression analysis. Inj Epidemiol. 2015;2(1):1–11.

    Article  Google Scholar 

  17. 17.

    Brito AL, Monteiro LD, Ramos Junior AN, Heukelbach J, Alencar CH. Temporal trends of leprosy in a Brazilian state capital in Northeast Brazil: epidemiology and analysis by joinpoints, 2001 to 2012. Braz J Epidemiol. 2016;19(1):194.

    Google Scholar 

  18. 18.

    Wang X, Fu Q, Li Z, Chen S, Liu Z, Nelson H, et al. Drug-resistant tuberculosis in Zhejiang province, China, 1999-2008. Emerging Infect Dis. 2012;18(3):496–8.

    Article  Google Scholar 

  19. 19.

    Li D, Wang J, Ji B, Cui J, Pan X, Fan C, et al. Persistently high prevalence of primary resistance and multidrug resistance of tuberculosis in Heilongjiang Province, China. BMC Infect Dis. 2016;16(1):516.

    Article  Google Scholar 

  20. 20.

    Li X, Lu W, Zu R, Zhu L, Yang H, Chen C, et al. Comparing risk factors for primary multidrug-resistant tuberculosis and primary drug-susceptible tuberculosis in Jiangsu Province, China: a matched-pairs case-control study. Am J Trop Med Hyg. 2015;92(2):280–5.

    Article  Google Scholar 

  21. 21.

    Raviglione M, Sulis G. Tuberculosis 2015: burden, challenges and strategy for control and elimination. Infect Dis Rep. 2016;8(2):33–7.

    Article  Google Scholar 

  22. 22.

    Karmaker H, Basar MA, Karim MR, et al. An epidemiological study of drug resistant tuberculosis cases: survey in the northern part of Bangladesh. Public Health Res. 2016;6(2):52–8.

    Google Scholar 

  23. 23.

    Samo Gudo P, Cuna Z, Coelho E, Maungate S, Borroni E, Miotto P, et al. Is multidrug-resistant tuberculosis on the rise in Mozambique? Results of a national drug resistance survey. Eur Respir J. 2011;38(1):222–4.

    CAS  Article  Google Scholar 

  24. 24.

    Mulisa G, Workneh T, Hordofa N, Suaudi M, Abebe G, Jarso G. Multidrug-resistant mycobacterium tuberculosis and associated risk factors in Oromia region of Ethiopia. Int J Infect Dis. 2015;39(C):57–61.

    Article  Google Scholar 

  25. 25.

    Liu Q, Zhu L, Shao Y, Song H, Li G, Zhou Y, et al. Rates and risk factors for drug resistance tuberculosis in Northeastern China. BMC Public Health. 2013;13(1):1171.

    Article  Google Scholar 

  26. 26.

    Global tuberculosis report 2017. Geneva: World Health Organization; 2017.

  27. 27.

    Li D, He W, Chen B, Lv P. Primary multidrug-resistant tuberculosis versus drug-sensitive tuberculosis in non-HIV-infected patients: comparisons of CT findings. PLoS One. 2017;12:e0176354.

    Article  Google Scholar 

  28. 28.

    Yeom JA, Jeong YJ, Jeon D, Kim KI, Kim CW, Park HK, et al. Imaging findings of primary multidrug-resistant tuberculosis: a comparison with findings of drug-sensitive tuberculosis. J Comput Assist Tomogr. 2009;33(6):956–60.

    Article  Google Scholar 

  29. 29.

    Fishman JE, Sais GJ, Schwartz DS, Otten J. Radiographic findings and patterns in multidrug-resistant tuberculosis. J Thorac Imaging. 1998;13(1):65–71.

    CAS  Article  Google Scholar 

  30. 30.

    Brandis G, Wrande M, Liljas L, Hughes D, Medicinska fakulteten, Medicinska och farmaceutiska vetenskapsområdet, et al. Fitness-compensatory mutations in rifampicin-resistant RNA polymerase. Mol Microbiol. 2012;85(1):142–51.

    CAS  Article  Google Scholar 

  31. 31.

    Cohen KA, Abeel T, Manson McGuire A, Desjardins CA, Munsamy V, Shea TP, et al. Evolution of extensively drug-resistant tuberculosis over four decades: whole genome sequencing and dating analysis of mycobacterium tuberculosis isolates from KwaZulu-Natal. PLoS Med. 2015;12(9):e1001880.

    Article  Google Scholar 

  32. 32.

    Zetola NM, Modongo C, Kip EC, Gross R, Bisson GP, Collman RG. Alcohol use and abuse among patients with multidrug-resistant tuberculosis in Botswana. Int J Tuberc Lung Dis. 2012;16(11):1529–34.

    CAS  Article  Google Scholar 

  33. 33.

    Pradhan N et al. Patterns of TB Drug-Resistance in a tertiary care facility in Pune, India. Clinical Microbiology: Open Access. 2013.

  34. 34.

    Marahatta SB, Kaewkungwal J, Ramasoota P, Singhasivanon P. Risk factors of multidrug resistant tuberculosis in central Nepal: a pilot study. Kathmandu Univ Med J (KUMJ). 2012;8(4):392–7.

    Article  Google Scholar 

  35. 35.

    Jia X, Chen J, Zhang S, Dai B, Long Q, Tang S. Implementing a “free” tuberculosis (TB) care policy under the integrated model in Jiangsu, China: practices and costs in the real world. Infect Dis Poverty. 2016;5(1):1.

    Article  Google Scholar 

  36. 36.

    Wang L, Zhang H, Ruan Y, Chin DP, Xia Y, Cheng S, et al. Tuberculosis prevalence in China, 1990–2010; a longitudinal analysis of national survey data. Lancet. 2014;383(9934):2057–64.

    Article  Google Scholar 

  37. 37.

    Wang L, Liu J, Chin DP. Progress in tuberculosis control and the evolving public-health system in China. Lancet. 2007;369(9562):691–6.

    Article  Google Scholar 

  38. 38.

    Ricks PM, Mavhunga F, Modi S, Indongo R, Zezai A, Lambert LA, et al. Characteristics of multidrug-resistant tuberculosis in Namibia. BMC Infect Dis. 2012;12(1):385.

    Article  Google Scholar 

  39. 39.

    Leung CC, Li T, Lam TH, Yew WW, Law WS, Tam CM, et al. Smoking and tuberculosis among the elderly in Hong Kong. Am J Respir Crit Care Med. 2004;170(9):1027–33.

    Article  Google Scholar 

  40. 40.

    Mesfin EA, Beyene D, Tesfaye A, Admasu A, Addise D, Amare M, et al. Drug-resistance patterns of mycobacterium tuberculosis strains and associated risk factors among multi drug-resistant tuberculosis suspected patients from Ethiopia. PLoS One. 2018;13(6):e0197737.

    Article  Google Scholar 

Download references


We thank Shandong Provincial Hospital, Shandong Provincial Chest Hospital, 13 municipal-level and 21 county-level local health departments for drug susceptibility data, sociodemographic and clinical data.


This work was supported by Department of Science & Technology of Shandong Province (CN)(No.2007GG30002033; No.2017GSF218052) and Jinan Science and Technology Bureau (CN) (No.201704100).

Author information




HCL and YFL conceived and designed the study. HCL, CBY, LLC, and LG directed its implementation including the data analysis and writing of the paper. WMS and YFL analyzed the data; YL, QYZ, JYL, TTX, SJL and NNT contributed materials/analytic tools; WMS, YFL and HCL wrote and revised the manuscript. All authors reviewed and approved the manuscript.

Corresponding author

Correspondence to Huai-chen Li.

Ethics declarations

Ethics approval and consent to participate

The protocols applied in this study were approved by the Ethics Committee of Shandong Provincial Hospital, affiliated with Shandong University (SPH) and and the Ethic Committee of Shandong Provincial Chest Hospital (SPCH), China. Before analysis, patient records were anonymized and deidentified.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.



Table 5 Changes in proportions of different primary drug-resistant mycobacterium tuberculosis subgroups, Shandong, China, 2004–2018a

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Song, Wm., Li, Yf., Ma, Xb. et al. Primary drug resistance of mycobacterium tuberculosis in Shandong, China, 2004–2018. Respir Res 20, 223 (2019).

Download citation


  • Tuberculosis
  • Primary drug resistance
  • Epidemiology
  • Risk factors
  • Joinpoint regression model