Transition bias influences the evolution of antibiotic resistance in Mycobacterium tuberculosis.

Transition bias, an overabundance of transitions relative to transversions, has been widely reported among studies of the rates and spectra of spontaneous mutations. However, demonstrating the role of transition bias in adaptive evolution remains challenging. In particular, it is unclear whether such biases direct the evolution of bacterial pathogens adapting to treatment. We addressed this challenge by analyzing adaptive antibiotic-resistance mutations in the major human pathogen Mycobacterium tuberculosis (MTB). We found strong evidence for transition bias in two independently curated data sets comprising 152 and 208 antibiotic-resistance mutations. This was true at the level of mutational paths (distinct adaptive DNA sequence changes) and events (individual instances of the adaptive DNA sequence changes) and across different genes and gene promoters conferring resistance to a diversity of antibiotics. It was also true for mutations that do not code for amino acid changes (in gene promoters and the 16S ribosomal RNA gene rrs) and for mutations that are synonymous to each other and are therefore likely to have similar fitness effects, suggesting that transition bias can be caused by a bias in mutation supply. These results point to a central role for transition bias in determining which mutations drive adaptive antibiotic resistance evolution in a key pathogen.

Expression Dysregulation as a Mediator of Fitness Costs in Antibiotic Resistance.

Antimicrobial resistance (AMR) poses a threat to global health and the economy. Rifampicin-resistant Mycobacterium tuberculosis accounts for a third of the global AMR burden. Gaining the upper hand on AMR requires a deeper understanding of the physiology of resistance. AMR often results in a fitness cost in the absence of drug. Identifying the molecular mechanisms underpinning this cost could help strengthen future treatment regimens. Here, we used a collection of M. tuberculosis strains that provide an evolutionary and phylogenetic snapshot of rifampicin resistance and subjected them to genome-wide transcriptomic and proteomic profiling to identify key perturbations of normal physiology. We found that the clinically most common rifampicin resistance-conferring mutation, RpoB Ser450Leu, imparts considerable gene expression changes, many of which are mitigated by the compensatory mutation in RpoC Leu516Pro. However, our data also provide evidence for pervasive epistasis-the same resistance mutation imposed a different fitness cost and functionally distinct changes to gene expression in genetically unrelated clinical strains. Finally, we report a likely posttranscriptional modulation of gene expression that is shared in most of the tested strains carrying RpoB Ser450Leu, resulting in an increased abundance of proteins involved in central carbon metabolism. These changes contribute to a more general trend in which the disruption of the composition of the proteome correlates with the fitness cost of the RpoB Ser450Leu mutation in different strains.

Whole-Genome Sequencing for Drug Resistance Profile Prediction in Mycobacterium tuberculosis.

Whole-genome sequencing allows rapid detection of drug-resistant Mycobacterium tuberculosis isolates. However, the availability of high-quality data linking quantitative phenotypic drug susceptibility testing (DST) and genomic data have thus far been limited. We determined drug resistance profiles of 176 genetically diverse clinical M. tuberculosis isolates from the Democratic Republic of the Congo, Ivory Coast, Peru, Thailand, and Switzerland by quantitative phenotypic DST for 11 antituberculous drugs using the BD Bactec MGIT 960 system and 7H10 agar dilution to generate a cross-validated phenotypic DST readout. We compared DST results with predicted drug resistance profiles inferred by whole-genome sequencing. Classification of strains by the two phenotypic DST methods into resistotype/wild-type populations was concordant in 73 to 99% of cases, depending on the drug. Our data suggest that the established critical concentration (5 mg/liter) for ethambutol resistance (MGIT 960 system) is too high and misclassifies strains as susceptible, unlike 7H10 agar dilution. Increased minimal inhibitory concentrations were explained by mutations identified by whole-genome sequencing. Using whole-genome sequences, we were able to predict quantitative drug resistance levels for the majority of drug resistance mutations. Predicting quantitative levels of drug resistance by whole-genome sequencing was partially limited due to incompletely understood drug resistance mechanisms. The overall sensitivity and specificity of whole-genome-based DST were 86.8% and 94.5%, respectively. Despite some limitations, whole-genome sequencing has the potential to infer resistance profiles without the need for time-consuming phenotypic methods.

Treemmer: a tool to reduce large phylogenetic datasets with minimal loss of diversity.

