China's tuberculosis epidemic stems from historical expansion of four strains of Mycobacterium tuberculosis

Genomic analysis of lineage-specific transmission of multidrug resistance tuberculosis in China

OBJECTIVES

We investigated the genetic diversities and lineage-specific transmission dynamics of multidrug-resistant tuberculosis (MDR-TB), with the goal of determining the potential factors driving the MDR epidemics in China.

METHODS

We curated a large nationwide Mycobacterium tuberculosis (M. tuberculosis) whole genome sequence data set, including 1313 MDR strains. We reconstructed the phylogeny and mapped the transmission networks of MDR-TB across China using Bayesian inference. To identify drug-resistance variants linked to enhanced transmissibility, we employed ordinary least-squares (OLS) regression analysis.

RESULT

The majority of MDR-TB strains in China belong to lineage 2.2.1. Transmission chain analysis has indicated that the repeated and frequent transmission of L2.2.1 plays a central role in the establishment of MDR epidemic in China, but no occurrence of a large predominant MDR outbreak was detected. Using OLS regression, the most common single nucleotide polymorphisms (SNPs) associated with resistance to isoniazid (katG_p.Ser315Thr and katG_p.Ser315Asn) and rifampicin (rpoB_p.Ser450Leu, rpoB_p.His445Tyr, rpoB_p.His445Arg, rpoB_p.His445Asp, and rpoB_p.His445Asn) were more likely to be found in L2 clustered strains. Several putative compensatory mutations in rpoA, rpoC, and katG were significantly associated with clustering. The eastern, central, and southern regions of China had a high level of connectivity for the migration of L2 MDR strains throughout the country. The skyline plot showed distinct population size expansion dynamics for MDR-TB lineages in China.

CONCLUSION

MDR-TB epidemic in China is predominantly driven by the spread of highly transmissible Beijing strains. A range of drug-resistance mutations of L2 MDR-TB strains displayed minimal fitness costs and may facilitate their transmission.

Genomic and spatial analysis reveal the transmission dynamics of tuberculosis in areas with high incidence of Zhejiang, China: A prospective cohort study

In the mountainous, rural regions of eastern China, tuberculosis (TB) remains a formidable challenge; however, the long-term molecular epidemiological surveillance in these regions is limited. This study aimed to investigate molecular and spatial epidemiology of TB in two mountainous, rural counties of Zhejiang Province, China, from 2015 to 2021, to elucidate the recent transmission and drug-resistance profiles. The predominant Lineage 2 (L2) Beijing family accounted for 80.1% of total 532 sequenced Mycobacterium tuberculosis (Mtb) strains, showing consistent prevalence over seven years. Gene mutations associated with drug resistance were identified in 19.4% (103/532) of strains, including 47 rifampicin or isoniazid-resistant strains, eight multi-drug-resistant (MDR) strains, and five pre-extensively drug-resistant (pre-XDR) strains. Genomic clustering revealed 53 distinct clusters with an overall transmission clustering rate of 23.9% (127/532). Patients with a history of retreatment and those infected with L2 strains had a higher risk of recent transmission. Spatial and epidemiological analysis unveiled significant transmission hotspots, especially in densely populated urban areas, involving various public places such as medical institutions, farmlands, markets, and cardrooms. The study emphasizes the pivotal role of Beijing strains and urban-based TB transmission in the western mountainous regions in Zhejiang, highlighting the urgent requirement for specific interventions to mitigate the impact of TB in these unique communities.

Prevalence and genetic basis of Mycobacterium tuberculosis resistance to pretomanid in China

BACKGROUND

Pretomanid is a key component of new regimens for the treatment of drug-resistant tuberculosis (TB) which are being rolled out globally. However, there is limited information on the prevalence of pre-existing resistance to the drug.

METHODS

To investigate pretomanid resistance rates in China and its underlying genetic basis, as well as to generate additional minimum inhibitory concentration (MIC) data for epidemiological cutoff (ECOFF)/breakpoint setting, we performed MIC determinations in the Mycobacterial Growth Indicator Tube™ (MGIT) system, followed by WGS analysis, on 475 Mycobacterium tuberculosis (MTB) isolated from Chinese TB patients between 2013 and 2020.

RESULTS

We observed a pretomanid MIC distribution with a 99% ECOFF equal to 0.5 mg/L. Of the 15 isolates with MIC values > 0.5 mg/L, one (MIC = 1 mg/L) was identified as MTB lineage 1 (L1), a genotype previously reported to be intrinsically less susceptible to pretomanid, two were borderline resistant (MIC = 2-4 mg/L) and the remaining 12 isolates were highly resistant (MIC ≥ 16 mg/L) to the drug. Five resistant isolates did not harbor mutations in the known pretomanid resistant genes.

CONCLUSIONS

Our results further support a breakpoint of 0.5 mg/L for a non-L1 MTB population, which is characteristic of China. Further, our data point to an unexpected high (14/475, 3%) pre-existing pretomanid resistance rate in the country, as well as to the existence of yet-to-be-discovered pretomanid resistance genes.

Iron-related gene mutations driving global Mycobacterium tuberculosis transmission revealed by whole-genome sequencing

BACKGROUND

Iron plays a crucial role in the growth of Mycobacterium tuberculosis (M. tuberculosis). However, the precise regulatory mechanism governing this system requires further elucidation. Additionally, limited studies have examined the impact of gene mutations related to iron on the transmission of M. tuberculosis globally. This research aims to investigate the correlation between mutations in iron-related genes and the worldwide transmission of M. tuberculosis.

RESULTS

A total of 13,532 isolates of M. tuberculosis were included in this study. Among them, 6,104 (45.11%) were identified as genomic clustered isolates, while 8,395 (62.04%) were classified as genomic clade isolates. Our results showed that a total of 12 single nucleotide polymorphisms (SNPs) showed a positive correlation with clustering, such as Rv1469 (ctpD, C758T), Rv3703c (etgB, G1122T), and Rv3743c (ctpJ, G676C). Additionally, seven SNPs, including Rv0104 (T167G, T478G), Rv0211 (pckA, A302C), Rv0283 (eccB3, C423T), Rv1436 (gap, G654T), ctpD C758T, and etgB C578A, demonstrated a positive correlation with transmission clades across different countries. Notably, our findings highlighted the positive association of Rv0104 T167G, pckA A302C, eccB3 C423T, ctpD C758T, and etgB C578A with transmission clades across diverse regions. Furthermore, our analysis identified 78 SNPs that exhibited significant associations with clade size.

CONCLUSIONS

Our study reveals the link between iron-related gene SNPs and M. tuberculosis transmission, offering insights into crucial factors influencing the pathogenicity of the disease. This research holds promise for targeted strategies in prevention and treatment, advancing research and interventions in this field.

Transmission dynamics and phylogeography of Mycobacterium tuberculosis in China based on whole-genome phylogenetic analysis

OBJECTIVES

This study aimed to describe the lineage-specific transmissibility and epidemiological migration of Mycobacterium tuberculosis in China.

METHODS

We curated a large set of whole-genome sequences from 3204 M. tuberculosis isolates, including thousands of newly sequenced genomes, and applied a series of metrics to compare the transmissibility of M. tuberculosis strains between lineages and sublineages. The countrywide transmission patterns of major lineages were explored.

RESULTS

We found that lineage 2 (L2) was the most prevalent lineage in China (85.7%), with the major sublineage 2.2.1 (80.9%), followed by lineage 4 (L4) (13.8%), which comprises major sublineages 4.2 (1.5%), 4.4 (6.2%) and 4.5 (5.8%). We showed evidence for frequent cross-regional spread and large cluster formation of L2.2.1 strains, whereas L4 strains were relatively geographically restricted in China. Next, we applied a series of genomic indices to evaluate M. tuberculosis strain transmissibility and uncovered higher transmissibility of L2.2.1 compared with the L2.2.2 and L4 sublineages. Phylogeographic analysis showed that southern, eastern, and northern China were highly connected regions for countrywide L2.2.1 strain spread.

CONCLUSIONS

The present study provides insights into the different transmission and migration patterns of the major M. tuberculosis lineages in China and highlights that transmissible L2.2.1 is a threat to tuberculosis control.

Performance evaluation of core genome multilocus sequence typing for genotyping of Mycobacterium tuberculosis strains in China: based on multicenter, population-based collection

PURPOSE

To evaluate the performance of core genome multilocus sequence typing (cgMLST) for genotyping Mycobacterium tuberculosis (M.tuberculosis) Strains in regions where the lineage 2 strains predominate.

METHODS

We compared clustering by whole-genome SNP typing with cgMLST clustering in the analysis of WGS data of 6240 strains from five regions of China. Using both the receiver operating characteristic (ROC) curve and epidemiological investigation to determine the optimal threshold for defining genomic clustering by cgMLST. The performance of cgMLST was evaluated by quantifying the sensitivity, specificity and concordance of clustering between two methods. Logistic regression was used to gauge the impact of strain genetic diversity and lineage on cgMLST clustering.

RESULTS

The optimal threshold for cgMLST to define genomic clustering was determined to be ≤ 10 allelic differences between strains. The overall sensitivity and specificity of cgMLST averaged 99.6% and 96.3%, respectively; the concordance of clustering between two methods averaged 97.1%. Concordance was significantly correlated with strain genetic diversity and was 3.99 times (95% CI, 2.94-5.42) higher in regions with high genetic diversity (π > 1.55 × 10-4) compared to regions with low genetic diversity. The difference missed statistical significance, while concordance for lineage 2 strains (96.8%) was less than that for lineage 4 strains (98.3%). CONCLUSION : cgMLST showed a discriminatory power comparable to whole-genome SNP typing and could be used to genotype clinical M.tuberculosis strains in different regions of China. The discriminative power of cgMLST was significantly correlated with strain genetic diversity and was slightly lower with strains from regions with low genetic diversity.

