Targeted resequencing of candidate genes using selector probes

Targeted genome enrichment is a powerful tool for making use of the massive throughput of novel DNA-sequencing instruments. We herein present a simple and scalable protocol for multiplex amplification of target regions based on the Selector technique. The updated version exhibits improved coverage and compatibility with next-generation-sequencing (NGS) library-construction procedures for shotgun sequencing with NGS platforms. To demonstrate the performance of the technique, all 501 exons from 28 genes frequently involved in cancer were enriched for and sequenced in specimens derived from cell lines and tumor biopsies. DNA from both fresh frozen and formalin-fixed paraffin-embedded biopsies were analyzed and 94% specificity and 98% coverage of the targeted region was achieved. Reproducibility between replicates was high (R² = 0, 98) and readily enabled detection of copy-number variations. The procedure can be carried out in <24 h and does not require any dedicated instrumentation.     

Comparison of Methods to Identify Pathogens and Associated Virulence Functional Genes in Biosolids from Two Different Wastewater Treatment Facilities in Canada

The use of treated municipal wastewater residues (biosolids) as fertilizers is an attractive, inexpensive option for growers and farmers. Various regulatory bodies typically employ indicator organisms (fecal coliforms, E. coli and Salmonella) to assess the adequacy and efficiency of the wastewater treatment process in reducing pathogen loads in the final product. Molecular detection approaches can offer some advantages over culture-based methods as they can simultaneously detect a wider microbial species range, including non-cultivable microorganisms. However, they cannot directly assess the viability of the pathogens. Here, we used bacterial enumeration methods together with molecular methods including qPCR, 16S rRNA and cpn60 gene amplicon sequencing and shotgun metagenomic sequencing to compare pre- and post-treatment biosolids from two Canadian wastewater treatment plants (WWTPs). Our results show that an anaerobic digestion WWTP was unsuccessful at reducing the live indicator organism load (coliforms, generic E. coli and Salmonella) below acceptable regulatory criteria, while biosolids from a dewatering/pelletization WWTP met these criteria. DNA from other pathogens was detected by the molecular methods, but these species were considered less abundant. Clostridium DNA increased significantly following anaerobic digestion treatments. In addition to pathogen DNA, genes related to virulence and antibiotic resistance were identified in treated biosolids. Shotgun metagenomics revealed the widest range of pathogen DNA and, among the approaches used here, was the only approach that could access functional gene information in treated biosolids. Overall, our results highlight the potential usefulness of amplicon sequencing and shotgun metagenomics as complementary screening methods that could be used in parallel with culture-based methods, although more detailed comparisons across a wider range of sites would be needed.     

An investigation of conventional microbial culture for the Naja atra bite wound,and the comparison between culture-based 16S Sanger sequencing and 16S metagenomics of the snake oropharyngeal bacterial microbiota

Naja atra is a major venomous snake found in Taiwan. The bite of this snake causes extensive wound necrosis or necrotizing soft tissue infection. Conventional microbial culture-based techniques may fail to identify potential human pathogens and render antibiotics ineffective in the management of wound infection. Therefore, we evaluated 16S Sanger sequencing and next-generation sequencing (NGS) to identify bacterial species in the oropharynx of N. atra. Using conventional microbial culture methods and the VITEK 2 system, we isolated nine species from snakebite wounds. On the basis of the 16S Sanger sequencing of bacterial clones from agar plates, we identified 18 bacterial species in the oropharynx of N. atra, including Morganella morganii, Proteus vulgaris, and Proteus mirabilis, which were also present in the infected bite wound. Using NGS of 16S metagenomics, we uncovered more than 286 bacterial species in the oropharynx of N. atra. In addition, the bacterial species identified using 16S Sanger sequencing accounted for only 2% of those identified through NGS of 16S metagenomics. The bacterial microbiota of the oropharynx of N. atra were modeled better using NGS of 16S metagenomics compared to microbial culture-based techniques. Stenotrophomonas maltophilia, Acinetobacter baumannii, and Proteus penneri were also identified in the NGS of 16S metagenomics. Understanding the bacterial microbiota that are native to the oropharynx of N. atra, in addition to the bite wound, may have additional therapeutic implications regarding empiric antibiotic selection for managing N. atra bites.     

