Optimizing sequencing protocols for leaderboard metagenomics by combining long and short reads

Dermal injury drives a skin to gut axis that disrupts the intestinal microbiome and intestinal immune homeostasis in mice

The composition of the microbial community in the intestine may influence the functions of distant organs such as the brain, lung, and skin. These microbes can promote disease or have beneficial functions, leading to the hypothesis that microbes in the gut explain the co-occurrence of intestinal and skin diseases. Here, we show that the reverse can occur, and that skin directly alters the gut microbiome. Disruption of the dermis by skin wounding or the digestion of dermal hyaluronan results in increased expression in the colon of the host defense genes Reg3 and Muc2, and skin wounding changes the composition and behavior of intestinal bacteria. Enhanced expression Reg3 and Muc2 is induced in vitro by exposure to hyaluronan released by these skin interventions. The change in the colon microbiome after skin wounding is functionally important as these bacteria penetrate the intestinal epithelium and enhance colitis from dextran sodium sulfate (DSS) as seen by the ability to rescue skin associated DSS colitis with oral antibiotics, in germ-free mice, and fecal microbiome transplantation to unwounded mice from mice with skin wounds. These observations provide direct evidence of a skin-gut axis by demonstrating that damage to the skin disrupts homeostasis in intestinal host defense and alters the gut microbiome.

Integration of polygenic and gut metagenomic risk prediction for common diseases

Multiomics has shown promise in noninvasive risk profiling and early detection of various common diseases. In the present study, in a prospective population-based cohort with ~18 years of e-health record follow-up, we investigated the incremental and combined value of genomic and gut metagenomic risk assessment compared with conventional risk factors for predicting incident coronary artery disease (CAD), type 2 diabetes (T2D), Alzheimer disease and prostate cancer. We found that polygenic risk scores (PRSs) improved prediction over conventional risk factors for all diseases. Gut microbiome scores improved predictive capacity over baseline age for CAD, T2D and prostate cancer. Integrated risk models of PRSs, gut microbiome scores and conventional risk factors achieved the highest predictive performance for all diseases studied compared with models based on conventional risk factors alone. The present study demonstrates that integrated PRSs and gut metagenomic risk models improve the predictive value over conventional risk factors for common chronic diseases.

A conserved interdomain microbial network underpins cadaver decomposition despite environmental variables

Microbial breakdown of organic matter is one of the most important processes on Earth, yet the controls of decomposition are poorly understood. Here we track 36 terrestrial human cadavers in three locations and show that a phylogenetically distinct, interdomain microbial network assembles during decomposition despite selection effects of location, climate and season. We generated a metagenome-assembled genome library from cadaver-associated soils and integrated it with metabolomics data to identify links between taxonomy and function. This universal network of microbial decomposers is characterized by cross-feeding to metabolize labile decomposition products. The key bacterial and fungal decomposers are rare across non-decomposition environments and appear unique to the breakdown of terrestrial decaying flesh, including humans, swine, mice and cattle, with insects as likely important vectors for dispersal. The observed lockstep of microbial interactions further underlies a robust microbial forensic tool with the potential to aid predictions of the time since death.

Occupational farm work activities influence workers' indoor home microbiome

BACKGROUND

Farm work entails a heterogeneous mixture of exposures that vary considerably across farms and farmers. Farm work is associated with various health outcomes, both adverse and beneficial. One mechanism by which farming exposures can impact health is through the microbiome, including the indoor home environment microbiome. It is unknown how individual occupational exposures shape the microbial composition in workers' homes.

OBJECTIVES

We investigated associations between farm work activities, including specific tasks and pesticide use, and the indoor microbiome in the homes of 468 male farmers.

METHODS

Participants were licensed pesticide applicators, mostly farmers, enrolled in the Agricultural Lung Health Study from 2008 to 2011. Vacuumed dust from participants' bedrooms underwent whole-genome shotgun sequencing for indoor microbiome assessment. Using questionnaire data, we evaluated 6 farm work tasks (processing of either hay, silage, animal feed, fertilizer, or soy/grains, and cleaning grain bins) and 19 pesticide ingredients currently used in the past year, plus 7 banned persistent pesticide ingredients ever used.

RESULTS

All 6 work tasks were associated with increased microbial diversity levels, with a positive dose-response for the total number of tasks performed (P = 0.001). All tasks were associated with altered microbial compositions (weighted UniFrac P = 0.001) and with higher abundance of specific microbes, including soil-based commensal microbes such as Haloterrigena. Among the 19 pesticides, current use of glyphosate and past use of lindane were associated with increased microbial diversity (P = 0.02-0.04). Ten currently used pesticides and all 7 banned pesticides were associated with altered microbial composition (P = 0.001-0.04). Six pesticides were associated with differential abundance of certain microbes.

DISCUSSION

Different farm activities and exposures can uniquely impact the dust microbiome inside homes. Our work suggests that changes to the home microbiome could serve as one pathway for how occupational exposures impact the health of workers and their cohabitating family members, offering possible future intervention targets.

Pangenome comparison of Bacteroides fragilis genomospecies unveil genetic diversity and ecological insights

Bacteroides fragilis is a Gram-negative commensal bacterium commonly found in the human colon that differentiates into two genomospecies termed division I and II. We leverage a comprehensive collection of 694 B. fragilis whole genome sequences and report differential gene abundance to further support the recent proposal that divisions I and II represent separate species. In division I strains, we identify an increased abundance of genes related to complex carbohydrate degradation, colonization, and host niche occupancy, confirming the role of division I strains as gut commensals. In contrast, division II strains display an increased prevalence of plant cell wall degradation genes and exhibit a distinct geographic distribution, primarily originating from Asian countries, suggesting dietary influences. Notably, division II strains have an increased abundance of genes linked to virulence, survival in toxic conditions, and antimicrobial resistance, consistent with a higher incidence of these strains in bloodstream infections. This study provides new evidence supporting a recent proposal for classifying divisions I and II B. fragilis strains as distinct species, and our comparative genomic analysis reveals their niche-specific roles.

