Molecular ecology of microbiomes in the wild: Common pitfalls, methodological advances and future directions.

The study of microbiomes across organisms and environments has become a prominent focus in molecular ecology. This perspective article explores common challenges, methodological advancements, and future directions in the field. Key research areas include understanding the drivers of microbiome community assembly, linking microbiome composition to host genetics, exploring microbial functions, transience and spatial partitioning, and disentangling non-bacterial components of the microbiome. Methodological advancements, such as quantifying absolute abundances, sequencing complete genomes, and utilizing novel statistical approaches, are also useful tools for understanding complex microbial diversity patterns. Our aims are to encourage robust practices in microbiome studies and inspire researchers to explore the next frontier of this rapidly changing field.

Coprophagy rapidly matures juvenile gut microbiota in a precocial bird.

Coprophagy is a behavior where animals consume feces, and has been observed across a wide range of species, including birds and mammals. The phenomenon is particularly prevalent in juveniles, but the reasons for this remain unclear. One hypothesis is that coprophagy enables offspring to acquire beneficial gut microbes that aid development. However, despite the potential importance of this behavior, studies investigating the effects in juveniles are rare. Here we experimentally test this idea by examining how ingestion of adult feces by ostrich chicks affects their gut microbiota development, growth, feeding behavior, pathogen abundance, and mortality. We conducted extensive longitudinal experiments for 8 weeks, repeated over 2 years. It involved 240 chicks, of which 128 were provided daily access to fresh fecal material from adults and 112 were simultaneously given a control treatment. Repeated measures, behavioral observations, and DNA metabarcoding of the microbial gut community, both prior to and over the course of the experiment, allowed us to evaluate multiple aspects of the behavior. The results show that coprophagy causes (a) marked shifts to the juvenile gut microbiota, including a major increase in diversity and rapid maturation of the microbial composition, (b) higher growth rates (fecal-supplemented chicks became 9.4% heavier at 8 weeks old), (c) changes to overall feeding behavior but no differences in feed intake, (d) lower abundance of a common gut pathogen (Clostridium colinum), and (e) lower mortality associated with gut disease. Together, our results suggest that the behavior of coprophagy in juveniles is highly beneficial and may have evolved to accelerate the development of gut microbiota.

Impact of ionizing radiation on the environmental microbiomes of Chornobyl wetlands.

Radioactive contamination has the potential to cause damage to DNA and other biomolecules. Anthropogenic sources of radioactive contamination include accidents in nuclear power plants, such as the one in Chornobyl in 1986 which caused long-term radioactive pollution. Studies on animals within radioactive zones have provided us with a greater understanding of how wildlife can persevere despite chronic radiation exposure. However, we still know very little about the effects of radiation on the microbial communities in the environment. We examined the impact of ionizing radiation and other environmental factors on the diversity and composition of environmental microbiomes in the wetlands of Chornobyl. We combined detailed field sampling along a gradient of radiation together with 16S rRNA high-throughput metabarcoding. While radiation did not affect the alpha diversity of the microbiomes in sediment, soil, or water, it had a significant effect on the beta diversity in all environment types, indicating that the microbial composition was affected by ionizing radiation. Specifically, we detected several microbial taxa that were more abundant in areas with high radiation levels within the Chornobyl Exclusion Zone, including bacteria and archaea known to be radioresistant. Our results reveal the existence of rich and diverse microbiomes in Chornobyl wetlands, with multiple taxonomic groups that are able to thrive despite the radioactive contamination. These results, together with additional field and laboratory-based approaches examining how microbes cope with ionizing radiation will help to forecast the functionality and re-naturalization dynamics of radiocontaminated environments.

Gene expression reveals immune response strategies of naïve Hawaiian honeycreepers experimentally infected with introduced avian malaria.

