A comparison of short-read, HiFi long-read, and hybrid strategies for genome-resolved metagenomics

Shotgun metagenomics enables the reconstruction of complex microbial communities at a high level of detail. Such an approach can be conducted using both short-read and long-read sequencing data, as well as a combination of both. To assess the pros and cons of these different approaches, we used 22 fecal DNA extracts collected weekly for 11 weeks from two respective lab mice to study seven performance metrics over four combinations of sequencing depth and technology: (i) 20 Gbp of Illumina short-read data, (ii) 40 Gbp of short-read data, (iii) 20 Gbp of PacBio HiFi long-read data, and (iv) 40 Gbp of hybrid (20 Gbp of short-read +20 Gbp of long-read) data. No strategy was best for all metrics; instead, each one excelled across different metrics. The long-read approach yielded the best assembly statistics, with the highest N50 and lowest number of contigs. The 40 Gbp short-read approach yielded the highest number of refined bins. Finally, the hybrid approach yielded the longest assemblies and the highest mapping rate to the bacterial genomes. Our results suggest that while long-read sequencing significantly improves the quality of reconstructed bacterial genomes, it is more expensive and requires deeper sequencing than short-read approaches to recover a comparable amount of reconstructed genomes. The most optimal strategy is study-specific and depends on how researchers assess the trade-off between the quantity and quality of recovered genomes.IMPORTANCEMice are an important model organism for understanding the gut microbiome. When studying these gut microbiomes using DNA techniques, researchers can choose from technologies that use short or long DNA reads. In this study, we perform an extensive benchmark between short- and long-read DNA sequencing for studying mice gut microbiomes. We find that no one approach was best for all metrics and provide information that can help guide researchers in planning their experiments.

The prevalence of controls in phyllosphere microbiome research: a methodological review

DNA contamination can critically confound microbiome studies. Here, we take a systematic approach to review the current literature and investigate the prevalence of contamination controls in phyllosphere microbiome research over the past decade. By utilising systematic review principles for this review, we were able to conduct a thorough investigation, screening 450 articles from three databases for eligibility and extracting data in a controlled and methodical manner. Worryingly, we observed a surprisingly low usage of both positive and negative contamination controls in phyllosphere research. As a result, we propose a set of minimum standards to combat the effects of contamination in future phyllosphere research.

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.

Mammals show distinct functional gut microbiome dynamics to identical series of environmental stressors

IMPORTANCE

In our manuscript, we report the first interspecific comparative study about the plasticity of the gut microbiota. We conducted a captivity experiment that exposed wild-captured mammals to a series of environmental challenges over 45 days. We characterized their gut microbial communities using genome-resolved metagenomics and modeled how the taxonomic, phylogenetic, and functional microbial dynamics varied across a series of disturbances in both species. Our results indicate that the intrinsic properties (e.g., diversity and functional redundancy) of microbial communities coupled with physiological attributes (e.g., thermal plasticity) of hosts shape the taxonomic, phylogenetic, and functional response of gut microbiomes to environmental stressors, which might influence their contribution to the acclimation and adaptation capacity of animal hosts.

Impact of microbial genome completeness on metagenomic functional inference

Inferring the functional capabilities of bacteria from metagenome-assembled genomes (MAGs) is becoming a central process in microbiology. Here we show that the completeness of genomes has a significant impact on the recovered functional signal, spanning all domains of metabolic functions. We identify factors that affect this relationship between genome completeness and function fullness, and provide baseline knowledge to guide efforts to correct for this overlooked bias in metagenomic functional inference.

Individuality and stability of the koala (Phascolarctos cinereus) faecal microbiota through time

