Disease-specific loss of microbial cross-feeding interactions in the human gut

Many gut microorganisms critical to human health rely on nutrients produced by each other for survival; however, these cross-feeding interactions are still challenging to quantify and remain poorly characterized. Here, we introduce a Metabolite Exchange Score (MES) to quantify those interactions. Using metabolic models of prokaryotic metagenome-assembled genomes from over 1600 individuals, MES allows us to identify and rank metabolic interactions that are significantly affected by a loss of cross-feeding partners in 10 out of 11 diseases. When applied to a Crohn's disease case-control study, our approach identifies a lack of species with the ability to consume hydrogen sulfide as the main distinguishing microbiome feature of disease. We propose that our conceptual framework will help prioritize in-depth analyses, experiments and clinical targets, and that targeting the restoration of microbial cross-feeding interactions is a promising mechanism-informed strategy to reconstruct a healthy gut ecosystem.

Bacterial clade-specific analysis identifies distinct epithelial responses in inflammatory bowel disease

Abnormal immune responses to the resident gut microbiome can drive inflammatory bowel disease (IBD). Here, we combine high-resolution, culture-based shotgun metagenomic sequencing and analysis with matched host transcriptomics across three intestinal sites (terminal ileum, cecum, rectum) from pediatric IBD (PIBD) patients (n = 58) and matched controls (n = 42) to investigate this relationship. Combining our site-specific approach with bacterial culturing, we establish a cohort-specific bacterial culture collection, comprising 6,620 isolates (170 distinct species, 32 putative novel), cultured from 286 mucosal biopsies. Phylogeny-based, clade-specific metagenomic analysis identifies key, functionally distinct Enterococcus clades associated with either IBD or health. Strain-specific in vitro validation demonstrates differences in cell cytotoxicity and inflammatory signaling in intestinal epithelial cells, consistent with the colonic mucosa-specific response measured in patients with IBD. This demonstrates the importance of strain-specific phenotypes and consideration of anatomical sites in exploring the dysregulated host-bacterial interactions in IBD.

expam-high-resolution analysis of metagenomes using distance trees

Summary

Shotgun metagenomic sequencing provides the capacity to understand microbial community structure and function at unprecedented resolution; however, the current analytical methods are constrained by a focus on taxonomic classifications that may obfuscate functional relationships. Here, we present expam, a tree-based, taxonomy agnostic tool for the identification of biologically relevant clades from shotgun metagenomic sequencing.

Availability and implementation

expam is an open-source Python application released under the GNU General Public Licence v3.0. expam installation instructions, source code and tutorials can be found at https://github.com/seansolari/expam.

Supplementary information

Supplementary data are available at Bioinformatics online.

Critical Assessment of Metagenome Interpretation: the second round of challenges

Evaluating metagenomic software is key for optimizing metagenome interpretation and focus of the Initiative for the Critical Assessment of Metagenome Interpretation (CAMI). The CAMI II challenge engaged the community to assess methods on realistic and complex datasets with long- and short-read sequences, created computationally from around 1,700 new and known genomes, as well as 600 new plasmids and viruses. Here we analyze 5,002 results by 76 program versions. Substantial improvements were seen in assembly, some due to long-read data. Related strains still were challenging for assembly and genome recovery through binning, as was assembly quality for the latter. Profilers markedly matured, with taxon profilers and binners excelling at higher bacterial ranks, but underperforming for viruses and Archaea. Clinical pathogen detection results revealed a need to improve reproducibility. Runtime and memory usage analyses identified efficient programs, including top performers with other metrics. The results identify challenges and guide researchers in selecting methods for analyses.

This study presents the results of the second round of the Critical Assessment of Metagenome Interpretation challenges (CAMI II), which is a community-driven effort for comprehensively benchmarking tools for metagenomics data analysis.

Characterization of the Gut Microbiome and Resistomes of Wild and Zoo-Captive Macaques