The unprecedented rise in the number of new and emerging infectious diseases in the last quarter century poses direct threats to human and wildlife health. The introduction to the Hawaiian archipelago of Plasmodium relictum and the mosquito vector that transmits the parasite has led to dramatic losses in endemic Hawaiian forest bird species. Understanding how mechanisms of disease immunity to avian malaria may evolve is critical as climate change facilitates increased disease transmission to high elevation habitats where malaria transmission has historically been low and the majority of the remaining extant Hawaiian forest bird species now reside. Here, we compare the transcriptomic profiles of highly susceptible Hawai'i 'amakihi (Chlorodrepanis virens) experimentally infected with P. relictum to those of uninfected control birds from a naïve high elevation population. We examined changes in gene expression profiles at different stages of infection to provide an in-depth characterization of the molecular pathways contributing to survival or mortality in these birds. We show that the timing and magnitude of the innate and adaptive immune response differed substantially between individuals that survived and those that succumbed to infection, and likely contributed to the observed variation in survival. These results lay the foundation for developing gene-based conservation strategies for Hawaiian honeycreepers by identifying candidate genes and cellular pathways involved in the pathogen response that correlate with a bird's ability to recover from malaria infection.

Microbiomes associated with avian malaria survival differ between susceptible Hawaiian honeycreepers and sympatric malaria-resistant introduced birds.

Of the estimated 55 Hawaiian honeycreepers (subfamily Carduelinae) only 17 species remain, nine of which the International Union for Conservation of Nature considers endangered. Among the most pressing threats to honeycreeper survival is avian malaria, caused by the introduced blood parasite Plasmodium relictum, which is increasing in distribution in Hawai'i as a result of climate change. Preventing further honeycreeper decline will require innovative conservation strategies that confront malaria from multiple angles. Research on mammals has revealed strong connections between gut microbiome composition and malaria susceptibility, illuminating a potential novel approach to malaria control through the manipulation of gut microbiota. One honeycreeper species, Hawai'i 'amakihi (Chlorodrepanis virens), persists in areas of high malaria prevalence, indicating they have acquired some level of immunity. To investigate if avian host-specific microbes may be associated with malaria survival, we characterized cloacal microbiomes and malaria infection for 174 'amakihi and 172 malaria-resistant warbling white-eyes (Zosterops japonicus) from Hawai'i Island using 16S rRNA gene metabarcoding and quantitative polymerase chain reaction. Neither microbial alpha nor beta diversity covaried with infection, but 149 microbes showed positive associations with malaria survivors. Among these were Escherichia and Lactobacillus spp., which appear to mitigate malaria severity in mammalian hosts, revealing promising candidates for future probiotic research for augmenting malaria immunity in sensitive endangered species.

Transcriptional response of individual Hawaiian Culex quinquefasciatus mosquitoes to the avian malaria parasite Plasmodium relictum.

BACKGROUND: Plasmodium parasites that cause bird malaria occur in all continents except Antarctica and are primarily transmitted by mosquitoes in the genus Culex. Culex quinquefasciatus, the mosquito vector of avian malaria in Hawai'i, became established in the islands in the 1820s. While the deadly effects of malaria on endemic bird species have been documented for many decades, vector-parasite interactions in avian malaria systems are relatively understudied. METHODS: To evaluate the gene expression response of mosquitoes exposed to a Plasmodium infection intensity known to occur naturally in Hawai'i, offspring of wild-collected Hawaiian Cx. quinquefasciatus were fed on a domestic canary infected with a fresh isolate of Plasmodium relictum GRW4 from a wild-caught Hawaiian honeycreeper. Control mosquitoes were fed on an uninfected canary. Transcriptomes of five infected and three uninfected individual mosquitoes were sequenced at each of three stages of the parasite life cycle: 24 h post feeding (hpf) during ookinete invasion; 5 days post feeding (dpf) when oocysts are developing; 10 dpf when sporozoites are released and invade the salivary glands. RESULTS: Differential gene expression analyses showed that during ookinete invasion (24 hpf), genes related to oxidoreductase activity and galactose catabolism had lower expression levels in infected mosquitoes compared to controls. Oocyst development (5 dpf) was associated with reduced expression of a gene with a predicted innate immune function. At 10 dpf, infected mosquitoes had reduced expression levels of a serine protease inhibitor, and further studies should assess its role as a Plasmodium agonist in C. quinquefasciatus. Overall, the differential gene expression response of Hawaiian Culex exposed to a Plasmodium infection intensity known to occur naturally in Hawai'i was low, but more pronounced during ookinete invasion. CONCLUSIONS: This is the first analysis of the transcriptional responses of vectors to malaria parasites in non-mammalian systems. Interestingly, few similarities were found between the response of Culex infected with a bird Plasmodium and those reported in Anopheles infected with human Plasmodium. The relatively small transcriptional changes observed in mosquito genes related to immune response and nutrient metabolism support conclusions of low fitness costs often documented in experimental challenges of Culex with avian Plasmodium.

