Gut microbiome biogeography in reindeer supersedes millennia of ecological and evolutionary separation.

Ruminants are dependent on their gut microbiomes for nutrient extraction from plant diets. However, knowledge about the composition, diversity, function, and spatial structure of gut microbiomes, especially in wild ruminants, is limited, largely because analysis has been restricted to faeces or the rumen. In two geographically separated reindeer subspecies, 16S rRNA gene amplicon sequencing revealed strong spatial structuring, and pronounced differences in microbial diversity of at least 33 phyla across the stomach, small intestine, and large intestine (including faeces). The main structural feature was the Bacteroidota to Firmicutes ratio, which declined from the stomach to the large intestine, likely reflecting functional adaptation. Metagenome shotgun sequencing also revealed highly significant structuring in the relative occurrence of carbohydrate-active enzymes (CAZymes). CAZymes were enriched in the rumen relative to the small and large intestines. Interestingly, taxonomic diversity was highest in the large intestine, suggesting an important and understudied role for this organ. Despite the two study populations being separated by an ocean and six millennia of evolutionary history, gut microbiome structuring was remarkably consistent. Our study suggests a strong selection for gut microbiome biogeography along the gastrointestinal tract in reindeer subspecies.

Changes in the fecal microbiota in dogs with acute hemorrhagic diarrhea during an outbreak in Norway.

BACKGROUND: A severe form of acute hemorrhagic diarrhea syndrome (AHDS) occurred in dogs in the Oslo region of Norway during autumn 2019. OBJECTIVES: To characterize the fecal microbiota of dogs with AHDS during the outbreak and compare it to that of healthy dogs from the same period and before the outbreak. ANIMALS: Dogs with AHDS (n = 50), dogs with nonhemorrhagic diarrhea (n = 3), and healthy dogs (n = 11) were sampled during the outbreak. In addition, 78 healthy dogs from the same region were sampled before the outbreak between 2017 and 2018. METHODS: Retrospective case-control study. The fecal microbiotas were characterized using 16S rRNA gene amplicon sequencing. RESULTS: Dogs with AHDS had significantly different microbiota composition (R2  = .07, P < .001) and decreased intestinal diversity relative to healthy dogs from the outbreak period (median, 2.7; range, 0.9-3.5 vs median, 3.2; range, 2.6-4.0; P < .001). The microbiota in dogs with AHDS was characterized by a decrease of Firmicutes and an outgrowth of Proteobacteria, with increased numbers of Clostridium perfringens and Providencia spp. Among the Providencia spp., 1 showed 100% sequence identity with a Providencia alcalifaciens strain that was cultivated and isolated from the same outbreak. No Providencia spp. was found in healthy dogs sampled before the outbreak. CONCLUSIONS AND CLINICAL IMPORTANCE: Dogs with AHDS had marked changes in fecal microbiota including increased numbers of Providencia spp. and C. perfringens, which may have contributed to the severity of this illness.

The Reproductive Microbiome: An Emerging Driver of Sexual Selection, Sexual Conflict, Mating Systems, and Reproductive Isolation.

All multicellular organisms host microbial communities in and on their bodies, and these microbiomes can have major influences on host biology. Most research has focussed on the oral, skin, and gut microbiomes, whereas relatively little is known about the reproductive microbiome. Here, we review empirical evidence to show that reproductive microbiomes can have significant effects on the reproductive function and performance of males and females. We then discuss the likely repercussions of these effects for evolutionary processes related to sexual selection and sexual conflict, as well as mating systems and reproductive isolation. We argue that knowledge of the reproductive microbiome is fundamental to our understanding of the evolutionary ecology of reproductive strategies and sexual dynamics of host organisms.

Individuality and convergence of the infant gut microbiota during the first year of life.

The human gut microbiota plays a vital role in health and disease, and microbial colonization is a key process in infant development. Here, we analyze 2684 fecal specimens from 12 infants during their first year of life, providing detailed insights into the human gut colonization process. Maturation of the gut microbial community shows strong temporal structure and specific developmental stages. At 2-4 months of age, there is a period of accelerated convergence concurrent with a bloom of Bifidobacterium, a genus associated with metabolism of oligosaccharides found in breast milk. The end of this period coincides with the introduction of solid food, a reduction in the relative abundance of Bifidobacterium, and an increase in several groups of Firmicutes. Our findings highlight the dynamic nature and individuality of the gut colonization process, and the need for high-frequency sampling over an extended period when designing and interpreting infant microbiome studies.

Multilevel social structure and diet shape the gut microbiota of the gelada monkey, the only grazing primate.