Ongoing evolution of the Mycobacterium tuberculosis lactate dehydrogenase reveals the pleiotropic effects of bacterial adaption to host pressure

The bacterial determinants that facilitate Mycobacterium tuberculosis (Mtb) adaptation to the human host environment are poorly characterized. We have sought to decipher the pressures facing the bacterium in vivo by assessing Mtb genes that are under positive selection in clinical isolates. One of the strongest targets of selection in the Mtb genome is lldD2, which encodes a quinone-dependent L-lactate dehydrogenase (LldD2) that catalyzes the oxidation of lactate to pyruvate. Lactate accumulation is a salient feature of the intracellular environment during infection and lldD2 is essential for Mtb growth in macrophages. We determined the extent of lldD2 variation across a set of global clinical isolates and defined how prevalent mutations modulate Mtb fitness. We show the stepwise nature of lldD2 evolution that occurs as a result of ongoing lldD2 selection in the background of ancestral lineage-defining mutations and demonstrate that the genetic evolution of lldD2 additively augments Mtb growth in lactate. Using quinone-dependent antibiotic susceptibility as a functional reporter, we also find that the evolved lldD2 mutations functionally increase the quinone-dependent activity of LldD2. Using 13C-lactate metabolic flux tracing, we find that lldD2 is necessary for robust incorporation of lactate into central carbon metabolism. In the absence of lldD2, label preferentially accumulates in dihydroxyacetone phosphate (DHAP) and glyceraldehyde-3-phosphate (G3P) and is associated with a discernible growth defect, providing experimental evidence for accrued lactate toxicity via the deleterious buildup of sugar phosphates. The evolved lldD2 variants increase lactate incorporation to pyruvate while altering triose phosphate flux, suggesting both an anaplerotic and detoxification benefit to lldD2 evolution. We further show that the mycobacterial cell is transcriptionally sensitive to the changes associated with altered lldD2 activity which affect the expression of genes involved in cell wall lipid metabolism and the ESX- 1 virulence system. Together, these data illustrate a multifunctional role of LldD2 that provides context for the selective advantage of lldD2 mutations in adapting to host stress.

The Recent Transmission and Associated Risk Factor of Mycobacterium tuberculosis in Golmud City, China

Background

Tuberculosis (TB) remains a severe public health problem globally, and it is essential to comprehend the transmission pattern to control tuberculosis. Herein, we evaluated the drug-resistant characteristics, recent transmission, and associated risk factors of TB in Golmud, Qinghai, China.

Methods

In this study, we performed a population-based study of patients diagnosed with TB in Golmud from 2013 to 2018. Drug-susceptibility testing and whole-genome sequencing were performed on 133 Mycobacterium tuberculosis strains. The genomic clustering rate was calculated to evaluate the level of recent transmission. Risk factors were identified by logistic regression analysis.

Results

Our results showed that 46.97% (62/132) of strains were phylogenetically clustered and formed into 23 transmission clusters, suggesting a high recent transmission of TB in the area. 12.78% (17/133) strains were multidrug-resistant/rifampicin tuberculosis (MDR/RR-TB), with a high drug-resistant burden. Based on drug resistance gene analysis, we found 23 strains belonging to genotype MDR/RR-TB, where some strains may have borderline mutations. Among these strains, 65.2% (15/23) were found within putative transmission clusters. Additionally, risk factor analysis showed that recent transmission of TB happened more in patients with Tibetan nationality or older age.

Conclusion

Overall our study indicates that the recent transmissions of MTB strains, especially genotypic MDR/RR strains, drive the tuberculosis epidemic in Golmud, which could contribute to developing effective TB prevention and control strategies.

Polyketide synthases mutation in tuberculosis transmission revealed by whole genomic sequence, China, 2011-2019

Introduction: Tuberculosis (TB) is an infectious disease caused by a bacterium called Mycobacterium tuberculosis (Mtb). Previous studies have primarily focused on the transmissibility of multidrug-resistant (MDR) or extensively drug-resistant (XDR) Mtb. However, variations in virulence across Mtb lineages may also account for differences in transmissibility. In Mtb, polyketide synthase (PKS) genes encode large multifunctional proteins which have been shown to be major mycobacterial virulence factors. Therefore, this study aimed to identify the role of PKS mutations in TB transmission and assess its risk and characteristics. Methods: Whole genome sequences (WGSs) data from 3,204 Mtb isolates was collected from 2011 to 2019 in China. Whole genome single nucleotide polymorphism (SNP) profiles were used for phylogenetic tree analysis. Putative transmission clusters (≤10 SNPs) were identified. To identify the role of PKS mutations in TB transmission, we compared SNPs in the PKS gene region between "clustered isolates" and "non-clustered isolates" in different lineages. Results: Cluster-associated mutations in ppsA, pks12, and pks13 were identified among different lineage isolates. They were statistically significant among clustered strains, indicating that they may enhance the transmissibility of Mtb. Conclusion: Overall, this study provides new insights into the function of PKS and its localization in M. tuberculosis. The study found that ppsA, pks12, and pks13 may contribute to disease progression and higher transmission of certain strains. We also discussed the prospective use of mutant ppsA, pks12, and pks13 genes as drug targets.

Whole genome sequencing of drug-resistant Mycobacterium tuberculosis isolates in Victoria, Australia

OBJECTIVES

Whole genome sequencing (WGS) can identify clusters, transmission patterns, and drug resistance mutations. This is important in low-burden settings such as Australia, as it can assist in efficient contact tracing and surveillance.

METHODS

We conducted a retrospective cohort study using WGS from 155 genomically defined drug-resistant Mycobacterium tuberculosis (DR-TB) isolates collected between 2018-2021 in Victoria, Australia. Bioinformatic analysis was performed to identify resistance-conferring mutations, lineages, clusters and understand how local sequences compared with international context.

RESULTS

Of the 155 sequences, 42% were identified as lineage 2 and 35% as lineage 1; 65.8% (102/155) were isoniazid mono-resistant, 8.4% were multi-drug resistant TB and 5.8% were pre-extensively drug-resistant / extensively drug-resistant TB. The most common mutations were observed in katG and fabG1 genes, especially at Ser315Thr and fabG1 -15 C>T for first-line drugs. Ser450Leu was the most frequent mutation in rpoB gene. Phylogenetic analysis confirmed that Victorian DR-TB were associated with importation events. There was little evidence of local transmission with only five isolate pairs.

CONCLUSION

Isoniazid-resistant TB is the commonest DR-TB in Victoria, and the mutation profile is similar to global circulating DR-TB. Most cases are diagnosed among migrants with limited transmission. This study highlights the value of WGS in identification of clusters and resistance-conferring mutations. This information is crucial in supporting disease mitigation and treatment strategies.

Association of lineage 4.2.2 of Mycobacterium tuberculosis with the 63-bp deletion variant of the mpt64 gene

IMPORTANCE

To date, rapid diagnostic methods based on the MPT64 antigen assay are increasingly utilized to differentiate between non-tuberculous mycobacteria and TB disease in clinical settings. Furthermore, numerous novel techniques based on the MPT64 release assay are continuously being developed and applied for the identification of both pulmonary and extrapulmonary TB. However, the diagnostic accuracy of the MPT64 antigen assay is influenced by the presence of 63 bp deletion variants within the mpt64 gene. To our knowledge, this is the first report on the association between the 63 bp deletion variant in mpt64 and Mycobacterium tuberculosis L4.2.2 globally, which highlights the need for the cautious utilization of MPT64-based testing in regions where L4.2.2 isolates are prevalent, such as China and Vietnam, and MPT64 negative results should be confirmed with another assay. In addition, further studies on vaccine development and immunology based on MPT64 should consider these isolates with 63 bp deletion variant.

Association between fatty acid metabolism gene mutations and Mycobacterium tuberculosis transmission revealed by whole genome sequencing

BACKGROUND

Fatty acid metabolism greatly promotes the virulence and pathogenicity of Mycobacterium tuberculosis (M.tb). However, the regulatory mechanism of fatty acid metabolism in M.tb remains to be elucidated, and limited evidence about the effects of gene mutations in fatty acid metabolism on the transmission of M.tb was reported.

RESULTS

Overall, a total of 3193 M.tb isolates were included in the study, of which 1596 (50%) were genomic clustered isolates. Most of the tuberculosis isolates belonged to lineage2(n = 2744,85.93%), followed by lineage4(n = 439,13.75%) and lineage3(n = 10,0.31%).Regression results showed that the mutations of gca (136,605, 317G > C, Arg106Pro; OR, 22.144; 95% CI, 2.591-189.272), ogt(1,477,346, 286G > C ,Gly96Arg; OR, 3.893; 95%CI, 1.432-10.583), and rpsA (1,834,776, 1235 C > T, Ala412Val; OR, 3.674; 95% CI, 1.217-11.091) were significantly associated with clustering; mutations in gca and rpsA were also significantly associated with clustering of lineage2. Mutation in arsA(3,001,498, 885 C > G, Thr295Thr; OR, 6.278; 95% CI, 2.508-15.711) was significantly associated with cross-regional clusters. We also found that 20 mutation sites were positively correlated with cluster size, while 11 fatty acid mutation sites were negatively correlated with cluster size.