Natural experiments and long-term monitoring are critical to understand and predict marine host-microbe ecology and evolution.

Marine multicellular organisms host a diverse collection of bacteria, archaea, microbial eukaryotes, and viruses that form their microbiome. Such host-associated microbes can significantly influence the host's physiological capacities; however, the identity and functional role(s) of key members of the microbiome ("core microbiome") in most marine hosts coexisting in natural settings remain obscure. Also unclear is how dynamic interactions between hosts and the immense standing pool of microbial genetic variation will affect marine ecosystems' capacity to adjust to environmental changes. Here, we argue that significantly advancing our understanding of how host-associated microbes shape marine hosts' plastic and adaptive responses to environmental change requires (i) recognizing that individual host-microbe systems do not exist in an ecological or evolutionary vacuum and (ii) expanding the field toward long-term, multidisciplinary research on entire communities of hosts and microbes. Natural experiments, such as time-calibrated geological events associated with well-characterized environmental gradients, provide unique ecological and evolutionary contexts to address this challenge. We focus here particularly on mutualistic interactions between hosts and microbes, but note that many of the same lessons and approaches would apply to other types of interactions.

Transcriptome assembly and differential gene expression of the invasive avian malaria parasite Plasmodium relictum in Hawai'i.

The malaria parasite Plasmodium relictum (lineage GRW4) was introduced less than a century ago to the native avifauna of Hawai'i, where it has since caused major declines of endemic bird populations. One of the native bird species that is frequently infected with GRW4 is the Hawai'i 'amakihi (Chlorodrepanis virens). To achieve a better understanding of the transcriptional activities of this virulent parasite, we performed a controlled challenge experiment of 15 'amakihi that were infected with GRW4. Blood samples containing malaria parasites were collected at two time points (intermediate and peak infection stages) from host individuals that were either experimentally infected by mosquitoes or inoculated with infected blood. We then used RNA sequencing to assemble a high-quality blood transcriptome of P. relictum GRW4, allowing us to quantify parasite expression levels inside individual birds. We found few significant differences (one to two transcripts) in GRW4 expression levels between host infection stages and between inoculation methods. However, 36 transcripts showed differential expression levels among all host individuals, indicating a potential presence of host-specific gene regulation across hosts. To reduce the extinction risk of the remaining native bird species in Hawai'i, genetic resources of the local Plasmodium lineage are needed to enable further molecular characterization of this parasite. Our newly built Hawaiian GRW4 transcriptome assembly, together with analyses of the parasite's transcriptional activities inside the blood of Hawai'i 'amakihi, can provide us with important knowledge on how to combat this deadly avian disease in the future.

Early-life gut dysbiosis linked to juvenile mortality in ostriches.