Comparative analysis of endophytic fungal communities in bamboo species Phyllostachys edulis, Bambusa rigida, and Pleioblastus amarus

Fungal endophytes in plant leaf mesophyll form mutually beneficial associations through carbon assimilation, synthesis of biologically active chemicals, and enhancement of aesthetic and nutritional value. Here, we compared community structure, diversity, and richness of endophytic fungi in the leaves of three bamboo species, including Phyllostachys edulis (MZ), Bambusa rigida (KZ), and Pleioblastus amarus (YT) via high-throughput Illumina sequencing. In total, 1070 operational taxonomic units (OTUs) were retrieved and classified into 7 phylum, 27 classes, 82 orders, 185 families, 310 genus, and 448 species. Dominant genera were Cladosporium, Trichomerium, Hannaella, Ascomycota, Sporobolomyces, Camptophora and Strelitziana. The highest fungal diversity was observed in Pleioblastus amarus, followed by Bambusa rigida, and Phyllostachys edulis. Comparatively, monopodial species Ph. edulis and sympodial B. rigida, mixed P. amarus revealed the highest richness of endophytic fungi. We retrieved a few biocontrol agents, Sarocladium and Paraconiothyrium, and unique Sporobolomyces, Camptophora, and Strelitziana genera. FUNGuild analysis revealed the surrounding environment (The annual average temperature is between 15 and 25 °C, and the relative humidity of the air is above 83% all year round) as a source of fungal accumulation in bamboo leaves and their pathogenic nature. Our results provide precise knowledge for better managing bamboo forests and pave the way for isolating secondary metabolites and potential bioactive compounds.

Maximizing the potential of high-throughput next-generation sequencing through precise normalization based on read count distribution

Next-generation sequencing technologies have enabled many advances across diverse areas of biology, with many benefiting from increased sample size. Although the cost of running next-generation sequencing instruments has dropped substantially over time, the cost of sample preparation methods has lagged behind. To counter this, researchers have adapted library miniaturization protocols and large sample pools to maximize the number of samples that can be prepared by a certain amount of reagents and sequenced in a single run. However, due to high variability of sample quality, over and underrepresentation of samples in a sequencing run has become a major issue in high-throughput sequencing. This leads to misinterpretation of results due to increased noise, and additional time and cost rerunning underrepresented samples. To overcome this problem, we present a normalization method that uses shallow iSeq sequencing to accurately inform pooling volumes based on read distribution. This method is superior to the widely used fluorometry methods, which cannot specifically target adapter-ligated molecules that contribute to sequencing output. Our normalization method not only quantifies adapter-ligated molecules but also allows normalization of feature space; for example, we can normalize to reads of interest such as non-ribosomal reads. As a result, this normalization method improves the efficiency of high-throughput next-generation sequencing by reducing noise and producing higher average reads per sample with more even sequencing depth. IMPORTANCE High-throughput next generation sequencing (NGS) has significantly contributed to the field of genomics; however, further improvements can maximize the potential of this important tool. Uneven sequencing of samples in a multiplexed run is a common issue that leads to unexpected extra costs or low-quality data. To mitigate this problem, we introduce a normalization method based on read counts rather than library concentration. This method allows for an even distribution of features of interest across samples, improving the statistical power of data sets and preventing the financial loss associated with resequencing libraries. This method optimizes NGS, which already has huge importance across many areas of biology.

Occupational Farm Work Activities Influence Workers' Indoor Home Microbiome

Background

Farm work entails a heterogeneous mixture of exposures that vary considerably across farms and farmers. Farm work is associated with various health outcomes, both adverse and beneficial. One mechanism by which farming exposures can impact health is through the microbiome, including the indoor built environment microbiome. It is unknown how individual occupational exposures shape the microbial composition in workers' homes.

Objectives

We investigated associations between farm work activities, including specific tasks and pesticide use, and the indoor microbiome in the homes of 468 male farmers.

Methods

Participants were licensed pesticide applicators, mostly farmers, enrolled in the Agricultural Lung Health Study from 2008-2011. Vacuumed dust from participants' bedrooms underwent whole-genome shotgun sequencing for indoor microbiome assessment. Using questionnaire data, we evaluated 6 farm work tasks (processing of either hay, silage, animal feed, fertilizer, or soy/grains, and cleaning grain bins) and 19 pesticide ingredients currently used in the past year, plus 7 persistent banned pesticide ingredients ever used.

Results

All 6 work tasks were associated with increased within-sample microbial diversity, with a positive dose-response for the sum of tasks (p=0.001). All tasks were associated with altered overall microbial compositions (weighted UniFrac p=0.001) and with higher abundance of specific microbes, including soil-based microbes such as Haloterrigena. Among the 19 pesticides, only current use of glyphosate and past use of lindane were associated with increased within-sample diversity (p=0.02-0.04). Ten currently used pesticides and all 7 banned pesticides were associated with altered microbial composition (p=0.001-0.04). Six pesticides were associated with differential abundance of certain microbes.

Discussion

Specific farm activities and exposures can impact the dust microbiome inside homes. Our work suggests that occupational farm exposures could impact the health of workers and their families through modifying the indoor environment, specifically the microbial composition of house dust, offering possible future intervention targets.

Metagenomics reveals novel microbial signatures of farm exposures in house dust

Indoor home dust microbial communities, important contributors to human health, are shaped by environmental factors, including farm-related exposures. Advanced metagenomic whole genome shotgun sequencing (WGS) improves detection and characterization of microbiota in the indoor built-environment dust microbiome, compared to conventional 16S rRNA amplicon sequencing (16S). We hypothesized that the improved characterization of indoor dust microbial communities by WGS will enhance detection of exposure-outcome associations. The objective of this study was to identify novel associations of environmental exposures with the dust microbiome from the homes of 781 farmers and farm spouses enrolled in the Agricultural Lung Health Study. We examined various farm-related exposures, including living on a farm, crop versus animal production, and type of animal production, as well as non-farm exposures, including home cleanliness and indoor pets. We assessed the association of the exposures on within-sample alpha diversity and between-sample beta diversity, and the differential abundance of specific microbes by exposure. Results were compared to previous findings using 16S. We found most farm exposures were significantly positively associated with both alpha and beta diversity. Many microbes exhibited differential abundance related to farm exposures, mainly in the phyla Actinobacteria, Bacteroidetes, Firmicutes, and Proteobacteria. The identification of novel differential taxa associated with farming at the genera level, including Rhodococcus, Bifidobacterium, Corynebacterium, and Pseudomonas, was a benefit of WGS compared to 16S. Our findings indicate that characterization of dust microbiota, an important component of the indoor environment relevant to human health, is heavily influenced by sequencing techniques. WGS is a powerful tool to survey the microbial community that provides novel insights on the impact of environmental exposures on indoor dust microbiota. These findings can inform the design of future studies in environmental health.