Rhesus macaques (Macaca mulatta) are the most widely distributed species of Old World monkey and are frequently used as animal models to study human health and disease. Their gastrointestinal microbial community likely plays a major role in their physiology, ecology and evolution. Herein, we compared the fecal microbiome and antibiotic resistance genes in 15 free-ranging and 81 zoo-captive rhesus macaques sampled from two zoos in China, using both 16S amplicon sequencing and whole genome shotgun DNA sequencing approaches. Our data revealed similar levels of microbial diversity/richness among the three groups, although the composition of each group differed significantly and were particularly marked between the two zoo-captive and one wild groups. Zoo-captive animals also demonstrated a greater abundance and diversity of antibiotic genes. Through whole genome shotgun sequencing we also identified a mammalian (simian) associated adenovirus. Overall, this study provides a comprehensive analysis of resistomes and microbiomes in zoo-captive and free-ranging monkeys, revealing that semi-captive wildlife might harbor a higher diversity of antimicrobial resistant genes.

The microbiome and host mucosal interactions in urinary tract diseases

The urinary tract consists of the bladder, ureters, and kidneys, and is an essential organ system for filtration and excretion of waste products and maintaining systemic homeostasis. In this capacity, the urinary tract is impacted by its interactions with other mucosal sites, including the genitourinary and gastrointestinal systems. Each of these sites harbors diverse ecosystems of microbes termed the microbiota, that regulates complex interactions with the local and systemic immune system. It remains unclear whether changes in the microbiota and associated metabolites may be a consequence or a driver of urinary tract diseases. Here, we review the current literature, investigating the impact of the microbiota on the urinary tract in homeostasis and disease including urinary stones, acute kidney injury, chronic kidney disease, and urinary tract infection. We propose new avenues for exploration of the urinary microbiome using emerging technology and discuss the potential of microbiome-based medicine for urinary tract conditions.

Key Technologies for Progressing Discovery of Microbiome-Based Medicines

A growing number of experimental and computational approaches are illuminating the “microbial dark matter” and uncovering the integral role of commensal microbes in human health. Through this work, it is now clear that the human microbiome presents great potential as a therapeutic target for a plethora of diseases, including inflammatory bowel disease, diabetes and obesity. The development of more efficacious and targeted treatments relies on identification of causal links between the microbiome and disease; with future progress dependent on effective links between state-of-the-art sequencing approaches, computational analyses and experimental assays. We argue determining causation is essential, which can be attained by generating hypotheses using multi-omic functional analyses and validating these hypotheses in complex, biologically relevant experimental models. In this review we discuss existing analysis and validation methods, and propose best-practice approaches required to enable the next phase of microbiome research.

Genomic adaptations to an endolithic lifestyle in the coral-associated alga Ostreobium

The green alga Ostreobium is an important coral holobiont member, playing key roles in skeletal decalcification and providing photosynthate to bleached corals that have lost their dinoflagellate endosymbionts. Ostreobium lives in the coral's skeleton, a low-light environment with variable pH and O₂ availability. We present the Ostreobium nuclear genome and a metatranscriptomic analysis of healthy and bleached corals to improve our understanding of Ostreobium's adaptations to its extreme environment and its roles as a coral holobiont member. The Ostreobium genome has 10,663 predicted protein-coding genes and shows adaptations for life in low and variable light conditions and other stressors in the endolithic environment. This alga presents a rich repertoire of light-harvesting complex proteins but lacks many genes for photoprotection and photoreceptors. It also has a large arsenal of genes for oxidative stress response. An expansion of extracellular peptidases suggests that Ostreobium may supplement its energy needs by feeding on the organic skeletal matrix, and a diverse set of fermentation pathways allows it to live in the anoxic skeleton at night. Ostreobium depends on other holobiont members for vitamin B12, and our metatranscriptomes identify potential bacterial sources. Metatranscriptomes showed Ostreobium becoming a dominant agent of photosynthesis in bleached corals and provided evidence for variable responses among coral samples and different Ostreobium genotypes. Our work provides a comprehensive understanding of the adaptations of Ostreobium to its extreme environment and an important genomic resource to improve our comprehension of coral holobiont resilience, bleaching, and recovery.