Gut microbiota studies often rely on a single sample taken per individual, representing a snapshot in time. However, we know that gut microbiota composition in many animals exhibits intra-individual variation over the course of days to months. Such temporal variations can be a confounding factor in studies seeking to compare the gut microbiota of different wild populations, or to assess the impact of medical/veterinary interventions. To date, little is known about the variability of the koala (Phascolarctos cinereus) gut microbiota through time. Here, we characterise the gut microbiota from faecal samples collected at eight timepoints over a month for a captive population of South Australian koalas (n individuals = 7), and monthly over 7 months for a wild population of New South Wales koalas (n individuals = 5). Using 16S rRNA gene sequencing, we found that microbial diversity was stable over the course of days to months. Each koala had a distinct faecal microbiota composition which in the captive koalas was stable across days. The wild koalas showed more variation across months, although each individual still maintained a distinct microbial composition. Per koala, an average of 57 (±16) amplicon sequence variants (ASVs) were detected across all time points; these ASVs accounted for an average of 97% (±1.9%) of the faecal microbial community per koala. The koala faecal microbiota exhibits stability over the course of days to months. Such knowledge will be useful for future studies comparing koala populations and developing microbiota interventions for this regionally endangered marsupial.

Sarcoptic mange changes bacterial and fungal microbiota of bare-nosed wombats (Vombatus ursinus)

Background

Sarcoptes scabiei is globally distributed and one of the most impactful mammalian ectoparasites. Sarcoptic mange, caused by infection with S. scabiei, causes disruption of the epidermis and its bacterial microbiota, but its effects on host fungal microbiota and on the microbiota of marsupials in general have not been studied. Here, we (i) examine bacterial and fungal microbiota changes associated with mange in wild bare-nosed wombats (BNWs) and (ii) evaluate whether opportunistic pathogens are potentiated by S. scabiei infection in this species.

Methods

Using Amplicon Sequencing of the 16S rRNA and ITS2 rDNA genes, we detected skin microbiota changes of the bare-nosed wombat (Vombatus ursinus). We compared the alpha and beta diversity among healthy, moderate, and severe disease states using ANOVA and PERMANOVA with nesting. Lastly, we identified taxa that differed between disease states using analysis of composition of microbes (ANCOM) testing.

Results

We detected significant changes in the microbial communities and diversity with mange in BNWs. Severely affected BNWs had lower amplicon sequence variant (ASV) richness compared to that of healthy individuals, and the microbial communities were significantly different between disease states with higher relative abundance of potentially pathogenic microbial taxa in mange-affected BNWs including Staphylococcus sciuri, Corynebacterium spp., Brevibacterium spp., Brachybacterium spp., and Pseudogymnascus spp. and Debaryomyces spp.

Conclusion

This study represents the first investigation of microbial changes in association with sarcoptic mange in a marsupial host, as well as the first investigation of fungal microbial changes on the skin of any host suffering from sarcoptic mange. Our results are broadly consistent with bacterial microbiota changes observed in humans, pigs, canids, and Iberian ibex, suggesting the epidermal microbial impacts of mange may be generalisable across host species. We recommend that future studies investigating skin microbiota changes include both bacterial and fungal data to gain a more complete picture of the effects of sarcoptic mange.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05452-y.

Characterising the Gut Microbiomes in Wild and Captive Short-Beaked Echidnas Reveals Diet-Associated Changes

The gut microbiome plays a vital role in health and wellbeing of animals, and an increasing number of studies are investigating microbiome changes in wild and managed populations to improve conservation and welfare. The short-beaked echidna (Tachyglossus aculeatus) is an iconic Australian species, the most widespread native mammal, and commonly held in zoos. Echidnas are cryptic animals, and much is still unknown about many aspects of their biology. Furthermore, some wild echidna populations are under threat, while echidnas held in captivity can have severe gastric health problems. Here, we used citizen science and zoos to collect echidna scats from across Australia to perform the largest gut microbiome study on any native Australian animal. Using 16S rRNA gene metabarcoding of scat samples, we characterised and compared the gut microbiomes of echidnas in wild (n = 159) and managed (n = 44) populations, which were fed four different diets. Wild echidna samples were highly variable, yet commonly dominated by soil and plant-fermenting bacteria, while echidnas in captivity were dominated by gut commensals and plant-fermenting bacteria, suggesting plant matter may play a significant role in echidna diet. This work demonstrates significant differences between zoo held and wild echidnas, as well as managed animals on different diets, revealing that diet is important in shaping the gut microbiomes in echidnas. This first analysis of echidna gut microbiome highlights extensive microbial diversity in wild echidnas and changes in microbiome composition in managed populations. This is a first step towards using microbiome analysis to better understand diet, gastrointestinal biology, and improve management in these iconic animals.

