High-resolution metagenomics targets specific functional types in complex microbial communities

Most microbes in the biosphere remain unculturable. Whole genome shotgun (WGS) sequencing of environmental DNA (metagenomics) can be used to study the genetic and metabolic properties of natural microbial communities. However, in communities of high complexity, metagenomics fails to link specific microbes to specific ecological functions. To overcome this limitation, we developed a method to target microbial subpopulations by labeling DNA through stable isotope probing (SIP), followed by WGS sequencing. Metagenome analysis of microbes from Lake Washington in Seattle that oxidize single-carbon (C1) compounds shows specific sequence enrichments in response to different C1 substrates, revealing the ecological roles of individual phylotypes. We also demonstrate the utility of our approach by extracting a nearly complete genome of a novel methylotroph, Methylotenera mobilis, reconstructing its metabolism and conducting genome-wide analyses. This high-resolution, targeted metagenomics approach may be applicable to a wide variety of ecosystems.     

Exploration of Genetic Variations through Single‐cell Whole‐genome Sequencing in the Model Ciliate Tetrahymena thermophila

Ciliates are unicellular eukaryotes with separate germline and somatic genomes and diverse life cycles, which make them a unique model to improve our understanding of population genetics through the detection of genetic variations. However, traditional sequencing methods cannot be directly applied to ciliates because the majority are uncultivated. Single‐cell whole‐genome sequencing (WGS) is a powerful tool for studying genetic variation in microbes, but no studies have been performed in ciliates. We compared the use of single‐cell WGS and bulk DNA WGS to detect genetic variation, specifically single nucleotide polymorphisms (SNPs), in the model ciliate Tetrahymena thermophila. Our analyses showed that (i) single‐cell WGS has excellent performance regarding mapping rate and genome coverage but lower sequencing uniformity compared with bulk DNA WGS due to amplification bias (which was reproducible); (ii) false‐positive SNP sites detected by single‐cell WGS tend to occur in genomic regions with particularly high sequencing depth and high rate of C:G to T:A base changes; (iii) SNPs detected in three or more cells should be reliable (an detection efficiency of 83.4–97.4% was obtained for combined data from three cells). This analytical method could be adapted to measure genetic variation in other ciliates and broaden research into ciliate population genetics.     

Microbiome of Pukzing Cave in India shows high antimicrobial activity against plant and animal pathogens

Pukzing cave, the largest cave of Mizoram, India was explored for bacterial diversity. Culture dependent method revealed 235 bacterial isolates using three different treatments. Identity of the microbial species was confirmed by 16S rDNA sequencing. The highest bacterial population was recovered from heat treatment (n = 97;41.2%) followed by normal (n = 79;33.6%) and cold treatment (n = 59;25.1%) indicating dominance of moderate thermophiles. Antimicrobial potential of isolates showed 20.4% isolates having antimicrobial ability against tested pathogens. Amplicon sequencing of PKSI, PKSII and NRP specific genes revealed presence of AMP genes in the microbial population. Six microbial pathogens were selected for screening as they are well known for different disease cause organism in various fields such as agriculture and human health. Cave environment harbors unique microbial flora and hypervariable region V4 is more informative. Higher activity of AMP assay against these microbes indicates that cave microbial communities could be potential source of future genomic resources.     

Intra-genomic rRNA gene variability of Nassellaria and Spumellaria (Rhizaria,Radiolaria) assessed by Sanger,MinION and Illumina sequencing

