Correction: The structural role of SARS-CoV-2 genetic background in the emergence and success of spike mutations: The case of the spike A222V mutation.

[This corrects the article DOI: 10.1371/journal.ppat.1010631.].

The structural role of SARS-CoV-2 genetic background in the emergence and success of spike mutations: The case of the spike A222V mutation.

The S:A222V point mutation, within the G clade, was characteristic of the 20E (EU1) SARS-CoV-2 variant identified in Spain in early summer 2020. This mutation has since reappeared in the Delta subvariant AY.4.2, raising questions about its specific effect on viral infection. We report combined serological, functional, structural and computational studies characterizing the impact of this mutation. Our results reveal that S:A222V promotes an increased RBD opening and slightly increases ACE2 binding as compared to the parent S:D614G clade. Finally, S:A222V does not reduce sera neutralization capacity, suggesting it does not affect vaccine effectiveness.

Pre-Columbian mycobacterial genomes reveal seals as a source of New World human tuberculosis.

Modern strains of Mycobacterium tuberculosis from the Americas are closely related to those from Europe, supporting the assumption that human tuberculosis was introduced post-contact. This notion, however, is incompatible with archaeological evidence of pre-contact tuberculosis in the New World. Comparative genomics of modern isolates suggests that M. tuberculosis attained its worldwide distribution following human dispersals out of Africa during the Pleistocene epoch, although this has yet to be confirmed with ancient calibration points. Here we present three 1,000-year-old mycobacterial genomes from Peruvian human skeletons, revealing that a member of the M. tuberculosis complex caused human disease before contact. The ancient strains are distinct from known human-adapted forms and are most closely related to those adapted to seals and sea lions. Two independent dating approaches suggest a most recent common ancestor for the M. tuberculosis complex less than 6,000 years ago, which supports a Holocene dispersal of the disease. Our results implicate sea mammals as having played a role in transmitting the disease to humans across the ocean.

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.

Acquired Resistance to Bedaquiline and Delamanid in Therapy for Tuberculosis.

Treatment of multidrug-resistant Mycobacterium tuberculosis is a challenge. This letter describes the emergence of resistance to new therapies, bedaquiline and delamanid.

Within Host Evolution Selects for a Dominant Genotype of Mycobacterium tuberculosis while T Cells Increase Pathogen Genetic Diversity.

Molecular epidemiological assessments, drug treatment optimization, and development of immunological interventions all depend on understanding pathogen adaptation and genetic variation, which differ for specific pathogens. Mycobacterium tuberculosis is an exceptionally successful human pathogen, yet beyond knowledge that this bacterium has low overall genomic variation but acquires drug resistance mutations, little is known of the factors that drive its population genomic characteristics. Here, we compared the genetic diversity of the bacteria that established infection to the bacterial populations obtained from infected tissues during murine M. tuberculosis pulmonary infection and human disseminated M. bovis BCG infection. We found that new mutations accumulate during in vitro culture, but that in vivo, purifying selection against new mutations dominates, indicating that M. tuberculosis follows a dominant lineage model of evolution. Comparing bacterial populations passaged in T cell-deficient and immunocompetent mice, we found that the presence of T cells is associated with an increase in the diversity of the M. tuberculosis genome. Together, our findings put M. tuberculosis genetic evolution in a new perspective and clarify the impact of T cells on sequence diversity of M. tuberculosis.

Consequences of genomic diversity in Mycobacterium tuberculosis.

The causative agent of human tuberculosis, Mycobacterium tuberculosis complex (MTBC), comprises seven phylogenetically distinct lineages associated with different geographical regions. Here we review the latest findings on the nature and amount of genomic diversity within and between MTBC lineages. We then review recent evidence for the effect of this genomic diversity on mycobacterial phenotypes measured experimentally and in clinical settings. We conclude that overall, the most geographically widespread Lineage 2 (includes Beijing) and Lineage 4 (also known as Euro-American) are more virulent than other lineages that are more geographically restricted. This increased virulence is associated with delayed or reduced pro-inflammatory host immune responses, greater severity of disease, and enhanced transmission. Future work should focus on the interaction between MTBC and human genetic diversity, as well as on the environmental factors that modulate these interactions.

Biological and Epidemiological Consequences of MTBC Diversity.

