Global expansion of Mycobacterium tuberculosis lineage 4 shaped by colonial migration and local adaptation

South Asian origin and global transmission history of Mycobacterium tuberculosis lineage 4

The origin and transmission history of Mycobacterium tuberculosis lineage 4 (MTB L4) has been extensively studied. However, different studies yielded different results; thus, the evolutionary history of MTB L4 still remains a subject of debate. Recently, a substantial amount of whole-genome sequencing data of MTB has become available, providing new data from diverse geographical locations worldwide. This study analyzed the most comprehensive global collection of MTB L4 genomes to date (n = 11,154), including new sequences from previously under-represented regions, to re-examine the bacterial evolutionary history. Our results suggest that the bacteria likely emerged and diversified in South Asia during the 12th and 13th centuries and then spread to various Old-World regions and the Americas between the 15th and 16th centuries, before undergoing extensive intercontinental transmissions starting from the 17th century onward. The effective population size of MTB L4 as a whole was estimated to expand steadily throughout its evolutionary history until the mid-20th century when a sharp drop occurred, coinciding with the introduction of antibiotics and significant improvements in human living conditions. Interestingly, this pattern was consistently observed across all major sub-lineages, indicating a broad impact of these factors on MTB L4 as a whole. Altogether, our results offer new refined insights into how anthropological changes might have shaped the bacterial evolutionary history and ultimately its global geographical distribution we observe today.IMPORTANCEContrary to previous studies, our analysis suggests that Mycobacterium tuberculosis lineage 4 (MTB L4) likely emerged and diversified in South Asia during the 12th and 13th centuries. It then spread to both the Americas and other Old-World regions between the 15th and 16th centuries, followed by extensive intercontinental transmissions beginning in the 17th century. These findings suggest that South Asia, one of the main crossroads of historical trade networks, might have played a pivotal role in promoting the early intercontinental spread of MTB L4. The effective population size of all major MTB L4 sub-lineages was inferred to increase steadily until the mid-20th century, after which a decline was observed. This decline coincides with the advent of antibiotics and improved living conditions, suggesting a wide impact of these factors on the entire bacterial population.

Ecology, global diversity and evolutionary mechanisms in the Mycobacterium tuberculosis complex

With the COVID-19 pandemic receding, tuberculosis (TB) is again the number one cause of human death to a single infectious agent. TB is caused by bacteria that belong to the Mycobacterium tuberculosis complex (MTBC). Recent advances in genome sequencing have provided new insights into the ecology and evolution of the MTBC. This includes the discovery of new phylogenetic lineages within the MTBC, a deeper understanding of the host tropism among the various animal-adapted lineages, enhanced knowledge on the evolutionary dynamics of antimicrobial resistance and transmission, as well as a better grasp of the within-host MTBC diversity. Moreover, advances in long-read sequencing are increasingly highlighting the relevance of structural genomic variation in the MTBC. These findings not only shed new light on the biology and epidemiology of TB, but also give rise to new questions and research avenues. The purpose of this Review is to summarize these new insights and discuss their implications for global TB control.

Phylogeographic History of Tomato Chlorosis Virus

Tomato chlorosis virus (ToCV), first reported in Florida, USA, in 1998, has since emerged in multiple regions worldwide, posing a significant threat to global tomato production. However, its origin, migration patterns, and evolutionary history remain poorly understood. In this study, we used Bayesian phylogeographic analysis of coat protein gene sequences from 155 ToCV isolates to reconstruct its phylogeographic history. Our results show that ToCV evolves at a rate of 6.24 × 10-4 subs/site/year (95% credibility interval: 4.35 × 10-4-8.28 × 10-4), with the most recent common ancestor dating back to 1882. The maximum clade credibility (MCC) tree revealed three major clades, with Clade 1-whose most recent common ancestor dates to approximately 1975-comprising over 90% of the isolates. Although the exact origin of ToCV remains uncertain, we identified five distinct migration pathways: one from Europe to the Americas, one from Europe to South Asia, one from the Middle East to East Asia, one from East Asia to mainland China, and one from mainland China to Europe. These findings underscore the complex global spread of ToCV and suggest that multiple geographic areas have contributed to its ongoing evolution and dissemination.

Nationwide phylogenomic surveillance of Mycobacterium tuberculosis in Mexico reveals pathogenic and drug resistant signatures of the prevailing L4 sublineage

BACKGROUND

Tuberculosis disease is a major global health concern. In Mexico, information regarding the genomic variants of Mycobacterium tuberculosis (MTB) prevailing in the country and the existence of specific biogeographical patterns remains extremely scarce.

OBJECTIVE

This study aimed to identify the genotypic patterns of MTB isolates in Mexico and determine the genes and specific single nucleotide polymorphisms involved in the evolution of these populations.

METHODS

Phylogenomic and pan-genomic analyses were performed using publicly available Mexican MTB genomes along with 33 newly sequenced genomes from Jalisco, considering a global context.

RESULTS

The L4 sublineages of MTB, such as L4.1.1 (X), L4.1.2 (H), and L4.3 (LAM), were the most prevalent in Mexico. We found exclusive mutations and gene clusters in a virulent sublineage L4.1.1.3 (X3), which is endemic to Mexico. These genes encoded three PE/PPE family proteins: a multidrug transporter, thioredoxin domain-containing protein, quinone-dependent l-lactate dehydrogenase, DUF1725 domain-containing protein, amidase, poly (A) polymerase, and six hypothetical/uncharacterised proteins. Additionally, the genes encode an ESX-1 secretion-associated protein and a deazaflavin-dependent nitroreductase (ddn).

CONCLUSION

X3 was distinguished from the rest of the sublineages by containing genes related to pathogenicity and virulence, as well as a gene linked to delamanid, an antibiotic for active multidrug-resistant tuberculosis. These findings provide valuable insight into the circulation and spread of MTB in Mexico.

Global epidemiology of Mycobacterium tuberculosis lineage 4 insights from Ecuadorian genomic data

Tuberculosis is a global public health concern, and understanding Mycobacterium tuberculosis transmission routes and genetic diversity of M. tuberculosis is crucial for outbreak control. This study aimed to explore the genomic epidemiology and genetic diversity of M. tuberculosis in Ecuador by analyzing 88 local isolates and 415 public genomes from 19 countries within the Euro-American lineage (L4). Our results revealed significant genomic diversity among the isolates, particularly in the genes related to protein processing, carbohydrate metabolism, lipid metabolism, and xenobiotic biodegradation and metabolism. The population structure analysis showed that sub-lineages 4.3.2/3 (35.4%), 4.1.2.1 (22.7%), 4.4.1 (12.7%), and 4.1.1. (10.7%) were the most prevalent. Phylogenetic and transmission network analyses suggest that these isolates circulating within Ecuador share genetic ties with isolates from other continents, implying historical and ongoing intercontinental transmission events. Our findings underscore the importance of integrating genomic data into public health strategies for tuberculosis control and suggest that enhanced genomic surveillance is essential for understanding and mitigating the global spread of M. tuberculosis. This study provides a comprehensive genomic framework for future epidemiological investigations and control measures targeting M. tuberculosis L4 in Ecuador.

Phylogeography and Evolutionary Dynamics of Tobacco Curly Shoot Virus

Tobacco curly shoot virus (TbCSV), a begomovirus, causes significant economic losses in tobacco and tomato crops across East, Southeast, and South Asia. Despite its agricultural importance, the evolutionary dynamics and emergence process of TbCSV remain poorly understood. This study analyzed the phylodynamics of TbCSV by examining its nucleotide sequences of the coat protein (CP) gene collected between 2000 and 2022. Using various combinations of priors, Bayes factor comparisons identified heterochronous datasets (3 × 100 million chains) generated from a strict molecular clock and Bayesian skyline tree priors as the most robust. The mean substitution rate of the CP gene was estimated at 6.50 × 10-4 substitutions/site/year (95% credibility interval: 4.74 × 10-4-8.50 × 10-4). TbCSV was inferred to have diverged around 1920 CE (95% credibility interval: 1887-1952), with its most probable origin in South Asia. These findings provide valuable insights for the phylogeography and evolutionary dynamics of TbCSV, and contribute to a broader understanding of begomovirus epidemiology.