Ribosomal RNA (rRNA) genes are known to be valuable markers for the barcoding of eukaryotic life and its phylogenetic classification at various taxonomic levels. The large-scale exploration of environmental microbial diversity through metabarcoding approaches has been focused mainly on the V4 and V9 regions of the 18S rRNA gene. The accurate interpretation of such environmental surveys is hampered by technical (e.g. PCR and sequencing errors) and biological biases (e.g. intra-genomic variability). Here we explored the intra-genomic diversity of Nassellaria and Spumellaria specimens (Radiolaria) by comparing Sanger sequencing with Illumina and Oxford Nanopore Technologies (MinION). Our analysis determined that intra-genomic variability of Nassellaria and Spumellaria is generally low, yet some Spumellaria specimens showed two different copies of the V4 with <97% similarity. Of the different sequencing methods, Illumina showed the highest number of contaminations (i.e. environmental DNA, cross-contamination, tag-jumping), revealed by its high sequencing depth; and MinION showed the highest sequencing rate error (~14%). Yet the long reads produced by MinION (~2900 bp) allowed accurate phylogenetic reconstruction studies. These results highlight the requirement for a careful interpretation of Illumina-based metabarcoding studies, in particular regarding low abundant amplicons, and open future perspectives towards full-length rDNA environmental metabarcoding surveys.     

Polyphasic approach of bacterial classification — An overview of recent advances

Classification of microorganisms on the basis of traditional microbiological methods (morphological, physiological and biochemical) creates a blurred image about their taxonomic status and thus needs further clarification. It should be based on a more pragmatic approach of deploying a number of methods for the complete characterization of microbes. Hence, the methods now employed for bacterial systematics include, the complete 16S rRNA gene sequencing and its comparative analysis by phylogenetic trees, DNA-DNA hybridization studies with related organisms, analyses of molecular markers and signature pattern(s), biochemical assays, physiological and morphological tests. Collectively these genotypic, chemotaxonomic and phenotypic methods for determining taxonomic position of microbes constitute what is known as the ‘polyphasic approach’ for bacterial systematics. This approach is currently the most popular choice for classifying bacteria and several microbes, which were previously placed under invalid taxa have now been resolved into new genera and species. This has been possible owing to rapid development in molecular biological techniques, automation of DNA sequencing coupled with advances in bioinformatic tools and access to sequence databases. Several DNA-based typing methods are known; these provide information for delineating bacteria into different genera and species and have the potential to resolve differences among the strains of a species. Therefore, newly isolated strains must be classified on the basis of the polyphasic approach. Also previously classified organisms, as and when required, can be reclassified on this ground in order to obtain information about their accurate position in the microbial world. Thus, current techniques enable microbiologists to decipher the natural phylogenetic relationships between microbes.     

Selective whole genome amplification and sequencing of Coxiella burnetii directly from environmental samples

Whole genome sequencing (WGS) is a widely available, inexpensive means of providing a wealth of information about an organism's diversity and evolution. However, WGS for many pathogenic bacteria remain limited because they are difficult, slow and/or dangerous to culture. To avoid culturing, metagenomic sequencing can be performed directly on samples, but the sequencing effort required to characterize low frequency organisms can be expensive. Recently developed methods for selective whole genome amplification (SWGA) can enrich target DNA to provide efficient sequencing. We amplified Coxiella burnetii (a bacterial select agent and human/livestock pathogen) from 3 three environmental samples that were overwhelmed with host DNA. The 68- to 147-fold enrichment of the bacterial sequences provided enough genome coverage for SNP analyses and phylogenetic placement. SWGA is a valuable tool for the study of difficult-to-culture organisms and has the potential to facilitate high-throughput population characterizations as well as targeted epidemiological or forensic investigations.     

Succession of Microbial Communities during Hot Composting as Detected by PCR–Single-Strand-Conformation Polymorphism-Based Genetic Profiles of Small-Subunit rRNA Genes