BACKGROUND: The gelada monkey (Theropithecus gelada), endemic to the Ethiopian highlands, is the only graminivorous primate, i.e., it feeds mainly on grasses and sedges. In spite of known dental, manual, and locomotor adaptations, the intestinal anatomy of geladas is similar to that of other primates. We currently lack a clear understanding of the adaptations in digestive physiology necessary for this species to subsist on a graminoid-based diet, but digestion in other graminivores, such as ruminants, relies heavily on the microbial community residing in the gastrointestinal (GI) system. Furthermore, geladas form complex, multilevel societies, making them a suitable system for investigating links between sociality and the GI microbiota. RESULTS: Here, we explore the gastrointestinal microbiota of gelada monkeys inhabiting an intact ecosystem and document how factors like multilevel social structure and seasonal changes in diet shape the GI microbiota. We compare the gelada GI microbiota to those of other primate species, reporting a gradient from geladas to herbivorous specialist monkeys to dietary generalist monkeys and lastly humans, the ultimate ecological generalists. We also compare the microbiotas of the gelada GI tract and the sheep rumen, finding that geladas are highly enriched for cellulolytic bacteria associated with ruminant digestion, relative to other primates. CONCLUSIONS: This study represents the first analysis of the gelada GI microbiota, providing insights into the adaptations underlying graminivory in a primate. Our results also highlight the role of social organization in structuring the GI microbiota within a society of wild animals.

Ecological plasticity in the gastrointestinal microbiomes of Ethiopian Chlorocebus monkeys.

Human activities can cause habitat degradation that may alter the types and quality of available food resources and thus influence the microbiomes of wild animal populations. Furthermore, seasonal shifts in food availability may cause adaptive responses in the gut microbiome to meet the need for different metabolic capabilities. Here, we demonstrate local-scale population structure in the gastrointestinal microbiotas of Chlorocebus monkeys, in southern Ethiopia, in response to varying degrees of human encroachment. We further provide evidence of adaptation to ecological conditions associated with the dry and wet seasons, and show seasonal effects to be more pronounced in areas with limited human activity. Finally, we report species-level microbiota differences between the endemic Ethiopian Bale monkey, an ecological specialist, and generalist Chlorocebus species from the same geographical region.

Linking Spatial Structure and Community-Level Biotic Interactions through Cooccurrence and Time Series Modeling of the Human Intestinal Microbiota.

The gastrointestinal (GI) microbiome is a densely populated ecosystem where dynamics are determined by interactions between microbial community members, as well as host factors. The spatial organization of this system is thought to be important in human health, yet this aspect of our resident microbiome is still poorly understood. In this study, we report significant spatial structure of the GI microbiota, and we identify general categories of spatial patterning in the distribution of microbial taxa along a healthy human GI tract. We further estimate the biotic interaction structure in the GI microbiota, both through time series and cooccurrence modeling of microbial community data derived from a large number of sequentially collected fecal samples. Comparison of these two approaches showed that species pairs involved in significant negative interactions had strong positive contemporaneous correlations and vice versa, while for species pairs without significant interactions, contemporaneous correlations were distributed around zero. We observed similar patterns when comparing these models to the spatial correlations between taxa identified in the adherent microbiota. This suggests that colocalization of microbial taxon pairs, and thus the spatial organization of the GI microbiota, is driven, at least in part, by direct or indirect biotic interactions. Thus, our study can provide a basis for an ecological interpretation of the biogeography of the human gut. IMPORTANCE The human gut microbiome is the subject of intense study due to its importance in health and disease. The majority of these studies have been based on the analysis of feces. However, little is known about how the microbial composition in fecal samples relates to the spatial distribution of microbial taxa along the gastrointestinal tract. By characterizing the microbial content both in intestinal tissue samples and in fecal samples obtained daily, we provide a conceptual framework for how the spatial structure relates to biotic interactions on the community level. We further describe general categories of spatial distribution patterns and identify taxa conforming to these categories. To our knowledge, this is the first study combining spatial and temporal analyses of the human gut microbiome. This type of analysis can be used for identifying candidate probiotics and designing strategies for clinical intervention.

A novel ultra high-throughput 16S rRNA gene amplicon sequencing library preparation method for the Illumina HiSeq platform.