Gut microbiota, circulating inflammatory markers and metabolites, and carotid artery atherosclerosis in HIV infection

BACKGROUND

Alterations in gut microbiota have been implicated in HIV infection and cardiovascular disease. However, how gut microbial alterations relate to host inflammation and metabolite profiles, and their relationships with atherosclerosis, have not been well-studied, especially in the context of HIV infection. Here, we examined associations of gut microbial species and functional components measured by shotgun metagenomics with carotid artery plaque assessed by B-mode carotid artery ultrasound in 320 women with or at high risk of HIV (65% HIV +) from the Women's Interagency HIV Study. We further integrated plaque-associated microbial features with serum proteomics (74 inflammatory markers measured by the proximity extension assay) and plasma metabolomics (378 metabolites measured by liquid chromatography tandem mass spectrometry) in relation to carotid artery plaque in up to 433 women.

RESULTS

Fusobacterium nucleatum, a potentially pathogenic bacteria, was positively associated with carotid artery plaque, while five microbial species (Roseburia hominis, Roseburia inulinivorans, Johnsonella ignava, Odoribacter splanchnicus, Clostridium saccharolyticum) were inversely associated with plaque. Results were consistent between women with and without HIV. Fusobacterium nucleatum was positively associated with several serum proteomic inflammatory markers (e.g., CXCL9), and the other plaque-related species were inversely associated with proteomic inflammatory markers (e.g., CX3CL1). These microbial-associated proteomic inflammatory markers were also positively associated with plaque. Associations between bacterial species (especially Fusobacterium nucleatum) and plaque were attenuated after further adjustment for proteomic inflammatory markers. Plaque-associated species were correlated with several plasma metabolites, including the microbial metabolite imidazole-propionate (ImP), which was positively associated with plaque and several pro-inflammatory markers. Further analysis identified additional bacterial species and bacterial hutH gene (encoding enzyme histidine ammonia-lyase in ImP production) associated with plasma ImP levels. A gut microbiota score based on these ImP-associated species was positively associated with plaque and several pro-inflammatory markers.

CONCLUSION

Among women living with or at risk of HIV, we identified several gut bacterial species and a microbial metabolite ImP associated with carotid artery atherosclerosis, which might be related to host immune activation and inflammation. Video Abstract.

Signatures of HIV and Major Depressive Disorder in the Plasma Microbiome

Inter-individual differences in the gut microbiome are linked to alterations in inflammation and blood-brain barrier permeability, which may increase the risk of depression in people with HIV (PWH). The microbiome profile of blood, which is considered by many to be typically sterile, remains largely unexplored. We aimed to characterize the blood plasma microbiome composition and assess its association with major depressive disorder (MDD) in PWH and people without HIV (PWoH). In this cross-sectional, observational cohort, we used shallow-shotgun metagenomic sequencing to characterize the plasma microbiome of 151 participants (84 PWH and 67 PWoH), all of whom underwent a comprehensive neuropsychiatric assessment. The microbial composition did not differ between PWH and PWoH or between participants with MDD and those without it. Using the songbird model, we computed the log ratio of the highest and lowest 30% of the ranked classes associated with HIV and MDD. We found that HIV infection and lifetime MDD were enriched in a set of differentially abundant inflammatory classes, such as Flavobacteria and Nitrospira. Our results suggest that the circulating plasma microbiome may increase the risk of MDD related to dysbiosis-induced inflammation in PWH. If confirmed, these findings may indicate new biological mechanisms that could be targeted to improve treatment of MDD in PWH.

Decoding dissolved information: environmental DNA sequencing at global scale to monitor a changing ocean

The use of environmental DNA (eDNA) technology for environmental monitoring is rapidly expanding, with applications for fisheries, coral reefs, harmful algal blooms, invasive and endangered species, and biodiversity monitoring. By enabling detection of species over space and time, eDNA fulfills a fundamental need of environmental surveys. Traditional surveys are expensive, require significant capital expenditure, and can be destructive; eDNA offers promise for cheaper, less invasive, and higher-resolution (i.e. genetic) assessments of environments and stocks. However, challenges in quantification, detection limits, biobanking capacity, reference databases, and data management and integration remain significant hurdles to efficient eDNA monitoring at global and decadal scale. Here, we consider the current state of eDNA technology and its suitability for the problems for which it is being used. We explore the current best practices, the logistical and social challenges that prevent eDNA from widespread adoption and benefit, and the emerging technologies that may address those challenges.

Metagenomics reveals novel microbial signatures of farm exposures in house dust