Tuberculosis is caused by different groups of bacteria belonging to the Mycobacterium tuberculosis complex (MTBC). The combined action of human factors, environmental conditions and bacterial virulence determine the extent and form of human disease. MTBC virulence is a composite of different clinical phenotypes such as transmission rate and disease severity among others. Clinical phenotypes are also influenced by cellular and immunological phenotypes. MTBC phenotypes are determined by the genotype, therefore finding genotypes responsible for clinical phenotypes would allow discovering MTBC virulence factors. Different MTBC strains display different cellular and clinical phenotypes. Strains from Lineage 5 and Lineage 6 are metabolically different, grow slower, and are less virulent. Also, at least certain groups of Lineage 2 and Lineage 4 strains are more virulent in terms of disease severity and human-to-human transmission. Because phenotypic differences are ultimately caused by genotypic differences, different genomic loci have been related to various cellular and clinical phenotypes. However, defining the impact of specific bacterial genomic loci on virulence when other bacterial determinants, human and environmental factors are also impacting the phenotype would contribute to a better knowledge of tuberculosis virulence and ultimately benefit tuberculosis control.

Standard Genotyping Overestimates Transmission of Mycobacterium tuberculosis among Immigrants in a Low-Incidence Country.

Immigrants from regions with a high incidence of tuberculosis (TB) are a risk group for TB in low-incidence countries such as Switzerland. In a previous analysis of a nationwide collection of 520 Mycobacterium tuberculosis isolates from 2000 to 2008, we identified 35 clusters comprising 90 patients based on standard genotyping (24-locus mycobacterial interspersed repetitive-unit-variable-number tandem-repeat [MIRU-VNTR] typing and spoligotyping). Here, we used whole-genome sequencing (WGS) to revisit these transmission clusters. Genome-based transmission clusters were defined as isolate pairs separated by ≤12 single nucleotide polymorphisms (SNPs). WGS confirmed 17/35 (49%) MIRU-VNTR typing clusters; the other 18 clusters contained pairs separated by >12 SNPs. Most transmission clusters (3/4) of Swiss-born patients were confirmed by WGS, as opposed to 25% (4/16) of the clusters involving only foreign-born patients. The overall clustering proportion was 17% (90 patients; 95% confidence interval [CI], 14 to 21%) by standard genotyping but only 8% (43 patients; 95% CI, 6 to 11%) by WGS. The clustering proportion was 17% (67/401; 95% CI, 13 to 21%) by standard genotyping and 7% (26/401; 95% CI, 4 to 9%) by WGS among foreign-born patients and 19% (23/119; 95% CI, 13 to 28%) and 14% (17/119; 95% CI, 9 to 22%), respectively, among Swiss-born patients. Using weighted logistic regression, we found weak evidence of an association between birth origin and transmission (adjusted odds ratio of 2.2 and 95% CI of 0.9 to 5.5 comparing Swiss-born patients to others). In conclusion, standard genotyping overestimated recent TB transmission in Switzerland compared to WGS, particularly among immigrants from regions with a high TB incidence, where genetically closely related strains often predominate. We recommend the use of WGS to identify transmission clusters in settings with a low incidence of TB.

Tracking a tuberculosis outbreak over 21 years: strain-specific single-nucleotide polymorphism typing combined with targeted whole-genome sequencing.

BACKGROUND: Whole-genome sequencing (WGS) is increasingly used in molecular-epidemiological investigations of bacterial pathogens, despite cost- and time-intensive analyses. We combined strain-specific single-nucleotide polymorphism (SNP) typing and targeted WGS to investigate a tuberculosis cluster spanning 21 years in Bern, Switzerland. METHODS: On the basis of genome sequences of 3 historical outbreak Mycobacterium tuberculosis isolates, we developed a strain-specific SNP-typing assay to identify further cases. We screened 1642 patient isolates and performed WGS on all identified cluster isolates. We extracted SNPs to construct genomic networks. Clinical and social data were retrospectively collected. RESULTS: We identified 68 patients associated with the outbreak strain. Most received a tuberculosis diagnosis in 1991-1995, but cases were observed until 2011. Two thirds were homeless and/or substance abusers. Targeted WGS revealed 133 variable SNP positions among outbreak isolates. Genomic network analyses suggested a single origin of the outbreak, with subsequent division into 3 subclusters. Isolates from patients with confirmed epidemiological links differed by 0-11 SNPs. CONCLUSIONS: Strain-specific SNP genotyping allowed rapid and inexpensive identification of M. tuberculosis outbreak isolates in a population-based strain collection. Subsequent targeted WGS provided detailed insights into transmission dynamics. This combined approach could be applied to track bacterial pathogens in real time and at high resolution.