BACKGROUND: Advances in sequencing technologies and bioinformatics have made the analysis of microbial communities almost routine. Nonetheless, the need remains to improve on the techniques used for gathering such data, including increasing throughput while lowering cost and benchmarking the techniques so that potential sources of bias can be better characterized. METHODS: We present a triple-index amplicon sequencing strategy to sequence large numbers of samples at significantly lower c ost and in a shorter timeframe compared to existing methods. The design employs a two-stage PCR protocol, incorpo rating three barcodes to each sample, with the possibility to add a fourth-index. It also includes heterogeneity spacers to overcome low complexity issues faced when sequencing amplicons on Illumina platforms. RESULTS: The library preparation method was extensively benchmarked through analysis of a mock community in order to assess biases introduced by sample indexing, number of PCR cycles, and template concentration. We further evaluated the method through re-sequencing of a standardized environmental sample. Finally, we evaluated our protocol on a set of fecal samples from a small cohort of healthy adults, demonstrating good performance in a realistic experimental setting. Between-sample variation was mainly related to batch effects, such as DNA extraction, while sample indexing was also a significant source of bias. PCR cycle number strongly influenced chimera formation and affected relative abundance estimates of species with high GC content. Libraries were sequenced using the Illumina HiSeq and MiSeq platforms to demonstrate that this protocol is highly scalable to sequence thousands of samples at a very low cost. CONCLUSIONS: Here, we provide the most comprehensive study of performance and bias inherent to a 16S rRNA gene amplicon sequencing method to date. Triple-indexing greatly reduces the number of long custom DNA oligos required for library preparation, while the inclusion of variable length heterogeneity spacers minimizes the need for PhiX spike-in. This design results in a significant cost reduction of highly multiplexed amplicon sequencing. The biases we characterize highlight the need for highly standardized protocols. Reassuringly, we find that the biological signal is a far stronger structuring factor than the various sources of bias.

Microsatellite Length Scoring by Single Molecule Real Time Sequencing - Effects of Sequence Structure and PCR Regime.

Microsatellites are DNA sequences consisting of repeated, short (1-6 bp) sequence motifs that are highly mutable by enzymatic slippage during replication. Due to their high intrinsic variability, microsatellites have important applications in population genetics, forensics, genome mapping, as well as cancer diagnostics and prognosis. The current analytical standard for microsatellites is based on length scoring by high precision electrophoresis, but due to increasing efficiency next-generation sequencing techniques may provide a viable alternative. Here, we evaluated single molecule real time (SMRT) sequencing, implemented in the PacBio series of sequencing apparatuses, as a means of microsatellite length scoring. To this end we carried out multiplexed SMRT sequencing of plasmid-carried artificial microsatellites of varying structure under different pre-sequencing PCR regimes. For each repeat structure, reads corresponding to the target length dominated. We found that pre-sequencing amplification had large effects on scoring accuracy and error distribution relative to controls, but that the effects of the number of amplification cycles were generally weak. In line with expectations enzymatic slippage decreased proportionally with microsatellite repeat unit length and increased with repetition number. Finally, we determined directional mutation trends, showing that PCR and SMRT sequencing introduced consistent but opposing error patterns in contraction and expansion of the microsatellites on the repeat motif and single nucleotide level.

Ecology of bacteria in the human gastrointestinal tract--identification of keystone and foundation taxa.

BACKGROUND: Determining ecological roles of community members and the impact of specific taxa on overall biodiversity in the gastrointestinal (GI) microbiota is of fundamental importance. A step towards a systems-level understanding of the GI microbiota is characterization of biotic interactions. Community time series analysis, an approach based on statistical analysis of changing population abundances within a single system over time, is needed in order to say with confidence that one population is affecting the dynamics of another. RESULTS: Here, we characterize biotic interaction structures and define ecological roles of major bacterial groups in four healthy individuals by analysing high-resolution, long-term (>180 days) GI bacterial community time series. Actinobacteria fit the description of a keystone taxon since they are relatively rare, but have a high degree of ecological connectedness, and are positively correlated with diversity both within and between individuals. Bacteriodetes were found to be a foundation taxon in that they are numerically dominant and interact extensively, in particular through positive interactions, with other taxa. Although community structure, diversity and biotic interaction patterns were specific to each individual, we observed a strong tendency towards more intense competition within than between phyla. This is in agreement with Darwin's limiting similarity hypothesis as well as a published biotic interaction model of the GI microbiota based on reverse ecology. Finally, we link temporal enterotype switching to a reciprocal positive interaction between two key genera. CONCLUSIONS: In this study, we identified ecological roles of key taxa in the human GI microbiota and compared our time series analysis results with those obtained through a reverse ecology approach, providing further evidence in favour of the limiting similarity hypothesis first put forth by Darwin. Larger longitudinal studies are warranted in order to evaluate the generality of basic ecological concepts as applied to the GI microbiota, but our results provide a starting point for achieving a more profound understanding of the GI microbiota as an ecological system.

