Contribution of neutral processes to the assembly of gut microbial communities in the zebrafish over host development

Despite their importance to host health and development, the communities of microorganisms associated with humans and other animals are characterized by a large degree of unexplained variation across individual hosts. The processes that drive such inter-individual variation are not well understood. To address this, we surveyed the microbial communities associated with the intestine of the zebrafish, Danio rerio, over developmental time. We compared our observations of community composition and distribution across hosts with that predicted by a neutral assembly model, which assumes that community assembly is driven solely by chance and dispersal. We found that as hosts develop from larvae to adults, the fit of the model to observed microbial distributions decreases, suggesting that the relative importance of non-neutral processes, such as microbe-microbe interactions, active dispersal, or selection by the host, increases as hosts mature. We also observed that taxa which depart in their distributions from the neutral prediction form ecologically distinct sub-groups, which are phylogenetically clustered with respect to the full metacommunity. These results demonstrate that neutral processes are sufficient to generate substantial variation in microbiota composition across individual hosts, and suggest that potentially unique or important taxa may be identified by their divergence from neutral distributions.     

Chronic Trichuris muris Infection Decreases Diversity of the Intestinal Microbiota and Concomitantly Increases the Abundance of Lactobacilli

The intestinal microbiota is vital for shaping the local intestinal environment as well as host immunity and metabolism. At the same time, epidemiological and experimental evidence suggest an important role for parasitic worm infections in maintaining the inflammatory and regulatory balance of the immune system. In line with this, the prevalence of persistent worm infections is inversely correlated with the incidence of immune-associated diseases, prompting the use of controlled parasite infections for therapeutic purposes. Despite this, the impact of parasite infection on the intestinal microbiota, as well as potential downstream effects on the immune system, remain largely unknown. We have assessed the influence of chronic infection with the large-intestinal nematode Trichuris muris, a close relative of the human pathogen Trichuris trichiura, on the composition of the murine intestinal microbiota by 16S ribosomal-RNA gene-based sequencing. Our results demonstrate that persistent T. muris infection dramatically affects the large-intestinal microbiota, most notably with a drop in the diversity of bacterial communities, as well as a marked increase in the relative abundance of the Lactobacillus genus. In parallel, chronic T. muris infection resulted in a significant shift in the balance between regulatory and inflammatory T cells in the intestinal adaptive immune system, in favour of inflammatory cells. Together, these data demonstrate that chronic parasite infection strongly influences the intestinal microbiota and the adaptive immune system. Our results illustrate the complex interactions between these factors in the intestinal tract, and contribute to furthering the understanding of this interplay, which is of crucial importance considering that 500 million people globally are suffering from these infections and their potential use for therapeutic purposes.     

Amphibian skin may select for rare environmental microbes

Host-microbe symbioses rely on the successful transmission or acquisition of symbionts in each new generation. Amphibians host a diverse cutaneous microbiota, and many of these symbionts appear to be mutualistic and may limit infection by the chytrid fungus, Batrachochytrium dendrobatidis, which has caused global amphibian population declines and extinctions in recent decades. Using bar-coded 454 pyrosequencing of the 16S rRNA gene, we addressed the question of symbiont transmission by examining variation in amphibian skin microbiota across species and sites and in direct relation to environmental microbes. Although acquisition of environmental microbes occurs in some host-symbiont systems, this has not been extensively examined in free-living vertebrate-microbe symbioses. Juvenile bullfrogs (Rana catesbeiana), adult red-spotted newts (Notophthalmus viridescens), pond water and pond substrate were sampled at a single pond to examine host-specificity and potential environmental transmission of microbiota. To assess population level variation in skin microbiota, adult newts from two additional sites were also sampled. Cohabiting bullfrogs and newts had distinct microbial communities, as did newts across the three sites. The microbial communities of amphibians and the environment were distinct; there was very little overlap in the amphibians'' core microbes and the most abundant environmental microbes, and the relative abundances of OTUs that were shared by amphibians and the environment were inversely related. These results suggest that, in a host species-specific manner, amphibian skin may select for microbes that are generally in low abundance in the environment.     

The microbiota of intertidal macroalgae Fucus distichus is site-specific and resistant to change following transplant

