Lipoarabinomannan modification as a source of phenotypic heterogeneity in host-adapted Mycobacterium abscessus isolates.

Mycobacterium abscessus is increasingly recognized as the causative agent of chronic pulmonary infections in humans. One of the genes found to be under strong evolutionary pressure during adaptation of M. abscessus to the human lung is embC which encodes an arabinosyltransferase required for the biosynthesis of the cell envelope lipoglycan, lipoarabinomannan (LAM). To assess the impact of patient-derived embC mutations on the physiology and virulence of M. abscessus, mutations were introduced in the isogenic background of M. abscessus ATCC 19977 and the resulting strains probed for phenotypic changes in a variety of in vitro and host cell-based assays relevant to infection. We show that patient-derived mutational variations in EmbC result in an unexpectedly large number of changes in the physiology of M. abscessus, and its interactions with innate immune cells. Not only did the mutants produce previously unknown forms of LAM with a truncated arabinan domain and 3-linked oligomannoside chains, they also displayed significantly altered cording, sliding motility, and biofilm-forming capacities. The mutants further differed from wild-type M. abscessus in their ability to replicate and induce inflammatory responses in human monocyte-derived macrophages and epithelial cells. The fact that different embC mutations were associated with distinct physiologic and pathogenic outcomes indicates that structural alterations in LAM caused by nonsynonymous nucleotide polymorphisms in embC may be a rapid, one-step, way for M. abscessus to generate broad-spectrum diversity beneficial to survival within the heterogeneous and constantly evolving environment of the infected human airway.

Human cytomegalovirus haplotype reconstruction reveals high diversity due to superinfection and evidence of within-host recombination.

Recent sequencing efforts have led to estimates of human cytomegalovirus (HCMV) genome-wide intrahost diversity that rival those of persistent RNA viruses [Renzette N, Bhattacharjee B, Jensen JD, Gibson L, Kowalik TF (2011) PLoS Pathog 7:e1001344]. Here, we deep sequence HCMV genomes recovered from single and longitudinally collected blood samples from immunocompromised children to show that the observations of high within-host HCMV nucleotide diversity are explained by the frequent occurrence of mixed infections caused by genetically distant strains. To confirm this finding, we reconstructed within-host viral haplotypes from short-read sequence data. We verify that within-host HCMV nucleotide diversity in unmixed infections is no greater than that of other DNA viruses analyzed by the same sequencing and bioinformatic methods and considerably less than that of human immunodeficiency and hepatitis C viruses. By resolving individual viral haplotypes within patients, we reconstruct the timing, likely origins, and natural history of superinfecting strains. We uncover evidence for within-host recombination between genetically distinct HCMV strains, observing the loss of the parental virus containing the nonrecombinant fragment. The data suggest selection for strains containing the recombinant fragment, generating testable hypotheses about HCMV evolution and pathogenesis. These results highlight that high HCMV diversity present in some samples is caused by coinfection with multiple distinct strains and provide reassurance that within the host diversity for single-strain HCMV infections is no greater than for other herpesviruses.

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.

Direct Whole-Genome Sequencing of Sputum Accurately Identifies Drug-Resistant Mycobacterium tuberculosis Faster than MGIT Culture Sequencing.

The current methods available to diagnose antimicrobial-resistant Mycobacterium tuberculosis infections require a positive culture or only test a limited number of resistance-associated mutations. A rapid accurate identification of antimicrobial resistance enables the prompt initiation of effective treatment. Here, we determine the utility of whole-genome sequencing (WGS) of M. tuberculosis directly from routinely obtained diagnostic sputum samples to provide a comprehensive resistance profile compared to that from mycobacterial growth indicator tube (MGIT) WGS. We sequenced M. tuberculosis from 43 sputum samples by targeted DNA enrichment using the Agilent SureSelectXT kit, and 43 MGIT positive samples from each participant. Thirty two (74%) sputum samples and 43 (100%) MGIT samples generated whole genomes. The times to antimicrobial resistance profiles and concordance were compared with Xpert MTB/RIF and phenotypic resistance testing from cultures of the same samples. Antibiotic susceptibility could be predicted from WGS of sputum within 5 days of sample receipt and up to 24 days earlier than WGS from MGIT culture and up to 31 days earlier than phenotypic testing. Direct sputum results could be reduced to 3 days with faster hybridization and if only regions encoding drug resistance are sequenced. We show that direct sputum sequencing has the potential to provide comprehensive resistance detection significantly faster than MGIT whole-genome sequencing or phenotypic testing of resistance from cultures in a clinical setting. This improved turnaround time enables prompt appropriate treatment with associated patient and health service benefits. Improvements in sample preparation are necessary to ensure comparable sensitivities and complete resistance profile predictions in all cases.