Genomic epidemiology of multidrug-resistant Mycobacterium tuberculosis during transcontinental spread.

The transcontinental spread of multidrug-resistant (MDR) tuberculosis is poorly characterized in molecular epidemiologic studies. We used genomic sequencing to understand the establishment and dispersion of MDR Mycobacterium tuberculosis within a group of immigrants to the United States. We used a genomic epidemiology approach to study a genotypically matched (by spoligotype, IS6110 restriction fragment length polymorphism, and mycobacterial interspersed repetitive units-variable number of tandem repeat signature) lineage 2/Beijing MDR strain implicated in an outbreak of tuberculosis among refugees in Thailand and consecutive cases within California. All 46 MDR M. tuberculosis genomes from both Thailand and California were highly related, with a median difference of 10 single-nucleotide polymorphisms (SNPs). The Wat Tham Krabok (WTK) strain is a new sequence type distinguished from all known Beijing strains by 55 SNPs and a genomic deletion (Rv1267c) associated with increased fitness. Sequence data revealed a highly prevalent MDR strain that included several closely related but distinct allelic variants within Thailand, rather than the occurrence of a single outbreak. In California, sequencing data supported multiple independent introductions of WTK with subsequent transmission and reactivation within the state, as well as a potential super spreader with a prolonged infectious period. Twenty-seven drug resistance-conferring mutations and 4 putative compensatory mutations were found within WTK strains. Genomic sequencing has substantial epidemiologic value in both low- and high-burden settings in understanding transmission chains of highly prevalent MDR strains.

Geographical and Temporal Structures of Legionella pneumophila Sequence Types in Comunitat Valenciana (Spain), 1998 to 2013.

Legionella pneumophila is an accidental human pathogen associated with aerosol formation in water-related sources. High recombination rates make Legionella populations genetically diverse, and nearly 2,000 different sequence types (STs) have been described to date for this environmental pathogen. The spatial distribution of STs is extremely heterogeneous, with some variants being present worldwide and others being detected at only a local scale. Similarly, some STs have been associated with disease outbreaks, such as ST578 or ST23. Spain is among the European countries with the highest incidences of reported legionellosis cases, and specifically, Comunitat Valenciana (CV) is the second most affected area in the country. In this work, we aimed at studying the overall diversity of Legionella pneumophila populations found in the period from 1998 to 2013 in 79 localities encompassing 23 regions within CV. To do so, we performed sequence-based typing (SBT) on 1,088 L. pneumophila strains detected in the area from both environmental and clinical sources. A comparison with the genetic structuring detected in a global data set that included 20 European and 7 non-European countries was performed. Our results reveal a level of diversity in CV that can be considered representative of the diversity found in other countries worldwide.

KvarQ: targeted and direct variant calling from fastq reads of bacterial genomes.

BACKGROUND: High-throughput DNA sequencing produces vast amounts of data, with millions of short reads that usually have to be mapped to a reference genome or newly assembled. Both reference-based mapping and de novo assembly are computationally intensive, generating large intermediary data files, and thus require bioinformatics skills that are often lacking in the laboratories producing the data. Moreover, many research and practical applications in microbiology require only a small fraction of the whole genome data. RESULTS: We developed KvarQ, a new tool that directly scans fastq files of bacterial genome sequences for known variants, such as single nucleotide polymorphisms (SNP), bypassing the need of mapping all sequencing reads to a reference genome and de novo assembly. Instead, KvarQ loads "testsuites" that define specific SNPs or short regions of interest in a reference genome, and directly synthesizes the relevant results based on the occurrence of these markers in the fastq files. KvarQ has a versatile command line interface and a graphical user interface. KvarQ currently ships with two "testsuites" for Mycobacterium tuberculosis, but new "testsuites" for other organisms can easily be created and distributed. In this article, we demonstrate how KvarQ can be used to successfully detect all main drug resistance mutations and phylogenetic markers in 880 bacterial whole genome sequences. The average scanning time per genome sequence was two minutes. The variant calls of a subset of these genomes were validated with a standard bioinformatics pipeline and revealed >99% congruency. CONCLUSION: KvarQ is a user-friendly tool that directly extracts relevant information from fastq files. This enables researchers and laboratory technicians with limited bioinformatics expertise to scan and analyze raw sequencing data in a matter of minutes. KvarQ is open-source, and pre-compiled packages with a graphical user interface are available at http://www.swisstph.ch/kvarq.

