Shotgun Sequencing Approaches Research Articles

Benchmarking a targeted 16S ribosomal RNA gene enrichment approach to reconstruct ancient microbial communities.

The taxonomic characterization of ancient microbiomes is a key step in the rapidly growing field of paleomicrobiology. While PCR amplification of the 16S ribosomal RNA (rRNA) gene is a widely used technique in modern microbiota studies, this method has systematic biases when applied to ancient microbial DNA. Shotgun metagenomic sequencing has proven to be the most effective method in reconstructing taxonomic profiles of ancient dental calculus samples. Nevertheless, shotgun sequencing approaches come with inherent limitations that could be addressed through hybridization enrichment capture. When employed together, shotgun sequencing and hybridization capture have the potential to enhance the characterization of ancient microbial communities. Here, we develop, test, and apply a hybridization enrichment capture technique to selectively target 16S rRNA gene fragments from the libraries of ancient dental calculus samples generated with shotgun techniques. We simulated data sets generated from hybridization enrichment capture, indicating that taxonomic identification of fragmented and damaged 16S rRNA gene sequences was feasible. Applying this enrichment approach to 15 previously published ancient calculus samples, we observed a 334-fold increase of ancient 16S rRNA gene fragments in the enriched samples when compared to unenriched libraries. Our results suggest that 16S hybridization capture is less prone to the effects of background contamination than 16S rRNA amplification, yielding a higher percentage of on-target recovery. While our enrichment technique detected low abundant and rare taxa within a given sample, these assignments may not achieve the same level of specificity as those achieved by unenriched methods.

Metagenomic exploration of Andaman region of the Indian Ocean

Ocean microbiome is crucial for global biogeochemical cycles and primary productivity. Despite numerous studies investigating the global ocean microbiomes, the microbiome composition of the Andaman region of the Indian Ocean remains largely unexplored. While this region harbors pristine biological diversity, the escalating anthropogenic activities along coastal habitats exert an influence on the microbial ecology and impact the aquatic ecosystems. We investigated the microbiome composition in the coastal waters of the Andaman Islands by 16S rRNA gene amplicon and metagenomic shotgun sequencing approaches and compared it with the Tara Oceans Consortium. In the coastal waters of the Andaman Islands, a significantly higher abundance and diversity of Synechococcus species was observed with a higher abundance of photosynthesis pigment-related genes to adapt to variable light conditions and nutrition. In contrast, Prochlorococcus species showed higher abundance in open ocean water samples of the Indian Ocean region, with a relatively limited functional diversity. A higher abundance of antibiotic-resistance genes was also noted in the coastal waters region. We also updated the ocean microbiome gene catalog with 93,172 unique genes from the Andaman coastal water microbiome. This study provides valuable insights into the Indian Ocean microbiome and supplements the global marine microbial ecosystem studies.

Microbiome of Ceca from Broiler Chicken Vaccinated or Not against Coccidiosis and Fed Berry Pomaces.

American cranberry (Vaccinium macrocarpon) and lowbush/wild blueberry (V. angustifolium) pomace are polyphenol-rich products having potentially beneficial effects in broiler chickens. This study investigated the cecal microbiome of broiler-vaccinated or non-vaccinated birds against coccidiosis. Birds in each of the two groups (vaccinated or non-vaccinated) were fed a basal non-supplemented diet (NC), a basal diet supplemented with bacitracin (BAC), American cranberry (CP), and lowbush blueberry (BP) pomace alone or in combination (CP + BP). At 21 days of age, cecal DNA samples were extracted and analyzed using both whole-metagenome shotgun sequencing and targeted-resistome sequencing approaches. Ceca from vaccinated birds showed a lower abundance of Lactobacillus and a higher abundance of Escherichia coli than non-vaccinated birds (p < 0.05). The highest and lowest abundance of L. crispatus and E. coli, respectively, were observed in birds fed CP, BP, and CP + BP compared to those from NC or BAC treatments (p < 0.05). Coccidiosis vaccination affected the abundance of virulence genes (VGs) related to adherence, flagella, iron utilization, and secretion system. Toxin-related genes were observed in vaccinated birds (p < 0.05) in general, with less prevalence in birds fed CP, BP, and CP + BP than NC and BAC (p < 0.05). More than 75 antimicrobial resistance genes (ARGs) detected by the shotgun metagenomics sequencing were impacted by vaccination. Ceca from birds fed CP, BP, and CP + BP showed the lowest (p < 0.05) abundances of ARGs related to multi-drug efflux pumps, modifying/hydrolyzing enzyme and target-mediated mutation, when compared to ceca from birds fed BAC. Targeted metagenomics showed that resistome from BP treatment was distant to other groups for antimicrobials, such as aminoglycosides (p < 0.05). Significant differences in the richness were observed between the vaccinated and non-vaccinated groups for aminoglycosides, β-lactams, lincosamides, and trimethoprim resistance genes (p < 0.05). Overall, this study demonstrated that dietary berry pomaces and coccidiosis vaccination significantly impacted cecal microbiota, virulome, resistome, and metabolic pathways in broiler chickens.