BACKGROUND: Imbalances in the gut microbial community (dysbiosis) of vertebrates have been associated with several gastrointestinal and autoimmune diseases. However, it is unclear which taxa are associated with gut dysbiosis, and if particular gut regions or specific time periods during ontogeny are more susceptible. We also know very little of this process in non-model organisms, despite an increasing realization of the general importance of gut microbiota for health. METHODS: Here, we examine the changes that occur in the microbiome during dysbiosis in different parts of the gastrointestinal tract in a long-lived bird with high juvenile mortality, the ostrich (Struthio camelus). We evaluated the 16S rRNA gene composition of the ileum, cecum, and colon of 68 individuals that died of suspected enterocolitis during the first 3 months of life (diseased individuals), and of 50 healthy individuals that were euthanized as age-matched controls. We combined these data with longitudinal environmental and fecal sampling to identify potential sources of pathogenic bacteria and to unravel at which stage of development dysbiosis-associated bacteria emerge. RESULTS: Diseased individuals had drastically lower microbial alpha diversity and differed substantially in their microbial beta diversity from control individuals in all three regions of the gastrointestinal tract. The clear relationship between low diversity and disease was consistent across all ages in the ileum, but decreased with age in the cecum and colon. Several taxa were associated with mortality (Enterobacteriaceae, Peptostreptococcaceae, Porphyromonadaceae, Clostridium), while others were associated with health (Lachnospiraceae, Ruminococcaceae, Erysipelotrichaceae, Turicibacter, Roseburia). Environmental samples showed no evidence of dysbiosis-associated bacteria being present in either the food, water, or soil substrate. Instead, the repeated fecal sampling showed that pathobionts were already present shortly after hatching and proliferated in individuals with low microbial diversity, resulting in high mortality several weeks later. CONCLUSIONS: Identifying the origins of pathobionts in neonates and the factors that subsequently influence the establishment of diverse gut microbiota may be key to understanding dysbiosis and host development. Video Abstract.

Host Transcriptional Responses to High- and Low-Virulent Avian Malaria Parasites.

The transcriptional response of hosts to genetically similar pathogens can vary substantially, with important implications for disease severity and host fitness. A low pathogen load can theoretically elicit both high and low host responses, as the outcome depends on both the effectiveness of the host at suppressing the pathogen and the ability of the pathogen to evade the immune system. Here, we investigate the transcriptional response of Eurasian siskins (Spinus spinus) to two closely related lineages of the malaria parasite Plasmodium relictum. Birds were infected with either the high-virulent lineage P. relictum SGS1, the low-virulent sister lineage P. relictum GRW4, or sham-injected (controls). Blood samples for RNA sequencing were collected at four time points during the course of infection, totaling 76 transcriptomes from 19 birds. Hosts infected with SGS1 experienced up to 87% parasitemia and major transcriptome shifts throughout the infection, and multiple genes showed strong correlation with parasitemia. In contrast, GRW4-infected hosts displayed low parasitemia (maximum 0.7%) with a minor transcriptional response. We furthermore demonstrate that the baseline gene expression levels of hosts prior to infection were irrelevant as immunocompetence markers, as they could not predict future pathogen load. This study shows that the magnitude of the host transcriptional response can differ markedly from related parasites with different virulence, and it enables a better understanding of the molecular interactions taking place between hosts and parasites.

Microbiome maturation during a unique developmental window.

Shortly after birth, mammals are colonized by a multitude of microbes derived from the mother and the environment. Studies in model organisms have demonstrated that the structure and composition of the gut microbiome of offspring steadily mature with increasing diversity during nursing and weaning (Sommer & Bäckhed, 2013). This period of microbiome assembly is critical for young mammals because the gut microbes they acquire will help train their immune system (Lathrop et al., 2011) with potential long-lasting effects on their health (Cox et al., 2014). In an article in this issue of Molecular Ecology, Stoffel et al. (2020) investigated the gut microbiota of northern elephant seals (Mirounga angustirostris) during a key developmental window. A month after giving birth, elephant seal mothers stop nursing their pups and return to the sea. As a consequence, their pups go from a diet of milk rich in fat to abruptly enter a post weaning fasting period which lasts for about two months while they remain with the colony. This particular life-history trait therefore offered the authors a unique and exciting opportunity to evaluate intrinsic factors contributing to gut microbiota development in a wild marine mammal.

