Commentary from PacBio users on their applications of SMRT Sequencing, including Ulf Gyllensten (Uppsala University), Tim Smith (USDA-ARS) and Bobby Sebra (Icahn School of Medicine)
Webinar: Complete genomes within reach – Closing bacterial genomes from the lakes of Minnesota to NYC hospitals
In this webinar, Ben Auch, Research Scientist, Innovation Lab, University of Minnesota Genomics Center, Cody Sheik, Assistant Professor of Biology, University of Minnesota Duluth, and Harm van Bakel, Assistant Professor…
Characterization of Reference Materials for Genetic Testing of CYP2D6 Alleles: A GeT-RM Collaborative Project.
Pharmacogenetic testing increasingly is available from clinical and research laboratories. However, only a limited number of quality control and other reference materials currently are available for the complex rearrangements and rare variants that occur in the CYP2D6 gene. To address this need, the Division of Laboratory Systems, CDC-based Genetic Testing Reference Material Coordination Program, in collaboration with members of the pharmacogenetic testing and research communities and the Coriell Cell Repositories (Camden, NJ), has characterized 179 DNA samples derived from Coriell cell lines. Testing included the recharacterization of 137 genomic DNAs that were genotyped in previous Genetic Testing Reference Material Coordination Program studies and 42 additional samples that had not been characterized previously. DNA samples were distributed to volunteer testing laboratories for genotyping using a variety of commercially available and laboratory-developed tests. These publicly available samples will support the quality-assurance and quality-control programs of clinical laboratories performing CYP2D6 testing.Published by Elsevier Inc.
A comparison of immunoglobulin IGHV, IGHD and IGHJ genes in wild-derived and classical inbred mouse strains.
The genomes of classical inbred mouse strains include genes derived from all three major subspecies of the house mouse, Mus musculus. We recently posited that genetic diversity in the immunoglobulin heavy chain (IGH) gene loci of C57BL/6 and BALB/c mice reflect differences in subspecies origin. To investigate this hypothesis, we conducted high-throughput sequencing of IGH gene rearrangements to document IGH variable (IGHV), joining (IGHJ), and diversity (IGHD) genes in four inbred wild-derived mouse strains (CAST/EiJ, LEWES/EiJ, MSM/MsJ, and PWD/PhJ), and a single disease model strain (NOD/ShiLtJ), collectively representing genetic backgrounds of several major mouse subspecies. A total of 341 germline IGHV sequences were inferred in the wild-derived strains, including 247 not curated in the International Immunogenetics Information System. In contrast, 83/84 inferred NOD IGHV genes had previously been observed in C57BL/6 mice. Variability among the strains examined was observed for only a single IGHJ gene, involving a description of a novel allele. In contrast, unexpected variation was found in the IGHD gene loci, with four previously unreported IGHD gene sequences being documented. Very few IGHV sequences of C57BL/6 and BALB/c mice were shared with strains representing major subspecies, suggesting that their IGH loci may be complex mosaics of genes of disparate origins. This suggests a similar level of diversity is likely present in the IGH loci of other classical inbred strains. This must now be documented if we are to properly understand inter-strain variation in models of antibody-mediated disease. This article is protected by copyright. All rights reserved.This article is protected by copyright. All rights reserved.
Complete genome screening of clinical MRSA isolates identifies lineage diversity and provides full resolution of transmission and outbreak events
Whole-genome sequencing (WGS) of Staphylococcus aureus is increasingly used as part of infection prevention practices, but most applications are focused on conserved core genomic regions due to limitations of short-read technologies. In this study we established a long-read technology-based WGS screening program of all first-episode MRSA blood infections at a major urban hospital. A survey of 132 MRSA genomes assembled from long reads revealed widespread gain/loss of accessory mobile genetic elements among established hospital- and community-associated lineages impacting >10% of each genome, and frequent megabase-scale inversions between endogenous prophages. We also characterized an outbreak of a CC5/ST105/USA100 clone among 3 adults and 18 infants in a neonatal intensive care unit (NICU) lasting 7 months. The pattern of changes among complete outbreak genomes provided full spatiotemporal resolution of its origins and progression, which was characterized by multiple sub-transmissions and likely precipitated by equipment sharing. Compared to other hospital strains, the outbreak strain carried distinct mutations and accessory genetic elements that impacted genes with roles in metabolism, resistance and persistence. This included a DNA-recognition domain recombination in the hsdS gene of a Type-I restriction-modification system that altered DNA methylation. RNA-Seq profiling showed that the (epi)genetic changes in the outbreak clone attenuated agr gene expression and upregulated genes involved in stress response and biofilm formation. Overall our findings demonstrate that long-read sequencing substantially improves our ability to characterize accessory genomic elements that impact MRSA virulence and persistence, and provides valuable information for infection control efforts.
Dissemination of multiple carbapenem resistance genes in an in vitro gut model simulating the human colon.
Carbapenemase-producing Enterobacteriaceae (CPE) pose a major global health risk. Mobile genetic elements account for much of the increasing CPE burden.To investigate CPE colonization and the impact of antibiotic exposure on subsequent resistance gene dissemination within the gut microbiota using a model to simulate the human colon.Gut models seeded with CPE-negative human faeces [screened with BioMérieux chromID® CARBA-SMART (Carba-Smart), Cepheid Xpert® Carba-R assay (XCR)] were inoculated with distinct carbapenemase-producing Klebsiella pneumoniae strains (KPC, NDM) and challenged with imipenem or piperacillin/tazobactam then meropenem. Resistant populations were enumerated daily on selective agars (Carba-Smart); CPE genes were confirmed by PCR (XCR, Check-Direct CPE Screen for BD MAX™). CPE gene dissemination was tracked using PacBio long-read sequencing.CPE populations increased during inoculation, plateauing at ~105?log10?cfu/mL in both models and persisting throughout the experiments (>65?days), with no evidence of CPE ‘washout’. After antibiotic administration, there was evidence of interspecies plasmid transfer of blaKPC-2 (111742?bp IncFII/IncR plasmid, 99% identity to pKpQIL-D2) and blaNDM-1 (~170?kb IncFIB/IncFII plasmid), and CPE populations rose from <0.01% to >45% of the total lactose-fermenting populations in the KPC model. Isolation of a blaNDM-1K. pneumoniae with one chromosomal single-nucleotide variant compared with the inoculated strain indicated clonal expansion within the model. Antibiotic administration exposed a previously undetected K. pneumoniae encoding blaOXA-232 (KPC model).CPE exposure can lead to colonization, clonal expansion and resistance gene transfer within intact human colonic microbiota. Furthermore, under antibiotic selective pressure, new resistant populations emerge, emphasizing the need to control exposure to antimicrobials. © The Author(s) 2019. Published by Oxford University Press on behalf of the British Society for Antimicrobial Chemotherapy. All rights reserved. For permissions, please email: firstname.lastname@example.org.
Single-molecule long-read sequencing datasets were generated for a son-father-mother trio of Han Chinese descent that is part of the Genome in a Bottle (GIAB) consortium portfolio. The dataset was generated using the Pacific Biosciences Sequel System. The son and each parent were sequenced to an average coverage of 60 and 30, respectively, with N50 subread lengths between 16 and 18?kb. Raw reads and reads aligned to both the GRCh37 and GRCh38 are available at the NCBI GIAB ftp site (ftp://ftp-trace.ncbi.nlm.nih.gov/giab/ftp/data/ChineseTrio/). The GRCh38 aligned read data are archived in NCBI SRA (SRX4739017, SRX4739121, and SRX4739122). This dataset is available for anyone to develop and evaluate long-read bioinformatics methods.
HIV-1 infection of primary CD4(+) T cells regulates the expression of specific HERV-K (HML-2) elements.
Endogenous retroviruses (ERVs) occupy extensive regions of the human genome. Although many of these retroviral elements have lost their ability to replicate, those whose insertion took place more recently, such as the HML-2 group of HERV-K elements, still retain intact open reading frames and the capacity to produce certain viral RNA and/or proteins. Transcription of these ERVs is, however, tightly regulated by dedicated epigenetic control mechanisms. Nonetheless, it has been reported that some pathologic states, such as viral infections and certain cancers, coincide with ERV expression suggesting transcriptional reawakening is possible. HML-2 elements are reportedly induced during HIV-1 infection, but the conserved nature of these elements has, until recently, rendered their expression profiling problematic.Here, we provide comprehensive HERV-K HML-2 expression profiles specific for productively HIV-1 infected primary human CD4(+) T cells. We combined enrichment of HIV-1 infected cells using a reporter virus expressing a surface reporter for gentle and efficient purification with long-read Single Molecule Real-Time sequencing. We show that three HML-2 proviruses, 6q25.1, 8q24.3, and 19q13.42 are up-regulated on average between 3- and 5-fold in HIV-1 infected CD4(+) T cells. One provirus, HML-2 12q24.33, in contrast, was repressed in the presence of active HIV replication.In conclusion, this report identifies the HERV-K HML-2 loci whose expression profiles differ upon HIV-1 infection in primary human CD4(+) T cells. These data will help pave the way for further studies on the influence of endogenous retroviruses on HIV-1 replication.Importance Endogenous retroviruses inhabit big portions of our genome. And although they are mainly inert some of the evolutionarily younger members maintain the ability to express both RNA as well as proteins. We have developed an approach using long-read SMRT sequencing that produces long reads, that provides us with ability to obtain detailed and accurate HERV-K HML-2 expression profiles. We have now applied this approach to study HERV-K expression in the presence and absence of productive HIV-1 infection of primary human CD4(+) T cells. In addition to using SMRT sequencing, our strategy also includes the magnetic selection of the infected cells so that levels of background expression due to uninfected cells are kept at a minimum. The results in this manuscript provide the blueprint for in-depth studies of the interactions of the authentic upregulated HERV-K HML-2 elements and HIV-1. Copyright © 2017 American Society for Microbiology.
Shotgun metagenomics methods enable characterization of microbial communities in human microbiome and environmental samples. Assembly of metagenome sequences does not output whole genomes, so computational binning methods have been developed to cluster sequences into genome ‘bins’. These methods exploit sequence composition, species abundance, or chromosome organization but cannot fully distinguish closely related species and strains. We present a binning method that incorporates bacterial DNA methylation signatures, which are detected using single-molecule real-time sequencing. Our method takes advantage of these endogenous epigenetic barcodes to resolve individual reads and assembled contigs into species- and strain-level bins. We validate our method using synthetic and real microbiome sequences. In addition to genome binning, we show that our method links plasmids and other mobile genetic elements to their host species in a real microbiome sample. Incorporation of DNA methylation information into shotgun metagenomics analyses will complement existing methods to enable more accurate sequence binning.
Improved OTU-picking using long-read 16S rRNA gene amplicon sequencing and generic hierarchical clustering
BACKGROUND: High-throughput bacterial 16S rRNA gene sequencing followed by clustering of short sequences into operational taxonomic units (OTUs) is widely used for microbiome profiling. However, clustering of short 16S rRNA gene reads into biologically meaningful OTUs is challenging, in part because nucleotide variation along the 16S rRNA gene is only partially captured by short reads. The recent emergence of long-read platforms, such as single-molecule real-time (SMRT) sequencing from Pacific Biosciences, offers the potential for improved taxonomic and phylogenetic profiling. Here, we evaluate the performance of long- and short-read 16S rRNA gene sequencing using simulated and experimental data, followed by OTU inference using computational pipelines based on heuristic and complete-linkage hierarchical clustering. RESULTS: In simulated data, long-read sequencing was shown to improve OTU quality and decrease variance. We then profiled 40 human gut microbiome samples using a combination of Illumina MiSeq and Blautia-specific SMRT sequencing, further supporting the notion that long reads can identify additional OTUs. We implemented a complete-linkage hierarchical clustering strategy using a flexible computational pipeline, tailored specifically for PacBio circular consensus sequencing (CCS) data that outperforms heuristic methods in most settings: https://github.com/oscar-franzen/oclust/. CONCLUSION: Our data demonstrate that long reads can improve OTU inference; however, the choice of clustering algorithm and associated clustering thresholds has significant impact on performance.
Little is known about the role of the microbiome in primary sclerosing cholangitis.To explore the mucosa-associated microbiota in primary sclerosing cholangitis (PSC) patients across different locations in the gut, and to compare it with inflammatory bowel disease (IBD)-only patients and healthy controls.Biopsies from the terminal ileum, right colon, and left colon were collected from patients and healthy controls undergoing colonoscopy. Microbiota profiling using bacterial 16S rRNA sequencing was performed on all biopsies.Forty-four patients were recruited: 20 with PSC (19 with PSC-IBD and one with PSC-only), 15 with IBD-only and nine healthy controls. The overall microbiome profile was similar throughout different locations in the gut. No differences in the global microbiome profile were found. However, we observed significant PSC-associated enrichment in Barnesiellaceae at the family level, and in Blautia and an unidentified Barnesiellaceae at the genus level. At the operational taxa unit level, most shifts in PSC were observed in Clostridiales and Bacteroidales orders, with approximately 86% of shifts occurring within the former order.The overall microbiota profile was similar across multiple locations in the gut from the same individual regardless of disease status. In this study, the mucosa associated-microbiota of patients with primary sclerosing cholangitis was characterised by enrichment of Blautia and Barnesiellaceae and by major shifts in operational taxa units within Clostridiales order.© 2016 John Wiley & Sons Ltd.
Iso-Seq analysis of Nepenthes ampullaria, Nepenthes rafflesiana and Nepenthes × hookeriana for hybridisation study in pitcher plants.
Tropical pitcher plants in the species-rich Nepenthaceae family of carnivorous plants possess unique pitcher organs. Hybridisation, natural or artificial, in this family is extensive resulting in pitchers with diverse features. The pitcher functions as a passive insect trap with digestive fluid for nutrient acquisition in nitrogen-poor habitats. This organ shows specialisation according to the dietary habit of different Nepenthes species. In this study, we performed the first single-molecule real-time isoform sequencing (Iso-Seq) analysis of full-length cDNA from Nepenthes ampullaria which can feed on leaf litter, compared to carnivorous Nepenthes rafflesiana, and their carnivorous hybrid Nepenthes × hookeriana. This allows the comparison of pitcher transcriptomes from the parents and the hybrid to understand how hybridisation could shape the evolution of dietary habit in Nepenthes. Raw reads have been deposited to SRA database with the accession numbers SRX2692198 (N. ampullaria), SRX2692197 (N. rafflesiana), and SRX2692196 (N. × hookeriana).
Methylation patterns can help determine the origin of metagenomic sequences.
Genetic basis of emerging vancomycin, linezolid, and daptomycin heteroresistance in a case of persistent Enterococcus faecium bacteremia.
Whole-genome sequencing was used to examine a persistent Enterococcus faecium bacteremia that acquired heteroresistance to three antibiotics in response to prolonged multidrug therapy. A comparison of the complete genomes before and after each change revealed the emergence of known resistance determinants for vancomycin and linezolid and suggested that a novel mutation in fabF, encoding a fatty acid synthase, was responsible for daptomycin nonsusceptibility. Plasmid recombination contributed to the progressive loss of vancomycin resistance after withdrawal of the drug. Copyright © 2018 Chacko et al.
During cefoxitin-based nasal screening, phenotypically categorized methicillin-resistant Staphylococcus aureus (MRSA) was isolated and tested negative for the presence of the mecA and mecC genes as well as for the SCCmec-orfX junction region. The isolate was found to carry a mecB gene previously described for Macrococcus caseolyticus but not for staphylococcal species. The gene is flanked by ß-lactam regulatory genes similar to mecR, mecI, and blaZ and is part of an 84.6-kb multidrug-resistance plasmid that harbors genes encoding additional resistances to aminoglycosides (aacA-aphD, aphA, and aadK) as well as macrolides (ermB) and tetracyclines (tetS). This further plasmidborne ß-lactam resistance mechanism harbors the putative risk of acceleration or reacceleration of MRSA spread, resulting in broad ineffectiveness of ß-lactams as a main therapeutic application against staphylococcal infections.