Three assays for in-solution enrichment of ancient human DNA at more than a million SNPs.

The strategy of in-solution enrichment for hundreds of thousands of single-nucleotide polymorphisms (SNPs) has been used to analyze >70% of individuals with genome-scale ancient DNA published to date. This approach makes it economical to study ancient samples with low proportions of human DNA and increases the rate of conversion of sampled remains into interpretable data. So far, nearly all such data have been generated using a set of bait sequences targeting about 1.24 million SNPs (the "1240k reagent"), but synthesis of the reagent has been cost-effective for only a few laboratories. In 2021, two companies, Daicel Arbor Biosciences and Twist Bioscience, made available assays that target the same core set of SNPs along with supplementary content. We test all three assays on a common set of 27 ancient DNA libraries and show that all three are effective at enriching many hundreds of thousands of SNPs. For all assays, one round of enrichment produces data that are as useful as two. In our testing, the "Twist Ancient DNA" assay produces the highest coverages, greatest uniformity on targeted positions, and almost no bias toward enriching one allele more than another relative to shotgun sequencing. We also identify hundreds of thousands of targeted SNPs for which there is minimal allelic bias when comparing 1240k data to either shotgun or Twist data. This facilitates coanalysis of the large data sets that have been generated using 1240k and Twist capture, as well as shotgun sequencing approaches.

The Application of Metagenomics to Study Microbial Communities and Develop Desirable Traits in Fermented Foods

The microbial communities present within fermented foods are diverse and dynamic, producing a variety of metabolites responsible for the fermentation processes, imparting characteristic organoleptic qualities and health-promoting traits, and maintaining microbiological safety of fermented foods. In this context, it is crucial to study these microbial communities to characterise fermented foods and the production processes involved. High Throughput Sequencing (HTS)-based methods such as metagenomics enable microbial community studies through amplicon and shotgun sequencing approaches. As the field constantly develops, sequencing technologies are becoming more accessible, affordable and accurate with a further shift from short read to long read sequencing being observed. Metagenomics is enjoying wide-spread application in fermented food studies and in recent years is also being employed in concert with synthetic biology techniques to help tackle problems with the large amounts of waste generated in the food sector. This review presents an introduction to current sequencing technologies and the benefits of their application in fermented foods.

Evaluation of Bacteria and Fungi DNA Abundance in Human Tissues.

Whereas targeted and shotgun sequencing approaches are both powerful in allowing the study of tissue-associated microbiota, the human: microorganism abundance ratios in tissues of interest will ultimately determine the most suitable sequencing approach. In addition, it is possible that the knowledge of the relative abundance of bacteria and fungi during a treatment course or in pathological conditions can be relevant in many medical conditions. Here, we present a qPCR-targeted approach to determine the absolute and relative amounts of bacteria and fungi and demonstrate their relative DNA abundance in nine different human tissue types for a total of 87 samples. In these tissues, fungi genomes are more abundant in stool and skin samples but have much lower levels in other tissues. Bacteria genomes prevail in stool, skin, oral swabs, saliva, and gastric fluids. These findings were confirmed by shotgun sequencing for stool and gastric fluids. This approach may contribute to a more comprehensive view of the human microbiota in targeted studies for assessing the abundance levels of microorganisms during disease treatment/progression and to indicate the most informative methods for studying microbial composition (shotgun versus targeted sequencing) for various samples types.