Indoor home dust microbial communities, important contributors to human health outcomes, are shaped by environmental factors, including farm-related exposures. Detection and characterization of microbiota are influenced by sequencing methodology; however, it is unknown if advanced metagenomic whole genome shotgun sequencing (WGS) can detect novel associations between environmental exposures and the indoor built-environment dust microbiome, compared to conventional 16S rRNA amplicon sequencing (16S). This study aimed to better depict indoor dust microbial communities using WGS to investigate novel associations with environmental risk factors from the homes of 781 farmers and farm spouses enrolled in the Agricultural Lung Health Study. We examined various farm-related exposures, including living on a farm, crop versus animal production, and type of animal production, as well as non-farm exposures, including home cleanliness and indoor pets. We assessed the association of the exposures on within-sample alpha diversity and between-sample beta diversity, and the differential abundance of specific microbes by exposure. Results were compared to previous findings using 16S. We found most farm exposures were significantly positively associated with both alpha and beta diversity. Many microbes exhibited differential abundance related to farm exposures, mainly in the phyla Actinobacteria, Bacteroidetes, Firmicutes , and Proteobacteria . The identification of novel differential taxa associated with farming at the genera level, including Rhodococcus, Bifidobacterium, Corynebacterium , and Pseudomonas , was a benefit of WGS compared to 16S. Our findings indicate that characterization of dust microbiota, an important component of the indoor environment relevant to human health, is heavily influenced by sequencing techniques. WGS is a powerful tool to survey the microbial community that provides novel insights on the impact of environmental exposures on indoor dust microbiota, and should be an important consideration in designing future studies in environmental health.

The gut microbiome is a significant risk factor for future chronic lung disease

BACKGROUND

The gut-lung axis is generally recognized, but there are few large studies of the gut microbiome and incident respiratory disease in adults.

OBJECTIVE

We sought to investigate the association and predictive capacity of the gut microbiome for incident asthma and chronic obstructive pulmonary disease (COPD).

METHODS

Shallow metagenomic sequencing was performed for stool samples from a prospective, population-based cohort (FINRISK02; N = 7115 adults) with linked national administrative health register-derived classifications for incident asthma and COPD up to 15 years after baseline. Generalized linear models and Cox regressions were used to assess associations of microbial taxa and diversity with disease occurrence. Predictive models were constructed using machine learning with extreme gradient boosting. Models considered taxa abundances individually and in combination with other risk factors, including sex, age, body mass index, and smoking status.

RESULTS

A total of 695 and 392 statistically significant associations were found between baseline taxonomic groups and incident asthma and COPD, respectively. Gradient boosting decision trees of baseline gut microbiome abundance predicted incident asthma and COPD in the validation data sets with mean area under the curves of 0.608 and 0.780, respectively. Cox analysis showed that the baseline gut microbiome achieved higher predictive performance than individual conventional risk factors, with C-indices of 0.623 for asthma and 0.817 for COPD. The integration of the gut microbiome and conventional risk factors further improved prediction capacities.

CONCLUSIONS

The gut microbiome is a significant risk factor for incident asthma and incident COPD and is largely independent of conventional risk factors.

A low-cost genomics workflow enables isolate screening and strain-level analyses within microbiomes

Earth's environments harbor complex consortia of microbes that affect processes ranging from host health to biogeochemical cycles. Understanding their evolution and function is limited by an inability to isolate genomes in a high-throughput manner. Here, we present a workflow for bacterial whole-genome sequencing using open-source labware and the OpenTrons robotics platform, reducing costs to approximately $10 per genome. We assess genomic diversity within 45 gut bacterial species from wild-living chimpanzees and bonobos. We quantify intraspecific genomic diversity and reveal divergence of homologous plasmids between hosts. This enables population genetic analyses of bacterial strains not currently possible with metagenomic data alone.

Towards standardized and reproducible research in skin microbiomes

Skin is a complex organ serving a critical role as a barrier and mediator of interactions between the human body and its environment. Recent studies have uncovered how resident microbial communities play a significant role in maintaining the normal healthy function of the skin and the immune system. In turn, numerous host-associated and environmental factors influence these communities' composition and diversity across the cutaneous surface. In addition, specific compositional changes in skin microbiota have also been connected to the development of several chronic diseases. The current era of microbiome research is characterized by its reliance on large data sets of nucleotide sequences produced with high-throughput sequencing of sample-extracted DNA. These approaches have yielded new insights into many previously uncharacterized microbial communities. Application of standardized practices in the study of skin microbial communities could help us understand their complex structures, functional capacities, and health associations and increase the reproducibility of the research. Here, we overview the current research in human skin microbiomes and outline challenges specific to their study. Furthermore, we provide perspectives on recent advances in methods, analytical tools and applications of skin microbiomes in medicine and forensics.

A comparison of six DNA extraction protocols for 16S, ITS and shotgun metagenomic sequencing of microbial communities

Microbial communities contain a broad phylogenetic diversity of organisms; however, the majority of methods center on describing bacteria and archaea. Fungi are important symbionts in many ecosystems and are potentially important members of the human microbiome, beyond those that can cause disease. To expand our analysis of microbial communities to include data from the fungal internal transcribed spacer (ITS) region, five candidate DNA extraction kits were compared against our standardized protocol for describing bacteria and archaea using 16S rRNA gene amplicon- and shotgun metagenomics sequencing. The results are presented considering a diverse panel of host-associated and environmental sample types and comparing the cost, processing time, well-to-well contamination, DNA yield, limit of detection and microbial community composition among protocols. Across all criteria, the MagMAX Microbiome kit was found to perform best. The PowerSoil Pro kit performed comparably but with increased cost per sample and overall processing time. The Zymo MagBead, NucleoMag Food and Norgen Stool kits were included.

Phylogeny-Aware Analysis of Metagenome Community Ecology Based on Matched Reference Genomes while Bypassing Taxonomy