Mixed infection by Legionella pneumophila in outbreak patients.

During the molecular epidemiological study of a legionellosis outbreak, we obtained sequence based typing (SBT) profiles from uncultured respiratory samples of 15 affected patients. We detected several distinct allelic profiles some of which were a mixture of alleles present in the more common profiles. Chromatograms from the sequences of one patient with mixed profile showed polymorphisms in several positions, which could result from the simultaneous presence of different Legionella variants in the sample. In order to test this possibility, we cloned PCR amplification products from six loci for two patients with a mixed profile and a patient with a pure profile. After obtaining around 20 sequences for each locus of three patients, we detected several variants in two of them and two variants in the third one. In summary, the three analyzed patients showed evidence of more than one Legionella variant during the acute infection. These results indicate that probably some patients were infected by more than one strain, which could be due to co-infection from the same environmental source or, alternatively, to independent infections in a very short period of time. Although our data cannot discriminate between these hypotheses, these results suggest that Legionella infection patterns can be more complex than previously assumed. None of the environmental samples analyzed during this outbreak was even similar to any of the clinical ones.

VIPERA: Viral Intra-Patient Evolution Reporting and Analysis.

Viral mutations within patients nurture the adaptive potential of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) during chronic infections, which are a potential source of variants of concern. However, there is no integrated framework for the evolutionary analysis of intra-patient SARS-CoV-2 serial samples. Herein, we describe Viral Intra-Patient Evolution Reporting and Analysis (VIPERA), a new software that integrates the evaluation of the intra-patient ancestry of SARS-CoV-2 sequences with the analysis of evolutionary trajectories of serial sequences from the same viral infection. We have validated it using positive and negative control datasets and have successfully applied it to a new case, which revealed population dynamics and evidence of adaptive evolution. VIPERA is available under a free software license at https://github.com/PathoGenOmics-Lab/VIPERA.

"Phylogenomic insights into brucellaceae: The Pseudochrobactrum algeriensis case".

The genus Pseudochrobactrum encompasses free-living bacteria phylogenetically close to Ochrobactrum opportunistic pathogens and to Brucella, facultative intracellular parasites causing brucellosis, a worldwide-extended and grave zoonosis. Recently, Pseudochrobactrum strains were isolated from Brucella natural hosts on Brucella selective media, potentially causing diagnostic confusions. Strikingly, P. algeriensis was isolated from cattle lymph nodes, organs that are inimical to bacteria. Here, we analyse P. algeriensis potential virulence factors in comparison with Ochrobactrum and Brucella. Consistent with genomic analyses, Western-Blot analyses confirmed that P. algeriensis lacks the ability to synthesize the N-formylperosamine O-polysaccharide characteristic of the lipopolysaccharide (LPS) of smooth Brucella core species. However, unlike other Pseudochrobactrum but similar to some early diverging brucellae, P. algeriensis carries genes potentially synthetizing a rhamnose-based O-polysaccharide LPS. Lipid A analysis by MALDI-TOF demonstrated that P. algeriensis LPS bears a lipid A with a reduced pathogen-associated molecular pattern, a trait shared with Ochrobactrum and Brucella that is essential to generate a highly stable outer membrane and to delay immune activation. Also, although not able to multiply intracellularly in macrophages, the analysis of P. algeriensis cell lipid envelope revealed the presence of large amounts of cationic aminolipids, which may account for the extremely high resistance of P. algeriensis to bactericidal peptides and could favor colonization of mucosae and transient survival in Brucella hosts. However, two traits critical in Brucella pathogenicity are either significantly different (T4SS [VirB]) or absent (erythritol catabolic pathway) in P. algeriensis. This work shows that, while diverging in other characteristics, lipidic envelope features relevant in Brucella pathogenicity are conserved in Brucellaceae. The constant presence of these features strongly suggests that reinforcement of the envelope integrity as an adaptive advantage in soil was maintained in Brucella because of the similarity of some environmental challenges, such as the action of cationic peptide antibiotics and host defense peptides. This information adds knowledge about the evolution of Brucellaceae, and also underlines the taxonomical differences of the three genera compared.

Novel Mycobacterium tuberculosis complex isolate from a wild chimpanzee.