Genomic analysis of the causative agents of coccidiosis in domestic chickens.

Global production of chickens has trebled in the past two decades and they are now the most important source of dietary animal protein worldwide. Chickens are subject to many infectious diseases that reduce their performance and productivity. Coccidiosis, caused by apicomplexan protozoa of the genus Eimeria, is one of the most important poultry diseases. Understanding the biology of Eimeria parasites underpins development of new drugs and vaccines needed to improve global food security. We have produced annotated genome sequences of all seven species of Eimeria that infect domestic chickens, which reveal the full extent of previously described repeat-rich and repeat-poor regions and show that these parasites possess the most repeat-rich proteomes ever described. Furthermore, while no other apicomplexan has been found to possess retrotransposons, Eimeria is home to a family of chromoviruses. Analysis of Eimeria genes involved in basic biology and host-parasite interaction highlights adaptations to a relatively simple developmental life cycle and a complex array of co-expressed surface proteins involved in host cell binding.

Phylogenomic exploration of the relationships between strains of Mycobacterium avium subspecies paratuberculosis.

BACKGROUND: Mycobacterium avium subspecies paratuberculosis (Map) is an infectious enteric pathogen that causes Johne's disease in livestock. Determining genetic diversity is prerequisite to understanding the epidemiology and biology of Map. We performed the first whole genome sequencing (WGS) of 141 global Map isolates that encompass the main molecular strain types currently reported. We investigated the phylogeny of the Map strains, the diversity of the genome and the limitations of commonly used genotyping methods. RESULTS: Single nucleotide polymorphism (SNP) and phylogenetic analyses confirmed two major lineages concordant with the former Type S and Type C designations. The Type I and Type III strain groups are subtypes of Type S, and Type B strains are a subtype of Type C and not restricted to Bison species. We found that the genome-wide SNPs detected provided greater resolution between isolates than currently employed genotyping methods. Furthermore, the SNP used for IS1311 typing is not informative, as it is likely to have occurred after Type S and C strains diverged and does not assign all strains to the correct lineage. Mycobacterial Interspersed Repetitive Unit-Variable Number Tandem Repeat (MIRU-VNTR) differentiates Type S from Type C but provides limited resolution between isolates within these lineages and the polymorphisms detected do not necessarily accurately reflect the phylogenetic relationships between strains. WGS of passaged strains and coalescent analysis of the collection revealed a very high level of genetic stability, with the substitution rate estimated to be less than 0.5 SNPs per genome per year. CONCLUSIONS: This study clarifies the phylogenetic relationships between the previously described Map strain groups, and highlights the limitations of current genotyping techniques. Map isolates exhibit restricted genetic diversity and a substitution rate consistent with a monomorphic pathogen. WGS provides the ultimate level of resolution for differentiation between strains. However, WGS alone will not be sufficient for tracing and tracking Map infections, yet importantly it can provide a phylogenetic context for affirming epidemiological connections.

Novel Single Nucleotide Polymorphism-Based Assay for Genotyping Mycobacterium avium subsp. paratuberculosis.