It is unclear how host-associated microbial communities will be affected by future environmental change. Characterizing how microbiota differ across sites with varying environmental conditions and assessing the stability of the microbiota in response to abiotic variation are critical steps towards predicting outcomes of environmental change. Intertidal organisms are valuable study systems because they experience extreme variation in environmental conditions on tractable timescales such as tide cycles and across small spatial gradients in the intertidal zone. Here we show a widespread intertidal macroalgae, Fucus distichus, hosts site-specific microbiota over small (meters to kilometres) spatial scales. We demonstrate stability of site-specific microbial associations by manipulating the host environment and microbial species pool with common garden and reciprocal transplant experiments. We hypothesized that F. distichus microbiota would readily shift to reflect the contemporary environment due to selective filtering by abiotic conditions and/or colonization by microbes from the new environment or nearby hosts. Instead, F. distichus microbiota was stable for days after transplantation in both the laboratory and field. Our findings expand the current understanding of microbiota dynamics on an intertidal foundation species. These results may also point to adaptations for withstanding short-term environmental variation, in hosts and/or microbes, facilitating stable host–microbial associations.     

Biphasic assembly of the murine intestinal microbiota during early development

The birth canal provides mammals with a primary maternal inoculum, which develops into distinctive body site-specific microbial communities post-natally. We characterized the distal gut microbiota from birth to weaning in mice. One-day-old mice had colonic microbiota that resembled maternal vaginal communities, but at days 3 and 9 of age there was a substantial loss of intestinal bacterial diversity and dominance of Lactobacillus. By weaning (21 days), diverse intestinal bacteria had established, including strict anaerobes. Our results are consistent with vertical transmission of maternal microbiota and demonstrate a nonlinear ecological succession involving an early drop in bacterial diversity and shift in dominance from Streptococcus to Lactobacillus, followed by an increase in diversity of anaerobes, after the introduction of solid food. Mammalian newborns are born highly susceptible to colonization, and lactation may control microbiome assembly during early development.     

The Colonization Dynamics of the Gut Microbiota in Tilapia Larvae

The gut microbiota of fish larvae evolves fast towards a complex community. Both host and environment affect the development of the gut microbiota; however, the relative importance of both is poorly understood. Determining specific changes in gut microbial populations in response to a change in an environmental factor is very complicated. Interactions between factors are difficult to separate and any response could be masked due to high inter-individual variation even for individuals that share a common environment. In this study we characterized and quantified the spatio-temporal variation in the gut microbiota of tilapia larvae, reared in recirculating aquaculture systems (RAS) or active suspension tanks (AS). Our results showed that variation in gut microbiota between replicate tanks was not significantly higher than within tank variation, suggesting that there is no tank effect on water and gut microbiota. However, when individuals were reared in replicate RAS, gut microbiota differed significantly. The highest variation was observed between individuals reared in different types of system (RAS vs. AS). Our data suggest that under experimental conditions in which the roles of deterministic and stochastic factors have not been precisely determined, compositional replication of the microbial communities of an ecosystem is not predictable.     

Evidence for a core gut microbiota in the zebrafish

Experimental analysis of gut microbial communities and their interactions with vertebrate hosts is conducted predominantly in domesticated animals that have been maintained in laboratory facilities for many generations. These animal models are useful for studying coevolved relationships between host and microbiota only if the microbial communities that occur in animals in lab facilities are representative of those that occur in nature. We performed 16S rRNA gene sequence-based comparisons of gut bacterial communities in zebrafish collected recently from their natural habitat and those reared for generations in lab facilities in different geographic locations. Patterns of gut microbiota structure in domesticated zebrafish varied across different lab facilities in correlation with historical connections between those facilities. However, gut microbiota membership in domesticated and recently caught zebrafish was strikingly similar, with a shared core gut microbiota. The zebrafish intestinal habitat therefore selects for specific bacterial taxa despite radical differences in host provenance and domestication status.     

Faecal Microbiota Composition in Adults Is Associated with the FUT2 Gene Determining the Secretor Status

The human intestine is colonised with highly diverse and individually defined microbiota, which likely has an impact on the host well-being. Drivers of the individual variation in the microbiota compositions are multifactorial and include environmental, host and dietary factors. We studied the impact of the host secretor status, encoded by fucosyltransferase 2 (FUT2) -gene, on the intestinal microbiota composition. Secretor status determines the expression of the ABH and Lewis histo-blood group antigens in the intestinal mucosa. The study population was comprised of 14 non-secretor (FUT2 rs601338 genotype AA) and 57 secretor (genotypes GG and AG) adult individuals of western European descent. Intestinal microbiota was analyzed by PCR-DGGE and for a subset of 12 non-secretor subjects and 12 secretor subjects additionally by the 16S rRNA gene pyrosequencing and the HITChip phylogenetic microarray analysis. All three methods showed distinct clustering of the intestinal microbiota and significant differences in abundances of several taxa representing dominant microbiota between the non-secretors and the secretors as well as between the FUT2 genotypes. In addition, the non-secretors had lower species richness than the secretors. The soft clustering of microbiota into enterotypes (ET) 1 and 3 showed that the non-secretors had a higher probability of belonging to ET1 and the secretors to ET3. Our study shows that secretor status and FUT2 polymorphism are associated with the composition of human intestinal microbiota, and appears thus to be one of the key drivers affecting the individual variation of human intestinal microbiota.     