A cultivation-independent technique for genetic profiling of PCR-amplified small-subunit rRNA genes (SSU rDNA) was chosen to characterize the diversity and succession of microbial communities during composting of an organic agricultural substrate. PCR amplifications were performed with DNA directly extracted from compost samples and with primers targeting either (i) the V4–V5 region of eubacterial 16S rRNA genes, (ii) the V3 region in the 16S rRNA genes of actinomycetes, or (iii) the V8–V9 region of fungal 18S rRNA genes. Homologous PCR products were converted to single-stranded DNA molecules by exonuclease digestion and were subsequently electrophoretically separated by their single-strand-conformation polymorphism (SSCP). Genetic profiles obtained by this technique showed a succession and increasing diversity of microbial populations with all primers. A total of 19 single products were isolated from the profiles by PCR reamplification and cloning. DNA sequencing of these molecular isolates showed similarities in the range of 92.3 to 100% to known gram-positive bacteria with a low or high G+C DNA content and to the SSU rDNA of γ-Proteobacteria. The amplified 18S rRNA gene sequences were related to the respective gene regions of Candida krusei and Candida tropicalis. Specific molecular isolates could be attributed to different composting stages. The diversity of cultivated bacteria isolated from samples taken at the end of the composting process was low. A total of 290 isolates were related to only 6 different species. Two or three of these species were also detectable in the SSCP community profiles. Our study indicates that community SSCP profiles can be highly useful for the monitoring of bacterial diversity and community successions in a biotechnologically relevant process.     

Microbial diversity and function in soil: from genes to ecosystems

Soils sustain an immense diversity of microbes, which, to a large extent, remains unexplored. A range of novel methods, most of which are based on rRNA and rDNA analyses, have uncovered part of the soil microbial diversity. The next step in the era of microbial ecology is to extract genomic, evolutionary and functional information from bacterial artificial chromosome libraries of the soil community genomes (the metagenome). Sophisticated analyses that apply molecular phylogenetics, DNA microarrays, functional genomics and in situ activity measurements will provide huge amounts of new data, potentially increasing our understanding of the structure and function of soil microbial ecosystems, and the interactions that occur within them. This review summarizes the recent progress in studies of soil microbial communities with focus on novel methods and approaches that provide new insight into the relationship between phylogenetic and functional diversity.     

Phylogenetic analysis of intestinal microflora indicates a novel Mycoplasma phylotype in farmed and wild salmon

All studies of the microbial community of the gastrointestinal tract of salmon to date have employed culture-based approaches, typically on pond- or tank-raised, freshwater animals. We present a phylogenetic survey of the bacterial populations present in the distal intestine of salmon from three different marine locations in Europe. This was accomplished through PCR amplification, cloning, and sequencing of partial 16S rDNA genes from microbial community DNA isolated from the contents of the GI tract distal to the pyloric ceca. Using this approach, the intestinal microbial communities of wild salmon from Scotland and pen-raised salmon from Scotland and Norway were compared. The predominating bacterial populations detected were Acinetobacter junii and a novel Mycoplasma phylotype. This Mycoplasma phylotype apparently comprised ~96% of the total microbes in the distal intestine of wild salmon. Substantial differences in intestinal microbial community composition and diversity were observed between the two groups of pen-raised salmon, which, in addition to geographical separation, were raised on different feeds. The microbial profiles found in this study were substantially different from those indicated in earlier culture-based studies for several species of fish, presumably because of the culture-independent techniques employed. Further, analysis of short-chain fatty acids in the digestive tract indicated that the decreasing redox gradient from proximal to distal reaches common to homeothermic animals was absent in salmon, and that the bacterial fermentation levels were much lower than are reported in homeothermic animals.     

Marine microbial community structure assessed from combined metagenomic analysis and ribosomal amplicon deep-sequencing

The microbial taxonomic composition of the three domains of life in two coastal plankton samples was assessed by random total community metagenomic sequencing and PCR-based rDNA amplicon deep-sequencing in order to compare the resulting diversity and investigate possible limitations and complementarities of each method. The various universal primer sets, used to amplify different hypervariable rDNA regions, revealed the same major high-level taxonomic groups in Bacteria and unicellular Eukaryota, and showed a scarce Archaea apparent richness. However, significant differences were found between the different primer sets (p-value < 0.05, with the Kolmogorov–Smirnov test), regarding both operational taxonomic unit (OTU) richness and relative abundance of the major high-level taxonomic groups detected. Based on the metagenomic approach, the phylum Bacteroidetes dominated the prokaryotic community, followed by Proteobacteria, while the detected eukaryotic unicellular taxa belonged to the groups of Alveolata, Fungi, Chlorophyta, Stramenopiles and Phaeophyceae. These groups were found to carry genes typically found in microbial communities, which are linked to DNA, RNA and protein metabolism and the synthesis of nucleotides, amino acids, carbohydrates and vitamins. Although our findings suggest that the total community metagenomic approach can provide a more comprehensive picture of the planktonic microbial community structure, a number of issues associated with this approach emerged. These issues include the still relatively high cost compared to amplicon sequencing, the possible low coverage of the full marine diversity, the insufficiency of databases for other gene markers than the small subunit gene, and the bias towards bacterial sequences because of their higher abundance relative to eukaryotes in marine environments.     

