Considering external information to improve the phylogenetic comparison of microbial communities: a new approach based on constrained Double Principal Coordinates Analysis (cDPCoA)

The use of next‐generation sequencing technologies is revolutionizing microbial ecology by allowing a deep phylogenetic coverage of tens to thousands of samples simultaneously. Double Principal Coordinates Analysis (DPCoA) is a multivariate method, developed in community ecology, able to integrate a distance matrix describing differences among species (e.g. phylogenetic distances) in the analysis of a species abundance matrix. This ordination technique has been used recently to describe microbial communities taking into account phylogenetic relatedness. In this work, we extend DPCoA to integrate the information of external variables measured on communities. The constrained Double Principal Coordinates Analysis (cDPCoA) is able to enforce a priori classifications to retrieve subtle differences and (or) remove the effect of confounding factors. We describe the main principles of this new approach and demonstrate its usefulness by providing application examples based on published 16S rRNA gene data sets.     

Benchmarking of 16S rRNA gene databases using known strain sequences

16S rRNA gene analysis is the most convenient and robust method for microbiome studies. Inaccurate taxonomic assignment of bacterial strains could have deleterious effects as all downstream analyses rely heavily on the accurate assessment of microbial taxonomy. The use of mock communities to check the reliability of the results has been suggested. However, often the mock communities used in most of the studies represent only a small fraction of taxa and are used mostly as validation of sequencing run to estimate sequencing artifacts. Moreover, a large number of databases and tools available for classification and taxonomic assignment of the 16S rRNA gene make it challenging to select the best-suited method for a particular dataset. In the present study, we used authentic and validly published 16S rRNA gene type strain sequences (full length, V3-V4 region) and analyzed them using a widely used QIIME pipeline along with different parameters of OTU clustering and QIIME compatible databases. Data Analysis Measures (DAM) revealed a high discrepancy in ratifying the taxonomy at different taxonomic hierarchies. Beta diversity analysis showed clear segregation of different DAMs. Limited differences were observed in reference data set analysis using partial (V3-V4) and full-length 16S rRNA gene sequences, which signify the reliability of partial 16S rRNA gene sequences in microbiome studies. Our analysis also highlights common discrepancies observed at various taxonomic levels using various methods and databases.     

Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms

DNA sequencing continues to decrease in cost with the Illumina HiSeq2000 generating up to 600 Gb of paired-end 100 base reads in a ten-day run. Here we present a protocol for community amplicon sequencing on the HiSeq2000 and MiSeq Illumina platforms, and apply that protocol to sequence 24 microbial communities from host-associated and free-living environments. A critical question as more sequencing platforms become available is whether biological conclusions derived on one platform are consistent with what would be derived on a different platform. We show that the protocol developed for these instruments successfully recaptures known biological results, and additionally that biological conclusions are consistent across sequencing platforms (the HiSeq2000 versus the MiSeq) and across the sequenced regions of amplicons.     

Phylogenetic stratigraphy in the Guerrero Negro hypersaline microbial mat

The microbial mats of Guerrero Negro (GN), Baja California Sur, Mexico historically were considered a simple environment, dominated by cyanobacteria and sulfate-reducing bacteria. Culture-independent rRNA community profiling instead revealed these microbial mats as among the most phylogenetically diverse environments known. A preliminary molecular survey of the GN mat based on only ∼1500 small subunit rRNA gene sequences discovered several new phylum-level groups in the bacterial phylogenetic domain and many previously undetected lower-level taxa. We determined an additional ∼119 000 nearly full-length sequences and 28 000 >200 nucleotide 454 reads from a 10-layer depth profile of the GN mat. With this unprecedented coverage of long sequences from one environment, we confirm the mat is phylogenetically stratified, presumably corresponding to light and geochemical gradients throughout the depth of the mat. Previous shotgun metagenomic data from the same depth profile show the same stratified pattern and suggest that metagenome properties may be predictable from rRNA gene sequences. We verify previously identified novel lineages and identify new phylogenetic diversity at lower taxonomic levels, for example, thousands of operational taxonomic units at the family-genus levels differ considerably from known sequences. The new sequences populate parts of the bacterial phylogenetic tree that previously were poorly described, but indicate that any comprehensive survey of GN diversity has only begun. Finally, we show that taxonomic conclusions are generally congruent between Sanger and 454 sequencing technologies, with the taxonomic resolution achieved dependent on the abundance of reference sequences in the relevant region of the rRNA tree of life.     

