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

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.Funding for this work was provided by The Wellcome Trust (investigator award no. 107032/Z/15/Z to R.A.F.), Fondation Botnar (Programme grant no. 6063) and the UK CF Trust (Innovation Hub award no. 001; Strategic Research Centre award no. 010). M.S., N.A.H. and R.M.D. acknowledge the Cystic Fibrosis Foundation for funding

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

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

Detection of Low Frequency Multi-Drug Resistance and Novel Putative Maribavir Resistance in Immunocompromised Pediatric Patients with Cytomegalovirus.

Human cytomegalovirus (HCMV) is a significant pathogen in immunocompromised individuals, with the potential to cause fatal pneumonitis and colitis, as well as increasing the risk of organ rejection in transplant patients. With the advent of new anti-HCMV drugs there is therefore considerable interest in using virus sequence data to monitor emerging resistance to antiviral drugs in HCMV viraemia and disease, including the identification of putative new mutations. We used target-enrichment to deep sequence HCMV DNA from 11 immunosuppressed pediatric patients receiving single or combination anti-HCMV treatment, serially sampled over 1-27 weeks. Changes in consensus sequence and resistance mutations were analyzed for three ORFs targeted by anti-HCMV drugs and the frequencies of drug resistance mutations monitored. Targeted-enriched sequencing of clinical material detected mutations occurring at frequencies of 2%. Seven patients showed no evidence of drug resistance mutations. Four patients developed drug resistance mutations a mean of 16 weeks after starting treatment. In two patients, multiple resistance mutations accumulated at frequencies of 20% or less, including putative maribavir and ganciclovir resistance mutations P522Q (UL54) and C480F (UL97). In one patient, resistance was detected 14 days earlier than by PCR. Phylogenetic analysis suggested recombination or superinfection in one patient. Deep sequencing of HCMV enriched from clinical samples excluded resistance in 7 of 11 subjects and identified resistance mutations earlier than conventional PCR-based resistance testing in 2 patients. Detection of multiple low level resistance mutations was associated with poor outcome

Two Generations of CubeSat Missions (CSSWE and CIRBE) to Take on the Challenges of Measuring Relativistic Electrons in the Earth’s Magnetosphere

The Colorado Student Space Weather Experiment (CSSWE) CubeSat, carrying the Relativistic Electron and Proton Telescope integrated little experiment (REPTile) to measure 0.5 to \u3e3.8 MeV electrons and 8-40 MeV protons, operated for over two years, 2012-2014, in low Earth orbit (LEO). There have been 25 peer-reviewed publications, including two in Nature, and five Ph.D. dissertations associated with CSSWE. Another 3U CubeSat mission: Colorado Inner Radiation Belt Electron Experiment (CIRBE), has been under development to address an unresolved science question: Where is the break point in terms of electron energy below which electrons can be transported into the inner belt from the outer belt but above which they cannot? This requires clean measurements of energetic electrons with fine energy resolution in an environment where all instruments are subject to the unforgiving penetration from highly energetic protons (tens of MeV to GeV). An advanced version of REPTile has been designed and built, REPTile-2. It has been integrated into the CIRBE bus, which has active attitude control, deployable solar panels, and a S-band radio, provided by Blue Canyon Technologies. CIRBE advances our science capabilities and has significantly improved performance vs. CSSWE and is ready to be launched into a LEO in early 2023

Stepwise pathogenic evolution of Mycobacterium abscessus.

Although almost all mycobacterial species are saprophytic environmental organisms, a few, such as Mycobacterium tuberculosis, have evolved to cause transmissible human infection. By analyzing the recent emergence and spread of the environmental organism M. abscessus through the global cystic fibrosis population, we have defined key, generalizable steps involved in the pathogenic evolution of mycobacteria. We show that epigenetic modifiers, acquired through horizontal gene transfer, cause saltational increases in the pathogenic potential of specific environmental clones. Allopatric parallel evolution during chronic lung infection then promotes rapid increases in virulence through mutations in a discrete gene network; these mutations enhance growth within macrophages but impair fomite survival. As a consequence, we observe constrained pathogenic evolution while person-to-person transmission remains indirect, but postulate accelerated pathogenic adaptation once direct transmission is possible, as observed for M. tuberculosis Our findings indicate how key interventions, such as early treatment and cross-infection control, might restrict the spread of existing mycobacterial pathogens and prevent new, emergent ones

Genomic epidemiology of a national outbreak of post-surgical Mycobacterium abscessus wound infections in Brazil.

An epidemic of post-surgical wound infections, caused by a non-tuberculous mycobacterium, has been on-going in Brazil. It has been unclear whether one or multiple lineages are responsible and whether their wide geographical distribution across Brazil is due to spread from a single point source or is the result of human-mediated transmission. 188 isolates, collected from nine Brazilian states, were whole genome sequenced and analysed using phylogenetic and comparative genomic approaches. The isolates from Brazil formed a single clade, which was estimated to have emerged in 2003. We observed temporal and geographic structure within the lineage that enabled us to infer the movement of sub-lineages across Brazil. The genome size of the Brazilian lineage was reduced relative to most strains in the three subspecies of Mycobacterium abscessus and contained a novel plasmid, pMAB02, in addition to the previously described pMAB01 plasmid. One lineage, which emerged just prior to the initial outbreak, is responsible for the epidemic of post-surgical wound infections in Brazil. Phylogenetic analysis indicates that multiple transmission events led to its spread. The presence of a novel plasmid and the reduced genome size suggest that the lineage has undergone adaptation to the surgical niche

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