Microbial biogeography of the wombat gastrointestinal tract

Most herbivorous mammals have symbiotic microbes living in their gastrointestinal tracts that help with harvesting energy from recalcitrant plant fibre. The bulk of research into these microorganisms has focused on samples collected from faeces, representing the distal region of the gastrointestinal (GI) tract. However, the GI tract in herbivorous mammals is typically long and complex, containing different regions with distinct physico-chemical properties that can structure resident microbial communities. Little work has been done to document GI microbial communities of herbivorous animals at these sites. In this study, we use 16S rRNA gene sequencing to characterize the microbial biogeography along the GI tract in two species of wombats. Specifically, we survey the microbes along four major gut regions (stomach, small intestine, proximal colon, distal colon) in a single bare-nosed wombat (Vombatus ursinus) and a single southern hairy-nosed wombat (Lasiorhinus latifrons). Our preliminary results show that GI microbial communities of wombats are structured by GI region. For both wombat individuals, we observed a trend of increasing microbial diversity from stomach to distal colon. The microbial composition in the first proximal colon region was more similar between wombat species than the corresponding distal colon region in the same species. We found several microbial genera that were differentially abundant between the first proximal colon (putative site for primary plant fermentation) and distal colon regions (which resemble faecal samples). Surprisingly, only 10.6% (98) and 18.8% (206) of amplicon sequence variants (ASVs) were shared between the first proximal colon region and the distal colon region for the bare-nosed and southern hairy-nosed wombat, respectively. These results suggest that microbial communities in the first proximal colon region-the putative site of primary plant fermentation in wombats-are distinct from the distal colon, and that faecal samples may have limitations in capturing the diversity of these communities. While faeces are still a valuable and effective means of characterising the distal colon microbiota, future work seeking to better understand how GI microbiota impact the energy economy of wombats (and potentially other hindgut-fermenting mammals) may need to take gut biogeography into account.

Paleo-diatom composition from Santa Barbara Basin deep-sea sediments: a comparison of 18S-V9 and diat-rbcL metabarcoding vs shotgun metagenomics

Sedimentary ancient DNA (sedaDNA) analyses are increasingly used to reconstruct marine ecosystems. The majority of marine sedaDNA studies use a metabarcoding approach (extraction and analysis of specific DNA fragments of a defined length), targeting short taxonomic marker genes. Promising examples are 18S-V9 rRNA (~121-130 base pairs, bp) and diat-rbcL (76 bp), targeting eukaryotes and diatoms, respectively. However, it remains unknown how 18S-V9 and diat-rbcL derived compositional profiles compare to metagenomic shotgun data, the preferred method for ancient DNA analyses as amplification biases are minimised. We extracted DNA from five Santa Barbara Basin sediment samples (up to ~11 000 years old) and applied both a metabarcoding (18S-V9 rRNA, diat-rbcL) and a metagenomic shotgun approach to (i) compare eukaryote, especially diatom, composition, and (ii) assess sequence length and database related biases. Eukaryote composition differed considerably between shotgun and metabarcoding data, which was related to differences in read lengths (~112 and ~161 bp, respectively), and overamplification of short reads in metabarcoding data. Diatom composition was influenced by reference bias that was exacerbated in metabarcoding data and characterised by increased representation of Chaetoceros, Thalassiosira and Pseudo-nitzschia. Our results are relevant to sedaDNA studies aiming to accurately characterise paleo-ecosystems from either metabarcoding or metagenomic data.