Comparative metagenomic,phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients

Terrestrial ecosystems are receiving elevated inputs of nitrogen (N) from anthropogenic sources and understanding how these increases in N availability affect soil microbial communities is critical for predicting the associated effects on belowground ecosystems. We used a suite of approaches to analyze the structure and functional characteristics of soil microbial communities from replicated plots in two long-term N fertilization experiments located in contrasting systems. Pyrosequencing-based analyses of 16S rRNA genes revealed no significant effects of N fertilization on bacterial diversity, but significant effects on community composition at both sites; copiotrophic taxa (including members of the Proteobacteria and Bacteroidetes phyla) typically increased in relative abundance in the high N plots, with oligotrophic taxa (mainly Acidobacteria) exhibiting the opposite pattern. Consistent with the phylogenetic shifts under N fertilization, shotgun metagenomic sequencing revealed increases in the relative abundances of genes associated with DNA/RNA replication, electron transport and protein metabolism, increases that could be resolved even with the shallow shotgun metagenomic sequencing conducted here (average of 75 000 reads per sample). We also observed shifts in the catabolic capabilities of the communities across the N gradients that were significantly correlated with the phylogenetic and metagenomic responses, indicating possible linkages between the structure and functioning of soil microbial communities. Overall, our results suggest that N fertilization may, directly or indirectly, induce a shift in the predominant microbial life-history strategies, favoring a more active, copiotrophic microbial community, a pattern that parallels the often observed replacement of K-selected with r-selected plant species with elevated N.     

Over 50,000 Metagenomically Assembled Draft Genomes for the Human Oral Microbiome Reveal New Taxa

The oral cavity of each person is home to hundreds of bacterial species. While taxa for oral diseases have been studied using culture-based characterization as well as amplicon sequencing,metagenomic and genomic information remains scarce compared to the fecal microbiome. Here,using metagenomic shotgun data for 3346 oral metagenomic samples together with 808 published samples, we obtain 56,213 metagenome-assembled genomes(MAGs), and more than 64% of the3589 species-level genome bins(SGBs) contai...     

The Validation and Clinical Implementation of BRCAplus: A Comprehensive High-Risk Breast Cancer Diagnostic Assay

Breast cancer is the most commonly diagnosed cancer in women, with 10% of disease attributed to hereditary factors. Although BRCA1 and BRCA2 account for a high percentage of hereditary cases, there are more than 25 susceptibility genes that differentially impact the risk for breast cancer. Traditionally, germline testing for breast cancer was performed by Sanger dideoxy terminator sequencing in a reflexive manner, beginning with BRCA1 and BRCA2. The introduction of next-generation sequencing (NGS) has enabled the simultaneous testing of all genes implicated in breast cancer resulting in diagnostic labs offering large, comprehensive gene panels. However, some physicians prefer to only test for those genes in which established surveillance and treatment protocol exists. The NGS based BRCAplus test utilizes a custom tiled PCR based target enrichment design and bioinformatics pipeline coupled with array comparative genomic hybridization (aCGH) to identify mutations in the six high-risk genes: BRCA1, BRCA2, PTEN, TP53, CDH1, and STK11. Validation of the assay with 250 previously characterized samples resulted in 100% detection of 3,025 known variants and analytical specificity of 99.99%. Analysis of the clinical performance of the first 3,000 BRCAplus samples referred for testing revealed an average coverage greater than 9,000X per target base pair resulting in excellent specificity and the sensitivity to detect low level mosaicism and allele-drop out. The unique design of the assay enabled the detection of pathogenic mutations missed by previous testing. With the abundance of NGS diagnostic tests being released, it is essential that clinicians understand the advantages and limitations of different test designs.     

VITCOMIC2: visualization tool for the phylogenetic composition of microbial communities based on 16S rRNA gene amplicons and metagenomic shotgun sequencing

Background

The 16S rRNA gene-based amplicon sequencing analysis is widely used to determine the taxonomic composition of microbial communities. Once the taxonomic composition of each community is obtained, evolutionary relationships among taxa are inferred by a phylogenetic tree. Thus, the combined representation of taxonomic composition and phylogenetic relationships among taxa is a powerful method for understanding microbial community structure; however, applying phylogenetic tree-based representation with information on the abundance of thousands or more taxa in each community is a difficult task. For this purpose, we previously developed the tool VITCOMIC (VIsualization tool for Taxonomic COmpositions of MIcrobial Community), which is based on the genome-sequenced microbes’ phylogenetic information. Here, we introduce VITCOMIC2, which incorporates substantive improvements over VITCOMIC that were necessary to address several issues associated with 16S rRNA gene-based analysis of microbial communities.