Tuberculosis (TB) is caused by gram-positive bacteria known as the Mycobacterium tuberculosis complex (MTBC). MTBC include several human-associated lineages and several variants adapted to domestic and, more rarely, wild animal species. We report an M. tuberculosis strain isolated from a wild chimpanzee in Côte d'Ivoire that was shown by comparative genomic and phylogenomic analyses to belong to a new lineage of MTBC, closer to the human-associated lineage 6 (also known as M. africanum West Africa 2) than to the other classical animal-associated MTBC strains. These results show that the general view of the genetic diversity of MTBC is limited and support the possibility that other MTBC variants exist, particularly in wild mammals in Africa. Exploring this diversity is crucial to the understanding of the biology and evolutionary history of this widespread infectious disease.

Quantifying nonvertical inheritance in the evolution of Legionella pneumophila.

The exchange of genetic material among bacterial strains and species is recognized as an important factor determining their evolutionary, population genetic, and epidemiological features. We present a detailed analysis of nonvertical inheritance in Legionella pneumophila, a human pathogen and facultative intracellular parasite of amoebas. We have analyzed the exchange of L. pneumophila genetic material with other bacteria at three different levels: population genetics, population genomics, and phylogenomics. At the population genetics level, we have analyzed 89 clinical and environmental isolates after sequencing six coding loci and three intergenic regions for a total of 3,923 bp. In the population genomics analysis, we have studied the roles of recombination and mutation in the common portion of the genome sequence of four L. pneumophila strains. In the phylogenomic analysis, we have studied the phylogenetic origin of 1,700 genes in the L. pneumophila pangenome. For this, we have considered 12 possible phylogenetic alternatives, derived from a reference tree obtained from 104 genes from 41 species, which have been tested under a rigorous statistical framework. The results obtained agree in assigning an important role to nonvertical inheritance in shaping the composition of the L. pneumophila genome and of the genetic variation in its populations. We have found a negative correlation between phylogenetic distance and likelihood of horizontal gene transfer. Phylogenetic proximity and increased chances resulting from sharing the ecological niche provided by the amoeba host have likely had a major influence on the rate of gene exchange in Legionella.

Staphylococcus prevails in the skin microbiota of long-term immunodeficient mice.

Host-commensal relationships in the skin are a complex system governed by variables related to the host, the bacteria and the environment. A disruption of this system may lead to new steady states, which, in turn, may lead to disease. We have studied one such disruption by characterizing the skin microbiota in healthy and immunodepressed (ID) mice. A detailed anatomopathological study failed to reveal any difference between the skin of healthy and ID mice. We sequenced the 16S rDNA V1-V2 gene region to saturation in 10 healthy and 10 ID 8 week-old mice, and found than all of the healthy and two of the ID mice had bacterial communities that were similar in composition to that of human skin, although, presumably because of the uniform raising conditions, less interindividual variation was found in mice. However, eight ID mice showed microbiota dominated by Staphylococcus epidermidis. Quantitative PCR amplification of 16S rDNA gene and of the Staphylococcus-specific TstaG region confirmed the previous results and indicated that the quantitative levels of Staphylococcus were similar in both groups while the total number of 16S copies was greater in the healthy mice. Thus, it is possible that, under long-term immunodeficiency, which removes the acquired but not the native immune system, S.epidermidis may inhibit the growth of other bacteria but does not cause a pathogenic state.

Drug resistance-conferring mutations in Mycobacterium tuberculosis from Madang, Papua New Guinea.

BACKGROUND: Monitoring drug resistance in Mycobacterium tuberculosis is essential to curb the spread of tuberculosis (TB). Unfortunately, drug susceptibility testing is currently not available in Papua New Guinea (PNG) and that impairs TB control in this country. We report for the first time M. tuberculosis mutations associated with resistance to first and second-line anti-TB drugs in Madang, PNG. A molecular cluster analysis was performed to identify M. tuberculosis transmission in that region. RESULTS: Phenotypic drug susceptibility tests showed 15.7% resistance to at least one drug and 5.2% multidrug resistant (MDR) TB. Rifampicin resistant strains had the rpoB mutations D516F, D516Y or S531L; Isoniazid resistant strains had the mutations katG S315T or inhA promoter C15T; Streptomycin resistant strains had the mutations rpsL K43R, K88Q, K88R), rrs A514C or gidB V77G. The molecular cluster analysis indicated evidence for transmission of resistant strain. CONCLUSIONS: We observed a substantial rate of MDR-TB in the Madang area of PNG associated with mutations in specific genes. A close monitoring of drug resistance is therefore urgently required, particularly in the presence of drug-resistant M. tuberculosis transmission. In the absence of phenotypic drug susceptibility testing in PNG, molecular assays for drug resistance monitoring would be of advantage.