Effectiveness of decontamination protocols when analyzing ancient DNA preserved in dental calculus

Ancient DNA analysis of human oral microbial communities within calcified dental plaque (calculus) has revealed key insights into human health, paleodemography, and cultural behaviors. However, contamination imposes a major concern for paleomicrobiological samples due to their low endogenous DNA content and exposure to environmental sources, calling into question some published results. Decontamination protocols (e.g. an ethylenediaminetetraacetic acid (EDTA) pre-digestion or ultraviolet radiation (UV) and 5% sodium hypochlorite immersion treatments) aim to minimize the exogenous content of the outer surface of ancient calculus samples prior to DNA extraction. While these protocols are widely used, no one has systematically compared them in ancient dental calculus. Here, we compare untreated dental calculus samples to samples from the same site treated with four previously published decontamination protocols: a UV only treatment; a 5% sodium hypochlorite immersion treatment; a pre-digestion in EDTA treatment; and a combined UV irradiation and 5% sodium hypochlorite immersion treatment. We examine their efficacy in ancient oral microbiota recovery by applying 16S rRNA gene amplicon and shotgun sequencing, identifying ancient oral microbiota, as well as soil and skin contaminant species. Overall, the EDTA pre-digestion and a combined UV irradiation and 5% sodium hypochlorite immersion treatment were both effective at reducing the proportion of environmental taxa and increasing oral taxa in comparison to untreated samples. This research highlights the importance of using decontamination procedures during ancient DNA analysis of dental calculus to reduce contaminant DNA.

Host reproductive cycle influences the pouch microbiota of wild southern hairy-nosed wombats (Lasiorhinus latifrons)

Background

Marsupials are born much earlier than placental mammals, with most crawling from the birth canal to the protective marsupium (pouch) to further their development. However, little is known about the microbiology of the pouch and how it changes throughout a marsupial’s reproductive cycle. Here, using stringent controls, we characterized the microbial composition of multiple body sites from 26 wild Southern Hairy-nosed Wombats (SHNWs), including pouch samples from animals at different reproductive stages.

Results

Using qPCR of the 16S rRNA gene we detected a microbial community in the SHNW pouch. We observed significant differences in microbial composition and diversity between the body sites tested, as well as between pouch samples from different reproductive stages. The pouches of reproductively active females had drastically lower microbial diversity (mean ASV richness 19 ± 8) compared to reproductively inactive females (mean ASV richness 941 ± 393) and were dominated by gram positive bacteria from the Actinobacteriota phylum (81.7–90.6%), with the dominant families classified as Brevibacteriaceae, Corynebacteriaceae, Microbacteriaceae, and Dietziaceae. Three of the five most abundant sequences identified in reproductively active pouches had closest matches to microbes previously isolated from tammar wallaby pouches.

Conclusions

This study represents the first contamination-controlled investigation into the marsupial pouch microbiota, and sets a rigorous framework for future pouch microbiota studies. Our results indicate that SHNW pouches contain communities of microorganisms that are substantially altered by the host reproductive cycle. We recommend further investigation into the roles that pouch microorganisms may play in marsupial reproductive health and joey survival.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-021-00074-8.

Signatures of landscape and captivity in the gut microbiota of Southern Hairy-nosed Wombats (Lasiorhinus latifrons)

BACKGROUND

Herbivorous mammals co-opt microbes to derive energy and nutrients from diets that are recalcitrant to host enzymes. Recent research has found that captive management-an important conservation tool for many species-can alter the gut microbiota of mammals. Such changes could negatively impact the ability of herbivorous mammals to derive energy from their native diets, and ultimately reduce host fitness. To date, nothing is known of how captivity influences the gut microbiota of the Southern Hairy-nosed Wombat (SHNW), a large herbivorous marsupial that inhabits South Australia. Here, using 16S rRNA gene sequencing, we characterized the faecal microbiota of SHNWs in captivity and from three wild populations, two from degraded habitats and one from an intact native grass habitat.