Next-generation sequencing reveals significant bacterial diversity of botrytized wine

While wine fermentation has long been known to involve complex microbial communities, the composition and role of bacteria other than a select set of lactic acid bacteria (LAB) has often been assumed either negligible or detrimental. This study served as a pilot study for using barcoded amplicon next-generation sequencing to profile bacterial community structure in wines and grape musts, comparing the taxonomic depth achieved by sequencing two different domains of prokaryotic 16S rDNA (V4 and V5). This study was designed to serve two goals: 1) to empirically determine the most taxonomically informative 16S rDNA target region for barcoded amplicon sequencing of wine, comparing V4 and V5 domains of bacterial 16S rDNA to terminal restriction fragment length polymorphism (TRFLP) of LAB communities; and 2) to explore the bacterial communities of wine fermentation to better understand the biodiversity of wine at a depth previously unattainable using other techniques. Analysis of amplicons from the V4 and V5 provided similar views of the bacterial communities of botrytized wine fermentations, revealing a broad diversity of low-abundance taxa not traditionally associated with wine, as well as atypical LAB communities initially detected by TRFLP. The V4 domain was determined as the more suitable read for wine ecology studies, as it provided greater taxonomic depth for profiling LAB communities. In addition, targeted enrichment was used to isolate two species of Alphaproteobacteria from a finished fermentation. Significant differences in diversity between inoculated and uninoculated samples suggest that Saccharomyces inoculation exerts selective pressure on bacterial diversity in these fermentations, most notably suppressing abundance of acetic acid bacteria. These results determine the bacterial diversity of botrytized wines to be far higher than previously realized, providing further insight into the fermentation dynamics of these wines, and demonstrate the utility of next-generation sequencing for wine ecology studies.     

Long insert whole genome sequencing for copy number variant and translocation detection

As next-generation sequencing continues to have an expanding presence in the clinic, the identification of the most cost-effective and robust strategy for identifying copy number changes and translocations in tumor genomes is needed. We hypothesized that performing shallow whole genome sequencing (WGS) of 900–1000-bp inserts (long insert WGS, LI-WGS) improves our ability to detect these events, compared with shallow WGS of 300–400-bp inserts. A priori analyses show that LI-WGS requires less sequencing compared with short insert WGS to achieve a target physical coverage, and that LI-WGS requires less sequence coverage to detect a heterozygous event with a power of 0.99. We thus developed an LI-WGS library preparation protocol based off of Illumina’s WGS library preparation protocol and illustrate the feasibility of performing LI-WGS. We additionally applied LI-WGS to three separate tumor/normal DNA pairs collected from patients diagnosed with different cancers to demonstrate our application of LI-WGS on actual patient samples for identification of somatic copy number alterations and translocations. With the evolution of sequencing technologies and bioinformatics analyses, we show that modifications to current approaches may improve our ability to interrogate cancer genomes.     

海洋微生物宏基因组工程进展与展望