Results

We developed VITCOMIC2 to provide (i) sequence identity searches against broad reference taxa including uncultured taxa; (ii) normalization of 16S rRNA gene copy number differences among taxa; (iii) rapid sequence identity searches by applying the graphics processing unit-based sequence identity search tool CLAST; (iv) accurate taxonomic composition inference and nearly full-length 16S rRNA gene sequence reconstructions for metagenomic shotgun sequencing; and (v) an interactive user interface for simultaneous representation of the taxonomic composition of microbial communities and phylogenetic relationships among taxa. We validated the accuracy of processes (ii) and (iv) by using metagenomic shotgun sequencing data from a mock microbial community.

Conclusions

The improvements incorporated into VITCOMIC2 enable users to acquire an intuitive understanding of microbial community composition based on the 16S rRNA gene sequence data obtained from both metagenomic shotgun and amplicon sequencing.

     

Metagenomic next-generation sequencing as a diagnostic tool for toxoplasmic encephalitis

Background

More than 100 different pathogens can cause encephalitis. Testing of all the neurological pathogens by conventional methods can be difficult. Metagenomic next-generation sequencing (NGS) could identify the infectious agents in a target-independent manner. The role of this novel method in clinical diagnostic microbiology still needs to be evaluated. In present study, we used metagenomic NGS to search for an infectious etiology in a human immunodeficiency virus (HIV)-infected patient with lethally diffuse brain lesions. Sequences mapping to Toxoplasma gondii were unexpectedly detected.

Case presentation

A 31-year-old HIV-infected patient presented to hospital in a critical ill condition with a Glasgow coma scale score of 3. Brain magnetic resonance imaging showed diffuse brain abnormalities with contrast enhancement. Metagenomic NGS was performed on DNA extract from 300 μL patient’s cerebrospinal fluid (CSF) with the BGISEQ-50 platform. The sequencing detection identified 65,357 sequence reads uniquely aligned to the Toxoplasma gondii genome. Presence of Toxoplasma gondii genome in CSF was further verified by Toxoplasma gondii-specific polymerase chain reaction and Sanger sequencing. Altogether, those results confirmed the diagnosis of toxoplasmic encephalitis.

Conclusions

This study suggests that metagenomic NGS may be a useful diagnostic tool for toxoplasmic encephalitis. As metagenomic NGS is able to identify all pathogens in a single run, it may be a promising strategy to explore the clinical causative pathogens in central nervous system infections with atypical features.

     

Targeted single molecule mutation detection with massively parallel sequencing

Next-generation sequencing (NGS) technologies have transformed genomic research and have the potential to revolutionize clinical medicine. However, the background error rates of sequencing instruments and limitations in targeted read coverage have precluded the detection of rare DNA sequence variants by NGS. Here we describe a method, termed CypherSeq, which combines double-stranded barcoding error correction and rolling circle amplification (RCA)-based target enrichment to vastly improve NGS-based rare variant detection. The CypherSeq methodology involves the ligation of sample DNA into circular vectors, which contain double-stranded barcodes for computational error correction and adapters for library preparation and sequencing. CypherSeq is capable of detecting rare mutations genome-wide as well as those within specific target genes via RCA-based enrichment. We demonstrate that CypherSeq is capable of correcting errors incurred during library preparation and sequencing to reproducibly detect mutations down to a frequency of 2.4 × 10⁻⁷ per base pair, and report the frequency and spectra of spontaneous and ethyl methanesulfonate-induced mutations across the Saccharomyces cerevisiae genome.     

Comparative Analysis of Functional Metagenomic Annotation and the Mappability of Short Reads