Variation in the intestinal microbiota at different developmental stages of Hynobius maoershanensis

Intestinal microbiota play an important role in the life of amphibians and its composition may vary by developmental stage. In this study, 16S rRNA high‐throughput sequencing was used to profile the intestinal microbiota of Hynobius maoershanensis, which exclusively inhabit the Maoer Mountain swamp at an altitude of approximately 2,000 m. We characterized the bacterial composition, structure, and function of the microbiota of H. maoershanensis at different developmental stages. The alpha diversity was not markedly different for the Simpson, Shannon, Ace, and Sobs indices of microbes. The beta diversity revealed that there were age‐related differences in the structure of the intestinal microbes of H. maoershanensis, specifically, at the phylum level. Bacteroidetes and Proteobacteria were the dominant bacteria present in the adult stage, and the relative abundance of Bacteroidetes was significantly higher compared with that of tadpoles. Firmicutes and Proteobacteria were the dominant phylum during the tadpole stage and their relative abundance was significantly higher compared with the adult period. Functional analysis revealed that the pathways associated with organismal systems and metabolism were significantly enriched in the adults, whereas human diseases, genetic information processing, and cellular processes were more enriched in the hindlimb bud stage. Human diseases and environmental information processing were more enriched in the forelimb bud stage at KEGG pathway level 1. Possibilities for the observed discrepancies include the adaptation to eating habits and the remodeling of the intestines during development. We speculated that H. maoershanensis adults may be more suitable to a high‐fiber diet, whereas the tadpoles are associated with a carnivorous diet. Our study provides evidence of variations in the intestinal microbiota during development in amphibians, highlighting the influence of historical developments on the intestinal microbiota and an increased understanding of the importance of physiological characteristics in shaping the intestinal microbiota of amphibians. These data will help us formulate more effective protection measures for H. maoershanensis.     

Bacterial community assembly and turnover within the intestines of developing zebrafish

Background

The majority of animal associated microorganisms are present in digestive tract communities. These intestinal communities arise from selective pressures of the gut habitats as well as host''s genotype are regarded as an extra ‘organ’ regulate functions that have not evolved wholly on the host. They are functionally essential in providing nourishment, regulating epithelial development, and influencing immunity in the vertebrate host. As vertebrates are born free of microorganisms, what is poorly understood is how intestinal bacterial communities assemble and develop in conjunction with the development of the host.

Methodology/Principal Findings

Set within an ecological framework, we investigated the bacterial community assembly and turnover within the intestinal habitats of developing zebrafish (from larvae to adult animals). Spatial and temporal species-richness relationships and Mantel and partial Mantel tests revealed that turnover was low and that richness and composition was best predicted by time and not intestinal volume (habitat size) or changes in food diet. We also observed that bacterial communities within the zebrafish intestines were deterministically assembled (reflected by the observed low turnover) switching to stochastic assembly in the later stages of zebrafish development.

Conclusions/Significance

This study is of importance as it provides a novel insight into how intestinal bacterial communities assemble in tandem with the host''s development (from early to adult stages). It is our hope that by studying intestinal microbiota of this vertebrate model with such or some more refined approaches in the future could well provide ecological insights for clinical benefit. In addition, this study also adds to our still fledgling knowledge of how spatial and temporal species-richness relationships are shaped and provides further mounting evidence that bacterial community assembly and dynamics are shaped by both deterministic and stochastic considerations.     

Spatiotemporal successions of shrimp gut microbial colonization: high consistency despite distinct species pool