据初步统计,生活于海洋环境包括大洋深处的微生物有100万种以上,构成了一个动态的遗传基因库,其中绝大多数微生物或者从来没有经过实验室培养,或者至今无法培养,因而其分类地位及其生态学功能尚未为人类所认识。随着16S rRNA序列分析与系统分类学的广泛应用,海洋微生物多样性研究领域已经发生了很可观的改变,这些变化极大的丰富了人们对的微生物多样性及其生态功能的认识和理解。这里结合笔者近十年来的工作实践,讨论近年来在海洋微生物资源开发利用方面的研究进展,提出一个带有自动化特征的宏基因组功能表达平台,探讨海洋微生物资源利用的新途径。可以预见在不久的将来,海洋环境宏基因组工程研究将在一定程度上使得传统未培养海洋微生物基因资源及其功能产物能够为人类所开发和利用。     

Rapid quantitative profiling of complex microbial populations

Diverse and complex microbial ecosystems are found in virtually every environment on earth, yet we know very little about their composition and ecology. Comprehensive identification and quantification of the constituents of these microbial communities—a ‘census’—is an essential foundation for understanding their biology. To address this problem, we developed, tested and optimized a DNA oligonucleotide microarray composed of 10 462 small subunit (SSU) ribosomal DNA (rDNA) probes (7167 unique sequences) selected to provide quantitative information on the taxonomic composition of diverse microbial populations. Using our optimized experimental approach, this microarray enabled detection and quantification of individual bacterial species present at fractional abundances of <0.1% in complex synthetic mixtures. The estimates of bacterial species abundance obtained using this microarray are similar to those obtained by phylogenetic analysis of SSU rDNA sequences from the same samples—the current ‘gold standard’ method for profiling microbial communities. Furthermore, probes designed to represent higher order taxonomic groups of bacterial species reliably detected microbes for which there were no species-specific probes. This simple, rapid microarray procedure can be used to explore and systematically characterize complex microbial communities, such as those found within the human body.     

Study of bacterial communities in Antarctic coastal waters by a combination of 16S rRNA and 16S rDNA sequencing

An ecological study on distribution of Antarctic bacterial communities was determined by 16S-based phylogenetic analyses of clone libraries derived from RNA and DNA extracted from two different marine areas and compared between each other. Superficial seawater samples were collected from four stations in Ross Sea, three of them located in Rod Bay and one in Evans Cove; for each station two clone libraries (16S rDNA and 16S rRNA) were prepared and evident divergences between DNA and RNA libraries of each site were obtained. Of all phylotypes 93.6% were found in RNA libraries; in contrast, only 31 phylotypes (70.5%) were retrieved from total microbial community (DNA libraries). DNA and RNA sequences related to gamma-Proteobacteria and Bacteroidetes groups, typical for Antarctic sea-ice bacterial communities, were detected in analysed sites. 16S rDNA and rRNA libraries derived from the two different areas were enriched by picophytoplanktonic 16S sequences of plastid and mitochondrion origins, reflecting that the algal blooms occurred during sampling (Antarctic summer 2003). The finding in Rod Bay libraries of high percentage of DNA clones apparently affiliated with beta-Proteobacteria typical for activated sludges and well water could be explained by the presence of a sewage depuration system at this site. Obtained results clearly demonstrate that combination of 16S rDNA and 16S rRNA gene sequencing is preferred approach to have a more reliable vision on the composition of microbial communities.     

Whole genome sequencing and future breeding of rice

Contemporary crop improvement relies on the genetic analysis of progeny derived from a cross between different lines with contrasting phenotypes. Such analysis allowed positioning of genes for agronomically important traits, enabling development of DNA makers for marker-assisted selection (MAS). So far the identification of loci for desirable traits have been carried out by linkage analysis using DNA markers. This process required the development of DNA markers that are distributed over the genome as well as the genotyping of each progeny. Due to recent development in next generation sequencing (NGS) technology, whole genome sequencing (WGS) is becoming easier and cheaper. Using NGS, we developed a new method called MutMap that allows rapid isolation of useful alleles from rice mutant lines. An important feature of MutMap is that it does not require marker development. We foresee that the era of genetic markers will be eventually eclipsed by that of WGS applied to all the individuals in the breeding processes.     

Pyrosequencing for Microbial Identification and Characterization