CCMetagen: comprehensive and accurate identification of eukaryotes and prokaryotes in metagenomic data

There is an increasing demand for accurate and fast metagenome classifiers that can not only identify bacteria, but all members of a microbial community. We used a recently developed concept in read mapping to develop a highly accurate metagenomic classification pipeline named CCMetagen. The pipeline substantially outperforms other commonly used software in identifying bacteria and fungi and can efficiently use the entire NCBI nucleotide collection as a reference to detect species with incomplete genome data from all biological kingdoms. CCMetagen is user-friendly, and the results can be easily integrated into microbial community analysis software for streamlined and automated microbiome studies.

The use of taxon-specific reference databases compromises metagenomic classification

A recent article in BMC Genomics describes a new bioinformatics tool, HumanMycobiomeScan, to classify fungal taxa in metagenomic samples. This tool was used to characterize the gut mycobiome of hunter-gatherers and Western populations, resulting in the identification of a range of fungal species in the vast majority of samples. In the HumanMycobiomeScan pipeline, sequence reads are mapped against a reference database containing fungal genome sequences only. We argue that using reference databases comprised of a single taxonomic group leads to an unacceptably high number of false-positives due to: (i) mapping to conserved genetic regions in reference genomes, and (ii) sequence contamination in the assembled reference genomes. To demonstrate this, we replaced the HumanMycobiomeScan's fungal reference database with one containing genome sequences of amphibians and reptiles and re-analysed their case study. The classification pipeline recovered all species present in the reference database, revealing turtles (Geoemydidae), bull frogs (Pyxicephalidae) and snakes (Colubridae) as the most abundant herpetological taxa in the human gut. We also re-analysed their case study using a kingdom-agnostic pipeline. This revealed that while the gut of hunter-gatherers and Western subjects may be colonized by a range of microbial eukaryotes, only three fungal families were retrieved. These results highlight the pitfalls of using taxon-specific reference databases for metagenome classification, even when they are comprised of curated whole genome data. We propose that databases containing all domains of life provide the most suitable option for metagenomic species profiling, especially when targeting microbial eukaryotes.

Meta-transcriptomics reveals a diverse antibiotic resistance gene pool in avian microbiomes

BACKGROUND

Antibiotic resistance is rendering common bacterial infections untreatable. Wildlife can incorporate and disperse antibiotic-resistant bacteria in the environment, such as water systems, which in turn serve as reservoirs of resistance genes for human pathogens. Anthropogenic activity may contribute to the spread of bacterial resistance cycling through natural environments, including through the release of human waste, as sewage treatment only partially removes antibiotic-resistant bacteria. However, empirical data supporting these effects are currently limited. Here we used bulk RNA-sequencing (meta-transcriptomics) to assess the diversity and expression levels of functionally viable resistance genes in the gut microbiome of birds with aquatic habits in diverse locations.

RESULTS

We found antibiotic resistance genes in birds from all localities, from penguins in Antarctica to ducks in a wastewater treatment plant in Australia. Comparative analysis revealed that birds feeding at the wastewater treatment plant carried the greatest resistance gene burden, including genes typically associated with multidrug resistance plasmids as the aac(6)-Ib-cr gene. Differences in resistance gene burden also reflected aspects of bird ecology, taxonomy, and microbial function. Notably, ducks, which feed by dabbling, carried a higher abundance and diversity of resistance genes than turnstones, avocets, and penguins, which usually prey on more pristine waters.

CONCLUSIONS

These transcriptome data suggest that human waste, even if it undergoes treatment, might contribute to the spread of antibiotic resistance genes to the wild. Differences in microbiome functioning across different bird lineages may also play a role in the antibiotic resistance burden carried by wild birds. In summary, we reveal the complex factors explaining the distribution of resistance genes and their exchange routes between humans and wildlife, and show that meta-transcriptomics is a valuable tool to access functional resistance genes in whole microbial communities.

Whole Genome Sequencing of Australian Candida glabrata Isolates Reveals Genetic Diversity and Novel Sequence Types