We introduce the operational genomic unit (OGU) method, a metagenome analysis strategy that directly exploits sequence alignment hits to individual reference genomes as the minimum unit for assessing the diversity of microbial communities and their relevance to environmental factors. This approach is independent of taxonomic classification, granting the possibility of maximal resolution of community composition, and organizes features into an accurate hierarchy using a phylogenomic tree. The outputs are suitable for contemporary analytical protocols for community ecology, differential abundance, and supervised learning while supporting phylogenetic methods, such as UniFrac and phylofactorization, that are seldom applied to shotgun metagenomics despite being prevalent in 16S rRNA gene amplicon studies. As demonstrated in two real-world case studies, the OGU method produces biologically meaningful patterns from microbiome data sets. Such patterns further remain detectable at very low metagenomic sequencing depths. Compared with taxonomic unit-based analyses implemented in currently adopted metagenomics tools, and the analysis of 16S rRNA gene amplicon sequence variants, this method shows superiority in informing biologically relevant insights, including stronger correlation with body environment and host sex on the Human Microbiome Project data set and more accurate prediction of human age by the gut microbiomes of Finnish individuals included in the FINRISK 2002 cohort. We provide Woltka, a bioinformatics tool to implement this method, with full integration with the QIIME 2 package and the Qiita web platform, to facilitate adoption of the OGU method in future metagenomics studies. IMPORTANCE Shotgun metagenomics is a powerful, yet computationally challenging, technique compared to 16S rRNA gene amplicon sequencing for decoding the composition and structure of microbial communities. Current analyses of metagenomic data are primarily based on taxonomic classification, which is limited in feature resolution. To solve these challenges, we introduce operational genomic units (OGUs), which are the individual reference genomes derived from sequence alignment results, without further assigning them taxonomy. The OGU method advances current read-based metagenomics in two dimensions: (i) providing maximal resolution of community composition and (ii) permitting use of phylogeny-aware tools. Our analysis of real-world data sets shows that it is advantageous over currently adopted metagenomic analysis methods and the finest-grained 16S rRNA analysis methods in predicting biological traits. We thus propose the adoption of OGUs as an effective practice in metagenomic studies.

The impact of maternal asthma on the preterm infants' gut metabolome and microbiome (MAP study)

Preterm infants are at a greater risk for the development of asthma and atopic disease, which can lead to lifelong negative health consequences. This may be due, in part, to alterations that occur in the gut microbiome and metabolome during their stay in the Neonatal Intensive Care Unit (NICU). To explore the differential roles of family history (i.e., predisposition due to maternal asthma diagnosis) and hospital-related environmental and clinical factors that alter microbial exposures early in life, we considered a unique cohort of preterm infants born ≤ 34 weeks gestational age from two local level III NICUs, as part of the MAP (Microbiome, Atopic disease, and Prematurity) Study. From MAP participants, we chose a sub-cohort of infants whose mothers had a history of asthma and matched gestational age and sex to infants of mothers without a history of asthma diagnosis (control). We performed a prospective, paired metagenomic and metabolomic analysis of stool and milk feed samples collected at birth, 2 weeks, and 6 weeks postnatal age. Although there were clinical factors associated with shifts in the diversity and composition of stool-associated bacterial communities, maternal asthma diagnosis did not play an observable role in shaping the infant gut microbiome during the study period. There were significant differences, however, in the metabolite profile between the maternal asthma and control groups at 6 weeks postnatal age. The most notable changes occurred in the linoleic acid spectral network, which plays a role in inflammatory and immune pathways, suggesting early metabolomic changes in the gut of preterm infants born to mothers with a history of asthma. Our pilot study suggests that a history of maternal asthma alters a preterm infants’ metabolomic pathways in the gut, as early as the first 6 weeks of life.

Gut Microbiome Composition Is Predictive of Incident Type 2 Diabetes in a Population Cohort of 5,572 Finnish Adults

OBJECTIVE

To examine the previously unknown long-term association between gut microbiome composition and incident type 2 diabetes in a representative population cohort.

RESEARCH DESIGN AND METHODS

We collected fecal samples from 5,572 Finns (mean age 48.7 years; 54.1% women) in 2002 who were followed up for incident type 2 diabetes until 31 December 2017. The samples were sequenced using shotgun metagenomics. We examined associations between gut microbiome composition and incident diabetes using multivariable-adjusted Cox regression models. We first used the eastern Finland subpopulation to obtain initial findings and validated these in the western Finland subpopulation.

RESULTS

Altogether, 432 cases of incident diabetes occurred over the median follow-up of 15.8 years. We detected four species and two clusters consistently associated with incident diabetes in the validation models. These four species were Clostridium citroniae (hazard ratio [HR] 1.21; 95% CI 1.04-1.42), C. bolteae (HR 1.20; 95% CI 1.04-1.39), Tyzzerella nexilis (HR 1.17; 95% CI 1.01-1.36), and Ruminococcus gnavus (HR 1.17; 95% CI 1.01-1.36). The positively associated clusters, cluster 1 (HR 1.18; 95% CI 1.02-1.38) and cluster 5 (HR 1.18; 95% CI 1.02-1.36), mostly consisted of these same species.

CONCLUSIONS

We observed robust species-level taxonomic features predictive of incident type 2 diabetes over long-term follow-up. These findings build on and extend previous mainly cross-sectional evidence and further support links between dietary habits, metabolic diseases, and type 2 diabetes that are modulated by the gut microbiome. The gut microbiome can potentially be used to improve disease prediction and uncover novel therapeutic targets for diabetes.

Disentangle genus microdiversity within a complex microbial community by using a multi-distance long-read binning method: example of Candidatus Accumulibacter

Complete genomes can be recovered from metagenomes by assembling and binning DNA sequences into metagenome assembled genomes (MAGs). Yet, the presence of microdiversity can hamper the assembly and binning processes, possibly yielding chimeric, highly fragmented and incomplete genomes. Here, the metagenomes of four samples of aerobic granular sludge bioreactors containing Candidatus (Ca.) Accumulibacter, a phosphate-accumulating organism of interest for wastewater treatment, were sequenced with both PacBio and Illumina. Different strategies of genome assembly and binning were investigated, including published protocols and a binning procedure adapted to the binning of long contigs (MuLoBiSC). Multiple criteria were considered to select the best strategy for Ca. Accumulibacter, whose multiple strains in every sample represent a challenging microdiversity. In this case, the best strategy relies on long-read only assembly and a custom binning procedure including MuLoBiSC in metaWRAP. Several high-quality Ca. Accumulibacter MAGs, including a novel species, were obtained independently from different samples. Comparative genomic analysis showed that MAGs retrieved in different samples harbour genomic rearrangements in addition to accumulation of point mutations. The microdiversity of Ca. Accumulibacter, likely driven by mobile genetic elements, causes major difficulties in recovering MAGs, but it is also a hallmark of the panmictic lifestyle of these bacteria.

Enhancing biological signals and detection rates in single-cell RNA-seq experiments with cDNA library equalization