BACKGROUND: Large sequence datasets are difficult to visualize and handle. Additionally, they often do not represent a random subset of the natural diversity, but the result of uncoordinated and convenience sampling. Consequently, they can suffer from redundancy and sampling biases. RESULTS: Here we present Treemmer, a simple tool to evaluate the redundancy of phylogenetic trees and reduce their complexity by eliminating leaves that contribute the least to the tree diversity. CONCLUSIONS: Treemmer can reduce the size of datasets with different phylogenetic structures and levels of redundancy while maintaining a sub-sample that is representative of the original diversity. Additionally, it is possible to fine-tune the behavior of Treemmer including any kind of meta-information, making Treemmer particularly useful for empirical studies.

The relative transmission fitness of multidrug-resistant Mycobacterium tuberculosis in a drug resistance hotspot.

Multidrug-resistant tuberculosis (MDR-TB) is among the most frequent causes of death due to antimicrobial resistance. Although only 3% of global TB cases are MDR, geographical hotspots with up to 40% of MDR-TB have been observed in countries of the former Soviet Union. While the quality of TB control and patient-related factors are known contributors to such hotspots, the role of the pathogen remains unclear. Here we show that in the country of Georgia, a known hotspot of MDR-TB, MDR Mycobacterium tuberculosis strains of lineage 4 (L4) transmit less than their drug-susceptible counterparts, whereas most MDR strains of L2 suffer no such defect. Our findings further indicate that the high transmission fitness of these L2 strains results from epistatic interactions between the rifampicin resistance-conferring mutation RpoB S450L, compensatory mutations in the RNA polymerase, and other pre-existing genetic features of L2/Beijing clones that circulate in Georgia. We conclude that the transmission fitness of MDR M. tuberculosis strains is heterogeneous, but can be as high as drug-susceptible forms, and that such highly drug-resistant and transmissible strains contribute to the emergence and maintenance of hotspots of MDR-TB. As these strains successfully overcome the metabolic burden of drug resistance, and given the ongoing rollout of new treatment regimens against MDR-TB, proper surveillance should be implemented to prevent these strains from acquiring resistance to the additional drugs.

An inactivated multivalent influenza A virus vaccine is broadly protective in mice and ferrets.

Influenza A viruses (IAVs) present major public health threats from annual seasonal epidemics and pandemics and from viruses adapted to a variety of animals including poultry, pigs, and horses. Vaccines that broadly protect against all such IAVs, so-called "universal" influenza vaccines, do not currently exist but are urgently needed. Here, we demonstrated that an inactivated, multivalent whole-virus vaccine, delivered intramuscularly or intranasally, was broadly protective against challenges with multiple IAV hemagglutinin and neuraminidase subtypes in both mice and ferrets. The vaccine is composed of four ß-propiolactone-inactivated low-pathogenicity avian IAV subtypes of H1N9, H3N8, H5N1, and H7N3. Vaccinated mice and ferrets demonstrated substantial protection against a variety of IAVs, including the 1918 H1N1 strain, the highly pathogenic avian H5N8 strain, and H7N9. We also observed protection against challenge with antigenically variable and heterosubtypic avian, swine, and human viruses. Compared to control animals, vaccinated mice and ferrets demonstrated marked reductions in viral titers, lung pathology, and host inflammatory responses. This vaccine approach indicates the feasibility of eliciting broad, heterosubtypic IAV protection and identifies a promising candidate for influenza vaccine clinical development.

Prisons as ecological drivers of fitness-compensated multidrug-resistant Mycobacterium tuberculosis.

Multidrug-resistant tuberculosis (MDR-TB) accounts for one third of the annual deaths due to antimicrobial resistance1. Drug resistance-conferring mutations frequently cause fitness costs in bacteria2-5. Experimental work indicates that these drug resistance-related fitness costs might be mitigated by compensatory mutations6-10. However, the clinical relevance of compensatory evolution remains poorly understood. Here we show that, in the country of Georgia, during a 6-year nationwide study, 63% of MDR-TB was due to patient-to-patient transmission. Compensatory mutations and patient incarceration were independently associated with transmission. Furthermore, compensatory mutations were overrepresented among isolates from incarcerated individuals that also frequently spilled over into the non-incarcerated population. As a result, up to 31% of MDR-TB in Georgia was directly or indirectly linked to prisons. We conclude that prisons fuel the epidemic of MDR-TB in Georgia by acting as ecological drivers of fitness-compensated strains with high transmission potential.

Lung epithelial and endothelial damage, loss of tissue repair, inhibition of fibrinolysis, and cellular senescence in fatal COVID-19.

Coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is characterized by respiratory distress, multiorgan dysfunction, and, in some cases, death. The pathological mechanisms underlying COVID-19 respiratory distress and the interplay with aggravating risk factors have not been fully defined. Lung autopsy samples from 18 patients with fatal COVID-19, with symptom onset-to-death times ranging from 3 to 47 days, and antemortem plasma samples from 6 of these cases were evaluated using deep sequencing of SARS-CoV-2 RNA, multiplex plasma protein measurements, and pulmonary gene expression and imaging analyses. Prominent histopathological features in this case series included progressive diffuse alveolar damage with excessive thrombosis and late-onset pulmonary tissue and vascular remodeling. Acute damage at the alveolar-capillary barrier was characterized by the loss of surfactant protein expression with injury to alveolar epithelial cells, endothelial cells, respiratory epithelial basal cells, and defective tissue repair processes. Other key findings included impaired clot fibrinolysis with increased concentrations of plasma and lung plasminogen activator inhibitor-1 and modulation of cellular senescence markers, including p21 and sirtuin-1, in both lung epithelial and endothelial cells. Together, these findings further define the molecular pathological features underlying the pulmonary response to SARS-CoV-2 infection and provide important insights into signaling pathways that may be amenable to therapeutic intervention.

Local adaptation in populations of Mycobacterium tuberculosis endemic to the Indian Ocean Rim.

Background: Lineage 1 (L1) and 3 (L3) are two lineages of the Mycobacterium tuberculosis complex (MTBC) causing tuberculosis (TB) in humans. L1 and L3 are prevalent around the rim of the Indian Ocean, the region that accounts for most of the world's new TB cases. Despite their relevance for this region, L1 and L3 remain understudied. Methods: We analyzed 2,938 L1 and 2,030 L3 whole genome sequences originating from 69 countries. We reconstructed the evolutionary history of these two lineages and identified genes under positive selection. Results: We found a strongly asymmetric pattern of migration from South Asia toward neighboring regions, highlighting the historical role of South Asia in the dispersion of L1 and L3. Moreover, we found that several genes were under positive selection, including genes involved in virulence and resistance to antibiotics. For L1 we identified signatures of local adaptation at the esxH locus, a gene coding for a secreted effector that targets the human endosomal sorting complex, and is included in several vaccine candidates. Conclusions: Our study highlights the importance of genetic diversity in the MTBC, and sheds new light on two of the most important MTBC lineages affecting humans.

Population Genomics of Francisella tularensis subsp. holarctica and its Implication on the Eco-Epidemiology of Tularemia in Switzerland.

Whole genome sequencing (WGS) methods provide new possibilities in the field of molecular epidemiology. This is particularly true for monomorphic organisms where the discriminatory power of traditional methods (e.g., restriction enzyme length polymorphism typing, multi locus sequence typing etc.) is inadequate to elucidate complex disease transmission patterns, as well as resolving the phylogeny at high resolution on a micro-geographic scale. In this study, we present insights into the population structure of Francisella tularensis subsp. holarctica, the causative agent of tularemia in Switzerland. A total of 59 Fth isolates were obtained from castor bean ticks (Ixodes ricinus), animals and humans and a high resolution phylogeny was inferred using WGS methods. The majority of the Fth population in Switzerland belongs to the west European B.11 clade and shows an extraordinary genetic diversity underlining the old evolutionary history of the pathogen in the alpine region. Moreover, a new B.11 subclade was identified which was not described so far. The combined analysis of the epidemiological data of human tularemia cases with the whole genome sequences of the 59 isolates provide evidence that ticks play a pivotal role in transmitting Fth to humans and other vertebrates in Switzerland. This is further underlined by the correlation of disease risk estimates with climatic and ecological factors influencing the survival of ticks.

Antimicrobial resistance in Mycobacterium tuberculosis: mechanistic and evolutionary perspectives.

Antibiotic-resistant Mycobacterium tuberculosis strains are threatening progress in containing the global tuberculosis epidemic. Mycobacterium tuberculosis is intrinsically resistant to many antibiotics, limiting the number of compounds available for treatment. This intrinsic resistance is due to a number of mechanisms including a thick, waxy, hydrophobic cell envelope and the presence of drug degrading and modifying enzymes. Resistance to the drugs which are active against M. tuberculosis is, in the absence of horizontally transferred resistance determinants, conferred by chromosomal mutations. These chromosomal mutations may confer drug resistance via modification or overexpression of the drug target, as well as by prevention of prodrug activation. Drug resistance mutations may have pleiotropic effects leading to a reduction in the bacterium's fitness, quantifiable e.g. by a reduction in the in vitro growth rate. Secondary so-called compensatory mutations, not involved in conferring resistance, can ameliorate the fitness cost by interacting epistatically with the resistance mutation. Although the genetic diversity of M. tuberculosis is low compared to other pathogenic bacteria, the strain genetic background has been demonstrated to influence multiple aspects in the evolution of drug resistance. The rate of resistance evolution and the fitness costs of drug resistance mutations may vary as a function of the genetic background.