CONCLUSION

Our research results suggested that mutations in genes related to fatty acid metabolism were related to the transmission of M.tb. This research could help in the future control of the transmission of M.tb.

Association between two-component systems gene mutation and Mycobacterium tuberculosis transmission revealed by whole genome sequencing

BACKGROUND

Two-component systems (TCSs) assume a pivotal function in Mycobacterium tuberculosis (M.tuberculosis) growth. However, the exact regulatory mechanism of this system needs to be elucidated, and only a few studies have investigated the effect of gene mutations within TCSs on M.tuberculosis transmission. This research explored the relationship between TCSs gene mutation and the global transmission of (M.tuberculosis).

RESULTS

A total of 13531 M.tuberculosis strains were enrolled in the study. Most of the M.tuberculosis strains belonged to lineage4 (n=6497,48.0%), followed by lineage2 (n=5136,38.0%). Our results showed that a total of 36 single nucleotide polymorphisms (SNPs) were positively correlated with clustering of lineage2, such as Rv0758 (phoR, C820G), Rv1747(T1102C), and Rv1057(C1168T). A total of 30 SNPs showed positive correlation with clustering of lineage4, such as phoR(C182A, C1184G, C662T, T758G), Rv3764c (tcrY, G1151T), and Rv1747 C20T. A total of 19 SNPs were positively correlated with cross-country transmission of lineage2, such as phoR A575C, Rv1028c (kdpD, G383T, G1246C), and Rv1057 G817T. A total of 41 SNPs were positively correlated with cross-country transmission of lineage4, such as phoR(T758G, T327G, C284G), kdpD(G1755A, G625C), Rv1057 C980T, and Rv1747 T373G.

CONCLUSIONS

Our study identified that SNPs in genes of two-component systems were related to the transmission of M. tuberculosis. This finding adds another layer of complexity to M. tuberculosis virulence and provides insight into future research that will help to elucidate a novel mechanism of M. tuberculosis pathogenicity.

Tuberculosis Severity Predictive Model Using Mtb Variants and Serum Biomarkers in a Colombian Cohort of APTB Patients

Currently, tuberculosis (TB) is a bacterial infection caused by Mycobacterium tuberculosis (Mtb) that primarily affects the lungs. The severity of active pulmonary TB (APTB) is an important determinant of transmission, morbidity, mortality, disease experience, and treatment outcomes. Several publications have shown a high prevalence of disabling complications in individuals who have had severe APTB. Furthermore, certain strains of Mtb were associated with more severe disease outcomes. The use of biomarkers to predict severe APTB patients who are candidates for host-directed therapies, due to the high risk of developing post-tuberculous lung disease (PTLD), has not yet been implemented in the management of TB patients. We followed 108 individuals with APTB for 6 months using clinical tools, flow cytometry, and whole-genome sequencing (WGS). The median age of the study population was 26.5 years, and the frequency of women was 53.7%. In this study, we aimed to identify biomarkers that could help us to recognize individuals with APTB and improve our understanding of the immunopathology in these individuals. In this study, we conducted a follow-up on the treatment progress of 121 cases of APTB. The follow-up process commenced at the time of diagnosis (T0), continued with a control visit at 2 months (T2), and culminated in an exit appointment at 6 months following the completion of medical treatment (T6). People classified with severe APTB showed significantly higher levels of IL-6 (14.7 pg/mL; p < 0.05) compared to those with mild APTB (7.7 pg/mL) at T0. The AUCs for the ROC curves and the Matthews correlation coefficient values (MCC) demonstrate correlations ranging from moderate to very strong. We conducted WGS on 88 clinical isolates of Mtb, and our analysis revealed a total of 325 genes with insertions and deletions (Indels) within their coding regions when compared to the Mtb H37Rv reference genome. The pattern of association was found between serum levels of CHIT1 and the presence of Indels in Mtb isolates from patients with severe APTB. A key finding in our study was the high levels of CHIT1 in severe APTB patients. We identified a biomarker profile (IL-6, IFN-γ, IL-33, and CHIT1) that allows us to identify individuals with severe APTB, as well as the identification of a panel of polymorphisms (125) in clinical isolates of Mtb from individuals with severe APTB. Integrating these findings into a predictive model of severity would show promise for the management of APTB patients in the future, to guide host-directed therapy and reduce the prevalence of PTLD.

Mycobacterial DnaQ is an Alternative Proofreader Ensuring DNA Replication Fidelity

Remove of mis-incorporated nucleotides ensures replicative fidelity. Although the ε-exonuclease DnaQ is a well-established proofreader in the model organism Escherichia coli, proofreading in mycobacteria relies on the polymerase and histidinol phosphatase (PHP) domain of replicative polymerase despite the presence of an alternative DnaQ homolog. Here, we show that depletion of DnaQ in Mycolicibacterium smegmatis results in increased mutation rate, leading to AT-biased mutagenesis and elevated insertions/deletions in homopolymer tract. We demonstrated that mycobacterial DnaQ binds to the b-clamp and functions synergistically with the PHP domain to correct replication errors. Further, we found that the mycobacterial DnaQ sustains replicative fidelity upon chromosome topological stress. Intriguingly, we showed that a naturally evolved DnaQ variant prevalent in clinical Mycobacterium tuberculosis isolates enables hypermutability and is associated with extensive drug resistance. These results collectively establish that the alternative DnaQ functions in proofreading, and thus reveal that mycobacteria deploy two proofreaders to maintain replicative fidelity.

Parallel signatures of Mycobacterium tuberculosis and human Y-chromosome phylogeography support the Two Layer model of East Asian population history

The Two Layer hypothesis is fast becoming the favoured narrative describing East Asian population history. Under this model, hunter-gatherer groups who initially peopled East Asia via a route south of the Himalayas were assimilated by agriculturalist migrants who arrived via a northern route across Eurasia. A lack of ancient samples from tropical East Asia limits the resolution of this model. We consider insight afforded by patterns of variation within the human pathogen Mycobacterium tuberculosis (Mtb) by analysing its phylogeographic signatures jointly with the human Y-chromosome. We demonstrate the Y-chromosome lineages enriched in the traditionally hunter-gatherer groups associated with East Asia's first layer of peopling to display deep roots, low long-term effective population size, and diversity patterns consistent with a southern entry route. These characteristics mirror those of the evolutionarily ancient Mtb lineage 1. The remaining East Asian Y-chromosome lineage is almost entirely absent from traditionally hunter-gatherer groups and displays spatial and temporal characteristics which are incompatible with a southern entry route, and which link it to the development of agriculture in modern-day China. These characteristics mirror those of the evolutionarily modern Mtb lineage 2. This model paves the way for novel host-pathogen coevolutionary research hypotheses in East Asia.

Ongoing evolution of the Mycobacterium tuberculosis lactate dehydrogenase reveals the pleiotropic effects of bacterial adaption to host pressure

The bacterial determinants that facilitate Mycobacterium tuberculosis (Mtb) adaptation to the human host environment are poorly characterized. We have sought to decipher the pressures facing the bacterium in vivo by assessing Mtb genes that are under positive selection in clinical isolates. One of the strongest targets of selection in the Mtb genome is lldD2 , which encodes a quinone-dependent L-lactate dehydrogenase (LldD2) that catalyzes the oxidation of lactate to pyruvate. Lactate accumulation is a salient feature of the intracellular environment during infection and lldD2 is essential for Mtb growth in macrophages. We determined the extent of lldD2 variation across a set of global clinical isolates and defined how prevalent mutations modulates Mtb fitness. We show the stepwise nature of lldD2 evolution that occurs as a result of ongoing lldD2 selection in the background of ancestral lineage defining mutations and demonstrate that the genetic evolution of lldD2 additively augments Mtb growth in lactate. Using quinone-dependent antibiotic susceptibility as a functional reporter, we also find that the evolved lldD2 mutations functionally increase the quinone-dependent activity of LldD2. Using 13 C-lactate metabolic flux tracing, we find that lldD2 is necessary for robust incorporation of lactate into central carbon metabolism. In the absence of lldD2 , label preferentially accumulates in methylglyoxal precursors dihydroxyacetone phosphate (DHAP) and glyceraldehyde-3-phosphate (G3P) and is associated with a discernible growth defect, providing experimental evidence for accumulated lactate toxicity via a methylglyoxal pathway that has been proposed previously. The evolved lldD2 variants increase lactate incorporation to pyruvate but also alter flux in the methylglyoxal pathway, suggesting both an anaplerotic and detoxification benefit to lldD2 evolution. We further show that the mycobacterial cell is transcriptionally sensitive to the changes associated with altered lldD2 activity which affect the expression of genes involved in cell wall lipid metabolism and the ESX-1 virulence system. Together, these data illustrate a multifunctional role of LldD2 that provide context for the selective advantage of lldD2 mutations in adapting to host stress.

Disruption of Aldehyde Dehydrogenase 2 protects against bacterial infection