To assess the functional capacities of microbial communities, including those inhabiting the human body, shotgun metagenomic reads are often aligned to a database of known genes. Such homology-based annotation practices critically rely on the assumption that short reads can map to orthologous genes of similar function. This assumption, however, and the various factors that impact short read annotation, have not been systematically evaluated. To address this challenge, we generated an extremely large database of simulated reads (totaling 15.9 Gb), spanning over 500,000 microbial genes and 170 curated genomes and including, for many genomes, every possible read of a given length. We annotated each read using common metagenomic protocols, fully characterizing the effect of read length, sequencing error, phylogeny, database coverage, and mapping parameters. We additionally rigorously quantified gene-, genome-, and protocol-specific annotation biases. Overall, our findings provide a first comprehensive evaluation of the capabilities and limitations of functional metagenomic annotation, providing crucial goal-specific best-practice guidelines to inform future metagenomic research.     

Metagenomic profiling of fresh water lakes at different altitudes in Pakistan

Freshwater lakes play a crucial role in sustainability of ecosystem in the favour of communities flourishing around it. In this study, differential metagenomic evaluation of three freshwater ecosystems in Pakistan was carried out. These included Lakes Keenjhar, Rawal and Saif-ul-Muluk situated at the altitudes 21, 527, 3224 m respectively. Former two lakes serve as water reservoirs for three major cities of Pakistan, while the high-altitude Lake Saif-ul-Muluk is located in a pristine environment at the meeting point of Himalayan and Karakoram ranges. Next generation sequencing (NGS) of metagenomic DNA of microbial communities was carried out using surface water samples from these lakes by Illumina Hi-Seq 2500 sequencing technology. Bioinformatics analysis identified Proteobacteria as the most dominant phylum (58–79%) followed by Planctomycetes (34%), Cyanobacteria (12%) and Bacteroidetes (15%). Lake Saif-ul-Muluk contained highest abundance of bacterial genera including Vibrio, Bordetella, Pseudomonas, Burkholderia and Escherichia; while Microcystis were found to be abundant in Lakes Keenjhar and Rawal. Non-parametric multi-dimensional scaling (NMDS) principle component analysis on genus level revealed that these ecosystems were divergent from each other and from metagenomes of lakes in the Americas. Characterization of virulence factors and antibiotics resistant genes in metagenomic contigs were also identified. Chemical profiling of Lakes Keenjhar and Rawal showed higher concentration of phosphorus which indicated eutrophication in these lakes might be due to Microcystis infestation.     

Microbial enrichment and multiplexed microfiltration for accelerated detection of Salmonella in spinach

To attain Salmonella detection thresholds in spinach suspensions using enrichment media requires at least 24 hr. Separation and concentration of selected microorganisms via microfiltration and microfugation reduce time for sample preparation, especially when working with large volumes of vegetable suspensions. This facilitates accelerated detection of Salmonella in spinach suspensions, and may contribute to effectively monitoring this pathogen before it reaches the consumer. We report a microfiltration-based protocol for accelerated sample preparation to concentrate and recover ≤1 colony forming unit (CFU) Salmonella/g pathogen-free spinach. Store-bought samples of spinach and a spinach plant subjected to two environmental conditions (temperature and light exposure) during its production were tested. The overall procedure involves extraction with buffer, a short enrichment step, prefiltration using a nylon filter, crossflow hollow fiber microfiltration, and retentate centrifugation to bring microbial cells to detection levels. Based on 1 CFU Salmonella/g frozen spinach, and a Poisson distribution statistical analyses with 99% probability, we calculated that 3 hr of incubation, when followed by microfiltration, is sufficient to reach the 2 log concentration required for Salmonella detection within 7 hr. Longer enrichment times (5 hr or more) is needed for concentrations lower than 1 CFU Salmonella/g of ready to eat spinach. The recovered microbial cells were identified and confirmed as Salmonella using both polymerase chain reaction (PCR) and plating methods. Different environmental conditions tested during production did not affect Salmonella viability; this demonstrated the broad adaptability of Salmonella and emphasized the need for methods that enable efficient monitoring of production for the presence of this pathogen.     

二代测序在脓毒血症患者病原学诊断中的应用