Typing of Mycobacterium avium subspecies paratuberculosis strains presents a challenge, since they are genetically monomorphic and traditional molecular techniques have limited discriminatory power. The recent advances and availability of whole-genome sequencing have extended possibilities for the characterization of Mycobacterium avium subspecies paratuberculosis, and whole-genome sequencing can provide a phylogenetic context to facilitate global epidemiology studies. In this study, we developed a single nucleotide polymorphism (SNP) assay based on PCR and restriction enzyme digestion or sequencing of the amplified product. The SNP analysis was performed using genome sequence data from 133 Mycobacterium avium subspecies paratuberculosis isolates with different genotypes from 8 different host species and 17 distinct geographic regions around the world. A total of 28,402 SNPs were identified among all of the isolates. The minimum number of SNPs required to distinguish between all of the 133 genomes was 93 and between only the type C isolates was 41. To reduce the number of SNPs and PCRs required, we adopted an approach based on sequential detection of SNPs and a decision tree. By the analysis of 14 SNPs Mycobacterium avium subspecies paratuberculosis isolates can be characterized within 14 phylogenetic groups with a higher discriminatory power than mycobacterial interspersed repetitive unit-variable number tandem repeat assay and other typing methods. Continuous updating of genome sequences is needed in order to better characterize new phylogenetic groups and SNP profiles. The novel SNP assay is a discriminative, simple, reproducible method and requires only basic laboratory equipment for the large-scale global typing of Mycobacterium avium subspecies paratuberculosis isolates.

Genomic Investigations unmask Mycoplasma amphoriforme, a new respiratory pathogen.

BACKGROUND: Mycoplasma amphoriforme has been associated with infection in patients with primary antibody deficiency (PAD). Little is known about the natural history of infection with this organism and its ability to be transmitted in the community. METHODS: The bacterial load was estimated in sequential sputum samples from 9 patients by quantitative polymerase chain reaction. The genomes of all available isolates, originating from patients in the United Kingdom, France, and Tunisia, were sequenced along with the type strain. Genomic data were assembled and annotated, and a high-resolution phylogenetic tree was constructed. RESULTS: By using high-resolution whole-genome sequencing (WGS) data, we show that patients can be chronically infected with M. amphoriforme manifesting as a relapsing-remitting bacterial load, interspersed by periods when the organism is undetectable. Importantly, we demonstrate transmission of strains within a clinical environment. Antibiotic resistance mutations accumulate in isolates taken from patients who received multiple courses of antibiotics. CONCLUSIONS: Mycoplasma amphoriforme isolates form a closely related species responsible for a chronic relapsing and remitting infection in PAD patients in the United Kingdom and from immunocompetent patients in other countries. We provide strong evidence of transmission between patients attending the same clinic, suggesting that screening and isolation may be necessary for susceptible patients. This work demonstrates the critical role that WGS can play in rapidly unraveling the biology of a novel pathogen.

Mycobacterium abscessus Complex Identification with Matrix-Assisted Laser Desorption Ionization-Time of Flight Mass Spectrometry.

We determined that the Vitek MS Plus matrix-assisted laser desorption ionization-time of flight mass spectrometry using research-use-only (RUO) v.4.12 and in vitro-diagnostic (IVD) v.3.0 databases accurately identified 41 Mycobacterium abscessus subsp. abscessus and 13 M. abscessus subsp. massiliense isolates identified by whole-genome sequencing to the species but not the subspecies level, from Middlebrook 7H11 and Burkholderia cepacia selective agars. Peak analysis revealed three peaks potentially able to differentiate between subspecies.

Rapid Whole-Genome Sequencing of Mycobacterium tuberculosis Isolates Directly from Clinical Samples.