Aquatic animals encounter suites of novel planktonic microbes during their development. Although hosts have been shown to exert strong selection on their gut microbiota from surrounding environment, to what extent and the generality that the gut microbiota and the underlying ecological processes are affected by biotic and abiotic variations are largely unclear. Here, these concerns were explored by coupling spatiotemporal data on gut and rearing water bacterial communities with environmental variables over shrimp life stages at spatially distant locations. Shrimp gut microbiotas significantly changed mirroring their development, as evidenced by gut bacterial signatures of shrimp life stage contributing 95.5% stratification accuracy. Shrimp sourced little (2.6%–15.8%) of their gut microbiota from their rearing water. This microbial resistance was reflected by weak compositional differences between shrimp farming spatially distinct locations where species pools were distinct. Consistently, the assembly of shrimp gut microbiota was not adequately explained by the rearing water variables and bacterial community, but rather by host-age-associated biotic features. The successions of shrimp gut microbiota were droved by replacement (βsim), rather than by nestedness (βnes), while those of bacterioplankton communities were equally governed by replacement and nestedness. Our study highlights how shrimp gut bacterial community assembly is coupled to their development, rearing species pool, and that the successional pattern of host-associated communities is differed from that of free-living bacteria.     

Variation in Taxonomic Composition of the Fecal Microbiota in an Inbred Mouse Strain across Individuals and Time

Genetics, diet, and other environmental exposures are thought to be major factors in the development and composition of the intestinal microbiota of animals. However, the relative contributions of these factors in adult animals, as well as variation with time in a variety of important settings, are still not fully understood. We studied a population of inbred, female mice fed the same diet and housed under the same conditions. We collected fecal samples from 46 individual mice over two weeks, sampling four of these mice for periods as long as 236 days for a total of 190 samples, and determined the phylogenetic composition of their microbial communities after analyzing 1,849,990 high-quality pyrosequencing reads of the 16S rRNA gene V3 region. Even under these controlled conditions, we found significant inter-individual variation in community composition, as well as variation within an individual over time, including increases in alpha diversity during the first 2 months of co-habitation. Some variation was explained by mouse membership in different cage and vendor shipment groups. The differences among individual mice from the same shipment group and cage were still significant. Overall, we found that 23% of the variation in intestinal microbiota composition was explained by changes within the fecal microbiota of a mouse over time, 12% was explained by persistent differences among individual mice, 14% by cage, and 18% by shipment group. Our findings suggest that the microbiota of controlled populations of inbred laboratory animals may not be as uniform as previously thought, that animal rearing and handling may account for some variation, and that as yet unidentified factors may explain additional components of variation in the composition of the microbiota within populations and individuals over time. These findings have implications for the design and interpretation of experiments involving laboratory animals.     

Dietary input of microbes and host genetic variation shape among-population differences in stickleback gut microbiota

To explain differences in gut microbial communities we must determine how processes regulating microbial community assembly (colonization, persistence) differ among hosts and affect microbiota composition. We surveyed the gut microbiota of threespine stickleback (Gasterosteus aculeatus) from 10 geographically clustered populations and sequenced environmental samples to track potential colonizing microbes and quantify the effects of host environment and genotype. Gut microbiota composition and diversity varied among populations. These among-population differences were associated with multiple covarying ecological variables: habitat type (lake, stream, estuary), lake geomorphology and food- (but not water-) associated microbiota. Fish genotype also covaried with gut microbiota composition; more genetically divergent populations exhibited more divergent gut microbiota. Our results suggest that population level differences in stickleback gut microbiota may depend more on internal sorting processes (host genotype) than on colonization processes (transient environmental effects).     

Aquatic animals promote antibiotic resistance gene dissemination in water via conjugation: Role of different regions within the zebra fish intestinal tract,and impact on fish intestinal microbiota

The aqueous environment is one of many reservoirs of antibiotic resistance genes (ARGs). Fish, as important aquatic animals which possess ideal intestinal niches for bacteria to grow and multiply, may ingest antibiotic resistance bacteria from aqueous environment. The fish gut would be a suitable environment for conjugal gene transfer including those encoding antibiotic resistance. However, little is known in relation to the impact of ingested ARGs or antibiotic resistance bacteria (ARB) on gut microbiota. Here, we applied the cultivation method, qPCR, nuclear molecular genetic marker and 16S rDNA amplicon sequencing technologies to develop a plasmid‐mediated ARG transfer model of zebrafish. Furthermore, we aimed to investigate the dissemination of ARGs in microbial communities of zebrafish guts after donors carrying self‐transferring plasmids that encode ARGs were introduced in aquaria. On average, 15% of faecal bacteria obtained ARGs through RP4‐mediated conjugal transfer. The hindgut was the most important intestinal region supporting ARG dissemination, with concentrations of donor and transconjugant cells almost 25 times higher than those of other intestinal segments. Furthermore, in the hindgut where conjugal transfer occurred most actively, there was remarkable upregulation of the mRNA expression of the RP4 plasmid regulatory genes, trbBp and trfAp. Exogenous bacteria seem to alter bacterial communities by increasing Escherichia and Bacteroides species, while decreasing Aeromonas compared with control groups. We identified the composition of transconjugants and abundance of both cultivable and uncultivable bacteria (the latter accounted for 90.4%–97.2% of total transconjugants). Our study suggests that aquatic animal guts contribute to the spread of ARGs in water environments.     