The ALDH2*2 (rs671) allele is one of the most common genetic mutations in humans, yet the positive evolutionary selective pressure to maintain this mutation is unknown, despite its association with adverse health outcomes. ALDH2 is responsible for the detoxification of metabolically produced aldehydes, including lipid-peroxidation end products derived from inflammation. Here, we demonstrate that host-derived aldehydes 4-hydroxynonenal (4HNE), malondialdehyde (MDA), and formaldehyde (FA), all of which are metabolized by ALDH2, are directly toxic to the bacterial pathogens Mycobacterium tuberculosis and Francisella tularensis at physiological levels. We find that Aldh2 expression in macrophages is decreased upon immune stimulation, and that bone marrow-derived macrophages from Aldh2 -/- mice contain elevated aldehydes relative to wild-type mice. Macrophages deficient for Aldh2 exhibited enhanced control of Francisella infection. Finally , mice lacking Aldh2 demonstrated increased resistance to pulmonary infection by M. tuberculosis , including in a hypersusceptible model of tuberculosis, and were also resistant to Francisella infection. We hypothesize that the absence of ALDH2 contributes to the host's ability to control infection by pathogens such as M. tuberculosis and F. tularensis , and that host-derived aldehydes act as antimicrobial factors during intracellular bacterial infections.

One sentence summary

Aldehydes produced by host cells contribute to the control of bacterial infections.

The Global Success of Mycobacterium tuberculosis Modern Beijing Family Is Driven by a Few Recently Emerged Strains

Strains of the Mycobacterium tuberculosis complex (MTBC) Beijing family aroused concern because they were often found in clusters and appeared to be exceptionally transmissible. However, it was later found that strains of the Beijing family were heterogeneous, and the transmission advantage was restricted to sublineage L2.3 or modern Beijing. In this study, we analyzed the previously published genome sequences of 7,896 L2.3 strains from 51 different countries. Using BEAST software to approximate the temporal emergence of L2.3, our calculations suggest that L2.3 initially emerged in northern East Asia during the early 15th century and subsequently diverged into six phylogenetic clades, identified as L2.3.1 through L2.3.6. Using terminal branch length and genomic clustering as proxies for transmissibility, we found that the six clades displayed distinct population dynamics, with the three recently emerged clades (L2.3.4 to L2.3.6) exhibiting significantly higher transmissibility than the older three clades (L2.3.1 to L2.3.3). Of the Beijing family strains isolated outside East Asia, 83.1% belonged to the clades L2.3.4 to L2.3.6, which were also associated with more cross-border transmission. This work reveals the heterogeneity in sublineage L2.3 and demonstrates that the global success of Beijing family strains is driven by the three recently emerged L2.3 clades. IMPORTANCE The recent population dynamics of the global tuberculosis epidemic are heavily shaped by Mycobacterium tuberculosis complex (MTBC) strains with enhanced transmissibility. The infamous Beijing family strain stands out because it has rapidly spread throughout the world. Identifying the strains responsible for the global expansion and tracing their evolution should help to understand the nature of high transmissibility and develop effective strategies to control transmission. In this study, we found that the L2.3 sublineage diversified into six phylogenetic clades (L2.3.1 to L2.3.6) with various transmission characteristics. Clades L2.3.4 to L2.3.6 exhibited significantly higher transmissibility than clades L2.3.1 to L2.3.3, which helps explain why more than 80% of Beijing family strains collected outside East Asia belong to these three clades. We conclude that the global success of L2.3 was not caused by the entire L2.3 sublineage but rather was due to the rapid expansion of L2.3.4 to L2.3.6. Tracking the transmission of L2.3.4 to L2.3.6 strains can help to formulate targeted TB prevention and control.

Whole Genomic Analysis Revealed High Genetic Diversity and Drug-Resistant Characteristics of Mycobacterium tuberculosis in Guangxi, China

Background

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), is a major public health issue in China. Nevertheless, the prevalence and drug resistance characteristics of isolates vary in different regions and provinces. In this study, we investigated the population structure, transmission dynamics and drug-resistant profiles of Mtb in Guangxi, located on the border of China.

Methods

From February 2016 to April 2017, 462 clinical M. tuberculosis isolates were selected from 5 locations in Guangxi. Drug-susceptibility testing was performed using 6 common anti-tuberculosis drugs. The genotypic drug resistance and transmission dynamics were analyzed by the whole genome sequence.

Results

Our data showed that the Mtb in Guangxi has high genetic diversity including Lineage 1 to Lineage 4, and mostly belong to Lineage 2 and Lineage 4. Novelty, 9.6% of Lineage 2 isolates were proto-Beijing genotype (L2.1), which is rare in China. About 12.6% of isolates were phylogenetically clustered and formed into 28 transmission clusters. We observed that the isolates with the high resistant rate of isoniazid (INH, 21.2%), followed by rifampicin (RIF, 13.2%), and 6.7%, 12.1%, 6.7% and 1.9% isolates were resistant to ethambutol (EMB), streptomycin (SM), ofloxacin (OFL) and kanamycin (KAN), respectively. Among these, 6.5% and 3.3% of isolates belong to MDR-TB and Pre-XDR, respectively, with a high drug-resistant burden. Genetic analysis identified the most frequently encountered mutations of INH, RIF, EMB, SM, OFL and KAN were katG_Ser315Thr (62.2%), rpoB_Ser450Leu (42.6%), embB_Met306Vol (45.2%), rpsL_Lys43Arg (53.6%), gyrA_Asp94Gly (29.0%) and rrs_A1401G (66.7%), respectively. Additionally, we discovered that isolates from border cities are more likely to be drug-resistant than isolates from non-border cities.

Conclusion

Our findings provide a deep analysis of the genomic population characteristics and drug-resistant of M. tuberculosis in Guangxi, which could contribute to developing effective TB prevention and control strategies.

Transmission of Mycobacterium tuberculosis in schools: a molecular epidemiological study using whole-genome sequencing in Guangzhou, China

Background

China is a country with a high burden of tuberculosis (TB). TB outbreaks are frequent in schools. Thus, understanding the transmission patterns is crucial for controlling TB.

Method

In this genomic epidemiological study, the conventional epidemiological survey data combined with whole-genome sequencing was used to assess the genotypic distribution and transmission characteristics of Mycobacterium tuberculosis strains isolated from patients with TB attending schools during 2015 to 2019 in Guangzhou, China.

Result

The TB incidence was mainly concentrated in regular secondary schools and technical and vocational schools. The incidence of drug resistance among the students was 16.30% (22/135). The phylogenetic tree showed that 79.26% (107/135) and 20.74% (28/135) of the strains belonged to lineage 2 (Beijing genotype) and lineage 4 (Euro-American genotype), respectively. Among the 135 isolates, five clusters with genomic distance within 12 single nucleotide polymorphisms were identified; these clusters included 10 strains, accounting for an overall clustering rate of 7.4% (10/135), which showed a much lower transmission index. The distance between the home or school address and the interval time of symptom onset or diagnosis indicated that campus dissemination and community dissemination may be existed both, and community dissemination is the main.

Conclusion and recommendation

TB cases in Guangzhou schools were mainly disseminated and predominantly originated from community transmission. Accordingly, surveillance needs to be strengthened to stop the spread of TB in schools.

Genomic Sequencing Profiles of Mycobacterium tuberculosis in Mandalay Region, Myanmar

This study aimed to characterize whole-genome sequencing (WGS) information of Mycobacterium tuberculosis (Mtb) in the Mandalay region of Myanmar. It was a cross-sectional study conducted with 151 Mtb isolates obtained from the fourth nationwide anti-tuberculosis (TB) drug-resistance survey. Frequency of lineages 1, 2, 3, and 4 were 55, 65, 9, and 22, respectively. The most common sublineage was L1.1.3.1 (n = 31). Respective multi-drug resistant tuberculosis (MDR-TB) frequencies were 1, 1, 0, and 0. Four clusters of 3 (L2), 2 (L4), 2 (L1), and 2 (L2) isolates defined by a 20-single-nucleotide variant (SNV) cutoff were detected. Simpson's index for sublineages was 0.0709. Such high diversity suggests that the area probably had imported Mtb from many geographical sources. Relatively few genetic clusters and MDR-TB suggest there is a chance the future control will succeed if it is carried out properly.

Deciphering the Tangible Spatio-Temporal Spread of a 25-Year Tuberculosis Outbreak Boosted by Social Determinants

Outbreak strains of Mycobacterium tuberculosis are promising candidates as targets in the search for intrinsic determinants of transmissibility, as they are responsible for many cases with sustained transmission; however, the use of low-resolution typing methods and restricted geographical investigations represent flaws in assessing the success of long-lived outbreak strains. We can now address the nature of outbreak strains by combining large genomic data sets and phylodynamic approaches. We retrospectively sequenced the whole genome of representative samples assigned to an outbreak circulating in the Canary Islands (the GC strain) since 1993, which accounts for ~20% of local tuberculosis cases. We selected a panel of specific single nucleotide polymorphism (SNP) markers for an in-silico search for additional outbreak-related sequences within publicly available tuberculosis genomic data. Using this information, we inferred the origin, spread, and epidemiological parameters of the GC strain. Our approach allowed us to accurately trace the historical and more recent dispersion of the GC strain. We provide evidence of a highly successful nature within the Canarian archipelago but limited expansion abroad. Estimation of epidemiological parameters from genomic data disagree with a distinctive biology of the GC strain. With the increasing availability of genomic data allowing for the accurate inference of strain spread and critical epidemiological parameters, we can now revisit the link between Mycobacterium tuberculosis genotypes and transmission, as is routinely carried out for SARS-CoV-2 variants of concern. We demonstrate that social determinants rather than intrinsically higher bacterial transmissibility better explain the success of the GC strain. Importantly, our approach can be used to trace and characterize strains of interest worldwide. IMPORTANCE Infectious disease outbreaks represent a significant problem for public health. Tracing outbreak expansion and understanding the main factors behind emergence and persistence remain critical to effective disease control. Our study allows researchers and public health authorities to use Whole-Genome Sequencing-based methods to trace outbreaks, and shows how available epidemiological information helps to evaluate the factors underpinning outbreak persistence. Taking advantage of all the freely available information placed in public repositories, researchers can accurately establish the expansion of an outbreak beyond original boundaries, and determine the potential risk of a strain to inform health authorities which, in turn, can define target strategies to mitigate expansion and persistence. Finally, we show the need to evaluate strain transmissibility in different geographic contexts to unequivocally associate spread to local or pathogenic factors, an important lesson taken from genomic surveillance of SARS-CoV-2.