Genomics of host-pathogen interactions: challenges and opportunities across ecological and spatiotemporal scales.

Evolutionary genomics has recently entered a new era in the study of host-pathogen interactions. A variety of novel genomic techniques has transformed the identification, detection and classification of both hosts and pathogens, allowing a greater resolution that helps decipher their underlying dynamics and provides novel insights into their environmental context. Nevertheless, many challenges to a general understanding of host-pathogen interactions remain, in particular in the synthesis and integration of concepts and findings across a variety of systems and different spatiotemporal and ecological scales. In this perspective we aim to highlight some of the commonalities and complexities across diverse studies of host-pathogen interactions, with a focus on ecological, spatiotemporal variation, and the choice of genomic methods used. We performed a quantitative review of recent literature to investigate links, patterns and potential tradeoffs between the complexity of genomic, ecological and spatiotemporal scales undertaken in individual host-pathogen studies. We found that the majority of studies used whole genome resolution to address their research objectives across a broad range of ecological scales, especially when focusing on the pathogen side of the interaction. Nevertheless, genomic studies conducted in a complex spatiotemporal context are currently rare in the literature. Because processes of host-pathogen interactions can be understood at multiple scales, from molecular-, cellular-, and physiological-scales to the levels of populations and ecosystems, we conclude that a major obstacle for synthesis across diverse host-pathogen systems is that data are collected on widely diverging scales with different degrees of resolution. This disparity not only hampers effective infrastructural organization of the data but also data granularity and accessibility. Comprehensive metadata deposited in association with genomic data in easily accessible databases will allow greater inference across systems in the future, especially when combined with open data standards and practices. The standardization and comparability of such data will facilitate early detection of emerging infectious diseases as well as studies of the impact of anthropogenic stressors, such as climate change, on disease dynamics in humans and wildlife.

Major shifts in gut microbiota during development and its relationship to growth in ostriches.

The development of gut microbiota during ontogeny is emerging as an important process influencing physiology, immunity and fitness in vertebrates. However, knowledge of how bacteria colonize the juvenile gut, how this is influenced by changes in the diversity of gut bacteria and to what extent this influences host fitness, particularly in nonmodel organisms, is lacking. Here we used 16S rRNA gene sequencing to describe the successional development of the faecal microbiome in ostriches (Struthio camelus, n = 66, repeatedly sampled) over the first 3 months of life and its relationship to growth. We found a gradual increase in microbial diversity with age that involved multiple colonization and extinction events and a major taxonomic shift in bacteria that coincided with the cessation of yolk absorption. Comparisons with the microbiota of adults (n = 5) revealed that the chicks became more similar in their microbial diversity and composition to adults as they aged. There was a five-fold difference in juvenile growth during development, and growth during the first week of age was strongly positively correlated with the abundance of the genus Bacteroides and negatively correlated with Akkermansia. After the first week, the abundances of six phylogenetically diverse families (Peptococcaceae, S24-7, Verrucomicrobiae, Anaeroplasmataceae, Streptococcaceae, Methanobacteriaceae) were associated with subsequent reductions in chick growth in an age-specific and transient manner. These results have broad implications for our understanding of the development of gut microbiota and its associations with animal growth.

Genomic Advances in Avian Malaria Research.

Haemosporidian parasites causing malaria-like diseases in birds are globally distributed and have been associated with reduced host fitness and mortality in susceptible bird species. This group of parasites has not only enabled a greater understanding of host specificity, virulence, and parasite dispersal, but has also been crucial in restructuring the evolutionary history of apicomplexans. Despite their importance, genomic resources of avian haemosporidians have proved difficult to obtain, and they have, as a result, been lagging behind the congeneric Plasmodium species infecting mammals. In this review, I discuss recent genomic advances in the field of avian malaria research, and outline outstanding questions that will become possible to investigate with the continued successful efforts to generate avian haemosporidian genomic data.