Phylodynamic analysis reveals disparate transmission dynamics of Mycobacterium tuberculosis-complex lineages in Botswana

Tuberculosis epidemics have traditionally been conceptualized as arising from a single uniform pathogen. However, Mycobacterium tuberculosis-complex (Mtbc), the pathogen causing tuberculosis in humans, encompasses multiple lineages exhibiting genetic and phenotypic diversity that may be responsible for heterogeneity in TB transmission. We analysed a population-based dataset of 1,354 Mtbc whole-genome sequences collected over four years in Botswana, a country with high HIV and tuberculosis burden. We identified Lineage 4 (L4) as the most prevalent (87.4%), followed by L1 (6.4%), L2 (5.3%), and L3 (0.9%). Within L4, multiple sublineages were identified, with L4.3.4 being the predominant sublineage. Phylodynamic analysis revealed L4.3.4 expanded steadily from late 1800s to early 2000s. Conversely, L1, L4.4, and L4.3.2 showed population trajectories closely aligned with the HIV epidemic. Meanwhile, L2 saw rapid expansion throughout most of the 20th century but declined sharply in early 1990s. Additionally, pairwise genome comparison of Mtbc highlighted differences in clustering proportions due to recent transmission at the sublineage level. These findings emphasize the diverse transmission dynamics of strains of different Mtbc lineages and highlight the potential for phylodynamic analysis of routine sequences to refine our understanding of lineage-specific behaviors.

A Genome-Focused Investigation Reveals the Emergence of a Mycobacterium tuberculosis Strain Related to Multidrug-Resistant Tuberculosis in the Amazon Region of Brazil

A previous study in Pará, Northern Brazil, described a strain of Mycobacterium tuberculosis with a unique genotype (SIT2517/T1) associated with multidrug-resistant tuberculosis (MDR-TB). To improve our understanding of MDR-TB transmission dynamics of these strains within this region, we performed phenotypic and genotypic drug susceptibility testing (pDST/gDST), 24-loci mycobacterial interspersed repetitive units (MIRU-VNTR) genotyping, whole-genome sequencing (WGS) and geo-epidemiology analysis. Of the 28 SIT2517/T1 isolates, 19 (67.9%) could be genotyped by 24-loci MIRU-VNTR and 15 by WGS. All belonged to sublineage 4.1.1.3, distinct from other representative Lineage 4 isolates identified in Brazil. The MDR phenotype determined by pDST was confirmed by gDST, the latter also demonstrating the presence of additional mutations conferring pre-extensively drug-resistance (pre-XDR). Discrepancies between gDST and pDST were observed for pyrazinamide and fluoroquinolones. Thirteen out of 15 isolates analyzed by WGS were clustered when applying a 12 single nucleotide polymorphisms (SNPs) cutoff. The SIT2517/T1 isolates were distributed across the metropolitan regions of Belém and Collares municipalities, showing no geographic clustering. WGS-transmission network analysis revealed a high likelihood of direct transmission and the formation of two closely linked transmission chains. This study highlights the need to implement TB genomic surveillance in the Brazilian Amazon region.

Delineating the evolutionary pathway to multidrug-resistant outbreaks of a Mycobacterium tuberculosis L4.1.2.1/Haarlem sublineage

OBJECTIVES

We sought to capture the evolutionary itinerary of the Mycobacterium tuberculosis L4.1.2.1/Haarlem sublineage in northern Tunisia, where it caused a major multidrug-resistant (MDR) tuberculosis outbreak in a context strictly negative for HIV infection.

METHODS

We combined whole genome sequencing and Bayesian approaches using a representative collection of drug-susceptible and drug-resistant L4.1.2.1/Haarlem clinical strains (n = 121) recovered from the outbreak region over 16 years.

RESULTS

In the absence of drug resistance, the L4.1.2.1/Haarlem sublineage showed a propensity for rapid transmission as witnessed by the high clustering (44.6%) and recent transmission rates (25%), as well as the reduced mean distance between genome pairs. The entire pool of L4.1.2.1/Haarlem MDR strains was found to be linked to either the aforementioned major outbreak (68 individuals, 2001-2016) or to a minor, newly uncovered outbreak (six cases, 2001-2011). Strikingly, the two outbreaks descended independently from a common ancestor that can be dated back to 1886.

CONCLUSIONS

Our data point to the intrinsic propensity for rapid transmission of the M. tuberculosis L4.1.2.1/Haarlem sublineage in northern Tunisia, linking the overall MDR tuberculosis epidemic to a single ancestor. These findings bring out the important role of the bacillus' genetic background in the emergence of successful MDR M. tuberculosis clones.

PE/PPE mutations in the transmission of Mycobacterium tuberculosis in China revealed by whole genome sequencing

OBJECTIVE

This study aims to examine the impact of PE/PPE gene mutations on the transmission of Mycobacterium tuberculosis (M. tuberculosis) in China.

METHODS

We collected the whole genome sequencing (WGS) data of 3202 M. tuberculosis isolates in China from 2007 to 2018 and investigated the clustering of strains from different lineages. To evaluate the potential role of PE/PPE gene mutations in the dissemination of the pathogen, we employed homoplastic analysis to detect homoplastic single nucleotide polymorphisms (SNPs) within these gene regions. Subsequently, logistic regression analysis was conducted to analyze the statistical association.

RESULTS

Based on nationwide M. tuberculosis WGS data, it has been observed that the majority of the M. tuberculosis burden in China is caused by lineage 2 strains, followed by lineage 4. Lineage 2 exhibited a higher number of transmission clusters, totaling 446 clusters, of which 77 were cross-regional clusters. Conversely, there were only 52 transmission clusters in lineage 4, of which 9 were cross-regional clusters. In the analysis of lineage 2 isolates, regression results showed that 4 specific gene mutations, PE4 (position 190,394; c.46G > A), PE_PGRS10 (839,194; c.744 A > G), PE16 (1,607,005; c.620T > G) and PE_PGRS44 (2,921,883; c.333 C > A), were significantly associated with the transmission of M. tuberculosis. Mutations of PE_PGRS10 (839,334; c.884 A > G), PE_PGRS11 (847,613; c.1455G > C), PE_PGRS47 (3,054,724; c.811 A > G) and PPE66 (4,189,930; c.303G > C) exhibited significant associations with the cross-regional clusters. A total of 13 mutation positions showed a positive correlation with clustering size, indicating a positive association. For lineage 4 strains, no mutations were found to enhance transmission, but 2 mutation sites were identified as risk factors for cross-regional clusters. These included PE_PGRS4 (338,100; c.974 A > G) and PPE13 (976,897; c.1307 A > C).

CONCLUSION

Our results indicate that some PE/PPE gene mutations can increase the risk of M. tuberculosis transmission, which might provide a basis for controlling the spread of tuberculosis.

Sustained transmission over two decades of a previously unrecognised MPT64 negative Mycobacterium tuberculosis strain in Queensland, Australia: a whole genome sequencing study

Background

MPT64 is a key protein used for Mycobacterium tuberculosis (MTB) complex strain identification. We describe protracted transmission of an MPT64 negative MTB strain in Queensland, Australia, and explore genomic factors related to its successful spread.

Methods

All MPT64 negative strains identified between 2002 and 2022 by the Queensland Mycobacteria Reference Laboratory, and an additional 2 isolates from New South Wales (NSW), were whole genome sequenced. Bayesian modelling and phylogeographical analyses were used to assess their evolutionary history and transmission dynamics. Protein structural modelling to understand the putative functional effects of the mutated gene coding for MPT64 protein was performed.

Findings