The rapid identification of antimicrobial resistance is essential for effective treatment of highly resistant Mycobacterium tuberculosis. Whole-genome sequencing provides comprehensive data on resistance mutations and strain typing for monitoring transmission, but unlike for conventional molecular tests, this has previously been achievable only from cultures of M. tuberculosis. Here we describe a method utilizing biotinylated RNA baits designed specifically for M. tuberculosis DNA to capture full M. tuberculosis genomes directly from infected sputum samples, allowing whole-genome sequencing without the requirement of culture. This was carried out on 24 smear-positive sputum samples, collected from the United Kingdom and Lithuania where a matched culture sample was available, and 2 samples that had failed to grow in culture. M. tuberculosis sequencing data were obtained directly from all 24 smear-positive culture-positive sputa, of which 20 were of high quality (>20× depth and >90% of the genome covered). Results were compared with those of conventional molecular and culture-based methods, and high levels of concordance between phenotypical resistance and predicted resistance based on genotype were observed. High-quality sequence data were obtained from one smear-positive culture-negative case. This study demonstrated for the first time the successful and accurate sequencing of M. tuberculosis genomes directly from uncultured sputa. Identification of known resistance mutations within a week of sample receipt offers the prospect for personalized rather than empirical treatment of drug-resistant tuberculosis, including the use of antimicrobial-sparing regimens, leading to improved outcomes.

Whole-genome sequencing to identify transmission of Mycobacterium abscessus between patients with cystic fibrosis: a retrospective cohort study.

BACKGROUND: Increasing numbers of individuals with cystic fibrosis are becoming infected with the multidrug-resistant non-tuberculous mycobacterium (NTM) Mycobacterium abscessus, which causes progressive lung damage and is extremely challenging to treat. How this organism is acquired is not currently known, but there is growing concern that person-to-person transmission could occur. We aimed to define the mechanisms of acquisition of M abscessus in individuals with cystic fibrosis. METHOD: Whole genome sequencing and antimicrobial susceptibility testing were done on 168 consecutive isolates of M abscessus from 31 patients attending an adult cystic fibrosis centre in the UK between 2007 and 2011. In parallel, we undertook detailed environmental testing for NTM and defined potential opportunities for transmission between patients both in and out of hospital using epidemiological data and social network analysis. FINDINGS: Phylogenetic analysis revealed two clustered outbreaks of near-identical isolates of the M abscessus subspecies massiliense (from 11 patients), differing by less than ten base pairs. This variation represents less diversity than that seen within isolates from a single individual, strongly indicating between-patient transmission. All patients within these clusters had numerous opportunities for within-hospital transmission from other individuals, while comprehensive environmental sampling, initiated during the outbreak, failed to detect any potential point source of NTM infection. The clusters of M abscessus subspecies massiliense showed evidence of transmission of mutations acquired during infection of an individual to other patients. Thus, isolates with constitutive resistance to amikacin and clarithromycin were isolated from several individuals never previously exposed to long-term macrolides or aminoglycosides, further indicating cross-infection. INTERPRETATION: Whole genome sequencing has revealed frequent transmission of multidrug resistant NTM between patients with cystic fibrosis despite conventional cross-infection measures. Although the exact transmission route is yet to be established, our epidemiological analysis suggests that it could be indirect. FUNDING: The Wellcome Trust, Papworth Hospital, NIHR Cambridge Biomedical Research Centre, UK Health Protection Agency, Medical Research Council, and the UKCRC Translational Infection Research Initiative.

Atypical at skew in Firmicute genomes results from selection and not from mutation.

The second parity rule states that, if there is no bias in mutation or selection, then within each strand of DNA complementary bases are present at approximately equal frequencies. In bacteria, however, there is commonly an excess of G (over C) and, to a lesser extent, T (over A) in the replicatory leading strand. The low G+C Firmicutes, such as Staphylococcus aureus, are unusual in displaying an excess of A over T on the leading strand. As mutation has been established as a major force in the generation of such skews across various bacterial taxa, this anomaly has been assumed to reflect unusual mutation biases in Firmicute genomes. Here we show that this is not the case and that mutation bias does not explain the atypical AT skew seen in S. aureus. First, recently arisen intergenic SNPs predict the classical replication-derived equilibrium enrichment of T relative to A, contrary to what is observed. Second, sites predicted to be under weak purifying selection display only weak AT skew. Third, AT skew is primarily associated with largely non-synonymous first and second codon sites and is seen with respect to their sense direction, not which replicating strand they lie on. The atypical AT skew we show to be a consequence of the strong bias for genes to be co-oriented with the replicating fork, coupled with the selective avoidance of both stop codons and costly amino acids, which tend to have T-rich codons. That intergenic sequence has more A than T, while at mutational equilibrium a preponderance of T is expected, points to a possible further unresolved selective source of skew.