Distribution and origins of Mycobacterium tuberculosis L4 in Southeast Asia

Molecular and genomic studies have revealed that Mycobacterium tuberculosis Lineage 4 (L4, Euro-American lineage) emerged in Europe before becoming distributed around the globe by trade routes, colonial migration and other historical connections. Although L4 accounts for tens or hundreds of thousands of tuberculosis (TB) cases in multiple Southeast Asian countries, phylogeographical studies have either focused on a single country or just included Southeast Asia as part of a global analysis. Therefore, we interrogated public genomic data to investigate the historical patterns underlying the distribution of L4 in Southeast Asia and surrounding countries. We downloaded 6037 genomes associated with 29 published studies, focusing on global analyses of L4 and Asian studies of M. tuberculosis. We identified 2256 L4 genomes including 968 from Asia. We show that 81 % of L4 in Thailand, 51 % of L4 in Vietnam and 9 % of L4 in Indonesia belong to sub-lineages of L4 that are rarely seen outside East and Southeast Asia (L4.2.2, L4.4.2 and L4.5). These sub-lineages have spread between East and Southeast Asian countries, with no recent European ancestor. Although there is considerable uncertainty about the exact direction and order of intra-Asian M. tuberculosis dispersal, due to differing sampling frames between countries, our analysis suggests that China may be the intermediate location between Europe and Southeast Asia for two of the three predominantly East and Southeast Asian L4 sub-lineages (L4.2.2 and L4.5). This new perspective on L4 in Southeast Asia raises the possibility of investigating host population-specific evolution and highlights the need for more structured sampling from Southeast Asian countries to provide more certainty of the historical and current routes of dispersal.

Tuberculosis treatment failure associated with evolution of antibiotic resilience

The widespread use of antibiotics has placed bacterial pathogens under intense pressure to evolve new survival mechanisms. Genomic analysis of 51,229 Mycobacterium tuberculosis (Mtb)clinical isolates has identified an essential transcriptional regulator, Rv1830, herein called resR for resilience regulator, as a frequent target of positive (adaptive) selection. resR mutants do not show canonical drug resistance or drug tolerance but instead shorten the post-antibiotic effect, meaning that they enable Mtb to resume growth after drug exposure substantially faster than wild-type strains. We refer to this phenotype as antibiotic resilience. ResR acts in a regulatory cascade with other transcription factors controlling cell growth and division, which are also under positive selection in clinical isolates of Mtb. Mutations of these genes are associated with treatment failure and the acquisition of canonical drug resistance.

An ancestral mycobacterial effector promotes dissemination of infection

The human pathogen Mycobacterium tuberculosis typically causes lung disease but can also disseminate to other tissues. We identified a M. tuberculosis (Mtb) outbreak presenting with unusually high rates of extrapulmonary dissemination and bone disease. We found that the causal strain carried an ancestral full-length version of the type VII-secreted effector EsxM rather than the truncated version present in other modern Mtb lineages. The ancestral EsxM variant exacerbated dissemination through enhancement of macrophage motility, increased egress of macrophages from established granulomas, and alterations in macrophage actin dynamics. Reconstitution of the ancestral version of EsxM in an attenuated modern strain of Mtb altered the migratory mode of infected macrophages, enhancing their motility. In a zebrafish model, full-length EsxM promoted bone disease. The presence of a derived nonsense variant in EsxM throughout the major Mtb lineages 2, 3, and 4 is consistent with a role for EsxM in regulating the extent of dissemination.

Need for speed: Key driver of host cell migration varies among mycobacteria

In this issue of Cell, Saelens et al. describe a new function for mycobacterial Type VII secretion systems: manipulation of host cell migration. They find that a substantial proportion of global TB cases arise from bacteria lacking this function, raising questions about its role in pathoadaptation of Mycobacterium tuberculosis.

Drug Degradation Caused by mce3R Mutations Confers Contezolid (MRX-I) Resistance in Mycobacterium tuberculosis

Contezolid (MRX-I), a safer antibiotic of the oxazolidinone class, is a promising new antibiotic with potent activity against Mycobacterium tuberculosis (MTB) both in vitro and in vivo. To identify resistance mechanisms of contezolid in MTB, we isolated several in vitro spontaneous contezolid-resistant MTB mutants, which exhibited 16-fold increases in the MIC of contezolid compared with the parent strain but were still unexpectedly susceptible to linezolid. Whole-genome sequencing revealed that most of the contezolid-resistant mutants bore mutations in the mce3R gene, which encodes a transcriptional repressor. The mutations in mce3R led to markedly increased expression of a monooxygenase encoding gene Rv1936. We then characterized Rv1936 as a putative flavin-dependent monooxygenase that catalyzes the degradation of contezolid into its inactive 2,3-dihydropyridin-4-one (DHPO) ring-opened metabolites, thereby conferring drug resistance. While contezolid is an attractive drug candidate with potent antimycobacterial activity and low toxicity, the occurrence of mutations in Mce3R should be considered when designing combination therapy using contezolid for treating tuberculosis.

High proportion of tuberculosis transmission among social contacts in rural China: a 12-year prospective population-based genomic epidemiological study

Tuberculosis (TB) is more prevalent in rural than urban areas in China, and delineating TB transmission patterns in rural populations could improve TB control. We conducted a prospective population-based study of culture-positive pulmonary TB patients diagnosed between July 1, 2009 and December 31, 2020 in two rural counties in China. Genomic clusters were defined with a threshold distance of 12-single-nucleotide-polymorphisms, based on whole-genome sequencing. Risk factors for clustering were identified by logistic regression. Transmission links were sought through epidemiological investigation of genomic-clustered patients. Of 1517 and 751 culture-positive pulmonary TB patients in Wusheng and Wuchang counties, respectively, 1289 and 699 strains were sequenced. Overall, 624 (31.4%, 624/1988) patients were grouped into 225 genomic clusters. Epidemiological links were confirmed in 41.8% (196/469) of clustered isolates, including family (32.7%, 64/196) and social contacts (67.3%, 132/196). Social contacts were generally with relatives, within the community or in shared aggregated settings outside the community, but the proportion of clustered contacts in each category differed between the two sites. The time interval between diagnosis of student cases and contacts was significantly shorter than family and social contacts, probably due to enhanced student contact screening. Transmission of multidrug-resistant (MDR) strains was likely responsible for 81.4% (83/102) of MDR-TB cases, with minimal acquisition of additional resistance mutations. A large proportion of TB transmission in rural China occurred among social contacts, suggesting that active screening and aggressive contact tracing could benefit TB control, but contact screening should be tailored to local patterns of social interactions.

Rapid adaptation of a complex trait during experimental evolution of Mycobacterium tuberculosis

Tuberculosis (TB), caused by Mycobacterium tuberculosis (M. tb), is a leading cause of death due to infectious disease. TB is not traditionally associated with biofilms, but M. tb biofilms are linked with drug and immune tolerance and there is increasing recognition of their contribution to the recalcitrance of TB infections. Here, we used M. tb experimental evolution to investigate this complex phenotype and identify candidate loci controlling biofilm formation. We identified novel candidate loci, adding to our understanding of the genetic architecture underlying M. tb biofilm development. Under selective pressure to grow as a biofilm, regulatory mutations rapidly swept to fixation and were associated with changes in multiple traits, including extracellular matrix production, cell size, and growth rate. Genetic and phenotypic paths to enhanced biofilm growth varied according to the genetic background of the parent strain, suggesting that epistatic interactions are important in M. tb adaptation to changing environments.

Tuberculosis reinfection and relapse in eastern China: A prospective study using whole-genome sequencing

OBJECTIVES

Tuberculosis recurrence after an initial successful treatment episode can occur from either reinfection or relapse. In a population-based sample and whole genome sequencing (WGS) in eastern China, we aimed to evaluate risk factors for tuberculosis recurrence, and assess the proportion of recurrence due to either reinfection or relapse.

METHODS

Successfully treated pulmonary tuberculosis patients with sputum culture positive results were recruited from five cities in Jiangsu Province from 2013-2015 and followed for two years for tuberculosis recurrence. Among patients developing a second tuberculosis episode, WGS was performed to distinguish relapse or reinfection through a distance threshold of 6-single-nucleotide polymorphisms (SNP). We analyzed risk factors for recurrence and epidemiological characteristics of different types of recurrent patients.

RESULTS

Of 1,897 successful treated tuberculosis patients, 7.4% (141/1879) developed recurrent tuberculosis. Compared with non-recurrent tuberculosis, patients were at higher risk of recurrence in older age (Adjusted Odds Ratio [AOR], 1.02 for each additional year; 95% CI, 1.01-1.03, P=0.003), patients previously treated for tuberculosis (AOR=2.22; 95% CI, 1.52-3.26, P<0.001), or with bilateral cavities (AOR, 1.56; 95% CI, 1.05-2.32, P=0.029). Among 27.0% (38/141) recurrent tuberculosis patients with successfully sequenced pairs, relapse was substantially more common than reinfection (71.1% versus 28.9%, P=0.014).