Forty-three MPT64 negative isolates were sequenced, belonging to a single MTB cluster of Lineage 4.1.1.1 strains. Combined with a UK dataset of the same lineage, molecular dating estimated 1990 (95% HPD 1987-1993) as the likely time of strain introduction into Australia. Although the strain has spread over a wide geographic area and new cases linked to the cluster continue to arise, phylodynamic analysis suggest the outbreak peaked around 2003. All MPT64 negative strains had a frame shift mutation (delAT, p.Val216fs) within the MPT64 gene, which confers two major structural rearrangements at the C-terminus of the protein.

Interpretation

This study uncovered the origins of an MPT64 negative MTB outbreak in Australia, providing a richer understanding of its biology and transmission dynamics, as well as guidance for clinical diagnosis and public health action. The potential spread of MPT64 negative strains undermines the diagnostic utility of the MPT64 immunochromatographic test.

Funding

This study was funded from an operational budget provided to the Queensland Mycobacterium Reference Laboratory by Pathology Queensland, Queensland Department of Health.

HIV co-infection is associated with reduced Mycobacterium tuberculosis transmissibility in sub-Saharan Africa

Persons living with HIV are known to be at increased risk of developing tuberculosis (TB) disease upon infection with Mycobacterium tuberculosis (Mtb). However, it has remained unclear how HIV co-infection affects subsequent Mtb transmission from these patients. Here, we customized a Bayesian phylodynamic framework to estimate the effects of HIV co-infection on the Mtb transmission dynamics from sequence data. We applied our model to four Mtb genomic datasets collected in sub-Saharan African countries with a generalized HIV epidemic. Our results confirm that HIV co-infection is a strong risk factor for developing active TB. Additionally, we demonstrate that HIV co-infection is associated with a reduced effective reproductive number for TB. Stratifying the population by CD4+ T-cell count yielded similar results, suggesting that, in this context, CD4+ T-cell count is not a better predictor of Mtb transmissibility than HIV infection status alone. Together, our genome-based analyses complement observational household contact studies, and more firmly establish the negative association between HIV co-infection and Mtb transmissibility.

Prevalence, drug resistance, and genotypic diversity of the RDRio subfamily of Mycobacterium tuberculosis in Ecuador: a retrospective analysis for years 2012-2016

Introduction

A major sublineage within the Mycobacterium tuberculosis (MTB) LAM family characterized by a new in-frame fusion gene Rv3346c/55c was discovered in Rio de Janeiro (Brazil) in 2007, called RDRio, associated to drug resistance. The few studies about prevalence of MTB RDRio strains in Latin America reported values ranging from 3% in Chile to 69.8% in Venezuela, although no information is available for countries like Ecuador.

Methods

A total of 814 MTB isolates from years 2012 to 2016 were screened by multiplex PCR for RDRio identification, followed by 24-loci MIRU-VNTR and spoligotyping.

Results

A total number of 17 MTB RDRio strains were identified, representing an overall prevalence of 2.09% among MTB strains in Ecuador. While 10.9% of the MTB isolates included in the study were multidrug resistance (MDR), 29.4% (5/17) of the RDRio strains were MDR.

Discussion

This is the first report of the prevalence of MTB RDRio in Ecuador, where a strong association with MDR was found, but also a very low prevalence compared to other countries in Latin America. It is important to improve molecular epidemiology tools as a part of MTB surveillance programs in Latin America to track the transmission of potentially dangerous MTB stains associated to MDR TB like MTB RDRio.

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.

Evolution and emergence of Mycobacterium tuberculosis

Tuberculosis (TB) remains one of the deadliest infectious diseases in human history, prevailing even in the 21st century. The causative agents of TB are represented by a group of closely related bacteria belonging to the Mycobacterium tuberculosis complex (MTBC), which can be subdivided into several lineages of human- and animal-adapted strains, thought to have shared a last common ancestor emerged by clonal expansion from a pool of recombinogenic Mycobacterium canettii-like tubercle bacilli. A better understanding of how MTBC populations evolved from less virulent mycobacteria may allow for discovering improved TB control strategies and future epidemiologic trends. In this review, we highlight new insights into the evolution of mycobacteria at the genus level, describing different milestones in the evolution of mycobacteria, with a focus on the genomic events that have likely enabled the emergence and the dominance of the MTBC. We also review the recent literature describing the various MTBC lineages and highlight their particularities and differences with a focus on host preferences and geographic distribution. Finally, we discuss on putative mechanisms driving the evolution of tubercle bacilli and mycobacteria in general, by taking the mycobacteria-specific distributive conjugal transfer as an example.

Detection of a historic reservoir of bedaquiline/clofazimine resistance-associated variants in Mycobacterium tuberculosis

BACKGROUND

Drug resistance in tuberculosis (TB) poses a major ongoing challenge to public health. The recent inclusion of bedaquiline into TB drug regimens has improved treatment outcomes, but this advance is threatened by the emergence of strains of Mycobacterium tuberculosis (Mtb) resistant to bedaquiline. Clinical bedaquiline resistance is most frequently conferred by off-target resistance-associated variants (RAVs) in the mmpR5 gene (Rv0678), the regulator of an efflux pump, which can also confer cross-resistance to clofazimine, another TB drug.

METHODS

We compiled a dataset of 3682 Mtb genomes, including 180 carrying variants in mmpR5, and its immediate background (i.e. mmpR5 promoter and adjacent mmpL5 gene), that have been associated to borderline (henceforth intermediate) or confirmed resistance to bedaquiline. We characterised the occurrence of all nonsynonymous mutations in mmpR5 in this dataset and estimated, using time-resolved phylogenetic methods, the age of their emergence.

RESULTS

We identified eight cases where RAVs were present in the genomes of strains collected prior to the use of bedaquiline in TB treatment regimes. Phylogenetic reconstruction points to multiple emergence events and circulation of RAVs in mmpR5, some estimated to predate the introduction of bedaquiline. However, epistatic interactions can complicate bedaquiline drug-susceptibility prediction from genetic sequence data. Indeed, in one clade, Ile67fs (a RAV when considered in isolation) was estimated to have emerged prior to the antibiotic era, together with a resistance reverting mmpL5 mutation.

CONCLUSIONS

The presence of a pre-existing reservoir of Mtb strains carrying bedaquiline RAVs prior to its clinical use augments the need for rapid drug susceptibility testing and individualised regimen selection to safeguard the use of bedaquiline in TB care and control.

Syringaldehyde Exhibits Antibacterial and Antioxidant Activities against Mycobacterium marinum Infection

Tuberculosis (TB) is caused by infection with Mycobacterium tuberculosis (Mtb), which has a unique resistance to many antimicrobial agents. TB has emerged as a significant worldwide health issue because of the rise of multidrug-resistant strains causing drug-resistant TB (DR-TB). As a result, the development of new drugs or effective strategies is crucial for patients with TB. Mycobacterium marinum (Mm) and Mtb are both species of mycobacteria. In zebrafish, Mm proliferates and forms chronic granulomatous infections, which are similar to Mtb infections in lung tissue. Syringaldehyde (SA) is a member of the phenolic aldehyde family found in various plants. Here, we investigated its antioxidative and antibacterial properties in Mm-infected cells and zebrafish. Our results demonstrated that SA inhibits Mm-infected pulmonary epithelial cells and inhibits the proliferation of Mm in Mm-infected zebrafish, suggesting that SA provides an antibacterial effect during Mm infection. Further study demonstrated that supplementation with SA inhibits the production of malondialdehyde (MDA) and reactive oxygen species (ROS) and increases the levels of reduced glutathione (GSH) in Mm-infection-induced macrophages. SA inhibits the levels of MDA in Mm-infected zebrafish, suggesting that SA exerts antioxidative effects in vivo. Additionally, we found that SA promotes the expression of NRF2/HO-1/NQO-1 and the activation of the AMPK-α1/AKT/GSK-3β signaling pathway. In summary, our data demonstrated that SA exerts antioxidative and antibacterial effects during Mm infection both in vivo and in vitro and that the antioxidative effects of SA may be due to the regulation of NRF2/HO-1/NQO-1 and the AMPK-α1/AKT/GSK-3β signaling pathway.

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.