Windows下16S rRNA基因扩增子测序数据分析的简易流程

微生物组数据分析需要掌握Linux系统操作,这对缺乏计算机知识的生物研究人员是一个很大的障碍。为此我们设计了一套在Windows的Linux子系统(WSL)下分析16S rRNA基因扩增子高通量测序数据的简易流程。本流程整合常用的开源软件VSEARCH与QIIME等,能对16S rRNA测序数据进行质量控制、OTU聚类、多样性分析及结果可视化呈现。以唾液微生物组分析为例,详细介绍从原始数据到多样性统计分析过程的参数和命令,及结果解读。教学实践证明,此流程易于学习,并有助于掌握微生物组的基本概念与方法。利用Windows系统最新的WSL功能,本流程方便Windows用户使用大量在Linux上运行的生物信息工具,有助于促进微生物组研究的发展。流程的安装程序与测序数据可从网址(http://www. ligene. cn/win16s/)免费下载使用。     

The influence of genetics,defensive chemistry and the fungal microbiome on disease outcome in whitebark pine trees

The invasive fungal pathogen Cronartium ribicola infects and kills whitebark pine (Pinus albicaulis) throughout western North America. Whitebark pine has been proposed for listing under the Endangered Species Act in the USA, and the loss of this species is predicted to have severe impacts on ecosystem composition and function in high‐elevation forests. Numerous fungal endophytes live inside whitebark pine tissues and may influence the severity of C. ribicola infection, either directly by inhibition of pathogen growth or indirectly by the induction of chemical defensive pathways in the tree. Terpenes, a form of chemical defence in pine trees, can also influence disease. In this study, we characterized fungal endophyte communities in whitebark pine seedlings before and after experimental inoculation with C. ribicola, monitored disease progression and compared fungal community composition in susceptible vs. resistant seedlings in a common garden. We analysed the terpene composition of these same seedlings. Seed family identity or maternal genetics influenced both terpenes and endophyte communities. Terpene and endophyte composition correlated with disease severity, and terpene concentrations differed in resistant vs. susceptible seedlings. These results suggest that the resistance to C. ribicola observed in natural whitebark pine populations is caused by the combined effects of genetics, endophytes and terpenes within needle tissue, in which initial interactions between microbes and hosts take place. Tree genotype, terpene and microbiome combinations associated with healthy trees could help to predict or reduce disease severity and improve outcomes of future tree breeding programmes.     

Using direct amplification and next-generation sequencing technology to explore foliar endophyte communities in experimentally inoculated western white pines

Fungal endophytes can influence survivability and disease severity of trees. Here we characterized the endophyte community in Pinus monticola (western white pine), an important species in the northwest USA, largely decimated by pathogenic fungi. We also assessed the ability to successfully inoculate seedlings with desirable endophytes, with the long-term goal of providing a protective microbiome and added defense from pathogens. P. monticola seedlings were inoculated in the field with potential pathogen antagonists and fungi isolated from healthy mature trees. Following inoculations direct amplification and next generation sequencing were used to characterize fungal endophyte communities, and explore interspecific competition, diversity, and co-occurrence patterns in needle tissues. Negative co-occurrence patterns between inoculated fungi and potential pathogens, as well as many other species, suggest early competitive interactions. Our study explores early endophyte community assemblage and shows that fungal inoculations may influence tree growth.