Comparing whole-genome shotgun sequencing and DNA metabarcoding approaches for species identification and quantification of pollen species mixtures.

Molecular identification of mixed‐species pollen samples has a range of applications in various fields of research. To date, such molecular identification has primarily been carried out via amplicon sequencing, but whole‐genome shotgun (WGS) sequencing of pollen DNA has potential advantages, including (1) more genetic information per sample and (2) the potential for better quantitative matching. In this study, we tested the performance of WGS sequencing methodology and publicly available reference sequences in identifying species and quantifying their relative abundance in pollen mock communities. Using mock communities previously analyzed with DNA metabarcoding, we sequenced approximately 200Mbp for each sample using Illumina HiSeq and MiSeq. Taxonomic identifications were based on the Kraken k‐mer identification method with reference libraries constructed from full‐genome and short read archive data from the NCBI database. We found WGS to be a reliable method for taxonomic identification of pollen with near 100% identification of species in mixtures but generating higher rates of false positives (reads not identified to the correct taxon at the required taxonomic level) relative to rbcL and ITS2 amplicon sequencing. For quantification of relative species abundance, WGS data provided a stronger correlation between pollen grain proportion and sequence read proportion, but diverged more from a 1:1 relationship, likely due to the higher rate of false positives. Currently, a limitation of WGS‐based pollen identification is the lack of representation of plant diversity in publicly available genome databases. As databases improve and costs drop, we expect that eventually genomics methods will become the methods of choice for species identification and quantification of mixed‐species pollen samples.

A primer on metagenomics and next‐generation sequencing in fish gut microbiome research

The fish gut is a complex ecosystem populated by an exceedingly diverse and assorted group of microbes, influenced by various factors such as habitat salinity, temperature, trophic level, taxonomy, feeding habits as well as host selective pressures within the gut. The gut microbiota is involved in several physiological functions of the host such as digestion, metabolism, reproduction, development, and immune response and is often heralded as an ‘extra organ’. Earlier studies employed culture-dependent techniques for the study of fish gut microflora. The establishment of great plate anomaly and technological strides in culture-independent techniques, particularly in the omics technologies such as metagenomics and next-generation sequencing techniques have led to a better understanding of the confederation of microorganisms with their respective hosts. This article reviews the applications of and prospects of metagenomics and amplicon and shotgun sequencing approaches in fish gut microbiota research.

Carotenoids: a missing link between Gut Microbiota and Obesity

Carotenoids are naturally occurring red, orange, and yellow pigments found mainly in fruits and vegetables. They may function as vitamin A precursors, have important antioxidant, immunological or metabolic functions. The overall absorption efficiency of carotenoids from whole foods is considerably lower in the small intestine. When carotenoids enter the colon together with indigestible dietary compounds, they can become a potential substrate for microbial fermentation. Carotenoids can be further released, or additional conversions can happen within colon epithelial cells. This study's overall objective is to determine the impact of carotenoids on obesity-associated gut dysfunction (microbial dysbiosis, impaired intestinal barrier function, and low-grade inflammation). Our central hypothesis is that carotenoids escaping the absorption in the small intestine can be transformed by the gut microbiota in the colon. First, we examined the molecular mechanisms behind the anti-inflammatory action of one of the major dietary carotenoid, β-carotene (BC), and its role in maintaining gut barrier integrity using in vitro intestinal epithelial cell model. We found out that BC could improve colonic inflammation by blocking inflammatory transcription factors and subsequently reducing pro-inflammatory cytokines in the colon. We have also seen that BC could improve gut barrier integrity within colon epithelial cells when cells challenged with endotoxins. Next, we determined the impact of microbial metabolism and fermentation on the carotenoids present in the colon using an ex vivo fecal fermentation model. We demonstrated that BC would be metabolized by the microbiome in human stool samples. Finally, we plan to examine the effect of carotenoid-rich diet on obesity-associated gut dysfunction, metabolic endotoxemia, and gut barrier in vivo using a murine model. We will characterize carotenoid-induced changes in the gut microbiota utilizing shotgun sequencing and multi-omic approaches, including metabolomics and shotgun metagenomics in the diet-induced obesity animal model.