脓毒血症是一种严重威胁生命的感染,精准、快速的病原学诊断可帮助临床医师优化抗菌药物的使用。目前,基于病原菌培养的方法仍是脓毒血症病原学诊断的主要手段,但具有耗时长、灵敏度低等不可忽视的缺点。近年来出现了一些不依赖培养的病原学诊断方法,其中基于聚合酶链反应(polymerase chain reaction,PCR)的方法已发展较为成熟。但PCR只能检测已知的特定病原体,临床定量PCR仅用于检测病毒及少数细菌,脓毒血症中的病原体PCR多仅为定性检测。目前,二代测序技术不断成熟并用于临床,成为病原学诊断的有力手段。与血培养等传统病原学检测方法相比,其具有快速、非选择性、可定量或半定量分析的优点。现阶段二代测序仍存在公认判读标准缺乏、测序结果与治疗关系不明确、耐药基因检测困难等不足,亦缺乏较大规模的二代测序与传统诊断方法比较验证的研究结果,尚有待更高级的循证医学证据支持。     

Shotgun Pyrosequencing Metagenomic Analyses of Dusts from Swine Confinement and Grain Facilities

Inhalation of agricultural dusts causes inflammatory reactions and symptoms such as headache, fever, and malaise, which can progress to chronic airway inflammation and associated diseases, e.g. asthma, chronic bronchitis, chronic obstructive pulmonary disease, and hypersensitivity pneumonitis. Although in many agricultural environments feed particles are the major constituent of these dusts, the inflammatory responses that they provoke are likely attributable to particle-associated bacteria, archaebacteria, fungi, and viruses. In this study, we performed shotgun pyrosequencing metagenomic analyses of DNA from dusts from swine confinement facilities or grain elevators, with comparisons to dusts from pet-free households. DNA sequence alignment showed that 19% or 62% of shotgun pyrosequencing metagenomic DNA sequence reads from swine facility or household dusts, respectively, were of swine or human origin, respectively. In contrast only 2% of such reads from grain elevator dust were of mammalian origin. These metagenomic shotgun reads of mammalian origin were excluded from our analyses of agricultural dust microbiota. The ten most prevalent bacterial taxa identified in swine facility compared to grain elevator or household dust were comprised of 75%, 16%, and 42% gram-positive organisms, respectively. Four of the top five swine facility dust genera were assignable (Clostridium, Lactobacillus, Ruminococcus, and Eubacterium, ranging from 4% to 19% relative abundance). The relative abundances of these four genera were lower in dust from grain elevators or pet-free households. These analyses also highlighted the predominance in swine facility dust of Firmicutes (70%) at the phylum level, Clostridia (44%) at the Class level, and Clostridiales at the Order level (41%). In summary, shotgun pyrosequencing metagenomic analyses of agricultural dusts show that they differ qualitatively and quantitatively at the level of microbial taxa present, and that the bioinformatic analyses used for such studies must be carefully designed to avoid the potential contribution of non-microbial DNA, e.g. from resident mammals.     

Comparison of Sample Preparation Methods Used for the Next-Generation Sequencing of Mycobacterium tuberculosis

The advent and widespread application of next-generation sequencing (NGS) technologies to the study of microbial genomes has led to a substantial increase in the number of studies in which whole genome sequencing (WGS) is applied to the analysis of microbial genomic epidemiology. However, microorganisms such as Mycobacterium tuberculosis (MTB) present unique problems for sequencing and downstream analysis based on their unique physiology and the composition of their genomes. In this study, we compare the quality of sequence data generated using the Nextera and TruSeq isolate preparation kits for library construction prior to Illumina sequencing-by-synthesis. Our results confirm that MTB NGS data quality is highly dependent on the purity of the DNA sample submitted for sequencing and its guanine-cytosine content (or GC-content). Our data additionally demonstrate that the choice of library preparation method plays an important role in mitigating downstream sequencing quality issues. Importantly for MTB, the Illumina TruSeq library preparation kit produces more uniform data quality than the Nextera XT method, regardless of the quality of the input DNA. Furthermore, specific genomic sequence motifs are commonly missed by the Nextera XT method, as are regions of especially high GC-content relative to the rest of the MTB genome. As coverage bias is highly undesirable, this study illustrates the importance of appropriate protocol selection when performing NGS studies in order to ensure that sound inferences can be made regarding mycobacterial genomes.     

Identification and Resolution of Microdiversity through Metagenomic Sequencing of Parallel Consortia