CONCLUSIONS

Endogenous relapse was significantly more common than exogenous reinfection in the first two years after treatment in eastern China. Prioritization of high-risk groups for recurrence, such as the elderly, with a previous tuberculosis diagnosis, or with bilateral cavities, may provide opportunities to reduce post-tuberculosis morbidity.

Prediction of drug resistance profile of multidrug-resistant Mycobacterium tuberculosis (MDR-MTB) isolates from newly diagnosed case by whole genome sequencing (WGS): a study from a high tuberculosis burden country

OBJECTIVES

Our aim was to assess the ability of the Whole-genome sequencing (WGS) in predicting drug resistance profile of multidrug-resistant mycobacterium tuberculosis (MDR-MTB) from newly diagnosed cases in China.

METHODS

We validated the Phenotypic drug Sensitivity Test (pDST) for 12 anti-tuberculosis drugs using the Bactec MGIT 960 system. We described the characteristics of the isolates enrolled and compared the pDST results with resistance profiles predicted by WGS.

RESULTS

The pDST showed that of the 43 isolates enrolled, 25.6% were sensitive to rifabutin (RFB); 97.7%、97.7%、93.0% and 93.0% were sensitive to cycloserine (Cs), amikacin/kanamycin (Ak/Km), para-aminosalicylic acid (Pas) and ethionamide Eto), respectively; 18.6% were resistant to fluoroquinolones (FQs) or second-line injections. Genotype DST determined by WGS of Ak/Km、Eto and RFP reached high consistency to 97.7% compared with pDST, followed by moxifloxacin (Mfx) 95.3%, levofloxaci (Lfx) and Pas 93%, streptomycin (Sm) 90.3%. The genotype DST of RFB and EMB showed low consistency with the pDST of 67.2 and 79.1%. WGS also detected 27.9% isolates of pyrazinamide(PZA)-related drug-resistant mutation. No mutations associated with linezolid (Lzd), bedaquiline (Bdq) and clofazimine (Cfz) were detectd.

CONCLUSIONS

WGS has the potential to infer resistance profiles without time-consuming phenotypic methods, which could be provide a basis to formulate reasonable treatment in high TB burden areas.

The aldehyde hypothesis: metabolic intermediates as antimicrobial effectors

There are many reactive intermediates found in metabolic pathways. Could these potentially toxic molecules be exploited for an organism's benefit? We propose that during certain microbial infections, the production of inherently reactive aldehydes by an infected host is a previously unappreciated innate immune defence mechanism. While there has been a significant focus on the effects of aldehydes on mammalian physiology, the idea that they might be exploited or purposefully induced to kill pathogens is new. Given that aldehydes are made as parts of metabolic programmes that accompany immune cell activation by the cytokine interferon-gamma (IFN-γ) during infections, we hypothesize that aldehydes are among the arsenal of IFN-γ-inducible effectors needed for pathogen control.

Distribution and identification of Mycobacterium tuberculosis lineage in Kashgar prefecture

Objectives

Kashgar prefecture is an important transportation and trade hub with a high incidence of tuberculosis. The following study analyzed the composition and differences in Mycobacterium tuberculosis (M.tb) lineage and specific tags to distinguish the lineage of the M.tb in Kashgar prefecture, thus providing a basis for the classification and diagnosis of tuberculosis in this area.

Methods

Whole-genome sequencing (WGS) of 161 M.tb clinical strains was performed. The phylogenetic tree was constructed using Maximum Likelihood (ML) based on single nucleotide polymorphisms (SNPs) and verified through principal component analysis (PCA). The composition structure of M.tb in different regions was analyzed by combining geographic information.

Results

M.tb clinical strains were composed of lineage 2 (73/161, 45.34%), lineage 3 (52/161, 32.30%) and lineage 4 (36/161, 22.36%). Moreover, the 3 lineages were subdivided into 11 sublineages, among which lineage 2 included lineage 2.2.2/Asia Ancestral 1 (9/73, 12.33%), lineage 2.2.1-Asia Ancestral 2 (9/73, 12.33%), lineage 2.2.1-Asia Ancestral 3 (18/73, 24.66%), and lineage 2.2.1-Modern Beijing (39/73, 53.42%). Lineage 3 included lineage 3.2 (14/52, 26.92%) and lineage 3.3 (38/52, 73.08%), while lineage 4 included lineage 4.1 (3/36, 8.33%), lineage 4.2 (2/36, 5.66%), lineage 4.4.2 (1/36, 2.78%), lineage 4.5 (28/36, 77.78%) and lineage 4.8 (2/36, 5.66%), all of which were consistent with the PCA results. One hundred thirty-six markers were proposed for discriminating known circulating strains. Reconstruction of a phylogenetic tree using the 136 SNPs resulted in a tree with the same number of delineated clades. Based on geographical location analysis, the composition of Lineage 2 in Kashgar prefecture (45.34%) was lower compared to other regions in China (54.35%-90.27%), while the composition of Lineage 3 (32.30%) was much higher than in other regions of China (0.92%-2.01%), but lower compared to the bordering Pakistan (70.40%).

Conclusion

Three lineages were identified in M.tb clinical strains from Kashgar prefecture, with 136 branch-specific SNP. Kashgar borders with countries that have a high incidence of tuberculosis, such as Pakistan and India, which results in a large difference between the M.tb lineage and sublineage distribution in this region and other provinces of China.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07307-4.

A Whole-Genome Sequencing-Based Study to Delineate the Risk and Characteristics of Tuberculosis Transmission in an Insular Population Over 10 Years in Shanghai

Background

Tuberculosis (TB) has remained a tough problem in China. This study aims to identify the risk of tuberculosis transmission and to assess its characteristics.

Methods

We performed a molecular epidemiological study for patients with culture-positive Mycobacterium tuberculosis (M. tuberculosis) in Shanghai, from 2009 to 2018. Demographic information was obtained from the Tuberculosis Information Management System. Whole-genome sequencing (WGS) was conducted with a threshold of 12 single-nucleotide polymorphisms (SNPs) to distinguish the genomic cluster. To analyze the characteristics of TB transmission, the contact investigation for clustered cases was performed.

Results

In total, 94 (27.25%) of the 345 enrolled patients were grouped into 42 genomic clusters, indicating local transmission of M. tuberculosis strains. Compared to a health system delay <14 days, patients with a health system delay ≥14 days [adjusted odds ratios (AOR) = 2.57, 95% confidence interval (CI): 1.34–4.95] were more likely to be clustered. Patients under 65 years old (AOR = 3.11, 95% CI: 1.76–5.49), residents (AOR = 2.43, 95% CI: 1.18–4.99), and Beijing genotype strains (AOR = 3.35, 95% CI: 1.32–8.53) were associated with increased risk of clustering. Interestingly, patients with resistance to isoniazid (AOR = 2.36, 95% CI: 1.15–4.88) had a higher risk of transmission. Sixteen confirmed/probable epidemiological links were identified. Local transmission of imported cases and household transmission were prominent.

Conclusion

Health system delay is a crucial factor for TB transmission. Patients with resistance to isoniazid should be priority targets for contact investigation to reduce transmission.

Connection between two historical tuberculosis outbreak sites in Japan, Honshu, by a new ancestral Mycobacterium tuberculosis L2 sublineage

By gathering 680 publicly available Sequence Read Archives from isolates of Mycobacterium tuberculosis complex (MTBC) including 190 belonging to the lineage 2 Beijing , and using an in-house bioinformatical pipeline, the TB-Annotator , that analyses more than 50 000 characters, we describe herein a new L2 sublineage from 20 isolates found in the Tochigi province, (Japan), that we designate as asia ancestral 5 (AAnc5). These isolates harbour a number of specific criteria (42 SNPs) and their intra-cluster pairwise distance suggests historical and not epidemiological transmission. These isolates harbour a mutation in rpoC , and do not fulfil, any of the modern Beijing lineage criteria, nor any of the other ancestral Beijing lineages described so far. Asia ancestral 5 isolates do not possess mutT2 58 and ogt 12 characteristics of modern Beijing , but possess ancestral Beijing SNPs characteristics. By looking into the literature, we found a reference isolate ID381, described in Kobe and Osaka belonging to the ‘G3’ group, sharing 36 out of the 42 specific SNPs found in AAnc5. We also assessed the intermediate position of the asia ancestral 4 (AAnc4) sublineage recently described in Thailand and propose an improved classification of the L2 that now includes AAnc4 and AAnc5. By increasing the recruitment into TB-Annotator to around 3000 genomes (including 642 belonging to L2), we confirmed our results and discovered additional historical ancestral L2 branches that remain to be investigated in more detail. We also present, in addition, some anthropological and historical data from Chinese and Japan history of tuberculosis, as well as from Korea, that could support our results on L2 evolution. This study shows that the reconstruction of the early history of tuberculosis in Asia is likely to reveal complex patterns since its emergence.

Whole genome sequencing of Mycobacterium tuberculosis for prediction of drug resistance