Evolution and host-specific adaptation of Pseudomonas aeruginosa.

The major human bacterial pathogen Pseudomonas aeruginosa causes multidrug-resistant infections in people with underlying immunodeficiencies or structural lung diseases such as cystic fibrosis (CF). We show that a few environmental isolates, driven by horizontal gene acquisition, have become dominant epidemic clones that have sequentially emerged and spread through global transmission networks over the past 200 years. These clones demonstrate varying intrinsic propensities for infecting CF or non-CF individuals (linked to specific transcriptional changes enabling survival within macrophages); have undergone multiple rounds of convergent, host-specific adaptation; and have eventually lost their ability to transmit between different patient groups. Our findings thus explain the pathogenic evolution of P. aeruginosa and highlight the importance of global surveillance and cross-infection prevention in averting the emergence of future epidemic clones.

Rapid whole-genome sequencing for investigation of a suspected tuberculosis outbreak.

Two Southeast Asian students attending the same school in the United Kingdom presented with pulmonary tuberculosis. An epidemiological investigation failed to link the two cases, and drug resistance profiles of the Mycobacterium tuberculosis isolates were discrepant. Whole-genome sequencing of the isolates found them to be genetically identical, suggesting a missed transmission event.

Inferring patient to patient transmission of Mycobacterium tuberculosis from whole genome sequencing data.

BACKGROUND: Mycobacterium tuberculosis is characterised by limited genomic diversity, which makes the application of whole genome sequencing particularly attractive for clinical and epidemiological investigation. However, in order to confidently infer transmission events, an accurate knowledge of the rate of change in the genome over relevant timescales is required. METHODS: We attempted to estimate a molecular clock by sequencing 199 isolates from epidemiologically linked tuberculosis cases, collected in the Netherlands spanning almost 16 years. RESULTS: Multiple analyses support an average mutation rate of ~0.3 SNPs per genome per year. However, all analyses revealed a very high degree of variation around this mean, making the confirmation of links proposed by epidemiology, and inference of novel links, difficult. Despite this, in some cases, the phylogenetic context of other strains provided evidence supporting the confident exclusion of previously inferred epidemiological links. CONCLUSIONS: This in-depth analysis of the molecular clock revealed that it is slow and variable over short time scales, which limits its usefulness in transmission studies. However, the superior resolution of whole genome sequencing can provide the phylogenetic context to allow the confident exclusion of possible transmission events previously inferred via traditional DNA fingerprinting techniques and epidemiological cluster investigation. Despite the slow generation of variation even at the whole genome level we conclude that the investigation of tuberculosis transmission will benefit greatly from routine whole genome sequencing.

Clinically relevant mutations in the PhoR sensor kinase of host-adapted Mycobacterium abscessus isolates impact response to acidic pH and virulence.

IMPORTANCE: Difficult-to-treat pulmonary infections caused by nontuberculous mycobacteria of the Mycobacterium abscessus group have been steadily increasing in the USA and globally. Owing to the relatively recent recognition of M. abscessus as a human pathogen, basic and translational research to address critical gaps in diagnosis, treatment, and prevention of diseases caused by this microorganism has been lagging behind that of the better-known mycobacterial pathogen, Mycobacterium tuberculosis. To begin unraveling the molecular mechanisms of pathogenicity of M. abscessus, we here focus on the study of a two-component regulator known as PhoPR which we found to be under strong evolutionary pressure during human lung infection. We show that PhoPR is activated at acidic pH and serves to regulate a defined set of genes involved in host adaptation. Accordingly, clinical isolates from chronically infected human lungs tend to hyperactivate this regulator enabling M. abscessus to escape macrophage killing.