From modern-day parasitology to paleoparasitology: the elusive past record and evolution of Cryptosporidium

Recent efforts have been made to review the state of the art on a variety of questions and targets in paleoparasitology, including protozoan taxa. Meanwhile, these efforts seemed to let aside Cryptosporidium, and we then intended to review its paleoparasitological record to assess its past distribution and favored detection methods, and eventually highlight needed research trajectories. This review shows that contrary to other parasites, most of the positive results came from South-American sites and coprolites rather than sediment samples, highlighting the need to test this kind of material, notably in Europe where many negative results were reported in the published literature from sediment samples. Moreover, aDNA-based detections are nearly absent from the paleoparasitological record of this parasite, though punctually shown successful. With their potential to address the evolutionary history of Cryptosporidium species, notably through their 18S rRNA tree, aDNA-based approaches should be encouraged in the future. In sum, and though the limits of currently used methods and materials remain unclear, this review highlights the potential role of coprolites and aDNA for the study of Cryptosporidium species in the past and how this history shaped their current diversity and distribution, notably among human populations but also farm animals.

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.

Genomic epidemiology analysis of drug-resistant Mycobacterium tuberculosis distributed in Mexico

Genomics has significantly revolutionized pathogen surveillance, particularly in epidemiological studies, the detection of drug-resistant strains, and disease control. Despite its potential, the representation of Latin American countries in the genomic catalogues of Mycobacterium tuberculosis (Mtb), the bacteria responsible for Tuberculosis (TB), remains limited. In this study, we present a whole genome sequencing (WGS)-based analysis of 85 Mtb clinical strains from 17 Mexican states, providing insights into local adaptations and drug resistance signatures in the region. Our results reveal that the Euro-American lineage (L4) accounts for 94% of our dataset, showing 4.1.2.1 (Haarlem, n = 32), and 4.1.1.3 (X-type, n = 34) sublineages as the most prevalent. We report the presence of the 4.1.1.3 sublineage, which is endemic to Mexico, in six additional locations beyond previous reports. Phenotypic drug resistance tests showed that 34 out of 85 Mtb samples were resistant, exhibiting a variety of resistance profiles to the first-line antibiotics tested. We observed high levels of discrepancy between phenotype and genotype associated with drug resistance in our dataset, including pyrazinamide-monoresistant Mtb strains lacking canonical variants of drug resistance. Expanding the Latin American Mtb genome databases will enhance our understanding of TB epidemiology and potentially provide new avenues for controlling the disease in the region.

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.

The Mycobacterium tuberculosis genome at 25 years: lessons and lingering questions

First achieved in 1998 by Cole et al., the complete genome sequence of Mycobacterium tuberculosis continues to provide an invaluable resource to understand tuberculosis (TB), the leading cause of global infectious disease mortality. At the 25-year anniversary of this accomplishment, we describe how insights gleaned from the M. tuberculosis genome have led to vital tools for TB research, epidemiology, and clinical practice. The increasing accessibility of whole-genome sequencing across research and clinical settings has improved our ability to predict antibacterial susceptibility, to track epidemics at the level of individual outbreaks and wider historical trends, to query the efficacy of the bacille Calmette-Guérin (BCG) vaccine, and to uncover targets for novel antitubercular therapeutics. Likewise, we discuss several recent efforts to extract further discoveries from this powerful resource.

Investigating the composition and recruitment of the mycobacterial ImuA'-ImuB-DnaE2 mutasome

A DNA damage-inducible mutagenic gene cassette has been implicated in the emergence of drug resistance in Mycobacterium tuberculosis during anti-tuberculosis (TB) chemotherapy. However, the molecular composition and operation of the encoded 'mycobacterial mutasome' - minimally comprising DnaE2 polymerase and ImuA' and ImuB accessory proteins - remain elusive. Following exposure of mycobacteria to DNA damaging agents, we observe that DnaE2 and ImuB co-localize with the DNA polymerase III β subunit (β clamp) in distinct intracellular foci. Notably, genetic inactivation of the mutasome in an imuBAAAAGG mutant containing a disrupted β clamp-binding motif abolishes ImuB-β clamp focus formation, a phenotype recapitulated pharmacologically by treating bacilli with griselimycin and in biochemical assays in which this β clamp-binding antibiotic collapses pre-formed ImuB-β clamp complexes. These observations establish the essentiality of the ImuB-β clamp interaction for mutagenic DNA repair in mycobacteria, identifying the mutasome as target for adjunctive therapeutics designed to protect anti-TB drugs against emerging resistance.

Towards a post-pandemic future for global pathogen genome sequencing

Pathogen genome sequencing has become a routine part of our response to active outbreaks of infectious disease and should be an important part of our preparations for future epidemics. In this Essay, we discuss the innovations that have enabled routine pathogen genome sequencing, as well as how genome sequences can be used to understand and control the spread of infectious disease. We also explore the impact of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) pandemic on the field of pathogen genomics and outline the challenges we must address to further improve the utility of pathogen genome sequencing in the future.

Phase variation as a major mechanism of adaptation in Mycobacterium tuberculosis complex

Phase variation induced by insertions and deletions (INDELs) in genomic homopolymeric tracts (HT) can silence and regulate genes in pathogenic bacteria, but this process is not characterized in MTBC (Mycobacterium tuberculosis complex) adaptation. We leverage 31,428 diverse clinical isolates to identify genomic regions including phase-variants under positive selection. Of 87,651 INDEL events that emerge repeatedly across the phylogeny, 12.4% are phase-variants within HTs (0.02% of the genome by length). We estimated the in-vitro frameshift rate in a neutral HT at 100× the neutral substitution rate at [Formula: see text] frameshifts/HT/year. Using neutral evolution simulations, we identified 4,098 substitutions and 45 phase-variants to be putatively adaptive to MTBC (P < 0.002). We experimentally confirm that a putatively adaptive phase-variant alters the expression of espA, a critical mediator of ESX-1-dependent virulence. Our evidence supports the hypothesis that phase variation in the ESX-1 system of MTBC can act as a toggle between antigenicity and survival in the host.

Genomic changes underpinning the emergence of a successful Mycobacterium tuberculosis Latin American and Mediterranean clonal complex

Introduction

The Latin American and Mediterranean sublineage (L4.3/LAM) is the most common generalist sublineage of Mycobacterium tuberculosis lineage 4 (L4), yet certain L4.3/LAM genotypes appear to be confined to particular geographic regions. This is typically the case of a L4.3/LAM clonal complex (CC), TUN4.3_CC1, which is the most preponderant in Tunisia (61.5% of L4.3/LAM).

Methods

Here, we used whole-genome sequencing data of 346 globally distributed L4 clinical strains, including 278 L4.3/LAM isolates, to reconstruct the evolutionary history of TUN4.3_CC1 and delineate critical genomic changes underpinning its success.

Results and Discussion

Phylogenomic coupled to phylogeographic analyses indicated that TUN4.3_CC1 has evolved locally, being confined mainly to North Africa. Maximum likelihood analyses using the site and branch-site models of the PAML package disclosed strong evidence of positive selection in the gene category "cell wall and cell processes" of TUN4.3_CC1. Collectively, the data indicate that TUN4.3_CC1 has inherited several mutations, which could have potentially contributed to its evolutionary success. Of particular interest are amino acid replacements at the esxK and eccC2 genes of the ESX/Type VII secretion system, which were found to be specific to TUN4.3_CC1, being common to almost all isolates. Because of its homoplastic nature, the esxK mutation could potentially have endowed TUN4.3_CC1 with a selective advantage. Moreover, we noticed the occurrence of additional, previously described homoplasic nonsense mutations in ponA1 and Rv0197. The mutation in the latter gene, a putative oxido-reductase, has previously been shown to be correlated with enhanced transmissibility in vivo. In sum, our findings unveiled several features underpinning the success of a locally evolved L4.3/LAM clonal complex, lending further support to the critical role of genes encoded by the ESX/type VII secretion system.