CaptureSeq: Hybridization-Based Enrichment of cpn60 Gene Fragments Reveals the Community Structures of Synthetic and Natural Microbial Ecosystems.

Background. The molecular profiling of complex microbial communities has become the basis for examining the relationship between the microbiome composition, structure and metabolic functions of those communities. Microbial community structure can be partially assessed with “universal” PCR targeting taxonomic or functional gene markers. Increasingly, shotgun metagenomic DNA sequencing is providing more quantitative insight into microbiomes. However, both amplicon-based and shotgun sequencing approaches have shortcomings that limit the ability to study microbiome dynamics. Methods. We present a novel, amplicon-free, hybridization-based method (CaptureSeq) for profiling complex microbial communities using probes based on the chaperonin-60 gene. Molecular profiles of a commercially available synthetic microbial community standard were compared using CaptureSeq, whole metagenome sequencing, and 16S universal target amplification. Profiles were also generated for natural ecosystems including antibiotic-amended soils, manure storage tanks, and an agricultural reservoir. Results. The CaptureSeq method generated a microbial profile that encompassed all of the bacteria and eukaryotes in the panel with greater reproducibility and more accurate representation of high G/C content microorganisms compared to 16S amplification. In the natural ecosystems, CaptureSeq provided a much greater depth of coverage and sensitivity of detection compared to shotgun sequencing without prior selection. The resulting community profiles provided quantitatively reliable information about all three domains of life (Bacteria, Archaea, and Eukarya) in the different ecosystems. The applications of CaptureSeq will facilitate accurate studies of host-microbiome interactions for environmental, crop, animal and human health. Conclusions: cpn60-based hybridization enriched for taxonomically informative DNA sequences from complex mixtures. In synthetic and natural microbial ecosystems, CaptureSeq provided sequences from prokaryotes and eukaryotes simultaneously, with quantitatively reliable read abundances. CaptureSeq provides an alternative to PCR amplification of taxonomic markers with deep community coverage while minimizing amplification biases.

Reduced metagenomic sequencing (RMS) approach to determine the gut-associated phageome in mother-child

The role of the human gut phageome (HGP) for a healthy gut microbiome is not well-established. This study aims to identify phages based on Reduced Metagenome Sequencing (RMS) fragments from an Indian mother and child cohort. For this study, fecal samples were collected from 17 mother-infant pairs at Nishanth Hospital, Tamil Nadu, India. RMS data analysis and shotgun sequencing approaches were used to assemble and identify the genome fragments. Out of the 156,926 RMS fragments, 434 were classified as bacteriophages by Kraken 2. Mapping of virus sequences in NCBI and de novo assembly with subsequent taxonomic assignment revealed 41 different phage species. The prevalence (>50%) of three bacteriophages was observed in mother and child; overall four phages were more prevalent in the mothers while one phage was more prevalent in the children. Even at the species level, mothers were found to have more diverse phage species than children. No significant association was observed for mother–child sharing of phages. This study highlights the prevalence of Caudovirales phages in healthy HGP and also the use of the RMS approach to study the phageome composition.

Systematic evaluation of supervised machine learning for sample origin prediction using metagenomic sequencing data

