Multi-amplicon microbiome data analysis pipelines for mixed orientation sequences using QIIME2: Assessing reference database, variable region and pre-processing bias in classification of mock bacterial community samples

Using short-read 16S rRNA sequencing of multiple variable regions to generate high-quality results to a species level

Introduction

Short-read amplicon sequencing studies have typically focused on 1-2 variable regions of the 16S rRNA gene. Species-level resolution is limited in these studies, as each variable region enables the characterisation of a different subsection of the microbiome. Although long-read sequencing techniques take advantage of all 9 variable regions by sequencing the entire 16S rRNA gene, they are substantially more expensive. This work assessed the feasibility of accurate species-level resolution and reproducibility using a relatively new sequencing kit and bioinformatics pipeline developed for short-read sequencing of multiple variable regions of the 16S rRNA gene. In addition, we evaluated the potential impact of different sample collection methods on our outcomes.

Methods

Using xGen™ 16S Amplicon Panel v2 kits, sequencing of all 9 variable regions of the 16S rRNA gene was carried out on an Illumina MiSeq platform. Mock cells and mock DNA for 8 bacterial species were included as extraction and sequencing controls respectively. Within-run and between-run replicate samples, and pairs of stool and rectal swabs collected at 0-5 weeks from the same participants, were incorporated. Observed relative abundances of each species were compared to theoretical abundances provided by ZymoBIOMICS. Paired Wilcoxon rank sum tests and distance-based intraclass correlation coefficients were used to statistically compare alpha and beta diversity measures, respectively, for pairs of replicates and stool/rectal swab sample pairs.

Results

Using multiple variable regions of the 16S ribosomal Ribonucleic Acid (rRNA) gene, we found that we could accurately identify taxa to a species level and obtain highly reproducible results at a species level. Yet, the microbial profiles of stool and rectal swab sample pairs differed substantially despite being collected concurrently from the same infants.

Conclusion

This protocol provides an effective means for studying infant gut microbial samples at a species level. However, sample collection approaches need to be accounted for in any downstream analysis.

A lipidome landscape of aging in mice

Understanding the molecular mechanisms of aging is crucial for enhancing healthy longevity. We conducted untargeted lipidomics across 13 biological samples from mice at various life stages (2, 12, 19 and 24 months) to explore the potential link between aging and lipid metabolism, considering sex (male or female) and microbiome (specific pathogen-free or germ-free) dependencies. By analyzing 2,704 molecules from 109 lipid subclasses, we characterized common and tissue-specific lipidome alterations associated with aging. For example, the levels of bis(monoacylglycero)phosphate containing polyunsaturated fatty acids increased in various organs during aging, whereas the levels of other phospholipids containing saturated and monounsaturated fatty acids decreased. In addition, we discovered age-dependent sulfonolipid accumulation, absent in germ-free mice, correlating with Alistipes abundance determined by 16S ribosomal RNA gene amplicon sequencing. In the male kidney, glycolipids such as galactosylceramides, galabiosylceramides (Gal2Cer), trihexosylceramides (Hex3Cer), and mono- and digalactosyldiacylglycerols were detected, with two lipid classes-Gal2Cer and Hex3Cer-being significantly enriched in aged mice. Integrated analysis of the kidney transcriptome revealed uridine diphosphate galactosyltransferase 8A (UGT8a), alkylglycerone phosphate synthase and fatty acyl-coenzyme A reductase 1 as potential enzymes responsible for the male-specific glycolipid biosynthesis in vivo, which would be relevant to sex dependency in kidney diseases. Inhibiting UGT8 reduced the levels of these glycolipids and the expression of inflammatory cytokines in the kidney. Our study provides a valuable resource for clarifying potential links between lipid metabolism and aging.

Mock community taxonomic classification performance of publicly available shotgun metagenomics pipelines

Shotgun metagenomic sequencing comprehensively samples the DNA of a microbial sample. Choosing the best bioinformatics processing package can be daunting due to the wide variety of tools available. Here, we assessed publicly available shotgun metagenomics processing packages/pipelines including bioBakery, Just a Microbiology System (JAMS), Whole metaGenome Sequence Assembly V2 (WGSA2), and Woltka using 19 publicly available mock community samples and a set of five constructed pathogenic gut microbiome samples. Also included is a workflow for labelling bacterial scientific names with NCBI taxonomy identifiers for better resolution in assessing results. The Aitchison distance, a sensitivity metric, and total False Positive Relative Abundance were used for accuracy assessments for all pipelines and mock samples. Overall, bioBakery4 performed the best with most of the accuracy metrics, while JAMS and WGSA2, had the highest sensitivities. Furthermore, bioBakery is commonly used and only requires a basic knowledge of command line usage. This work provides an unbiased assessment of shotgun metagenomics packages and presents results assessing the performance of the packages using mock community sequence data.

Mock microbial community meta-analysis using different trimming of amplicon read lengths

Trimming of sequencing reads is a pre-processing step that aims to discard sequence segments such as primers, adapters and low quality nucleotides that will interfere with clustering and classification steps. We evaluated the impact of trimming length of paired-end 16S and 18S rRNA amplicon reads on the ability to reconstruct the taxonomic composition and relative abundances of communities with a known composition in both even and uneven proportions. We found that maximizing read retention maximizes recall but reduces precision by increasing false positives. The presence of expected taxa was accurately predicted across broad trim length ranges but recovering original relative proportions remains a difficult challenge. We show that parameters that maximize taxonomic recovery do not simultaneously maximize relative abundance accuracy. Trim length represents one of several experimental parameters that have non-uniform impact across microbial clades, making it a difficult parameter to optimize. This study offers insights, guidelines, and helps researchers assess the significance of their decisions when trimming raw reads in a microbiome analysis based on overlapping or non-overlapping paired-end amplicons.

Effects of blister blight disease on endophytic microbial diversity and community structure in tea (Camellia sinensis) leaves

In this study, metagenomic sequencing technology was employed to analyze the ITS1 region sequence of the ITS rDNA gene of endophytic fungi and 16S sequence of endophytic bacteria in tea leaves with varying degrees of infection by tea blister blight disease as well as healthy tea leaves. Subsequently, a comparative analysis was conducted on the endophytic microbial diversity and the community structure in tea leaves. The findings of this investigation reveal a shift in the dominant endophytic fungal genera from Ascomycota to Basidiomycota as the disease progressed. Furthermore, a negative correlation was observed between Exobasidium and Talaromyce, with Talaromyce exhibiting potential as an antagonist against the disease. Meanwhile, our findings reveal that Proteobacteria, Firmicutes, and Actinobacteria were the three most abundant bacteria phyla in tea leaves. As the disease progressed, there was an increase in the relative abundance of Actinobacteria, while Variovorax, Sphingomonas, and Pseudomonas were found to have higher abundance in later stages. The diversity analysis results indicated that the endophytic microbial diversity and the community structure in tea leaves in the diseased group were lower than those in the healthy control group. In general, blister blight disease altered the community structure of endophytic microorganisms in tea leaves, resulting in a few species with high abundance. The study lays a foundation for investigating the pathogenic mechanism of tea blister disease and establishing a theoretical basis for controlling diseases in tea trees.