Considerable effort has been devoted to refining experimental protocols to reduce levels of technical variability and artifacts in single-cell RNA-sequencing data (scRNA-seq). We here present evidence that equalizing the concentration of cDNA libraries prior to pooling, a step not consistently performed in single-cell experiments, improves gene detection rates, enhances biological signals, and reduces technical artifacts in scRNA-seq data. To evaluate the effect of equalization on various protocols, we developed Scaffold, a simulation framework that models each step of an scRNA-seq experiment. Numerical experiments demonstrate that equalization reduces variation in sequencing depth and gene-specific expression variability. We then performed a set of experiments in vitro with and without the equalization step and found that equalization increases the number of genes that are detected in every cell by 17-31%, improves discovery of biologically relevant genes, and reduces nuisance signals associated with cell cycle. Further support is provided in an analysis of publicly available data.

Standardized multi-omics of Earth's microbiomes reveals microbial and metabolite diversity

Despite advances in sequencing, lack of standardization makes comparisons across studies challenging and hampers insights into the structure and function of microbial communities across multiple habitats on a planetary scale. Here we present a multi-omics analysis of a diverse set of 880 microbial community samples collected for the Earth Microbiome Project. We include amplicon (16S, 18S, ITS) and shotgun metagenomic sequence data, and untargeted metabolomics data (liquid chromatography-tandem mass spectrometry and gas chromatography mass spectrometry). We used standardized protocols and analytical methods to characterize microbial communities, focusing on relationships and co-occurrences of microbially related metabolites and microbial taxa across environments, thus allowing us to explore diversity at extraordinary scale. In addition to a reference database for metagenomic and metabolomic data, we provide a framework for incorporating additional studies, enabling the expansion of existing knowledge in the form of an evolving community resource. We demonstrate the utility of this database by testing the hypothesis that every microbe and metabolite is everywhere but the environment selects. Our results show that metabolite diversity exhibits turnover and nestedness related to both microbial communities and the environment, whereas the relative abundances of microbially related metabolites vary and co-occur with specific microbial consortia in a habitat-specific manner. We additionally show the power of certain chemistry, in particular terpenoids, in distinguishing Earth's environments (for example, terrestrial plant surfaces and soils, freshwater and marine animal stool), as well as that of certain microbes including Conexibacter woesei (terrestrial soils), Haloquadratum walsbyi (marine deposits) and Pantoea dispersa (terrestrial plant detritus). This Resource provides insight into the taxa and metabolites within microbial communities from diverse habitats across Earth, informing both microbial and chemical ecology, and provides a foundation and methods for multi-omics microbiome studies of hosts and the environment.

Post-weaning shifts in microbiome composition and metabolism revealed by over 25 000 pig gut metagenome-assembled genomes

Using a previously described metagenomics dataset of 27 billion reads, we reconstructed over 50 000 metagenome-assembled genomes (MAGs) of organisms resident in the porcine gut, 46.5 % of which were classified as >70 % complete with a <10 % contamination rate, and 24.4 % were nearly complete genomes. Here, we describe the generation and analysis of those MAGs using time-series samples. The gut microbial communities of piglets appear to follow a highly structured developmental programme in the weeks following weaning, and this development is robust to treatments including an intramuscular antibiotic treatment and two probiotic treatments. The high resolution we obtained allowed us to identify specific taxonomic ‘signatures’ that characterize the gut microbial development immediately after weaning. Additionally, we characterized the carbohydrate repertoire of the organisms resident in the porcine gut. We tracked the abundance shifts of 294 carbohydrate active enzymes, and identified the species and higher-level taxonomic groups carrying each of these enzymes in their MAGs. This knowledge can contribute to the design of probiotics and prebiotic interventions as a means to modify the piglet gut microbiome.

The Potential Role of Clinical Metagenomics in Infectious Diseases: Therapeutic Perspectives

Clinical metagenomics (CMg) is the process of sequencing nucleic acid of clinical samples to obtain clinically relevant information such as the identification of microorganisms and their susceptibility to antimicrobials. Over the last decades, sequencing and bioinformatic solutions supporting CMg have much evolved and an increasing number of case reports and series covering various infectious diseases have been published. Metagenomics is a new approach to infectious disease diagnosis that is currently being developed and is certainly one of the most promising for the coming years. However, most CMg studies are retrospective, and few address the potential impact CMg could have on patient management, including initiation, adaptation, or cessation of antimicrobials. In this narrative review, we have discussed the potential role of CMg in bacteriology, virology, mycology, and parasitology. Several reports and case-series confirm that CMg is an innovative tool with which one can (i) identify more microorganisms than with conventional methods in a single test, (ii) obtain results within hours, and (iii) tailor the antimicrobial regimen of patients. However, the cost-efficiency of CMg and its real impact on patient management are still to be determined.

Taxonomic signatures of cause-specific mortality risk in human gut microbiome

The collection of fecal material and developments in sequencing technologies have enabled standardised and non-invasive gut microbiome profiling. Microbiome composition from several large cohorts have been cross-sectionally linked to various lifestyle factors and diseases. In spite of these advances, prospective associations between microbiome composition and health have remained uncharacterised due to the lack of sufficiently large and representative population cohorts with comprehensive follow-up data. Here, we analyse the long-term association between gut microbiome variation and mortality in a well-phenotyped and representative population cohort from Finland (n = 7211). We report robust taxonomic and functional microbiome signatures related to the Enterobacteriaceae family that are associated with mortality risk during a 15-year follow-up. Our results extend previous cross-sectional studies, and help to establish the basis for examining long-term associations between human gut microbiome composition, incident outcomes, and general health status.

Spinal Cord Injury Changes the Structure and Functional Potential of Gut Bacterial and Viral Communities