Whole-genome sequencing (WGS) has shown tremendous potential in rapid diagnosis of drug-resistant tuberculosis (TB). In the current study, we performed WGS on drug-resistant Mycobacterium tuberculosis isolates obtained from Shanghai (n = 137) and Russia (n = 78). We aimed to characterise the underlying and high-frequency novel drug-resistance-conferring mutations, and also create valuable combinations of resistance mutations with high predictive sensitivity to predict multidrug- and extensively drug-resistant tuberculosis (MDR/XDR-TB) phenotype using a bootstrap method. Most strains belonged to L2.2, L4.2, L4.4, L4.5 and L4.8 lineages. We found that WGS could predict 82.07% of phenotypically drug-resistant domestic strains. The prediction sensitivity for rifampicin (RIF), isoniazid (INH), ethambutol (EMB), streptomycin (STR), ofloxacin (OFL), amikacin (AMK) and capreomycin (CAP) was 79.71%, 86.30%, 76.47%, 88.37%, 83.33%, 70.00% and 70.00%, respectively. The mutation combination with the highest sensitivity for MDR prediction was rpoB S450L + rpoB H445A/P + katG S315T + inhA I21T + inhA S94A, with a sensitivity of 92.17% (0.8615, 0.9646), and the mutation combination with highest sensitivity for XDR prediction was rpoB S450L + katG S315T + gyrA D94G + rrs A1401G, with a sensitivity of 92.86% (0.8158, 0.9796). The molecular information presented here will be of particular value for the rapid clinical detection of MDR- and XDR-TB isolates through laboratory diagnosis.

Revised nomenclature and SNP barcode for Mycobacterium tuberculosis lineage 2

Mycobacterium tuberculosis (Mtb) lineage 2 (L2) strains are present globally, contributing to a widespread tuberculosis (TB) burden, particularly in Asia where both prevalence of TB and numbers of drug resistant TB are highest. The increasing availability of whole-genome sequencing (WGS) data worldwide provides an opportunity to improve our understanding of the global genetic diversity of Mtb L2 and its association with the disease epidemiology and pathogenesis. However, existing L2 sublineage classification schemes leave >20 % of the Modern Beijing isolates unclassified. Here, we present a revised SNP-based classification scheme of L2 in a genomic framework based on phylogenetic analysis of >4000 L2 isolates from 34 countries in Asia, Eastern Europe, Oceania and Africa. Our scheme consists of over 30 genotypes, many of which have not been described before. In particular, we propose six main genotypes of Modern Beijing strains, denoted L2.2.M1–L2.2.M6. We also provide SNP markers for genotyping L2 strains from WGS data. This fine-scale genotyping scheme, which can classify >98 % of the studied isolates, serves as a basis for more effective monitoring and reporting of transmission and outbreaks, as well as improving genotype-phenotype associations such as disease severity and drug resistance. This article contains data hosted by Microreact.

Tuberculosis Outbreak in an Educational Institution in Henan Province, China

On June 17, 2018, a case of pulmonary tuberculosis (TB) was reported among students at a senior high school in Luoning, China. The outbreak encompassed a total of 23 cases along with TB screening in the whole school by means of PPD and chest X-ray. By the end of September 2018, the entire 9 cases cultured positive had epidemiological association. All of the 9 Mycobacterium tuberculosis (Mtb) isolates available were sensitive to all drugs tested and had similar spoligotyping and 15 loci mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) profile. Whole-genome sequencing (WGS) of the Mtb isolates revealed 20 variable nucleotide positions within 8 cases, indicating a clonal outbreak. The index case, which was first identified and diagnosed, is separated from the cluster by a minimum number of 95 distinct SNPs. Minimum distance spanning tree (MST) indicted that the 8 cases were indeed part of a single transmission chain. It was concluded that this is an epidemic situation of TB outbreak exposed by the aggrieved index case at school, which was caused by the veiled infectious case wherein a student was suffering from TB and attending school simultaneously.

Genetic composition and evolution of the prevalent Mycobacterium tuberculosis lineages 2 and 4 in the Chinese and Zhejiang Province populations

Background

There are seven human-adaptation lineages of Mycobacterium tuberculosis (Mtb). Tuberculosis (TB) dissemination is strongly influenced by human movements and host genetics. The detailed lineage distribution evolution of Mtb in Zhejiang Province is unknown. We aim to determine how different sub-lineages are transmitted and distributed within China and Zhejiang Province.

Methods

We analysed whole-genome sequencing data for a worldwide collection of 1154 isolates and a provincial collection of 1296 isolates, constructed the best-scoring maximum likelihood phylogenetic tree. Bayesian evolutionary analysis was used to calculate the latest common ancestor of lineages 2 and 4. The antigenic diversity of human T cell epitopes was evaluated by calculating the pairwise dN/dS ratios.

Results

Of the Zhejiang isolates, 964 (74.38%) belonged to lineage 2 and 332 (25.62%) belonged to lineage 4. The distributions of the sub-lineages varied across the geographic regions of Zhejiang Province. L2.2 is the most ancient sub-lineage in Zhejiang, first appearing approximately 6897 years ago (95% highest posterior density interval (HDI): 6513–7298). L4.4 is the most modern sub-lineage, first appearing approximately 2217 years ago (95% HDI: 1864–2581). The dN/dS ratios showed that the epitope and non-epitope regions of lineage 2 strains were significantly (P < 0.001) more conserved than those of lineage 4.

Conclusions

An increase in the frequency of lineage 4 may reflect its successful transmission over the last 20 years. The recent common ancestors of the sub-lineages and their transmission routes are relevant to the entry of humans into China and Zhejiang Province. Diversity in T cell epitopes may prevent Mycobacterium tuberculosis from being recognized by the immune system.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-021-00673-7.

A genome epidemiological study of mycobacterium tuberculosis in subpopulations with high and low incidence rate in Guangxi, South China

BACKGROUND

Tuberculosis (TB) is caused by a bacterium called Mycobacterium tuberculosis (Mtb). China is the third in top 8 high TB burden countries and Guangxi is one of the high incidence areas in South China. Determine bacterial factors that affected TB incidence rate is a step toward Ending the TB epidemic.

RESULTS

Genomes of M. tuberculosis cultures from a relatively high and low incidence region in Guangxi have been sequenced. 347 of 358(96.9%) were identified as M. tuberculosis. All the strains belong to Lineage 2 and Lineage 4, except for one in Lineage 1. We found that the genetic structure of the M. tuberculosis population in each county varies enormously. Low incidence rate regions have a lower prevalence of Beijing genotypes than other regions. Four isolates which harbored mutT4-48 also had mutT2-58 mutations. It is suggested that strains from the ancestors of modern Beijing lineage is circulating in Guangxi. Strains of modern Beijing lineage (OR=2.04) were more likely to acquire drug resistances than Lineage 4. Most of the lineage differentiation SNPs are related to cell wall biosynthetic pathways.

CONCLUSIONS

These results provided a higher resolution to better understand the history of transmission of M. tuberculosis from/to South China. And the incidence rate of tuberculosis might be affected by bacterial population structure shaped by demographic history. Our findings also support the hypothesis that Modern Beijing lineage originated in South China.

Molecular epidemiology and drug-resistance of tuberculosis in Luodian revealed by whole genome sequencing

In this study, we aimed to investigate the molecular epidemiology and drug-resistance profiles of tuberculosis (TB) in Luodian, an area with highest TB incidence and limited healthcare resources in Guizhou, China. The passive case finding strategy was used to identify suspected pulmonary TB with symptoms, and individuals with positive Mycobacterium tuberculosis (MTB) culture were enrolled from May 22, 2018 to April 21, 2019. All the 107 cases except three came from nine towns, including 55.1% from Longping and Bianyang. The phylogeny tree showed that 53.3% of strains were Lineage 2 (Beijing genotype), while 46.7% were Lineage 4 (Euro-American genotype). Among Lineage 2 strains, 66.7% were of "modern" Beijing type. Seven clusters with genomic distance within 12 SNPs were identified. The clusters included 14 strains, accounting for a clustering rate of 13.1%. The distance separating the clustered cases was between 2.1 and 71.0 km (Km), with an average paired distance of 21.8 Km (interquartile range, 2.8-38.0 Km). Based on the gene mutations associated with drug-resistance, we predicted that 4.8% of strains were resistant to isoniazid, 3.7% to rifampicin, and 3.7% to streptomycin; only one strain (0.9%) had multidrug resistance (MDR). This study found low drug-resistance rates in Luodian, and the sub-lineage of the "modern" Beijing branch has recent expansion in Luodian. This work may also serve as a genomic baseline to assess the evolution and spread of MTB in Guizhou.

Multiple Merger Genealogies in Outbreaks of Mycobacterium tuberculosis

The Kingman coalescent and its developments are often considered among the most important advances in population genetics of the last decades. Demographic inference based on coalescent theory has been used to reconstruct the population dynamics and evolutionary history of several species, including Mycobacterium tuberculosis (MTB), an important human pathogen causing tuberculosis. One key assumption of the Kingman coalescent is that the number of descendants of different individuals does not vary strongly, and violating this assumption could lead to severe biases caused by model misspecification. Individual lineages of MTB are expected to vary strongly in reproductive success because 1) MTB is potentially under constant selection due to the pressure of the host immune system and of antibiotic treatment, 2) MTB undergoes repeated population bottlenecks when it transmits from one host to the next, and 3) some hosts show much higher transmission rates compared with the average (superspreaders).

Here, we used an approximate Bayesian computation approach to test whether multiple-merger coalescents (MMC), a class of models that allow for large variation in reproductive success among lineages, are more appropriate models to study MTB populations. We considered 11 publicly available whole-genome sequence data sets sampled from local MTB populations and outbreaks and found that MMC had a better fit compared with the Kingman coalescent for 10 of the 11 data sets. These results indicate that the null model for analyzing MTB outbreaks should be reassessed and that past findings based on the Kingman coalescent need to be revisited.