BackgroundThe advent of metagenomic sequencing provides microbial abundance patterns that can be leveraged for sample origin prediction. Supervised machine learning classification approaches have been reported to predict sample origin accurately when the origin has been previously sampled. Using metagenomic datasets provided by the 2019 CAMDA challenge, we evaluated the influence of variable technical, analytical and machine learning approaches for result interpretation and novel source prediction.ResultsComparison between 16S rRNA amplicon and shotgun sequencing approaches as well as metagenomic analytical tools showed differences in normalized microbial abundance, especially for organisms present at low abundance. Shotgun sequence data analyzed using Kraken2 and Bracken, for taxonomic annotation, had higher detection sensitivity. As classification models are limited to labeling pre-trained origins, we took an alternative approach using Lasso-regularized multivariate regression to predict geographic coordinates for comparison. In both models, the prediction errors were much higher in Leave-1-city-out than in 10-fold cross validation, of which the former realistically forecasted the increased difficulty in accurately predicting samples from new origins. This challenge was further confirmed when applying the model to a set of samples obtained from new origins. Overall, the prediction performance of the regression and classification models, as measured by mean squared error, were comparable on mystery samples. Due to higher prediction error rates for samples from new origins, we provided an additional strategy based on prediction ambiguity to infer whether a sample is from a new origin. Lastly, we report increased prediction error when data from different sequencing protocols were included as training data.ConclusionsHerein, we highlight the capacity of predicting sample origin accurately with pre-trained origins and the challenge of predicting new origins through both regression and classification models. Overall, this work provides a summary of the impact of sequencing technique, protocol, taxonomic analytical approaches, and machine learning approaches on the use of metagenomics for prediction of sample origin.

Genomic Adaptations and Evolutionary History of the Extinct Scimitar-Toothed Cat, Homotherium latidens

SummaryHomotherium was a genus of large-bodied scimitar-toothed cats, morphologically distinct from any extant felid species, that went extinct at the end of the Pleistocene [1, 2, 3, 4]. They possessed large, saber-form serrated canine teeth, powerful forelimbs, a sloping back, and an enlarged optic bulb, all of which were key characteristics for predation on Pleistocene megafauna [5]. Previous mitochondrial DNA phylogenies suggested that it was a highly divergent sister lineage to all extant cat species [6, 7, 8]. However, mitochondrial phylogenies can be misled by hybridization [9], incomplete lineage sorting (ILS), or sex-biased dispersal patterns [10], which might be especially relevant for Homotherium since widespread mito-nuclear discrepancies have been uncovered in modern cats [10]. To examine the evolutionary history of Homotherium, we generated a ∼7x nuclear genome and a ∼38x exome from H. latidens using shotgun and target-capture sequencing approaches. Phylogenetic analyses reveal Homotherium as highly divergent (∼22.5 Ma) from living cat species, with no detectable signs of gene flow. Comparative genomic analyses found signatures of positive selection in several genes, including those involved in vision, cognitive function, and energy consumption, putatively consistent with diurnal activity, well-developed social behavior, and cursorial hunting [5]. Finally, we uncover relatively high levels of genetic diversity, suggesting that Homotherium may have been more abundant than the limited fossil record suggests [3, 4, 11, 12, 13, 14]. Our findings complement and extend previous inferences from both the fossil record and initial molecular studies, enhancing our understanding of the evolution and ecology of this remarkable lineage.

Depicting Temporal, Functional, and Phylogenetic Patterns in Estuarine Diazotrophic Communities from Environmental DNA and RNA.

Seasonally nitrogen-limited and phosphorus-replete temperate coastal waters generally host dense and diverse diazotrophic communities. Despite numerous studies in marine systems, little is known about diazotrophs and their functioning in oligohaline estuarine environments. Here we applied a combination of nifH transcript and metagenomic shotgun sequencing approaches to investigate temporal shifts in taxonomic composition and nifH activity of size-fractionated diazotrophic communities in a shallow and mostly freshwater coastal lagoon. Patterns in active nifH phylotypes exhibited a clear seasonal succession, which reflected their different tolerances to temperature change and nitrogen (N) availability. Thus, in spring, heterotrophic diazotrophs (Proteobacteria) dominated the nifH phylotypes, while increasing water temperature and depletion of inorganic N fostered heterocystous Cyanobacteria in summer. Metagenomic data demonstrated four main N-cycling pathways and three of them with a clear seasonal pattern: denitrification (spring) → N2 fixation (summer) → assimilative NO3- reduction (fall), with NH4+ uptake into cells occurring across all seasons. Although a substantial denitrification signal was observed in spring, it could have originated from the re-suspended benthic rather than planktonic community. Our results contribute to a better understanding of the realized genetic potential of pelagic N2 fixation and its seasonal dynamics in oligohaline estuarine ecosystems, which are natural coastal biogeochemical reactors.