Transmission Dynamics of a Mycobacterium tuberculosis Complex Outbreak in an Indigenous Population in the Colombian Amazon Region

Whole genome sequencing (WGS) has become the main tool for studying the transmission of Mycobacterium tuberculosis complex (MTBC) strains; however, the clonal expansion of one strain often limits its application in local MTBC outbreaks. The use of an alternative reference genome and the inclusion of repetitive regions in the analysis could potentially increase the resolution, but the added value has not yet been defined. Here, we leveraged short and long WGS read data of a previously reported MTBC outbreak in the Colombian Amazon Region to analyze possible transmission chains among 74 patients in the indigenous setting of Puerto Nariño (March to October 2016). In total, 90.5% (67/74) of the patients were infected with one distinct MTBC strain belonging to lineage 4.3.3. Employing a reference genome from an outbreak strain and highly confident single nucleotide polymorphisms (SNPs) in repetitive genomic regions, e.g., the proline-glutamic acid/proline-proline-glutamic-acid (PE/PPE) gene family, increased the phylogenetic resolution compared to a classical H37Rv reference mapping approach. Specifically, the number of differentiating SNPs increased from 890 to 1,094, which resulted in a more granular transmission network as judged by an increasing number of individual nodes in a maximum parsimony tree, i.e., 5 versus 9 nodes. We also found in 29.9% (20/67) of the outbreak isolates, heterogenous alleles at phylogenetically informative sites, suggesting that these patients are infected with more than one clone. In conclusion, customized SNP calling thresholds and employment of a local reference genome for a mapping approach can improve the phylogenetic resolution in highly clonal MTBC populations and help elucidate within-host MTBC diversity. IMPORTANCE The Colombian Amazon around Puerto Nariño has a high tuberculosis burden with a prevalence of 1,267/100,000 people in 2016. Recently, an outbreak of Mycobacterium tuberculosis complex (MTBC) bacteria among the indigenous populations was identified with classical MTBC genotyping methods. Here, we employed a whole-genome sequencing-based outbreak investigation in order to improve the phylogenetic resolution and gain new insights into the transmission dynamics in this remote Colombian Amazon Region. The inclusion of well-supported single nucleotide polymorphisms in repetitive regions and a de novo-assembled local reference genome provided a more granular picture of the circulating outbreak strain and revealed new transmission chains. Multiple patients from different settlements were possibly infected with at least two different clones in this high-incidence setting. Thus, our results have the potential to improve molecular surveillance studies in other high-burden settings, especially regions with few clonal multidrug-resistant (MDR) MTBC lineages/clades.

Large genomic deletions delineate Mycobacterium tuberculosis L4 sublineages in South American countries

Mycobacterium tuberculosis (Mtb) is still one of the primary pathogens of humans causing tuberculosis (TB) disease. Mtb embraces nine well-defined phylogenetic lineages with biological and geographical disparities. The lineage L4 is the most globally widespread of all lineages and was introduced to America with European colonization. Taking advantage of many genome projects available in public repositories, we undertake an evolutionary and comparative genomic analysis of 522 L4 Latin American Mtb genomes. Initially, we performed careful quality control of public read datasets and applied several thresholds to filter out low-quality data. Using a genome de novo assembly strategy and phylogenomic methods, we spotted novel south American clades that have not been revealed yet. Additionally, we describe genomic deletion profiles of these strains from an evolutionary perspective and report Mycobacterium tuberculosis L4 sublineages signature-like gene deletions, some of the novel. One is a specific deletion of 6.5 kbp that is only present in sublineage 4.1.2.1. This deletion affects a complex group of 10 genes with putative products annotated, among others, as a lipoprotein, transmembrane protein, and toxin/antitoxin system proteins. The second novel deletion spans for 4.9 kbp and specific of a particular clade of the 4.8 sublineage and affects 7 genes. The last novel deletion affects 4 genes, extends for 4.8 kbp., and is specific to some strains within the 4.1.2.1 sublineage that are present in Colombia, Peru and Brasil.

Evolution and epidemic success of Mycobacterium tuberculosis in eastern China: evidence from a prospective study

BACKGROUND

Lineage distribution of Mycobacterium tuberculosis (Mtb) isolates is strongly associated with geographically distinct human populations, and its transmission can be further impacted by the bacterial genome. However, the epidemic success of Mtb isolates at an individual level was unknown in eastern China. Knowledge regarding the emergence and transmission of Mtb isolates as well as relevant factors may offer a new solution to curb the spread of the disease. Thus, this study aims to reveal the evolution and epidemic success of Mtb isolates in eastern China.

RESULTS

Of initial 1040 isolates, 997 were retained after removing duplicates and those with insufficient sequencing depth. Of the final samples, 733 (73.52%) were from Zhejiang Province, and 264 (26.48%) were from Shanghai City. Lineage 2 and lineage 4 accounted for 80.44% and 19.56%, with common ancestors dating around 7017 years ago and 6882 years ago, respectively. Sub-lineage L2.2 (80.34%) contributed the majority of total isolates, followed by L4.4 (8.93%) and L4.5 (8.43%). Additionally, 51 (5.12%) isolates were identified to be multidrug-resistant (MDR), of which 21 (29.17%) were pre-extensively drug-resistant (pre-XDR). One clade harboring katG S315T mutation may date back to 65 years ago and subsequently acquired mutations conferring resistance to another five antibiotic drugs. The prevalence of compensatory mutation was the highest in pre-XDR isolates (76.19%), followed by MDR isolates (47.06%) and other drug-resistant isolates (20.60%). Time-scaled haplotypic density analyses suggested comparable success indices between lineage 2 and lineage 4 (P = 0.306), and drug resistance did not significantly promote the transmission of Mtb isolates (P = 0.340). But for pre-XDR isolates, we found a higher success index in those with compensatory mutations (P = 0.025). Mutations under positive selection were found in genes associated with resistance to second-line injectables (whiB6) and drug tolerance (prpR) in both lineage 2 and lineage 4.

CONCLUSIONS

Our study demonstrates the population expansion of lineage 2 and lineage 4 in eastern China, with comparable transmission capacity, while accumulation of resistance mutations does not necessarily facilitate the success of Mtb isolates. Compensatory mutations usually accompany drug resistance and significantly contribute to the epidemiological transmission of pre-XDR strains. Prospective molecular surveillance is required to further monitor the emergence and spread of pre-XDR/XDR strains in eastern China.

Strain diversity and gene mutations associated with presumptive multidrug-resistant Mycobacterium tuberculosis complex isolates in Northwest Ethiopia

OBJECTIVES

In this study, we assessed the genetic diversity and gene mutations that confer resistance to rifampicin (RIF), isoniazid (INH), fluoroquinolone (FQ), and second-line injectable (SLI) drugs in RIF-resistant (RR)/multidrug-resistant tuberculosis (MDR-TB) isolates in Northwest Ethiopia.

METHODS

Spoligotyping was used to assign isolates to TB lineages (Ls), and Hain line probe assays were used to detect resistance to RIF, INH, and FQs, and SLIs.

RESULTS

Among 130 analyzed strains, 68.5% were RR, and four major Mycobacterium tuberculosis complex lineages (L1, L3, L4, and L7) were identified with a predominance of the Euro-American L4 (72, 54.7%), while L7 genotypes were less common (3, 2.3%). Overall, the L4-T3-ETH (41, 32.0%), L3-CAS1-Delhi (29, 22.7%), and L3-CAS1-Killi (19, 14.8%) families were most common. Line probe analysis showed that among rpoB mutants, 65.2% were S450L, while 87.8% of katG mutants were S315T. Only three isolates showed mutation (c-15t) at the inhA gene, and no double mutation with katG and inhA genes was found. Six strains, two each of L1, L3, and L4, were resistant to FQs, having gyrA mutations (D94G, S91P), of which three isolates had additional resistance to SLI (rrs A1401G or C1402T mutations) including one isolate with low-level kanamycin (KAN) resistance.

CONCLUSIONS