Emerging data indicate that gut dysbiosis contributes to many human diseases, including several comorbidities that develop after traumatic spinal cord injury (SCI). To date, all analyses of SCI-induced gut dysbiosis have used 16S rRNA amplicon sequencing. This technique has several limitations, including being susceptible to taxonomic "blind spots," primer bias, and an inability to profile microbiota functions or identify viruses. Here, SCI-induced gut dysbiosis was assessed by applying genome- and gene-resolved metagenomic analysis of murine stool samples collected 21 days after an experimental SCI at the 4th thoracic spine (T4) or 10th thoracic spine (T10) spinal level. These distinct injuries partially (T10) or completely (T4) abolish sympathetic tone in the gut. Among bacteria, 105 medium- to high-quality metagenome-assembled genomes (MAGs) were recovered, with most (n = 96) representing new bacterial species. Read mapping revealed that after SCI, the relative abundance of beneficial commensals (Lactobacillus johnsonii and CAG-1031 spp.) decreased, while potentially pathogenic bacteria (Weissella cibaria, Lactococcus lactis _A, Bacteroides thetaiotaomicron) increased. Functionally, microbial genes encoding proteins for tryptophan, vitamin B6, and folate biosynthesis, essential pathways for central nervous system function, were reduced after SCI. Among viruses, 1,028 mostly novel viral populations were recovered, expanding known murine gut viral species sequence space ∼3-fold compared to that of public databases. Phages of beneficial commensal hosts (CAG-1031, Lactobacillus, and Turicibacter) decreased, while phages of pathogenic hosts (Weissella, Lactococcus, and class Clostridia) increased after SCI. Although the microbiomes and viromes were changed in all SCI mice, some of these changes varied as a function of spinal injury level, implicating loss of sympathetic tone as a mechanism underlying gut dysbiosis.IMPORTANCE To our knowledge, this is the first article to apply metagenomics to characterize changes in gut microbial population dynamics caused by a clinically relevant model of central nervous system (CNS) trauma. It also utilizes the most current approaches in genome-resolved metagenomics and viromics to maximize the biological inferences that can be made from these data. Overall, this article highlights the importance of autonomic nervous system regulation of a distal organ (gut) and its microbiome inhabitants after traumatic spinal cord injury (SCI). By providing information on taxonomy, function, and viruses, metagenomic data may better predict how SCI-induced gut dysbiosis influences systemic and neurological outcomes after SCI.

Evaluation of the Effect of Storage Methods on Fecal, Saliva, and Skin Microbiome Composition

Expanding our knowledge of microbial communities across diverse environments includes collecting samples in places far from the laboratory. Identifying cost-effective preservatives that will enable room temperature storage of microbial communities for sequencing analysis is crucial to enabling microbiome analyses across diverse populations.

As the number of human microbiome studies expand, it is increasingly important to identify cost-effective, practical preservatives that allow for room temperature sample storage. Here, we reanalyzed 16S rRNA gene amplicon sequencing data from a large sample storage study published in 2016 and performed shotgun metagenomic sequencing on remnant DNA from this experiment. Both results support the initial findings that 95% ethanol, a nontoxic, cost-effective preservative, is effective at preserving samples at room temperature for weeks. We expanded on this analysis by collecting a new set of fecal, saliva, and skin samples to determine the optimal ratio of 95% ethanol to sample. We identified optimal collection protocols for fecal samples (storing a fecal swab in 95% ethanol) and saliva samples (storing unstimulated saliva in 95% ethanol at a ratio of 1:2). Storing skin swabs in 95% ethanol reduced microbial biomass and disrupted community composition, highlighting the difficulties of low biomass sample preservation. The results from this study identify practical solutions for large-scale analyses of fecal and oral microbial communities.

IMPORTANCE Expanding our knowledge of microbial communities across diverse environments includes collecting samples in places far from the laboratory. Identifying cost-effective preservatives that will enable room temperature storage of microbial communities for sequencing analysis is crucial to enabling microbiome analyses across diverse populations. Here, we validate findings that 95% ethanol efficiently preserves microbial composition at room temperature for weeks. We also identified the optimal ratio of 95% ethanol to sample for stool and saliva to preserve both microbial load and composition. These results provide rationale for an accessible, nontoxic, cost-effective solution that will enable crowdsourcing microbiome studies, such as The Microsetta Initiative, and lower the barrier for collecting diverse samples.

Recovery of complete genomes and non-chromosomal replicons from activated sludge enrichment microbial communities with long read metagenome sequencing

New long read sequencing technologies offer huge potential for effective recovery of complete, closed genomes from complex microbial communities. Using long read data (ONT MinION) obtained from an ensemble of activated sludge enrichment bioreactors we recover 22 closed or complete genomes of community members, including several species known to play key functional roles in wastewater bioprocesses, specifically microbes known to exhibit the polyphosphate- and glycogen-accumulating organism phenotypes (namely Candidatus Accumulibacter and Dechloromonas, and Micropruina, Defluviicoccus and Candidatus Contendobacter, respectively), and filamentous bacteria (Thiothrix) associated with the formation and stability of activated sludge flocs. Additionally we demonstrate the recovery of close to 100 circularised plasmids, phages and small microbial genomes from these microbial communities using long read assembled sequence. We describe methods for validating long read assembled genomes using their counterpart short read metagenome-assembled genomes, and assess the influence of different correction procedures on genome quality and predicted gene quality. Our findings establish the feasibility of performing long read metagenome-assembled genome recovery for both chromosomal and non-chromosomal replicons, and demonstrate the value of parallel sampling of moderately complex enrichment communities to obtaining high quality reference genomes of key functional species relevant for wastewater bioprocesses.

A comparison of DNA/RNA extraction protocols for high-throughput sequencing of microbial communities

One goal of microbial ecology researchers is to capture the maximum amount of information from all organisms in a sample. The recent COVID-19 pandemic, caused by the RNA virus SARS-CoV-2, has highlighted a gap in traditional DNA-based protocols, including the high-throughput methods the authors previously established as field standards. To enable simultaneous SARS-CoV-2 and microbial community profiling, the authors compared the relative performance of two total nucleic acid extraction protocols with the authors' previously benchmarked protocol. The authors included a diverse panel of environmental and host-associated sample types, including body sites commonly swabbed for COVID-19 testing. Here the authors present results comparing the cost, processing time, DNA and RNA yield, microbial community composition, limit of detection and well-to-well contamination between these protocols.