RESULTS

We found that captive SHNWs had gut microbiota that were compositionally different and less diverse compared to wild SHNWs. There were major differences in gut microbiota community membership between captive and wild animals, both in statistically significant changes in relative abundance of microbes, and in the presence/absence of microbes. We also observed differences in microbial composition between wild populations, with the largest difference associated with native vs. degraded habitat.

CONCLUSIONS

These results suggest that captivity has a major impact on the gut microbiota of SHNWs, and that different wild populations harbour distinct microbial compositions. Such findings warrant further work to determine what impacts these changes have on the fitness of SHNWs, and whether they could be manipulated to improve future management of the species.

Investigating the demographic history of Japan using ancient oral microbiota

While microbial communities in the human body (microbiota) are now commonly associated with health and disease in industrialised populations, we know very little about how these communities co-evolved and changed with humans throughout history and deep prehistory. We can now examine these communities by sequencing ancient DNA preserved within calcified dental plaque (calculus), providing insights into the origins of disease and their links to human history. Here, we examine ancient DNA preserved within dental calculus samples and their associations with two major cultural periods in Japan: the Jomon period hunter–gatherers approximately 3000 years before present (BP) and the Edo period agriculturalists 400–150 BP. We investigate how human oral microbiomes have changed in Japan through time and explore the presence of microorganisms associated with oral diseases (e.g. periodontal disease, dental caries) in ancient Japanese populations. Finally, we explore oral microbial strain diversity and its potential links to ancient demography in ancient Japan by performing phylogenomic analysis of a widely conserved oral species—Anaerolineaceae oral taxon 439. This research represents, to our knowledge, the first study of ancient oral microbiomes from Japan and demonstrates that the analysis of ancient dental calculus can provide key information about the origin of non-infectious disease and its deep roots with human demography.

This article is part of the theme issue ‘Insights into health and disease from ancient biomolecules’.

An optimized method for the extraction of ancient eukaryote DNA from marine sediments

Marine sedimentary ancient DNA (sedaDNA) provides a powerful means to reconstruct marine palaeo-communities across the food web. However, currently there are few optimized sedaDNA extraction protocols available to maximize the yield of small DNA fragments typical of ancient DNA (aDNA) across a broad diversity of eukaryotes. We compared seven combinations of sedaDNA extraction treatments and sequencing library preparations using marine sediments collected at a water depth of 104 m off Maria Island, Tasmania, in 2018. These seven methods contrasted frozen versus refrigerated sediment, bead-beating induced cell lysis versus ethylenediaminetetraacetic acid (EDTA) incubation, DNA binding in silica spin columns versus in silica-solution, diluted versus undiluted DNA in shotgun library preparations to test potential inhibition issues during amplification steps, and size-selection of low molecular-weight (LMW) DNA to increase the extraction efficiency of sedaDNA. Maximum efficiency was obtained from frozen sediments subjected to a combination of EDTA incubation and bead-beating, DNA binding in silica-solution, and undiluted DNA in shotgun libraries, across 45 marine eukaryotic taxa. We present an optimized extraction protocol integrating these steps, with an optional post-library LMW size-selection step to retain DNA fragments of ≤500 base pairs. We also describe a stringent bioinformatic filtering approach for metagenomic data and provide a comprehensive list of contaminants as a reference for future sedaDNA studies. The new extraction and data-processing protocol should improve quantitative paleo-monitoring of eukaryotes from marine sediments, as well as other studies relying on the detection of highly fragmented and degraded eukaryote DNA in sediments.

Microbiome applications for pathology: challenges of low microbial biomass samples during diagnostic testing