This study showed a predominance of L4-T3-ETH, L3-CAS1-Delhi, and L3-CAS1-Killi families, with a high rate of rpoB_S450L and katG_S315T mutations and a low proportion of gyrA and rrs mutations. L7 was less frequently observed in this study. Further investigations are, therefore, needed to understand L7 and other lineages with undefined mutations.

Genomic evidence for adaptation to tuberculosis in the Andes before European contact

Most studies focusing on human high-altitude adaptation in the Andean highlands have thus far been focused on Peruvian populations. We present high-coverage whole genomes from Indigenous people living in the Ecuadorian highlands and perform multi-method scans to detect positive natural selection. We identified regions of the genome that show signals of strong selection to both cardiovascular and hypoxia pathways, which are distinct from those uncovered in Peruvian populations. However, the strongest signals of selection were related to regions of the genome that are involved in immune function related to tuberculosis. Given our estimated timing of this selection event, the Indigenous people of Ecuador may have adapted to Mycobacterium tuberculosis thousands of years before the arrival of Europeans. Furthermore, we detect a population collapse that coincides with the arrival of Europeans, which is more severe than other regions of the Andes, suggesting differing effects of contact across high-altitude populations.

Reintroduction of highly pathogenic avian influenza A H7N9 virus in southwestern China

Highly pathogenic (HP) avian influenza A H7N9 virus has emerged in China since 2016. In recent years, it has been most prevalent in northern China. However, several strains of HP H7N9 reappeared in southwestern China (Yunnan Province) in 2021. As a result, we are wondering if these viruses have re-emerged in situ or been reintroduced. Here, we present phylogenetic evidence that the HP H7N9 viruses isolated in Yunnan emigrated from northern to southwestern China in 2020. The northern subregion of China has become a novel epicenter in HP H7N9 dissemination. Meanwhile, a cleavage motif re-emerged due to the T341I mutation, implying a parallel evolution. This cross-region transmission, which originated in non-adjacent provinces and traveled a great geographic distance in an unknown way, indicates that HP H7N9 dissemination did not halt in 2020, even under the shadow of the COVID-19 pandemic. Additional surveillance studies in poultry are required to determine the HP H7N9 virus's geographic distribution and spread.

Epidemiology and Drug Resistance Patterns of Mycobacterium tuberculosis in High-Burden Area in Western Siberia, Russia

Russia is a high-burden area for multidrug-resistant tuberculosis (MDR-TB). Here, we studied the epidemiological situation and drug resistance patterns of Mycobacterium tuberculosis in the Omsk region in Western Siberia. M. tuberculosis isolates (n = 851) were recovered from newly diagnosed TB patients in 2021. The isolates were tested by bacteriological and molecular methods, and long-term epidemiological data were analyzed. The TB incidence dec, this is not variablereased from 93.9 in 2012 to 48.1 in 2021, per 100,000 population, but the primary MDR-TB rate increased from 19.2% to 26.4%. The destructive forms of tuberculosis accounted for 37.8% of all cases, while 35.5% of patients were smear-positive. Of all isolates tested, 55.2% were culture-positive, of which 94.5% were further tested for phenotypic drug resistance and associated mutations. More than half (53.4%) of isolates were drug-resistant, 13.9% were monoresistant and 67.9% were MDR. Among MDR isolates, 40.4% were pre-XDR, and 19.2% were XDR. The spectrum of drug resistance included second-line drugs (new-generation fluoroquinolones, linezolid), which significantly increase the risk of an adverse outcome in patients. In conclusion, our results highlight the critical importance of monitoring drug resistance in circulating M. tuberculosis strains emerging due to ineffective treatment and active transmission.

Impact and mitigation of sampling bias to determine viral spread: Evaluating discrete phylogeography through CTMC modeling and structured coalescent model approximations

Bayesian phylogeographic inference is a powerful tool in molecular epidemiological studies, which enables reconstruction of the origin and subsequent geographic spread of pathogens. Such inference is, however, potentially affected by geographic sampling bias. Here, we investigated the impact of sampling bias on the spatiotemporal reconstruction of viral epidemics using Bayesian discrete phylogeographic models and explored different operational strategies to mitigate this impact. We considered the continuous-time Markov chain (CTMC) model and two structured coalescent approximations (Bayesian structured coalescent approximation [BASTA] and marginal approximation of the structured coalescent [MASCOT]). For each approach, we compared the estimated and simulated spatiotemporal histories in biased and unbiased conditions based on the simulated epidemics of rabies virus (RABV) in dogs in Morocco. While the reconstructed spatiotemporal histories were impacted by sampling bias for the three approaches, BASTA and MASCOT reconstructions were also biased when employing unbiased samples. Increasing the number of analyzed genomes led to more robust estimates at low sampling bias for the CTMC model. Alternative sampling strategies that maximize the spatiotemporal coverage greatly improved the inference at intermediate sampling bias for the CTMC model, and to a lesser extent, for BASTA and MASCOT. In contrast, allowing for time-varying population sizes in MASCOT resulted in robust inference. We further applied these approaches to two empirical datasets: a RABV dataset from the Philippines and a SARS-CoV-2 dataset describing its early spread across the world. In conclusion, sampling biases are ubiquitous in phylogeographic analyses but may be accommodated by increasing the sample size, balancing spatial and temporal composition in the samples, and informing structured coalescent models with reliable case count data.

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.

Phylogeography and transmission of Mycobacterium tuberculosis spanning prisons and surrounding communities in Paraguay

Recent rises in incident tuberculosis (TB) cases in Paraguay and the increasing concentration of TB within prisons highlight the urgency of targeting strategies to interrupt transmission and prevent new infections. However, whether specific cities or carceral institutions play a disproportionate role in transmission remains unknown. We conducted prospective genomic surveillance, sequencing 471 Mycobacterium tuberculosis complex genomes, from inside and outside prisons in Paraguay's two largest urban areas, Asunción and Ciudad del Este, from 2016 to 2021. We found genomic evidence of frequent recent transmission within prisons and transmission linkages spanning prisons and surrounding populations. We identified a signal of frequent M. tuberculosis spread between urban areas and marked recent population size expansion of the three largest genomic transmission clusters. Together, our findings highlight the urgency of strengthening TB control programs to reduce transmission risk within prisons in Paraguay, where incidence was 70 times that outside prisons in 2021.

Evolutionary history and spread of the Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage, Tunisia

BACKGROUND

To infer the origin and spread of the Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage in a Mediterranean country, Tunisia, where it predominates.

METHODS

We combined Bayesian (STRUCTURE) and maximum likelihood (MIGRAINE) estimation approaches based on a global 24-loci mycobacterial interspersed repetitive units-variable numbers of tandem repeats (MIRU-VNTR24) genotyping dataset consisting of 1573 L4.3/LAM clinical strains from four continents, including 252 isolates originating from Tunisia.

RESULTS

Phylogenetic analyses coupled with Bayesian estimations suggested that the most predominant L4.3/LAM subpopulation in Tunisia (65.07%), which is dominated by a single clonal complex, TUN4.3_CC1 (94.51%), has evolved from an ancestral pool that is restricted to Europe and Africa, contrasting with the remaining L4.3/LAM subpopulations whose ancestry was traced all over the word. Maximum likelihood analysis revealed that TUN4.3_CC1 has been undergoing a demographic expansion since 131 years ago (CI95%: 90.7-205), thus explaining its preponderance relative to the second most predominant CC, TUN4.3_CC2, whose population was found under contraction.

CONCLUSIONS

The preponderance of L4.3/LAM in Tunisia stems from a 130-year expansion process of a locally evolved clone.

Evolutionary history and global spatiotemporal pattern of alfalfa mosaic virus