More for less: predicting and maximizing genomic variant discovery via Bayesian nonparametrics

While the cost of sequencing genomes has decreased dramatically in recent years, this expense often remains nontrivial. Under a fixed budget, scientists face a natural trade-off between quantity and quality: spending resources to sequence a greater number of genomes or spending resources to sequence genomes with increased accuracy. Our goal is to find the optimal allocation of resources between quantity and quality. Optimizing resource allocation promises to reveal as many new variations in the genome as possible. We introduce a Bayesian nonparametric methodology to predict the number of new variants in a follow-up study based on a pilot study. When experimental conditions are kept constant between the pilot and follow-up, we find that our prediction is competitive with the best existing methods. Unlike current methods, though, our new method allows practitioners to change experimental conditions between the pilot and the follow-up. We demonstrate how this distinction allows our method to be used for more realistic predictions and for optimal allocation of a fixed budget between quality and quantity. We validate our method on cancer and human genomics data.

Links between gut microbiome composition and fatty liver disease in a large population sample

Fatty liver disease is the most common liver disease in the world. Its connection with the gut microbiome has been known for at least 80 y, but this association remains mostly unstudied in the general population because of underdiagnosis and small sample sizes. To address this knowledge gap, we studied the link between the Fatty Liver Index (FLI), a well-established proxy for fatty liver disease, and gut microbiome composition in a representative, ethnically homogeneous population sample of 6,269 Finnish participants. We based our models on biometric covariates and gut microbiome compositions from shallow metagenome sequencing. Our classification models could discriminate between individuals with a high FLI (≥60, indicates likely liver steatosis) and low FLI (<60) in internal cross-region validation, consisting of 30% of the data not used in model training, with an average AUC of 0.75 and AUPRC of 0.56 (baseline at 0.30). In addition to age and sex, our models included differences in 11 microbial groups from class Clostridia, mostly belonging to orders Lachnospirales and Oscillospirales. Our models were also predictive of the high FLI group in a different Finnish cohort, consisting of 258 participants, with an average AUC of 0.77 and AUPRC of 0.51 (baseline at 0.21). Pathway analysis of representative genomes of the positively FLI-associated taxa in (NCBI) Clostridium subclusters IV and XIVa indicated the presence of, e.g., ethanol fermentation pathways. These results support several findings from smaller case-control studies, such as the role of endogenous ethanol producers in the development of the fatty liver.

Translocation of Viable Gut Microbiota to Mesenteric Adipose Drives Formation of Creeping Fat in Humans

A mysterious feature of Crohn's disease (CD) is the extra-intestinal manifestation of "creeping fat" (CrF), defined as expansion of mesenteric adipose tissue around the inflamed and fibrotic intestine. In the current study, we explore whether microbial translocation in CD serves as a central cue for CrF development. We discovered a subset of mucosal-associated gut bacteria that consistently translocated and remained viable in CrF in CD ileal surgical resections, and identified Clostridium innocuum as a signature of this consortium with strain variation between mucosal and adipose isolates, suggesting preference for lipid-rich environments. Single-cell RNA sequencing characterized CrF as both pro-fibrotic and pro-adipogenic with a rich milieu of activated immune cells responding to microbial stimuli, which we confirm in gnotobiotic mice colonized with C. innocuum. Ex vivo validation of expression patterns suggests C. innocuum stimulates tissue remodeling via M2 macrophages, leading to an adipose tissue barrier that serves to prevent systemic dissemination of bacteria.

Bacterial modification of the host glycosaminoglycan heparan sulfate modulates SARS-CoV-2 infectivity

The human microbiota has a close relationship with human disease and it remodels components of the glycocalyx including heparan sulfate (HS). Studies of the severe acute respiratory syndrome coronavirus (SARS-CoV-2) spike protein receptor binding domain suggest that infection requires binding to HS and angiotensin converting enzyme 2 (ACE2) in a codependent manner. Here, we show that commensal host bacterial communities can modify HS and thereby modulate SARS-CoV-2 spike protein binding and that these communities change with host age and sex. Common human-associated commensal bacteria whose genomes encode HS-modifying enzymes were identified. The prevalence of these bacteria and the expression of key microbial glycosidases in bronchoalveolar lavage fluid (BALF) was lower in adult COVID-19 patients than in healthy controls. The presence of HS-modifying bacteria decreased with age in two large survey datasets, FINRISK 2002 and American Gut, revealing one possible mechanism for the observed increase in COVID-19 susceptibility with age. In vitro , bacterial glycosidases from unpurified culture media supernatants fully blocked SARS-CoV-2 spike binding to human H1299 protein lung adenocarcinoma cells. HS-modifying bacteria in human microbial communities may regulate viral adhesion, and loss of these commensals could predispose individuals to infection. Understanding the impact of shifts in microbial community composition and bacterial lyases on SARS-CoV-2 infection may lead to new therapeutics and diagnosis of susceptibility.

Microbiome analyses of blood and tissues suggest cancer diagnostic approach

Systematic characterization of the cancer microbiome provides the opportunity to develop techniques that exploit non-human, microorganism-derived molecules in the diagnosis of a major human disease. Following recent demonstrations that some types of cancer show substantial microbial contributions1-10, we re-examined whole-genome and whole-transcriptome sequencing studies in The Cancer Genome Atlas11 (TCGA) of 33 types of cancer from treatment-naive patients (a total of 18,116 samples) for microbial reads, and found unique microbial signatures in tissue and blood within and between most major types of cancer. These TCGA blood signatures remained predictive when applied to patients with stage Ia-IIc cancer and cancers lacking any genomic alterations currently measured on two commercial-grade cell-free tumour DNA platforms, despite the use of very stringent decontamination analyses that discarded up to 92.3% of total sequence data. In addition, we could discriminate among samples from healthy, cancer-free individuals (n = 69) and those from patients with multiple types of cancer (prostate, lung, and melanoma; 100 samples in total) solely using plasma-derived, cell-free microbial nucleic acids. This potential microbiome-based oncology diagnostic tool warrants further exploration.