Modelling spatial patterns in host-associated microbial communities

Microbial communities exhibit spatial structure at different scales, due to constant interactions with their environment and dispersal limitation. While this spatial structure is often considered in studies focusing on free-living environmental communities, it has received less attention in the context of host-associated microbial communities or microbiota. The wider adoption of methods accounting for spatial variation in these communities will help to address open questions in basic microbial ecology as well as realize the full potential of microbiome-aided medicine. Here, we first overview known factors affecting the composition of microbiota across diverse host types and at different scales, with a focus on the human gut as one of the most actively studied microbiota. We outline a number of topical open questions in the field related to spatial variation and patterns. We then review the existing methodology for the spatial modelling of microbiota. We suggest that methodology from related fields, such as systems biology and macro-organismal ecology, could be adapted to obtain more accurate models of spatial structure. We further posit that methodological developments in the spatial modelling and analysis of microbiota could in turn broadly benefit theoretical and applied ecology and contribute to the development of novel industrial and clinical applications.     

Mining zebrafish microbiota reveals key community-level resistance against fish pathogen infection

The long-known resistance to pathogens provided by host-associated microbiota fostered the notion that adding protective bacteria could prevent or attenuate infection. However, the identification of endogenous or exogenous bacteria conferring such protection is often hindered by the complexity of host microbial communities. Here, we used zebrafish and the fish pathogen Flavobacterium columnare as a model system to study the determinants of microbiota-associated colonization resistance. We compared infection susceptibility in germ-free, conventional and reconventionalized larvae and showed that a consortium of 10 culturable bacterial species are sufficient to protect zebrafish. Whereas survival to F. columnare infection does not rely on host innate immunity, we used antibiotic dysbiosis to alter zebrafish microbiota composition, leading to the identification of two different protection strategies. We first identified that the bacterium Chryseobacterium massiliae individually protects both larvae and adult zebrafish. We also showed that an assembly of 9 endogenous zebrafish species that do not otherwise protect individually confer a community-level resistance to infection. Our study therefore provides a rational approach to identify key endogenous protecting bacteria and promising candidates to engineer resilient microbial communities. It also shows how direct experimental analysis of colonization resistance in low-complexity in vivo models can reveal unsuspected ecological strategies at play in microbiota-based protection against pathogens.Subject terms: Microbiome, Microbial ecology     

Neonatal environment exerts a sustained influence on the development of the intestinal microbiota and metabolic phenotype

The postnatal environment, including factors such as weaning and acquisition of the gut microbiota, has been causally linked to the development of later immunological diseases such as allergy and autoimmunity, and has also been associated with a predisposition to metabolic disorders. We show that the very early-life environment influences the development of both the gut microbiota and host metabolic phenotype in a porcine model of human infants. Farm piglets were nursed by their mothers for 1 day, before removal to highly controlled, individual isolators where they received formula milk until weaning at 21 days. The experiment was repeated, to create two batches, which differed only in minor environmental fluctuations during the first day. At day 1 after birth, metabolic profiling of serum by ¹H nuclear magnetic resonance spectroscopy demonstrated significant, systemic, inter-batch variation which persisted until weaning. However, the urinary metabolic profiles demonstrated that significant inter-batch effects on 3-hydroxyisovalerate, trimethylamine-N-oxide and mannitol persisted beyond weaning to at least 35 days. Batch effects were linked to significant differences in the composition of colonic microbiota at 35 days, determined by 16 S pyrosequencing. Different weaning diets modulated both the microbiota and metabolic phenotype independently of the persistent batch effects. We demonstrate that the environment during the first day of life influences development of the microbiota and metabolic phenotype and thus should be taken into account when interrogating experimental outcomes. In addition, we suggest that intervention at this early time could provide ‘metabolic rescue'' for at-risk infants who have undergone aberrant patterns of initial intestinal colonisation.     