The human microbiome can play key roles in disease, and diagnostic testing will soon have the ability to examine these roles in the context of clinical applications. Currently, most diagnostic testing in pathology applications focuses on a small number of disease‐causing microbes and dismisses the whole microbial community that causes or is modulated by disease. Microbiome modifications have already provided clinically relevant insights in gut and oral diseases, such as irritable bowel disease, but there are currently limitations when clinically examining microbiomes outside of these body sites. This is critical, as the majority of microbial samples used in pathology originate from body sites that contain low concentrations of microbial DNA, including skin, tissue, blood, and urine. These samples, also known as low microbial biomass samples, are difficult to examine without careful consideration and precautions to mitigate contamination and biases. Here, we present the limitations when analysing low microbial biomass samples using current protocols and techniques and highlight the advantages that microbiome testing can offer diagnostics in the future, if the proper precautions are implemented. Specifically, we discuss the sources of contamination and biases that may result in false assessments for these sample types. Finally, we provide recommendations to mitigate contamination and biases from low microbial biomass samples during diagnostic testing, which will be especially important to effectively diagnose and treat patients using microbiome analyses.

A new home for microbes

Modern microorganisms growing in fossils provide major challenges for researchers trying to detect ancient molecules in the same fossils.

Laboratory contamination over time during low-biomass sample analysis

Bacteria are not only ubiquitous on earth but can also be incredibly diverse within clean laboratories and reagents. The presence of both living and dead bacteria in laboratory environments and reagents is especially problematic when examining samples with low endogenous content (e.g., skin swabs, tissue biopsies, ice, water, degraded forensic samples or ancient material), where contaminants can outnumber endogenous microorganisms within samples. The contribution of contaminants within high‐throughput studies remains poorly understood because of the relatively low number of contaminant surveys. Here, we examined 144 negative control samples (extraction blank and no‐template amplification controls) collected in both typical molecular laboratories and an ultraclean ancient DNA laboratory over 5 years to characterize long‐term contaminant diversity. We additionally compared the contaminant content within a home‐made silica‐based extraction method, commonly used to analyse low endogenous content samples, with a widely used commercial DNA extraction kit. The contaminant taxonomic profile of the ultraclean ancient DNA laboratory was unique compared to modern molecular biology laboratories, and changed over time according to researcher, month and season. The commercial kit also contained higher microbial diversity and several human‐associated taxa in comparison to the home‐made silica extraction protocol. We recommend a minimum of two strategies to reduce the impacts of laboratory contaminants within low‐biomass metagenomic studies: (a) extraction blank controls should be included and sequenced with every batch of extractions and (b) the contributions of laboratory contamination should be assessed and reported in each high‐throughput metagenomic study.

Assessing alignment-based taxonomic classification of ancient microbial DNA

The field of palaeomicrobiology-the study of ancient microorganisms-is rapidly growing due to recent methodological and technological advancements. It is now possible to obtain vast quantities of DNA data from ancient specimens in a high-throughput manner and use this information to investigate the dynamics and evolution of past microbial communities. However, we still know very little about how the characteristics of ancient DNA influence our ability to accurately assign microbial taxonomies (i.e. identify species) within ancient metagenomic samples. Here, we use both simulated and published metagenomic data sets to investigate how ancient DNA characteristics affect alignment-based taxonomic classification. We find that nucleotide-to-nucleotide, rather than nucleotide-to-protein, alignments are preferable when assigning taxonomies to short DNA fragment lengths routinely identified within ancient specimens (<60 bp). We determine that deamination (a form of ancient DNA damage) and random sequence substitutions corresponding to ∼100,000 years of genomic divergence minimally impact alignment-based classification. We also test four different reference databases and find that database choice can significantly bias the results of alignment-based taxonomic classification in ancient metagenomic studies. Finally, we perform a reanalysis of previously published ancient dental calculus data, increasing the number of microbial DNA sequences assigned taxonomically by an average of 64.2-fold and identifying microbial species previously unidentified in the original study. Overall, this study enhances our understanding of how ancient DNA characteristics influence alignment-based taxonomic classification of ancient microorganisms and provides recommendations for future palaeomicrobiological studies.