Biotic interactions and temporal dynamics of the human gastrointestinal microbiota.

The human gastrointestinal (GI) microbiota is important to human health and imbalances or shifts in the gut microbial community have been linked to many diseases. Most studies of the GI microbiota only capture snapshots of this dynamic community at one or a few time points. Although this is valuable in terms of providing knowledge of community composition and variability between individuals, it does not provide the foundation for going beyond descriptive studies and toward truly predictive ecological models. In order to achieve this goal, we need longitudinal data of appropriate temporal and taxonomic resolution, so that established time series analysis tools for identifying and quantifying putative interactions among community members can be used. Here, we present new analyses of existing data to illustrate the potential usefulness of this approach. We discuss challenges related to sampling and data processing, as well as analytical approaches and considerations for future studies of the GI microbiota and other complex microbial systems.

Large-scale genetic structuring of a widely distributed carnivore--the Eurasian lynx (Lynx lynx).

Over the last decades the phylogeography and genetic structure of a multitude of species inhabiting Europe and North America have been described. The flora and fauna of the vast landmasses of north-eastern Eurasia are still largely unexplored in this respect. The Eurasian lynx is a large felid that is relatively abundant over much of the Russian sub-continent and the adjoining countries. Analyzing 148 museum specimens collected throughout its range over the last 150 years we have described the large-scale genetic structuring in this highly mobile species. We have investigated the spatial genetic patterns using mitochondrial DNA sequences (D-loop and cytochrome b) and 11 microsatellite loci, and describe three phylogenetic clades and a clear structuring along an east-west gradient. The most likely scenario is that the contemporary Eurasian lynx populations originated in central Asia and that parts of Europe were inhabited by lynx during the Pleistocene. After the Last Glacial Maximum (LGM) range expansions lead to colonization of north-western Siberia and Scandinavia from the Caucasus and north-eastern Siberia from a refugium further east. No evidence of a Berinigan refugium could be detected in our data. We observed restricted gene flow and suggest that future studies of the Eurasian lynx explore to what extent the contemporary population structure may be explained by ecological variables.

Context-dependent competition in a model gut bacterial community.

Understanding the ecological processes that generate complex community structures may provide insight into the establishment and maintenance of a normal microbial community in the human gastrointestinal tract, yet very little is known about how biotic interactions influence community dynamics in this system. Here, we use natural strains of Escherichia coli and a simplified model microbiota to demonstrate that the colonization process on the strain level can be context dependent, in the sense that the outcome of intra-specific competition may be determined by the composition of the background community. These results are consistent with previous models for competition between organisms where one competitor has adapted to low resource environments whereas the other is optimized for rapid reproduction when resources are abundant. The genomic profiles of E. coli strains representing these differing ecological strategies provide clues for deciphering the genetic underpinnings of niche adaptation within a single species. Our findings extend the role of ecological theory in understanding microbial systems and the conceptual toolbox for describing microbial community dynamics. There are few, if any, concrete examples of context-dependent competition on a single trophic level. However, this phenomenon can have potentially dramatic effects on which bacteria will successfully establish and persist in the gastrointestinal system, and the principle should be equally applicable to other microbial ecosystems.

Comparisons of infant Escherichia coli isolates link genomic profiles with adaptation to the ecological niche.

BACKGROUND: Despite being one of the most intensely studied model organisms, many questions still remain about the evolutionary biology and ecology of Escherichia coli. An important step toward achieving a more complete understanding of E.coli biology entails elucidating relationships between gene content and adaptation to the ecological niche. RESULTS: Here, we present genome comparisons of 16 E.coli strains that represent commensals and pathogens isolated from infants during a specific time period in Trondheim, Norway. Using differential gene content, we characterized enrichment profiles of the collection of strains relating to phylogeny, early vs. late colonization, pathogenicity and growth rate. We found clear gene content distinctions relating to the various grouping criteria. We also found that different categories of strains use different genetic elements for similar biological processes. The sequenced genomes included two pairs of strains where each pair was isolated from the same infant at different time points. One pair, in which the strains were isolated four months apart, showed maintenance of an early colonizer genome profile but also gene content and codon usage changes toward the late colonizer profile. Lastly, we placed our sequenced isolates into a broader genomic context by comparing them with 25 published E.coli genomes that represent a variety of pathotypes and commensal strains. This analysis demonstrated the importance of geography in shaping strain level gene content profiles. CONCLUSIONS: Our results indicate a general pattern where alternative genetic pathways lead toward a consistent ecological role for E.coli as a species. Within this framework however, we saw selection shaping the coding repertoire of E.coli strains toward distinct ecotypes with different phenotypic properties.