Alfalfa mosaic virus (AMV) is an important plant virus causing considerable economic loss to alfalfa production. Knowledge of the evolutionary and demographic history of the pathogen is limited but essential to the development of effective and sustainable pathogen management schemes. In this study, we performed worldwide phylodynamic analyses of AMV based on 154 nucleotide sequences of the coat protein gene, sampled from 1985 to 2020, to understand the epidemiology of this pathogen. Bayesian phylogenetic reconstruction estimates that the crown group of AMV dates back to 1840 (95% credibility interval, 1687-1955). We revealed that AMV continuously evolves at a rate of 4.14 × 10^-4 substitutions/site/year (95% credibility interval, 1.04 × 10^-4 - 6.68 × 10^-4). Our phylogeographic analyses identified multiple migration links between Europe and other regions, implying that Europe played a key role in spreading the virus worldwide. Further analyses showed that the clustering pattern of AMV isolates is significantly correlated to geographic regions, indicating that geography-driven adaptation may be a factor that affects the evolution of AMV. Our findings may be potentially used in the development of effective control strategies for AMV.

The Evolution and Global Spatiotemporal Dynamics of Senecavirus A

Recurrent outbreaks of senecavirus A (SVA)-associated vesicular disease have led to a large number of infected pigs being culled and has caused considerable economic losses to the swine industry. Although SVA was discovered 2 decades ago, knowledge about the evolutionary and transmission histories of SVA remains unclear. Herein, we performed an integrated analysis of the recombination, phylogeny, selection, and spatiotemporal dynamics of SVA. Phylogenetic analysis demonstrated that SVA diverged into two main branches, clade I (pre-2007 strains) and clade II (post-2007 strains). Importantly, analysis of selective strength showed that clade II was evolving under relaxed selection compared with clade I. Positive selection analysis identified 27 positive selective sites, most of which are located on the outer surface of capsid protomer or on the important functional domains of nonstructure proteins. Bayesian phylodynamics suggested that the estimated time to the most recent common ancestor of SVA was around 1986, and the estimated substitution rate of SVA was 3.3522 × 10-3 nucleotide substitutions/site/year. Demographic history analysis revealed that the effective population size of SVA has experienced a gradually increasing trend with slight fluctuation until 2017 followed by a sharp decline. Notably, Bayesian phylogeographic analysis inferred that Brazil might be the source of SVA's global transmission since 2015. In summary, these data illustrated that the ongoing evolution of SVA drove the lineage-specific innovation and potentially phenotypically important variation. Our study sheds new light on the fundamental understanding of SVA evolution and spread history. IMPORTANCE Recurrent outbreaks and global epidemics of senecavirus A-associated vesicular disease have caused heavy economic losses and have threatened the development of the pig industry. However, the question of where the virus comes from has been one of the biggest puzzles due to the stealthy nature of the virus. Consequently, tracing the source, evolution, and transmission pattern of SVA is a very challenging task. Based on the most comprehensive analysis, we revealed the origin time, rapid evolution, epidemic dynamics, and selection of SVA. We observed two main genetic branches, clade I (pre-2007 strains) and clade II (post-2007 strains), and described the epidemiological patterns of SVA in different countries. We also first identified Brazil as the source of SVA's global transmission since 2015. Findings in this study provide important implications for the control and prevention of the virus.

Comparative genomics of drug-resistant strains of Mycobacterium tuberculosis in Ecuador

BACKGROUND

Tuberculosis is a serious infectious disease affecting millions of people. In spite of efforts to reduce the disease, increasing antibiotic resistance has contributed to persist in the top 10 causes of death worldwide. In fact, the increased cases of multi (MDR) and extreme drug resistance (XDR) worldwide remains the main challenge for tuberculosis control. Whole genome sequencing is a powerful tool for predicting drug resistance-related variants, studying lineages, tracking transmission, and defining outbreaks. This study presents the identification and characterization of resistant clinical isolates of Mycobacterium tuberculosis including a phylogenetic and molecular resistance profile study by sequencing the complete genome of 24 strains from different provinces of Ecuador.

RESULTS

Genomic sequencing was used to identify the variants causing resistance. A total of 15/21 isolates were identified as MDR, 4/21 as pre-XDR and 2/21 as XDR, with three isolates discarded due to low quality; the main sub-lineage was LAM (61.9%) and Haarlem (19%) but clades X, T and S were identified. Of the six pre-XDR and XDR strains, it is noteworthy that five come from females; four come from the LAM sub-lineage and two correspond to the X-class sub-lineage. A core genome of 3,750 genes, distributed in 295 subsystems, was determined. Among these, 64 proteins related to virulence and implicated in the pathogenicity of M. tuberculosis and 66 possible pharmacological targets stand out. Most variants result in nonsynonymous amino acid changes and the most frequent genotypes were identified as conferring resistance to rifampicin, isoniazid, ethambutol, para-aminosalicylic acid and streptomycin. However, an increase in the resistance to fluoroquinolones was detected.

CONCLUSION

This work shows for the first time the variability of circulating resistant strains between men and women in Ecuador, highlighting the usefulness of genomic sequencing for the identification of emerging resistance. In this regard, we found an increase in fluoroquinolone resistance. Further sampling effort is needed to determine the total variability and associations with the metadata obtained to generate better health policies.

Biocultural perspectives of infectious diseases and demographic evolution: Tuberculosis and its comorbidities through history

Anthropologists recognize the importance of conceptualizing health in the context of the mutually evolving nature of biology and culture through the biocultural approach, but biocultural anthropological perspectives of infectious diseases and their impacts on humans (and vice versa) through time are relatively underrepresented. Tuberculosis (TB) has been a constant companion of humans for thousands of years and has heavily influenced population health in almost every phase of cultural and demographic evolution. TB in human populations has been dramatically influenced by behavior, demographic and epidemiological shifts, and other comorbidities through history. This paper critically discusses TB and some of its major comorbidities through history within a biocultural framework to show how transitions in human demography and culture affected the disease-scape of TB. In doing so, I address the potential synthesis of biocultural and epidemiological transition theory to better comprehend the mutual evolution of infectious diseases and humans.

Bayesian Phylogeographic Inference Suggests Japan as the Center for the Origin and Dissemination of Rice Stripe Virus

Rice stripe virus (RSV) is one of the most important viral pathogens of rice in East Asia. The origin and dispersal of RSV remain poorly understood, but an emerging hypothesis suggests that: (i) RSV originates from Yunnan, a southwest province of China; and (ii) some places of eastern China have acted as a center for the international dissemination of RSV. This hypothesis, however, has never been tested rigorously. Using a data set comprising more than 200 time-stamped coat protein gene sequences of RSV from Japan, China and South Korea, we reconstructed the phylogeographic history of RSV with Bayesian phylogeographic inference. Unexpectedly, the results did not support the abovementioned hypothesis. Instead, they suggested that RSV originates from Japan and Japan has been the major center for the dissemination of RSV in the past decades. Based on these data and the temporal dynamics of RSV reported recently by another group, we proposed a new hypothesis to explain the origin and dispersal of RSV. This new hypothesis may be valuable for further studies aiming to clarify the epidemiology of RSV. It may also be useful in designing management strategies against this devastating virus.

Evolution of Mycobacterium tuberculosis drug resistance in the genomic era

Mycobacterium tuberculosis has acquired drug resistance to all drugs that have been used against it, including those only recently introduced into clinical practice. Compared to other bacteria, it has a well conserved genome due to its role as an obligate human pathogen that has adapted to a niche over five to ten thousand years. These features facilitate reconstruction and dating of M. tuberculosis phylogenies, giving key insights into how resistance has been acquired and spread globally. Resistance to each new drug has occurred within five to ten years of clinical use and has occurred even more rapidly with recently introduced drugs. In most cases, resistance-conferring mutations come with a fitness cost, but this can be overcome by compensatory mutations which restore fitness to that of wild-type bacteria. It is likely that M. tuberculosis acquires drug resistance while maintaining limited genomic variability due the generation of low frequency within-host variation, combined with ongoing purifying selection causing loss of variants without a clear fitness advantage. However, variants that do confer an advantage, such as drug resistance, can increase in prevalence amongst all bacteria within a host and become the dominant clone. These resistant strains can then be transmitted leading to primary drug resistant infection in a new host. As many countries move towards genomic methods for diagnosis of M. tuberculosis infection and drug resistance, it is important to be aware of the implications for the evolution of resistance. Currently, understanding of resistance-conferring mutations is incomplete, and some targeted genetic diagnostics create their own selective pressures. We discuss an example where a rifampicin resistance-conferring mutation which was not routinely covered by standard testing became dominant. Finally, resistance to new drugs such as bedaquiline and delamanid is caused by individually rare mutations occurring across a large mutational genomic target that have been detected over a short time, and do not provide statistical power for genotype-phenotype correlation – in contrast to longer-established drugs that form the backbone of drug-sensitive antituberculosis therapy. Therefore, we need a different approach to identify resistance-conferring mutations of new drugs before their resistance becomes widespread, abrogating their usefulness.