Pyrosequencing is a versatile technique that facilitates microbial genome sequencing that can be used to identify bacterial species, discriminate bacterial strains and detect genetic mutations that confer resistance to anti-microbial agents. The advantages of pyrosequencing for microbiology applications include rapid and reliable high-throughput screening and accurate identification of microbes and microbial genome mutations. Pyrosequencing involves sequencing of DNA by synthesizing the complementary strand a single base at a time, while determining the specific nucleotide being incorporated during the synthesis reaction. The reaction occurs on immobilized single stranded template DNA where the four deoxyribonucleotides (dNTP) are added sequentially and the unincorporated dNTPs are enzymatically degraded before addition of the next dNTP to the synthesis reaction. Detection of the specific base incorporated into the template is monitored by generation of chemiluminescent signals. The order of dNTPs that produce the chemiluminescent signals determines the DNA sequence of the template. The real-time sequencing capability of pyrosequencing technology enables rapid microbial identification in a single assay. In addition, the pyrosequencing instrument, can analyze the full genetic diversity of anti-microbial drug resistance, including typing of SNPs, point mutations, insertions, and deletions, as well as quantification of multiple gene copies that may occur in some anti-microbial resistance patterns.     

Analysis of the Cystic Fibrosis Lung Microbiota via Serial Illumina Sequencing of Bacterial 16S rRNA Hypervariable Regions

The characterization of bacterial communities using DNA sequencing has revolutionized our ability to study microbes in nature and discover the ways in which microbial communities affect ecosystem functioning and human health. Here we describe Serial Illumina Sequencing (SI-Seq): a method for deep sequencing of the bacterial 16S rRNA gene using next-generation sequencing technology. SI-Seq serially sequences portions of the V5, V6 and V7 hypervariable regions from barcoded 16S rRNA amplicons using an Illumina short-read genome analyzer. SI-Seq obtains taxonomic resolution similar to 454 pyrosequencing for a fraction of the cost, and can produce hundreds of thousands of reads per sample even with very high multiplexing. We validated SI-Seq using single species and mock community controls, and via a comparison to cystic fibrosis lung microbiota sequenced using 454 FLX Titanium. Our control runs show that SI-Seq has a dynamic range of at least five orders of magnitude, can classify >96% of sequences to the genus level, and performs just as well as 454 and paired-end Illumina methods in estimation of standard microbial ecology diversity measurements. We illustrate the utility of SI-Seq in a pilot sample of central airway secretion samples from cystic fibrosis patients.     

Microbial diversity of topographical gradient profiles in Fushan forest soils of Taiwan

To evaluate the microbial diversity of Fushan forest soils, the variation of soil properties, microbial populations, and soil DNA with soil depth in three sites of different altitude were analyzed. Microbial population, moisture content, total organic carbon (C_org), and total nitrogen (N_tot) decreased with increasing soil depth. The valley site had the lowest microbial populations among the three tested sites due to the low organic matter content. Bacterial population was the highest among the microbial populations. The ratios of cellulolytic microbes to the total bacteria in organic layers were high, implying their roles in the carbon cycle. The microbial biomass carbon (C_mic) and nitrogen (N_mic) contents ranged from 130.5 to 564.1 μg g⁻¹ and from 16.7 to 95.4 μg g⁻¹, respectively. The valley had the lowest C_mic and N_mic. The organic layer had the highest C_mic and N_mic and decreased with soil depth. Analysis using denaturing gradient gel electrophoresis (DGGE) of polymerase chain reaction (PCR) amplicons of 16S rDNA showed that the bacterial diversity of the three sites were very similar to each other in the major bands, and the variation was in the minor bands. However, the patterns in PCR-DGGE profile through gradient horizons were different, indicating the prevalence of specific microbes at different horizons. These results suggest that the microbial diversity in the deeper horizons is not simply the diluted analogs of the surface soils and that some microbes dominate only in the deeper horizons. Topography influenced the quantity and diversity of microbial populations.