Candida glabrata is a pathogen with reduced susceptibility to azoles and echinocandins. Analysis by traditional multilocus sequence typing (MLST) has recognized an increasing number of sequence types (STs), which vary with geography. Little is known about STs of C. glabrata in Australia. Here, we utilized whole genome sequencing (WGS) to study the genetic diversity of 51 Australian C. glabrata isolates and sought associations between STs over two time periods (2002-2004, 2010-2017), and with susceptibility to fluconazole by principal component analysis (PCA). Antifungal susceptibility was determined using Sensititre YeastOne^TM Y010 methodology and WGS performed on the NextSeq 500 platform (Illumina) with in silico MLST STs inferred by WGS data. Single nucleotide polymorphisms (SNPs) in genes linked to echinocandin, azole and 5-fluorocytosine resistance were analyzed. Of 51 isolates, WGS identified 18 distinct STs including four novel STs (ST123, ST124, ST126, and ST127). Four STs accounted for 49% of isolates (ST3, 15.7%; ST83, 13.7%; ST7, 9.8%; ST26, 9.8%). Split-tree network analysis resolved isolates to terminal branches; many of these comprised multiple isolates from disparate geographic settings but four branches contained Australian isolates only. ST3 isolates were common in Europe, United States and now Australia, whilst ST8 and ST19, relatively frequent in the United States, were rare/absent amongst our isolates. There was no association between ST distribution (genomic similarity) and the two time periods or with fluconazole susceptibility. WGS identified mutations in the FKS1 (S629P) and FKS2 (S663P) genes in three, and one, echinocandin-resistant isolate(s), respectively. Both mutations confer phenotypic drug resistance. Twenty-five percent (13/51) of isolates were fluconazole-resistant (MIC ≥ 64 μg/ml) of which 9 (18%) had non wild-type MICs to voriconazole and posaconazole. Multiple SNPs were present in genes linked to azole resistance such as CgPDR1 and CgCDR1, as well as several in MSH2; however, SNPs occurred in both azole-susceptible and azole-resistant isolates. Although no particular SNP in these genes was definitively associated with resistance, azole-resistant/non-wild type isolates had a propensity to harbor SNPs resulting in amino acid substitutions in Pdr1 beyond the first 250 amino acid positions. The presence of SNPs may be markers of STs. Our study shows the value of WGS for high-resolution sequence typing of C. glabrata, discovery of novel STs and potential to monitor trends in genetic diversity. WGS assessment for echinocandin resistance augments phenotypic susceptibility testing.

Large Diversity of Nonstandard Genes and Dynamic Evolution of Chloroplast Genomes in Siphonous Green Algae (Bryopsidales, Chlorophyta)

Chloroplast genomes have undergone tremendous alterations through the evolutionary history of the green algae (Chloroplastida). This study focuses on the evolution of chloroplast genomes in the siphonous green algae (order Bryopsidales). We present five new chloroplast genomes, which along with existing sequences, yield a data set representing all but one families of the order. Using comparative phylogenetic methods, we investigated the evolutionary dynamics of genomic features in the order. Our results show extensive variation in chloroplast genome architecture and intron content. Variation in genome size is accounted for by the amount of intergenic space and freestanding open reading frames that do not show significant homology to standard plastid genes. We show the diversity of these nonstandard genes based on their conserved protein domains, which are often associated with mobile functions (reverse transcriptase/intron maturase, integrases, phage- or plasmid-DNA primases, transposases, integrases, ligases). Investigation of the introns showed proliferation of group II introns in the early evolution of the order and their subsequent loss in the core Halimedineae, possibly through RT-mediated intron loss.

Highly structured prokaryote communities exist within the skeleton of coral colonies