The Spatial Diffusion of Cherry Leaf Roll Virus Revealed by a Bayesian Phylodynamic Analysis

Cherry leaf roll virus (CLRV) is an important plant pathogen that causes severe and detrimental effects on cherry and other fruit plants. Despite recent progress in plant pathology, molecular biology, and population genetics of CLRV, the spatiotemporal spread of this virus remains poorly studied. In this study, we employed a Bayesian phylodynamics framework to investigate the spatial diffusion patterns of CLRV by analyzing the coat protein gene sequences of 81 viral isolates collected from five different countries. Consistent with the trade of cherry, our Bayesian phylodynamic analyses pointed to viral origins in New Zealand and identified multiple migration pathways between Germany and other countries, suggesting that Germany has played an important role in CLRV transmission. The results of our study will be useful in developing sustainable management strategies to control this pathogen.

Transcontinental spread and evolution of Mycobacterium tuberculosis W148 European/Russian clade toward extensively drug resistant tuberculosis

Transmission-driven multi-/extensively drug resistant (M/XDR) tuberculosis (TB) is the largest single contributor to human mortality due to antimicrobial resistance. A few major clades of the Mycobacterium tuberculosis complex belonging to lineage 2, responsible for high prevalence of MDR-TB in Eurasia, show outstanding transnational distributions. Here, we determined factors underlying the emergence and epidemic spread of the W148 clade by genome sequencing and Bayesian demogenetic analyses of 720 isolates from 23 countries. We dated a common ancestor around 1963 and identified two successive epidemic expansions in the late 1980s and late 1990s, coinciding with major socio-economic changes in the post-Soviet Era. These population expansions favored accumulation of resistance mutations to up to 11 anti-TB drugs, with MDR evolving toward additional resistances to fluoroquinolones and second-line injectable drugs within 20 years on average. Timescaled haplotypic density analysis revealed that widespread acquisition of compensatory mutations was associated with transmission success of XDR strains. Virtually all W148 strains harbored a hypervirulence-associated ppe38 gene locus, and incipient recurrent emergence of prpR mutation-mediated drug tolerance was detected. The outstanding genetic arsenal of this geographically widespread M/XDR strain clade represents a “perfect storm” that jeopardizes the successful introduction of new anti-M/XDR-TB antibiotic regimens.

An outbreak of the multidrug-resistant Mycobacterium tuberculosis lineage W148 has spread widely across Russia, Central Asia and Europe. Here, the authors use whole genome sequences of ~700 isolates of this lineage collected over ~20 years to analyze its spread, evolution of drug resistance, and impact of compensatory mutations.

Genetic Diversity of Mycobacterium tuberculosis Isolates From an Amerindian Population in Chiapas, México

This is the first report of the genetic diversity of the Mycobacterium tuberculosis complex isolates found in a Mexican-Amerindian setting. In this study, we analyzed isolates collected from the Highlands region of Chiapas, Mexico, by using spoligotyping and whole-genome sequencing analyses. Seventy-three M. tuberculosis isolates were analyzed initially by spoligotyping; no new spoligotypes were identified. Nineteen percent of the isolates were identified as SIT53 (T1) (n = 14), followed by SIT42 (14%, n = 10, LAM9) and SIT119 (11%; n = 8, X1). SIT53, SIT42, and orphan isolates (16.4%, n = 12) constituted about 50% of the isolates studied and were subjected to whole-genome sequencing (WGS) analysis. Most SIT53 (10/12) isolates belonged to the Euro-American sub-lineage 4.8. Most SIT42 isolates (4/7) as .well as most orphan isolates (5/8) belonged to the lineage 4.3.3 LAM group. By comparing the single-nucleotide polymorphism (SNP) patterns of the SIT53 isolates, we found one clone (&lt;7 SNPs) and four clustered isolates (&lt;15 SNPs). In isolates from the SIT42 and orphan groups, we did not find any clones or clusters. This work demonstrates the success of sub-lineage 4.8 to predominate in Mexico and confirms the dominion of sub-lineage 4.3.3 in Central and South America.

Whole genome sequencing analysis to evaluate the influence of T2DM on polymorphisms associated with drug resistance in M. tuberculosis

Background

Type 2 diabetes mellitus (T2DM) has been associated with treatment failure, and the development of drug resistance in tuberculosis (TB). Also, whole-genome sequencing has provided a better understanding and allowed the growth of knowledge about polymorphisms in genes associated with drug resistance. Considering the above, this study analyzes genome sequences to evaluate the influence of type 2 diabetes mellitus in the development of mutations related to tuberculosis drug resistance. M. tuberculosis isolates from individuals with (n = 74), and without (n = 74) type 2 diabetes mellitus was recovered from online repositories, and further analyzed.

Results

The results showed the presence of 431 SNPs with similar proportions between diabetics, and non-diabetics individuals (48% vs. 52%), but with no significant relationship. A greater number of mutations associated with rifampicin resistance was observed in the T2DM-TB individuals (23.2% vs. 16%), and the exclusive presence of rpoBQ432L, rpoBQ432P, rpoBS441L, and rpoBH445L variants. While these variants are not private to T2DM-TB cases they are globally rare highlighting a potential role of T2DM. The phylogenetic analysis showed 12 sublineages, being 4.1.1.3, and 4.1.2.1 the most prevalent in T2DM-TB individuals but not differing from those most prevalent in their geographic location. Four clonal complexes were found, however, no significant relationship with T2DM was observed. Samples size and potential sampling biases prevented us to look for significant associations.

Conclusions

The occurrence of globally rare rifampicin variants identified only in isolates from individuals with T2DM could be due to the hyperglycemic environment within the host. Therefore, further studies about the dynamics of SNPs’ generation associated with antibiotic resistance in patients with diabetes mellitus are necessary.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08709-z.

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.

The host-pathogen-environment triad: Lessons learned through the study of the multidrug-resistant Mycobacterium tuberculosis M strain

Multidrug-resistant tuberculosis is one of the major obstacles that face the tuberculosis eradication efforts. Drug-resistant Mycobacterium tuberculosis clones were initially disregarded as a public health threat, because they were assumed to have paid a high fitness cost in exchange of resistance acquisition. However, some genotypes manage to overcome the impact of drug-resistance conferring mutations, retain transmissibility and cause large outbreaks. In Argentina, the HIV-AIDS epidemics fuelled the expansion of the so-called M strain in the early 1990s, which is responsible for the largest recorded multidrug-resistant tuberculosis cluster of Latin America. The aim of this work is to review the knowledge gathered after nearly three decades of multidisciplinary research on epidemiological, microbiological and immunological aspects of this highly successful strain. Collectively, our results indicate that the successful transmission of the M strain could be ascribed to its unaltered virulence, low Th1/Th17 response, a low fitness cost imposed by the resistance conferring mutations and a high resistance to host-related stress. In the early 2000s, the incident cases due to the M strain steadily declined and stabilized in the latest years. Improvements in the management, diagnosis and treatment of multidrug-resistant tuberculosis along with societal factors such as the low domestic and international mobility of the patients affected by this strain probably contributed to the outbreak containment. This stresses the importance of sustaining the public health interventions to avoid the resurgence of this conspicuous multidrug-resistant strain.