Web of ecological interactions in an experimental gut microbiota.

The dynamics of all ecosystems are dictated by intrinsic, density-dependent mechanisms and by density-independent environmental forcing. In spite of the importance of the gastrointestinal microbiota in health and disease, the ecology of this system remains largely unknown. Here, we take an ecological approach to gut microbial community analysis, with statistical modelling of time series data from chemostats. This approach removes effects of host forcing, allowing us to describe a network of intrinsic interactions determining the dynamic structure of an experimental gut microbiota. Surprisingly, the main colonization pattern in this simplified model system resembled that of the human infant gut, suggesting a potentially important role of density-dependent interactions in the early gut microbiota. Knowledge of ecological structures in microbial systems may provide us with a means of controlling such systems by modifying the strength and nature of interactions among microbes and between the microbes and their environment.

Convergent temporal dynamics of the human infant gut microbiota.

Temporal dynamics of the human gut microbiota is of fundamental importance for the development of proper gut function and maturation of the immune system. Here we present a description of infant gut ecological dynamics using a combination of nonlinear data modeling and simulations of the early infant gut colonization processes. Principal component analysis of infant microbiota 16S rRNA gene microarray data showed that the main directions of variation were defined by three phylum-specific probes targeting Bacteroides, Proteobacteria and Firmicutes. Nonlinear regression analysis identified several dynamic interactions between these three phyla. Simulations of the early phylum-level colonization process showed the relatively rapid establishment of an equilibrium community after an unstable initial phase. In general, varying the initial composition of phyla in the simulations had little bearing on the final equilibrium. The dynamic interaction model was found to maintain its predictive ability for Proteobacteria and Firmicutes well into the simulation, whereas Bacteroides densities tended to be underestimated, possibly due to host top-down selection for Bacteroides. In accordance with our model, initial perturbation of the microbiota by different mode of delivery (vaginal and C-section) did not affect the later phylum composition in the infants investigated. Considering the predictive ability and convergence of our phylum-level model, we now propose that deterministic bacterial-bacterial interactions are more important for shaping the human infant gut microbiota than previously anticipated.

Characterizing mixed microbial population dynamics using time-series analysis.

Due to a general shortage of temporal population data, dynamic structures in microbial communities remain largely unexplored. Knowledge of community dynamics is, however, essential for understanding the mechanisms by which microbes interact. Here, we have used a computational approach for quantification of bacteria in multispecies populations, generating data for time-series modeling. Moreover, we have used online FR-IR spectroscopy to monitor the main metabolic processes. The approach enabled us to provide a functional description of the parameters governing the population dynamics in a three-species model bacterial community, demonstrating density-dependent regulation, interspecies competition and even a case of cooperation between two species. Since the field of microbial ecology has yet to embrace many of the concepts and methods developed for the study of ecology of higher plants and animals, the realization that microbial systems can be analyzed within the same conceptual framework as other ecosystems is of fundamental importance.

Co-infection dynamics of a major food-borne zoonotic pathogen in chicken.

A major bottleneck in understanding zoonotic pathogens has been the analysis of pathogen co-infection dynamics. We have addressed this challenge using a novel direct sequencing approach for pathogen quantification in mixed infections. The major zoonotic food-borne pathogen Campylobacter jejuni, with an important reservoir in the gastrointestinal (GI) tract of chickens, was used as a model. We investigated the co-colonisation dynamics of seven C. jejuni strains in a chicken GI infection trial. The seven strains were isolated from an epidemiological study showing multiple strain infections at the farm level. We analysed time-series data, following the Campylobacter colonisation, as well as the dominant background flora of chickens. Data were collected from the infection at day 16 until the last sampling point at day 36. Chickens with two different background floras were studied, mature (treated with Broilact, which is a product consisting of bacteria from the intestinal flora of healthy hens) and spontaneous. The two treatments resulted in completely different background floras, yet similar Campylobacter colonisation patterns were detected in both groups. This suggests that it is the chicken host and not the background flora that is important in determining the Campylobacter colonisation pattern. Our results showed that mainly two of the seven C. jejuni strains dominated the Campylobacter flora in the chickens, with a shift of the dominating strain during the infection period. We propose a model in which multiple C. jejuni strains can colonise a single host, with the dominant strains being replaced as a consequence of strain-specific immune responses. This model represents a new understanding of C. jejuni epidemiology, with future implications for the development of novel intervention strategies.