Using full chloroplast genomes of \u2018red\u2019 and \u2018yellow\u2019 Bixa orellana (achiote) for kmer based identification and phylogenetic inference

BackgroundFull chloroplast genomes provide high resolution taxonomic discrimination between closely related plant species and are quickly replacing single and multi-locus barcoding regions as reference materials of choice for DNA based taxonomic annotation of plants. Bixa orellana, commonly known as “achiote” and “annatto” is a plant used for both human and animal foods and was thus identified for full chloroplast sequencing for the Center for Veterinary Medicine (CVM) Complete Chloroplast Animal Feed database. This work was conducted in collaboration with the Instituto de Medicina Tradicional (IMET) in Iquitos, Peru. There is a wide range of color variation in pods of Bixa orellana for which genetic loci that distinguish phenotypes have not yet been identified. Here we apply whole chloroplast genome sequencing of “red” and “yellow” individuals of Bixa orellana to provide high quality reference genomes to support kmer database development for use identifying this plant from complex mixtures using shotgun data. Additionally, we describe chloroplast gene content, synteny and phylogeny, and identify an indel and snp that may be associated with seed pod color.ResultsFully assembled chloroplast genomes were produced for both red and yellow Bixa orellana accessions (158,918 and 158,823 bp respectively). Synteny and gene content was identical to the only other previously reported full chloroplast genome of Bixa orellana (NC_041550). We observed a 17 base pair deletion at position 58,399–58,415 in both accessions, relative to NC_041550 and a 6 bp deletion at position 75,531–75,526 and a snp at position 86,493 in red Bixa orellana.ConclusionsOur data provide high quality reference genomes of individuals of red and yellow Bixa orellana to support kmer based identity markers for use with shotgun sequencing approaches for rapid, precise identification of Bixa orellana from complex mixtures. Kmer based phylogeny of full chloroplast genomes supports monophylly of Bixaceae consistent with alignment based approaches. A potentially discriminatory indel and snp were identified that may be correlated with the red phenotype.

Development of microsatellite markers for the threatened species Coleocephalocereus purpureus (Cactaceae) using next-generation sequencing.

Ten microsatellite loci were developed and validated for the endangered cactus species Coleocephalocereus purpureus. The markers were obtained from sequences generated by whole genome shotgun sequencing approaches. A testing group of 36 specimens of the main population were genotyped and all described markers presented suitable outcomes to population genetic studies, showing polymorphic status for C. purpureus testing group with clean and reproducible amplification. No evidence for scoring errors, null alleles or linkage disequilibrium was detected. Number of alleles per locus ranged from 3 to 6 and expected heterozygosity ranged from 0.78 to 0.99. These new microsatellite loci are suitable to be used in future diversity and structure population studies of C. purpureus.

Antimicrobial activity and carbohydrate metabolism in the bacterial metagenome of the soil-living invertebrate Folsomia candida

The microbiome associated with an animal’s gut and other organs is considered an integral part of its ecological functions and adaptive capacity. To better understand how microbial communities influence activities and capacities of the host, we need more information on the functions that are encoded in a microbiome. Until now, the information about soil invertebrate microbiomes is mostly based on taxonomic characterization, achieved through culturing and amplicon sequencing. Using shotgun sequencing and various bioinformatics approaches we explored functions in the bacterial metagenome associated with the soil invertebrate Folsomia candida, an established model organism in soil ecology with a fully sequenced, high-quality genome assembly. Our metagenome analysis revealed a remarkable diversity of genes associated with antimicrobial activity and carbohydrate metabolism. The microbiome also contains several homologs to F. candida genes that were previously identified as candidates for horizontal gene transfer (HGT). We suggest that the carbohydrate- and antimicrobial-related functions encoded by Folsomia’s metagenome play a role in the digestion of recalcitrant soil-born polysaccharides and the defense against pathogens, thereby significantly contributing to the adaptation of these animals to life in the soil. Furthermore, the transfer of genes from the microbiome may constitute an important source of new functions for the springtail.

Sequencing Library Preparation from Degraded Samples for Non-illumina Sequencing Platforms.

