Reference materials for clinical applications of human genome sequencing
The Genome in a Bottle Consortium is developing the reference materials, reference methods , and reference data n
The Genome in a Bottle Consortium is developing the reference materials, reference methods , and reference data n
Purpose: Clinical laboratories, research laboratories and technology developers all need DNA samples with reliably known genotypes in order to help validate and improve their methods. The Genome in a Bottle Consortium (genomeinabottle.org) has been developing Reference Materials with high-accuracy whole genome sequences to support these efforts.Methodology: Our pilot reference material is based on Coriell sample NA12878 and was released in May 2015 as NIST RM 8398 (tinyurl.com/giabpilot). To minimize bias and improve accuracy, 11 whole-genome and 3 exome data sets produced using 5 different technologies were integrated using a systematic arbitration method [1]. The Genome in a Bottle Analysis Group is adapting these methods and developing new methods to characterize 2 families, one Asian and one Ashkenazi Jewish from the Personal Genome Project, which are consented for public release of sequencing and phenotype data. We have generated a larger and even more diverse data set on these samples, including high-depth Illumina paired-end and mate-pair, Complete Genomics, and Ion Torrent short-read data, as well as Moleculo, 10X, Oxford Nanopore, PacBio, and BioNano Genomics long-read data. We are analyzing these data to provide an accurate assessment of not just small variants but also large structural variants (SVs) in both “easy” regions of the genome and in some “hard” repetitive regions. We have also made all of the input data sources publicly available for download, analysis, and publication.Results: Our arbitration method produced a reference data set of 2,787,291 single nucleotide variants (SNVs), 365,135 indels, 2744 SVs, and 2.2 billion homozygous reference calls for our pilot genome. We found that our call set is highly sensitive and specific in comparison to independent reference data sets. We have also generated preliminary assemblies and structural variant calls for the next 2 trios from long read data and are currently integrating and validating these.Discussion: We combined the strengths of each of our input datasets to develop a comprehensive and accurate benchmark call set. In the short time it has been available, over 20 published or submitted papers have used our data. Many challenges exist in comparing to our benchmark calls, and thus we have worked with the Global Alliance for Genomics and Health to develop standardized methods, performance metrics, and software to assist in its use.[1] Zook et al, Nat Biotech. 2014.
At AGBT 2017, the Broad Institute’s Daniel Neafsey reported a large collaborative effort to sequence the mosquito that carries Zika virus. The team is using long-read PacBio sequencing to produce…
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.
The antimicrobial resistance (AMR) crisis represents a serious threat to public health and has resulted in concentrated efforts to accelerate development of rapid molecular diagnostics for AMR. In combination with publicly-available web-based AMR databases, whole genome sequencing (WGS) offers the capacity for rapid detection of antibiotic resistance genes. Here we studied the concordance between WGS-based resistance prediction and phenotypic susceptibility testing results for methicillin-resistant Staphylococcus aureus (MRSA) and vancomycin resistant Enterococcus (VRE) clinical isolates using publicly-available tools and databases.Clinical isolates prospectively collected at the University of Pittsburgh Medical Center between December 2016 and December 2017 underwent WGS. Antibiotic resistance gene content was assessed from assembled genomes by BLASTn search of online databases. Concordance between WGS-predicted resistance profile and phenotypic susceptibility as well as sensitivity, specificity, positive and negative predictive values (NPV, PPV) were calculated for each antibiotic/organism combination, using the phenotypic results as the gold standard.Phenotypic susceptibility testing and WGS results were available for 1242 isolate/antibiotic combinations. Overall concordance was 99.3% with a sensitivity, specificity, PPV, NPV of 98.7% (95% CI, 97.2-99.5%), 99.6% (95 % CI, 98.8-99.9%), 99.3% (95% CI, 98.0-99.8%), 99.2% (95% CI, 98.3-99.7%), respectively. Additional identification of point mutations in housekeeping genes increased the concordance to 99.4% and the sensitivity to 99.3% (95% CI, 98.2-99.8%) and NPV to 99.4% (95% CI, 98.4-99.8%).WGS can be used as a reliable predicator of phenotypic resistance for both MRSA and VRE using readily-available online tools.Copyright © 2019. Published by Elsevier Ltd.
The evolution and global transmission of antimicrobial resistance has been well documented in Gram-negative bacteria and healthcare-associated epidemic pathogens, often emerging from regions with heavy antimicrobial use. However, the degree to which similar processes occur with Gram-positive bacteria in the community setting is less well understood. Here, we trace the recent origins and global spread of a multidrug resistant, community-associated Staphylococcus aureus lineage from the Indian subcontinent, the Bengal Bay clone (ST772). We generated whole genome sequence data of 340 isolates from 14 countries, including the first isolates from Bangladesh and India, to reconstruct the evolutionary history and genomic epidemiology of the lineage. Our data shows that the clone emerged on the Indian subcontinent in the early 1970s and disseminated rapidly in the 1990s. Short-term outbreaks in community and healthcare settings occurred following intercontinental transmission, typically associated with travel and family contacts on the subcontinent, but ongoing endemic transmission was uncommon. Acquisition of a multidrug resistance integrated plasmid was instrumental in the divergence of a single dominant and globally disseminated clade in the early 1990s. Phenotypic data on biofilm, growth and toxicity point to antimicrobial resistance as the driving force in the evolution of ST772. The Bengal Bay clone therefore combines the multidrug resistance of traditional healthcare-associated clones with the epidemiological transmission of community-associated MRSA. Our study demonstrates the importance of whole genome sequencing for tracking the evolution of emerging and resistant pathogens. It provides a critical framework for ongoing surveillance of the clone on the Indian subcontinent and elsewhere.Importance The Bengal Bay clone (ST772) is a community-acquired and multidrug-resistant Staphylococcus aureus lineage first isolated from Bangladesh and India in 2004. In this study, we show that the Bengal Bay clone emerged from a virulent progenitor circulating on the Indian subcontinent. Its subsequent global transmission was associated with travel or family contact in the region. ST772 progressively acquired specific resistance elements at limited cost to its fitness and continues to be exported globally resulting in small-scale community and healthcare outbreaks. The Bengal Bay clone therefore combines the virulence potential and epidemiology of community-associated clones with the multidrug-resistance of healthcare-associated S. aureus lineages. This study demonstrates the importance of whole genome sequencing for the surveillance of highly antibiotic resistant pathogens, which may emerge in the community setting of regions with poor antibiotic stewardship and rapidly spread into hospitals and communities across the world.
Microbacterium sp. strain SGAir0570 was isolated from air samples collected in Singapore. Its genome was assembled using single-molecule real-time sequencing and MiSeq short reads. It has one chromosome with a length of 3.38?Mb and one 59.2-kb plasmid. It contains 3,170 protein-coding genes, 48 tRNAs, and 6 rRNAs.Copyright © 2019 Kalsi et al.
The Gram-negative pathogenic spirochetal bacteria Leptospira spp. cause leptospirosis in humans and livestock animals. Leptospira kmetyi strain LS 001/16 was isolated from a soil sample associated with a leptospirosis patient in Kelantan, which is among the states in Malaysia with a high reported number of disease cases. Here, we report the complete genome sequence of Leptospira kmetyi strain LS 001/16. Copyright © 2019 Yusof et al.
Several emerging pathogens have arisen as a result of selection pressures exerted by modern health care. Klebsiella quasipneumoniae was recently defined as a new species, yet its prevalence, niche, and propensity to acquire antimicrobial resistance genes are not fully described. We have been tracking inter- and intraspecies transmission of the Klebsiella pneumoniae carbapenemase (KPC) gene, blaKPC, between bacteria isolated from a single institution. We applied a combination of Illumina and PacBio whole-genome sequencing to identify and compare K. quasipneumoniae from patients and the hospital environment over 10- and 5-year periods, respectively. There were 32 blaKPC-positive K. quasipneumoniae isolates, all of which were identified as K. pneumoniae in the clinical microbiology laboratory, from 8 patients and 11 sink drains, with evidence for seven separate blaKPC plasmid acquisitions. Analysis of a single subclade of K. quasipneumoniae subsp. quasipneumoniae (n?=?23 isolates) from three patients and six rooms demonstrated seeding of a sink by a patient, subsequent persistence of the strain in the hospital environment, and then possible transmission to another patient. Longitudinal analysis of this strain demonstrated the acquisition of two unique blaKPC plasmids and then subsequent within-strain genetic rearrangement through transposition and homologous recombination. Our analysis highlights the apparent molecular propensity of K. quasipneumoniae to persist in the environment as well as acquire carbapenemase plasmids from other species and enabled an assessment of the genetic rearrangements which may facilitate horizontal transmission of carbapenemases. Copyright © 2019 Mathers et al.
Four Enterobacteriaceae clinical isolates bearing mcr-1 gene-harboring plasmids were characterized. All isolates demonstrated the ability to transfer colistin resistance to Escherichia coli; plasmids were stable in conjugants after multiple passages on nonselective media. mcr-1 was located on an IncX4 (n?=?3) or IncN (n?=?1) plasmid. The IncN plasmid harbored 13 additional antimicrobial resistance genes. Results indicate that the mcr-1-bearing plasmids in this study were highly transferable in vitro and stable in the recipients.This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.
New Delhi metallo-beta-lactamases (NDMs) are an uncommon but emerging cause of carbapenem resistance in the United States. Genomic factors promoting their domestic spread remain poorly characterized. A prospective genomic surveillance program among Boston-area hospitals identified multiple new occurrences of NDM-carrying strains of Escherichia coli and Enterobacter cloacae complex in inpatient and outpatient settings, representing the first occurrences of NDM-mediated resistance since initiating genomic surveillance in 2011. Cases included domestic patients with no international exposures. PacBio sequencing of isolates identified strain characteristics, resistance genes, and the complement of mobile vectors mediating spread. Analyses revealed a common 3,114-bp region containing the blaNDM gene, with carriage of this conserved region among unique strains by diverse transposon and plasmid backbones. Functional studies revealed a broad capacity for blaNDM transmission by conjugation, transposition, and complex interplasmid recombination events. NDMs represent a rapidly spreading form of drug resistance that can occur in inpatient and outpatient settings and in patients without international exposures. In contrast to Tn4401-based spread of Klebsiella pneumoniae carbapenemases (KPCs), diverse transposable elements mobilize NDM enzymes, commonly with other resistance genes, enabling naive strains to acquire multi- and extensively drug-resistant profiles with single transposition or plasmid conjugation events. Genomic surveillance provides effective means to rapidly identify these gene-level drivers of resistance and mobilization in order to inform clinical decisions to prevent further spread.Copyright © 2019 American Society for Microbiology.
Targeted PCR amplification and high-throughput sequencing (amplicon sequencing) of 16S rRNA gene fragments is widely used to profile microbial communities. New long-read sequencing technologies can sequence the entire 16S rRNA gene, but higher error rates have limited their attractiveness when accuracy is important. Here we present a high-throughput amplicon sequencing methodology based on PacBio circular consensus sequencing and the DADA2 sample inference method that measures the full-length 16S rRNA gene with single-nucleotide resolution and a near-zero error rate. In two artificial communities of known composition, our method recovered the full complement of full-length 16S sequence variants from expected community members without residual errors. The measured abundances of intra-genomic sequence variants were in the integral ratios expected from the genuine allelic variants within a genome. The full-length 16S gene sequences recovered by our approach allowed Escherichia coli strains to be correctly classified to the O157:H7 and K12 sub-species clades. In human fecal samples, our method showed strong technical replication and was able to recover the full complement of 16S rRNA alleles in several E. coli strains. There are likely many applications beyond microbial profiling for which high-throughput amplicon sequencing of complete genes with single-nucleotide resolution will be of use. © The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.
We report the complete genome sequence for Streptococcus pyogenes strain 1838 (type emm3) isolated from a patient with toxic shock syndrome. The strain lacked the speK- and sla-encoding prophage frequently encountered among emm3 strains and possessed an Arg66His mutation in CovR of the 2-component virulence regulatory system CovRS.
We characterized 170 complete genome assemblies from clinical Bordetella pertussis isolates representing geographic and temporal diversity in the United States. These data capture genotypic shifts, including increased pertactin deficiency, occurring amid the current pertussis disease resurgence and provide a foundation for needed research to direct future public health control strategies.
Historically, the major cause of meningococcal epidemics in the meningitis belt of sub-Saharan Africa has been Neisseria meningitidis serogroup A (NmA), but the incidence has been substantially reduced since the introduction of a serogroup A conjugate vaccine starting in 2010. We performed whole-genome sequencing on isolates collected post-2010 to assess their phylogenetic relationships and inter-country transmission.A total of 716 invasive meningococcal isolates collected between 2011 and 2016 from 11 meningitis belt countries were whole-genome sequenced for molecular characterization by the three WHO Collaborating Centers for Meningitis.We identified three previously-reported clonal complexes (CC): CC11 (n?=?434), CC181 (n?=?62) and CC5 (n?=?90) primarily associated with NmW, NmX, and NmA, respectively, and an emerging CC10217 (n?=?126) associated with NmC. CC11 expanded throughout the meningitis belt independent of the 2000 Hajj outbreak strain, with isolates from Central African countries forming a distinct sub-lineage within this expansion. Two major sub-lineages were identified for CC181 isolates, one mainly expanding in West African countries and the other found in Chad. CC10217 isolates from the large outbreaks in Nigeria and Niger were more closely related than those from the few cases in Mali and Burkina Faso.Whole-genome based phylogenies revealed geographically distinct strain circulation as well as inter-country transmission events. Our results stress the importance of continued meningococcal molecular surveillance in the region, as well as the development of an affordable vaccine targeting these strains. FUND: Meningitis Research Foundation; CDC’s Office of Advanced Molecular Detection; GAVI, the Vaccine Alliance. Copyright © 2019. Published by Elsevier B.V.
If you have a question, need to check the status of an order, or are interested in purchasing an instrument, we're here to help.