Inhibitory bacteria reduce fungi on early life stages of endangered Colorado boreal toads (Anaxyrus boreas)

Increasingly, host-associated microbiota are recognized to mediate pathogen establishment, providing new ecological perspectives on health and disease. Amphibian skin-associated microbiota interact with the fungal pathogen, Batrachochytrium dendrobatidis (Bd), but little is known about microbial turnover during host development and associations with host immune function. We surveyed skin microbiota of Colorado''s endangered boreal toads (Anaxyrus boreas), sampling 181 toads across four life stages (tadpoles, metamorphs, subadults and adults). Our goals were to (1) understand variation in microbial community structure among individuals and sites, (2) characterize shifts in communities during development and (3) examine the prevalence and abundance of known Bd-inhibitory bacteria. We used high-throughput 16S and 18S rRNA gene sequencing (Illumina MiSeq) to characterize bacteria and microeukaryotes, respectively. Life stage had the largest effect on the toad skin microbial community, and site and Bd presence also contributed. Proteobacteria dominated tadpole microbial communities, but were later replaced by Actinobacteria. Microeukaryotes on tadpoles were dominated by the classes Alveolata and Stramenopiles, while fungal groups replaced these groups after metamorphosis. Using a novel database of Bd-inhibitory bacteria, we found fewer Bd-inhibitory bacteria in post-metamorphic stages correlated with increased skin fungi, suggesting that bacteria have a strong role in early developmental stages and reduce skin-associated fungi.     

Ontogeny and nutritional control of adipogenesis in zebrafish (Danio rerio)

The global obesity epidemic demands an improved understanding of the developmental and environmental factors regulating fat storage. Adipocytes serve as major sites of fat storage and as regulators of energy balance and inflammation. The optical transparency of developing zebrafish provides new opportunities to investigate mechanisms governing adipocyte biology, however zebrafish adipocytes remain uncharacterized. We have developed methods for visualizing zebrafish adipocytes in vivo by labeling neutral lipid droplets with Nile Red. Our results establish that neutral lipid droplets first accumulate in visceral adipocytes during larval stages and increase in number and distribution as zebrafish grow. We show that the cellular anatomy of zebrafish adipocytes is similar to mammalian white adipocytes and identify peroxisome-proliferator activated receptor γ and fatty acid binding protein 11a as markers of the zebrafish adipocyte lineage. By monitoring adipocyte development prior to neutral lipid deposition, we find that the first visceral preadipocytes appear in association with the pancreas shortly after initiation of exogenous nutrition. Zebrafish reared in the absence of food fail to form visceral preadipocytes, indicating that exogenous nutrition is required for adipocyte development. These results reveal homologies between zebrafish and mammalian adipocytes and establish the zebrafish as a new model for adipocyte research.     

Host genetics exerts lifelong effects upon hindgut microbiota and its association with bovine growth and immunity

The gut microbiota is a complex ecological community that plays multiple critical roles within a host. Known intrinsic and extrinsic factors affect gut microbiota structure, but the influence of host genetics is understudied. To investigate the role of host genetics upon the gut microbiota structure, we performed a longitudinal study in which we evaluated the hindgut microbiota and its association with animal growth and immunity across life. We evaluated three different growth stages in an Angus-Brahman multibreed population with a graduated spectrum of genetic variation, raised under variable environmental conditions and diets. We found the gut microbiota structure was changed significantly during growth when preweaning, and fattening calves experienced large variations in diet and environmental changes. However, regardless of the growth stage, we found gut microbiota is significantly influenced by breed composition throughout life. Host genetics explained the relative abundances of 52.2%, 40.0%, and 37.3% of core bacterial taxa at the genus level in preweaning, postweaning, and fattening calves, respectively. Sutterella, Oscillospira, and Roseburia were consistently associated with breed composition at these three growth stages. Especially, butyrate-producing bacteria, Roseburia and Oscillospira, were associated with nine single-nucleotide polymorphisms (SNPs) located in genes involved in the regulation of host immunity and metabolism in the hindgut. Furthermore, minor allele frequency analysis found breed-associated SNPs in the short-chain fatty acids (SCFAs) receptor genes that promote anti-inflammation and enhance intestinal epithelial barrier functions. Our findings provide evidence of dynamic and lifelong host genetic effects upon gut microbiota, regardless of growth stages. We propose that diet, environmental changes, and genetic components may explain observed variation in critical hindgut microbiota throughout life.Subject terms: Microbiome, Agricultural genetics