Efficient methods for building genomic sequencing libraries from degraded DNA have been in place for Illumina sequencing platforms for some years now, but such methods are still lacking for other sequencing platforms. Here, we provide a protocol for building genomic libraries from degraded DNA (archival or ancient sample material) for sequencing on the Ion Torrent™ high-throughput sequencing platforms. In addition to a reduction in time and cost in comparison to commercial kits, this protocol removes purification steps prior to library amplification, an important consideration for work involving historical samples. Libraries prepared using this method are appropriate for either shotgun sequencing or enrichment-based downstream approaches.

Combining 16S rRNA gene variable regions enables high-resolution microbial community profiling

BackgroundMost of our knowledge about the remarkable microbial diversity on Earth comes from sequencing the 16S rRNA gene. The use of next-generation sequencing methods has increased sample number and sequencing depth, but the read length of the most widely used sequencing platforms today is quite short, requiring the researcher to choose a subset of the gene to sequence (typically 16–33% of the total length). Thus, many bacteria may share the same amplified region, and the resolution of profiling is inherently limited. Platforms that offer ultra-long read lengths, whole genome shotgun sequencing approaches, and computational frameworks formerly suggested by us and by others all allow different ways to circumvent this problem yet suffer various shortcomings. There is a need for a simple and low-cost 16S rRNA gene-based profiling approach that harnesses the short read length to provide a much larger coverage of the gene to allow for high resolution, even in harsh conditions of low bacterial biomass and fragmented DNA.ResultsThis manuscript suggests Short MUltiple Regions Framework (SMURF), a method to combine sequencing results from different PCR-amplified regions to provide one coherent profiling. The de facto amplicon length is the total length of all amplified regions, thus providing much higher resolution compared to current techniques. Computationally, the method solves a convex optimization problem that allows extremely fast reconstruction and requires only moderate memory. We demonstrate the increase in resolution by in silico simulations and by profiling two mock mixtures and real-world biological samples. Reanalyzing a mock mixture from the Human Microbiome Project achieved about twofold improvement in resolution when combing two independent regions. Using a custom set of six primer pairs spanning about 1200 bp (80%) of the 16S rRNA gene, we were able to achieve ~ 100-fold improvement in resolution compared to a single region, over a mock mixture of common human gut bacterial isolates. Finally, the profiling of a Drosophila melanogaster microbiome using the set of six primer pairs provided a ~ 100-fold increase in resolution and thus enabling efficient downstream analysis.ConclusionsSMURF enables the identification of near full-length 16S rRNA gene sequences in microbial communities, having resolution superior compared to current techniques. It may be applied to standard sample preparation protocols with very little modifications. SMURF also paves the way to high-resolution profiling of low-biomass and fragmented DNA, e.g., in the case of formalin-fixed and paraffin-embedded samples, fossil-derived DNA, or DNA exposed to other degrading conditions. The approach is not restricted to combining amplicons of the 16S rRNA gene and may be applied to any set of amplicons, e.g., in multilocus sequence typing (MLST).

Microbial diversity within the airway microbiome in chronic pediatric lung diseases

The study of the airway microbiome in children is an area of emerging research, especially in relation to the role microbial diversity may play in acute and chronic inflammation. Three such pediatric airway diseases include cystic fibrosis, asthma, and chronic lung disease of prematurity. In cystic fibrosis, the presence of Pseudomonas spp. is associated with decreased microbial diversity. Decreasing microbial diversity is also associated with poor lung function. In asthma, early viral infections appear to drive changes in bacterial diversity which may be associated with asthma risk. Premature infants with Ureaplasma spp. are at higher risk for chronic lung disease due to inflammation. Microbiome changes due to prematurity also appear to affect the inflammatory response to viral infections post-natally. Importantly, microbial diversity can be measured using metataxonomic (e.g., 16S rRNA sequencing) and metagenomic (e.g., shotgun sequencing) approaches. A metagenomics approach may be preferable as it can provide further granularity of the sample composition, identifying the bacterial species or strain, information on additional microbial components, including fungal and viral components, information about functional genomics of the microbiome, and information about antimicrobial resistance mutations. Future studies of pediatric airway diseases incorporating these techniques may provide evidence for new treatment approaches for these vulnerable patient populations.