Insights into Avian Incomplete Dosage Compensation: Sex-Biased Gene Expression Coevolves with Sex Chromosome Degeneration in the Common Whitethroat.

Non-recombining sex chromosomes (Y and W) accumulate deleterious mutations and degenerate. This poses a problem for the heterogametic sex (XY males; ZW females) because a single functional gene copy often implies less gene expression and a potential imbalance of crucial expression networks. Mammals counteract this by dosage compensation, resulting in equal sex chromosome expression in males and females, whereas birds show incomplete dosage compensation with significantly lower expression in females (ZW). Here, we study the evolution of Z and W sequence divergence and sex-specific gene expression in the common whitethroat (Sylvia communis), a species within the Sylvioidea clade where a neo-sex chromosome has been formed by a fusion between an autosome and the ancestral sex chromosome. In line with data from other birds, females had lower expression than males at the majority of sex-linked genes. Results from the neo-sex chromosome region showed that W gametologs have diverged functionally to a higher extent than their Z counterparts, and that the female-to-male expression ratio correlated negatively with the degree of functional divergence of these gametologs. We find it most likely that sex-linked genes are being suppressed in females as a response to W chromosome degradation, rather than that these genes experience relaxed selection, and thus diverge more, by having low female expression. Overall, our data of this unique avian neo-sex chromosome system suggest that incomplete dosage compensation evolves, at least partly, through gradual accumulation of deleterious mutations at the W chromosome and declining female gene expression.

Plasmodium parasites of birds have the most AT-rich genes of eukaryotes.

The genomic architecture of organisms, including nucleotide composition, can be highly variable, even among closely-related species. To better understand the causes leading to structural variation in genomes, information on distinct and diverse genomic features is needed. Malaria parasites are known for encompassing a wide range of genomic GC-content and it has long been thought that Plasmodium falciparum, the virulent malaria parasite of humans, has the most AT-biased eukaryotic genome. Here, I perform comparative genomic analyses of the most AT-rich eukaryotes sequenced to date, and show that the avian malaria parasites Plasmodium gallinaceum, P. ashfordi, and P. relictum have the most extreme coding sequences in terms of AT-bias. Their mean GC-content is 21.21, 21.22 and 21.60 %, respectively, which is considerably lower than the transcriptome of P. falciparum (23.79 %) and other eukaryotes. This information enables a better understanding of genome evolution and raises the question of how certain organisms are able to prosper despite severe compositional constraints.

Measuring the gut microbiome in birds: Comparison of faecal and cloacal sampling.

The gut microbiomes of birds and other animals are increasingly being studied in ecological and evolutionary contexts. Numerous studies on birds and reptiles have made inferences about gut microbiota using cloacal sampling; however, it is not known whether the bacterial community of the cloaca provides an accurate representation of the gut microbiome. We examined the accuracy with which cloacal swabs and faecal samples measure the microbiota in three different parts of the gastrointestinal tract (ileum, caecum, and colon) using a case study on juvenile ostriches, Struthio camelus, and high-throughput 16S rRNA sequencing. We found that faeces were significantly better than cloacal swabs in representing the bacterial community of the colon. Cloacal samples had a higher abundance of Gammaproteobacteria and fewer Clostridia relative to the gut and faecal samples. However, both faecal and cloacal samples were poor representatives of the microbial communities in the caecum and ileum. Furthermore, the accuracy of each sampling method in measuring the abundance of different bacterial taxa was highly variable: Bacteroidetes was the most highly correlated phylum between all three gut sections and both methods, whereas Actinobacteria, for example, was only strongly correlated between faecal and colon samples. Based on our results, we recommend sampling faeces, whenever possible, as this sample type provides the most accurate assessment of the colon microbiome. The fact that neither sampling technique accurately portrayed the bacterial community of the ileum nor the caecum illustrates the difficulty in noninvasively monitoring gut bacteria located further up in the gastrointestinal tract. These results have important implications for the interpretation of avian gut microbiome studies.