Physiological performance, disease and bleaching prevalence are often patchy within individual coral colonies. These responses are largely influenced by coral-associated microbes, but how the coral microbiome changes over small spatial scales has never been quantified before. We performed a high-resolution quantification of the spatial scale of microbial species turnover (β-diversity) within skeletons of boulder-forming Porites corals. We found very strong prokaryotic species turnover across spatial scales ranging from 4 mm to 2 m within individual colonies, possibly resulting from dispersal limitation and microbial interactions. The microalgal community was more homogeneously distributed, which is likely due to these photosymbionts actively boring through limestone. Our findings highlight unprecedented levels of intra-colony heterogeneity in the skeletal microbiome, which has implications for the experimental design of coral microbiome studies and for our understanding of coral resilience.

Phylogenetic position of the coral symbiont Ostreobium (Ulvophyceae) inferred from chloroplast genome data

The green algal genus Ostreobium is an important symbiont of corals, playing roles in reef decalcification and providing photosynthates to the coral during bleaching events. A chloroplast genome of a cultured strain of Ostreobium was available, but low taxon sampling and Ostreobium's early-branching nature left doubt about its phylogenetic position. Here, we generate and describe chloroplast genomes from four Ostreobium strains as well as Avrainvillea mazei and Neomeris sp., strategically sampled early-branching lineages in the Bryopsidales and Dasycladales respectively. At 80,584 bp, the chloroplast genome of Ostreobium sp. HV05042 is the most compact yet found in the Ulvophyceae. The Avrainvillea chloroplast genome is ~94 kbp and contains introns in infA and cysT that have nearly complete sequence identity except for an open reading frame (ORF) in infA that is not present in cysT. In line with other bryopsidalean species, it also contains regions with possibly bacteria-derived ORFs. The Neomeris data did not assemble into a canonical circular chloroplast genome but a large number of contigs containing fragments of chloroplast genes and showing evidence of long introns and intergenic regions, and the Neomeris chloroplast genome size was estimated to exceed 1.87 Mb. Chloroplast phylogenomics and 18S nrDNA data showed strong support for the Ostreobium lineage being sister to the remaining Bryopsidales. There were differences in branch support when outgroups were varied, but the overall support for the placement of Ostreobium was strong. These results permitted us to validate two suborders and introduce a third, the Ostreobineae.

Evolutionary Dynamics of Chloroplast Genomes in Low Light: A Case Study of the Endolithic Green Alga Ostreobium quekettii

Some photosynthetic organisms live in extremely low light environments. Light limitation is associated with selective forces as well as reduced exposure to mutagens, and over evolutionary timescales it can leave a footprint on species’ genomes. Here, we present the chloroplast genomes of four green algae (Bryopsidales, Ulvophyceae), including the endolithic (limestone-boring) alga Ostreobium quekettii, which is a low light specialist. We use phylogenetic models and comparative genomic tools to investigate whether the chloroplast genome of Ostreobium corresponds to our expectations of how low light would affect genome evolution. Ostreobium has the smallest and most gene-dense chloroplast genome among Ulvophyceae reported to date, matching our expectation that light limitation would impose resource constraints reflected in the chloroplast genome architecture. Rates of molecular evolution are significantly slower along the phylogenetic branch leading to Ostreobium, in agreement with the expected effects of low light and energy levels on molecular evolution. We expected the ability of Ostreobium to perform photosynthesis in very low light to be associated with positive selection in genes related to the photosynthetic machinery, but instead, we observed that these genes may be under stronger purifying selection. Besides shedding light on the genome dynamics associated with a low light lifestyle, this study helps to resolve the role of environmental factors in shaping the diversity of genome architectures observed in nature.

Ecological niche models of invasive seaweeds

Ecological niche models (ENMs) are commonly used to calculate habitat suitability from species' occurrence and macroecological data. In invasive species biology, ENMs can be applied to anticipate whether invasive species are likely to establish in an area, to identify critical routes and arrival points, to build risk maps and to predict the extent of potential spread following an introduction. Most studies using ENMs focus on terrestrial organisms and applications in the marine realm are still relatively rare. Here, we review some common methods to build ENMs and their application in seaweed invasion biology. We summarize methods and concepts involved in the development of niche models, show examples of how they have been applied in studies on algae and discuss the application of ENMs in invasive algae research and to predict effects of climate change on seaweed distributions.