PheWAS and cross-disorder analysis reveal genetic architecture, pleiotropic loci and phenotypic correlations across 11 autoimmune disorders.

Introduction: Autoimmune disorders (ADs) are a group of about 80 disorders that occur when self-attacking autoantibodies are produced due to failure in the self-tolerance mechanisms. ADs are polygenic disorders and associations with genes both in the human leukocyte antigen (HLA) region and outside of it have been described. Previous studies have shown that they are highly comorbid with shared genetic risk factors, while epidemiological studies revealed associations between various lifestyle and health-related phenotypes and ADs. Methods: Here, for the first time, we performed a comparative polygenic risk score (PRS) - Phenome Wide Association Study (PheWAS) for 11 different ADs (Juvenile Idiopathic Arthritis, Primary Sclerosing Cholangitis, Celiac Disease, Multiple Sclerosis, Rheumatoid Arthritis, Psoriasis, Myasthenia Gravis, Type 1 Diabetes, Systemic Lupus Erythematosus, Vitiligo Late Onset, Vitiligo Early Onset) and 3,254 phenotypes available in the UK Biobank that include a wide range of socio-demographic, lifestyle and health-related outcomes. Additionally, we investigated the genetic relationships of the studied ADs, calculating their genetic correlation and conducting cross-disorder GWAS meta-analyses for the observed AD clusters. Results: In total, we identified 508 phenotypes significantly associated with at least one AD PRS. 272 phenotypes were significantly associated after excluding variants in the HLA region from the PRS estimation. Through genetic correlation and genetic factor analyses, we identified four genetic factors that run across studied ADs. Cross-trait meta-analyses within each factor revealed pleiotropic genome-wide significant loci. Discussion: Overall, our study confirms the association of different factors with genetic susceptibility for ADs and reveals novel observations that need to be further explored.

Mycobacterium abscessus pathogenesis identified by phenogenomic analyses.

The medical and scientific response to emerging and established pathogens is often severely hampered by ignorance of the genetic determinants of virulence, drug resistance and clinical outcomes that could be used to identify therapeutic drug targets and forecast patient trajectories. Taking the newly emergent multidrug-resistant bacteria Mycobacterium abscessus as an example, we show that combining high-dimensional phenotyping with whole-genome sequencing in a phenogenomic analysis can rapidly reveal actionable systems-level insights into bacterial pathobiology. Through phenotyping of 331 clinical isolates, we discovered three distinct clusters of isolates, each with different virulence traits and associated with a different clinical outcome. We combined genome-wide association studies with proteome-wide computational structural modelling to define likely causal variants, and employed direct coupling analysis to identify co-evolving, and therefore potentially epistatic, gene networks. We then used in vivo CRISPR-based silencing to validate our findings and discover clinically relevant M. abscessus virulence factors including a secretion system, thus illustrating how phenogenomics can reveal critical pathways within emerging pathogenic bacteria.

Dissemination of Mycobacterium abscessus via global transmission networks.

Mycobacterium abscessus, a multidrug-resistant nontuberculous mycobacterium, has emerged as a major pathogen affecting people with cystic fibrosis (CF). Although originally thought to be acquired independently from the environment, most individuals are infected with one of several dominant circulating clones (DCCs), indicating the presence of global transmission networks of M. abscessus. How and when these clones emerged and spread globally is unclear. Here, we use evolutionary analyses of isolates from individuals both with and without CF to reconstruct the population history, spatiotemporal spread and recent transmission networks of the DCCs. We demonstrate synchronous expansion of six unrelated DCCs in the 1960s, a period associated with major changes in CF care and survival. Each of these clones has spread globally as a result of rare intercontinental transmission events. We show that the DCCs, but not environmentally acquired isolates, exhibit a specific smoking-associated mutational signature and that current transmission networks include individuals both with and without CF. We therefore propose that the DCCs initially emerged in non-CF populations but were then amplified and spread through the CF community. While individuals with CF are probably the most permissive host, non-CF individuals continue to play a key role in transmission networks and may facilitate long-distance transmission.