Local adaptation of Mycobacterium tuberculosis on the Tibetan Plateau

During its global dispersal, Mycobacterium tuberculosis (Mtb) has encountered varied geographic environments and host populations. Although local adaptation seems to be a plausible model for describing long-term host-pathogen interactions, genetic evidence for this model is lacking. Here, we analyzed 576 whole-genome sequences of Mtb strains sampled from different regions of high-altitude Tibet. Our results show that, after sequential introduction of a few ancestral strains, the Tibetan Mtb population diversified locally while maintaining strict separation from the Mtb populations on the lower altitude plain regions of China. The current population structure and estimated past population dynamics suggest that the modern Beijing sublineage strains, which expanded over most of China and other global regions, did not show an expansion advantage in Tibet. The mutations in the Tibetan strains showed a higher proportion of A > G/T > C transitions than strains from the plain regions, and genes encoding DNA repair enzymes showed evidence of positive selection. Moreover, the long-term Tibetan exclusive selection for truncating mutations in the thiol-oxidoreductase encoding sseA gene suggests that Mtb was subjected to local selective pressures associated with oxidative stress. Collectively, the population genomics of Mtb strains in the relatively isolated population of Tibet provides genetic evidence that Mtb has adapted to local environments.

A population-based genomic epidemiological study of the source of tuberculosis infections in an emerging city: Shenzhen, China

Background

Tuberculosis (TB) in emerging cities is often a disease of recent immigrants, and understanding this epidemiology is crucial for designing effective control and prevention strategies.

Methods

We conducted a retrospective population-based genomic epidemiological study of culture-positive pulmonary TB patients diagnosed between June 2014 and November 2017 in the Bao'an District of Shenzhen, a Chinese city with dramatic recent growth. After whole genome sequencing, transmission clusters were defined as strains differing by no more than 12 SNPs.

Findings

Of 1696 culture-positive TB patients, 93•8% (1591/1696) were migrants, with 51•6% (821/1591) employed in housekeeping or unemployed. Of the 1460 migrants with known residence time, 47•7% (697/1460) developed TB within two years after arriving in Bao'an. Only 12•2% (207/1696) of Bao'an isolates were in genomic clusters, indicating that recent transmission was not the primary cause of TB in Bao'an. The isolates’ median terminal branch length was 56 SNPs, more than could have accumulated since the arrival of the migrants in Bao'an. The migrants’ isolates had genotypic distributions similar to those in their home provinces. One strain isolated in Bao'an belonged to a clade circulating in the patient's home province, providing further evidence that the strains were brought to Bao'an with the migrants.

Interpretation

TB in the Bao'an District is principally caused by reactivation of infections acquired by migrants in their home provinces. Nearly half developed TB within two years after arriving in Bao'an, suggesting a need for increased TB screening of migrants, especially housekeeping workers and the unemployed.

Funding

Sanming Project of Medicine in Shenzhen; National Science and Technology Major Project of China; Natural Science Foundation of China.

Genomic evidence supporting the clonal expansion of extensively drug-resistant tuberculosis bacteria belonging to a rare protoBeijing genotype

Tuberculosis disease (TB), caused by Mycobacterium tuberculosis, is a major public health issue in Thailand. The high prevalence of modern Beijing (Lineage 2.2.1) strains has been associated with multi- and extensively drug-resistant infections (MDR-, XDR-TB), complicating disease control. The impact of rarer proto-Beijing (L2.1) strains is less clear. In our study of thirty-seven L2.1 clinical isolates spanning thirteen years, we found a high prevalence of XDR-TB cases (32.4%). With ≤ 12 pairwise SNP distances, 43.2% of L2.1 patients belong to MDR-TB or XDR-TB transmission clusters suggesting a high level of clonal expansion across four Thai provinces. All XDR-TB (100%) were likely due to transmission rather than inadequate treatment. We found a 47 mutation signature and a partial deletion of the fadD14 gene in the circulating XDR-TB cluster, which can be used for surveillance of this rare and resilient M. tuberculosis strain-type that is causing increasing health burden. We also detected three novel deletion positions, a deletion of 1285 bp within desA3 (Rv3230c), large deletions in the plcB, plcA, and ppe38 gene which may play a role in the virulence, pathogenesis or evolution of the L2.1 strain-type.

Mutations in dnaA and a cryptic interaction site increase drug resistance in Mycobacterium tuberculosis

Genomic dissection of antibiotic resistance in bacterial pathogens has largely focused on genetic changes conferring growth above a single critical concentration of drug. However, reduced susceptibility to antibiotics-even below this breakpoint-is associated with poor treatment outcomes in the clinic, including in tuberculosis. Clinical strains of Mycobacterium tuberculosis exhibit extensive quantitative variation in antibiotic susceptibility but the genetic basis behind this spectrum of drug susceptibility remains ill-defined. Through a genome wide association study, we show that non-synonymous mutations in dnaA, which encodes an essential and highly conserved regulator of DNA replication, are associated with drug resistance in clinical M. tuberculosis strains. We demonstrate that these dnaA mutations specifically enhance M. tuberculosis survival during isoniazid treatment via reduced expression of katG, the activator of isoniazid. To identify DnaA interactors relevant to this phenotype, we perform the first genome-wide biochemical mapping of DnaA binding sites in mycobacteria which reveals a DnaA interaction site that is the target of recurrent mutation in clinical strains. Reconstructing clinically prevalent mutations in this DnaA interaction site reproduces the phenotypes of dnaA mutants, suggesting that clinical strains of M. tuberculosis have evolved mutations in a previously uncharacterized DnaA pathway that quantitatively increases resistance to the key first-line antibiotic isoniazid. Discovering genetic mechanisms that reduce drug susceptibility and support the evolution of high-level drug resistance will guide development of biomarkers capable of prospectively identifying patients at risk of treatment failure in the clinic.

The population genomics of within-host Mycobacterium tuberculosis

Recent progress in genomic sequencing from patient samples has allowed for the first detailed insight into the within-host genetic diversity of Mycobacterium tuberculosis (M.TB), revealing remarkably low levels of variation. While this has often been attributed to low mutation rates, other factors have been described, including resistance evolution (i.e., selective sweeps), widespread purifying and background selection, and, more recently, progeny skew. Here we review recent findings pertaining to the processes governing the evolutionary dynamics of M.TB, discuss their implications for improving our understanding of this important human pathogen, and make recommendations for future work. Significantly, this emerging evolutionary framework involving the joint estimation of demographic, selective, and reproductive processes is forming a new paradigm for the study of within-host pathogen evolution that will be widely applicable across organisms.

Insertion and deletion evolution reflects antibiotics selection pressure in a Mycobacterium tuberculosis outbreak

In genome evolution, genetic variants are the source of diversity, which natural selection acts upon. Treatment of human tuberculosis (TB) induces a strong selection pressure for the emergence of antibiotic resistance-conferring variants in the infecting Mycobacterium tuberculosis (MTB) strains. MTB evolution in response to treatment has been intensively studied and mainly attributed to point substitutions. However, the frequency and contribution of insertions and deletions (indels) to MTB genome evolution remains poorly understood. Here, we analyzed a multi-drug resistant MTB outbreak for the presence of high-quality indels and substitutions. We find that indels are significantly enriched in genes conferring antibiotic resistance. Furthermore, we show that indels are inherited during the outbreak and follow a molecular clock with an evolutionary rate of 5.37e-9 indels/site/year, which is 23 times lower than the substitution rate. Inherited indels may co-occur with substitutions in genes along related biological pathways; examples are iron storage and resistance to second-line antibiotics. This suggests that epistatic interactions between indels and substitutions affect antibiotic resistance and compensatory evolution in MTB.

Genomic variant-identification methods may alter Mycobacterium tuberculosis transmission inferences

Pathogen genomic data are increasingly used to characterize global and local transmission patterns of important human pathogens and to inform public health interventions. Yet, there is no current consensus on how to measure genomic variation. To test the effect of the variant-identification approach on transmission inferences for Mycobacterium tuberculosis, we conducted an experiment in which five genomic epidemiology groups applied variant-identification pipelines to the same outbreak sequence data. We compared the variants identified by each group in addition to transmission and phylogenetic inferences made with each variant set. To measure the performance of commonly used variant-identification tools, we simulated an outbreak. We compared the performance of three mapping algorithms, five variant callers and two variant filters in recovering true outbreak variants. Finally, we investigated the effect of applying increasingly stringent filters on transmission inferences and phylogenies. We found that variant-calling approaches used by different groups do not recover consistent sets of variants, which can lead to conflicting transmission inferences. Further, performance in recovering true variation varied widely across approaches. While no single variant-identification approach outperforms others in both recovering true genome-wide and outbreak-level variation, variant-identification algorithms calibrated upon real sequence data or that incorporate local reassembly outperform others in recovering true pairwise differences between isolates. The choice of variant filters contributed to extensive differences across pipelines, and applying increasingly stringent filters rapidly eroded the accuracy of transmission inferences and quality of phylogenies reconstructed from outbreak variation. Commonly used approaches to identify M. tuberculosis genomic variation have variable performance, particularly when predicting potential transmission links from pairwise genetic distances. Phylogenetic reconstruction may be improved by less stringent variant filtering. Approaches that improve variant identification in repetitive, hypervariable regions, such as long-read assemblies, may improve transmission inference.