To gain a predictive understanding of the interspecies interactions within microbial communities that govern community function, the genomic complement of every member population must be determined. Although metagenomic sequencing has enabled the de novo reconstruction of some microbial genomes from environmental communities, microdiversity confounds current genome reconstruction techniques. To overcome this issue, we performed short-read metagenomic sequencing on parallel consortia, defined as consortia cultivated under the same conditions from the same natural community with overlapping species composition. The differences in species abundance between the two consortia allowed reconstruction of near-complete (at an estimated >85% of gene complement) genome sequences for 17 of the 20 detected member species. Two Halomonas spp. indistinguishable by amplicon analysis were found to be present within the community. In addition, comparison of metagenomic reads against the consensus scaffolds revealed within-species variation for one of the Halomonas populations, one of the Rhodobacteraceae populations, and the Rhizobiales population. Genomic comparison of these representative instances of inter- and intraspecies microdiversity suggests differences in functional potential that may result in the expression of distinct roles in the community. In addition, isolation and complete genome sequence determination of six member species allowed an investigation into the sensitivity and specificity of genome reconstruction processes, demonstrating robustness across a wide range of sequence coverage (9× to 2,700×) within the metagenomic data set.     

Metagenomic analysis of historical herbarium specimens reveals a postmortem microbial community

Advances in DNA extraction and next‐generation sequencing have made a vast number of historical herbarium specimens available for genomic investigation. These specimens contain not only genomic information from the individual plants themselves, but also from associated microorganisms such as bacteria and fungi. These microorganisms may have colonized the living plant (e.g., pathogens or host‐associated commensal taxa) or may result from postmortem colonization that may include decomposition processes or contamination during sample handling. Here we characterize the metagenomic profile from shotgun sequencing data from herbarium specimens of two widespread plant species (Ambrosia artemisiifolia and Arabidopsis thaliana) collected up to 180 years ago. We used blast searching in combination with megan and were able to infer the metagenomic community even from the oldest herbarium sample. Through comparison with contemporary plant collections, we identify three microbial species that are nearly exclusive to herbarium specimens, including the fungus Alternaria alternata, which can comprise up to 7% of the total sequencing reads. This species probably colonizes the herbarium specimens during preparation for mounting or during storage. By removing the probable contaminating taxa, we observe a temporal shift in the metagenomic composition of the invasive weed Am. artemisiifolia. Our findings demonstrate that it is generally possible to use herbarium specimens for metagenomic analyses, but that the results should be treated with caution, as some of the identified species may be herbarium contaminants rather than representing the natural metagenomic community of the host plant.     

Quantitative and Sensitive Detection of GNAS Mutations Causing McCune-Albright Syndrome with Next Generation Sequencing

Somatic activating GNAS mutations cause McCune-Albright syndrome (MAS). Owing to low mutation abundance, mutant-specific enrichment procedures, such as the peptide nucleic acid (PNA) method, are required to detect mutations in peripheral blood. Next generation sequencing (NGS) can analyze millions of PCR amplicons independently, thus it is expected to detect low-abundance GNAS mutations quantitatively. In the present study, we aimed to develop an NGS-based method to detect low-abundance somatic GNAS mutations. PCR amplicons encompassing exons 8 and 9 of GNAS, in which most activating mutations occur, were sequenced on the MiSeq instrument. As expected, our NGS-based method could sequence the GNAS locus with very high read depth (approximately 100,000) and low error rate. A serial dilution study with use of cloned mutant and wildtype DNA samples showed a linear correlation between dilution and measured mutation abundance, indicating the reliability of quantification of the mutation. Using the serially diluted samples, the detection limits of three mutation detection methods (the PNA method, NGS, and combinatory use of PNA and NGS [PNA-NGS]) were determined. The lowest detectable mutation abundance was 1% for the PNA method, 0.03% for NGS and 0.01% for PNA-NGS. Finally, we analyzed 16 MAS patient-derived leukocytic DNA samples with the three methods, and compared the mutation detection rate of them. Mutation detection rate of the PNA method, NGS and PNA-NGS in 16 patient-derived peripheral blood samples were 56%, 63% and 75%, respectively. In conclusion, NGS can detect somatic activating GNAS mutations quantitatively and sensitively from peripheral blood samples. At present, the PNA-NGS method is likely the most sensitive method to detect low-abundance GNAS mutation.