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.

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.

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.

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.

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.

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.

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.

De novo synthesis of thiamine (vitamin B1) is the ancestral state in Plasmodium parasites - evidence from avian haemosporidians.

Parasites often have reduced genomes as their own genes become redundant when utilizing their host as a source of metabolites, thus losing their own de novo production of metabolites. Primate malaria parasites can synthesize vitamin B1 (thiamine) de novo but rodent malaria and other genome-sequenced apicomplexans cannot, as the three essential genes responsible for this pathway are absent in their genomes. The unique presence of functional thiamine synthesis genes in primate malaria parasites and their sequence similarities to bacterial orthologues, have led to speculations that this pathway was horizontally acquired from bacteria. Here we show that the genes essential for the de novo synthesis of thiamine are found also in avian Plasmodium species. Importantly, they are also present in species phylogenetically basal to all mammalian and avian Plasmodium parasites, i.e. Haemoproteus. Furthermore, we found that these genes are expressed during the blood stage of the avian malaria infection, indicating that this metabolic pathway is actively transcribed. We conclude that the ability to synthesize thiamine is widespread among haemosporidians, with a recent loss in the rodent malaria species.

Strong Maternal Effects on Gene Expression in Arabidopsis lyrata Hybrids.

Hybridization between populations or species can have pronounced fitness consequences. Yet little is known about how hybridization affects gene regulation. Three main models have been put forward to explain gene expression patterns in hybrids: additive, dominance, or parental effects. Here, we use high throughput RNA-sequencing to examine the extent to which hybrid gene expression follows predictions by each of the three models. We performed a reciprocal crossing experiment between two differentiated populations of the perennial herb Arabidopsis lyrata and sequenced RNA in rosette leaves of 12-week-old plants grown in greenhouse conditions. The two parental populations had highly differentiated gene expression patterns. In hybrids, a majority of genes showed intermediate expression relative to that of their parental populations (i.e., additive effects), but expression was frequently more similar to the maternal than to their paternal population (i.e., maternal effects). Allele-specific expression analyses showed that in the vast majority of cases, genes with pronounced maternal effect expressed both the maternal and the paternal allele. Maternal effects on hybrid gene expression have rarely been documented previously and our study suggests it could be more common than previously assumed. Whether the maternal effect on gene expression persists to later life-stages, and whether the variation in gene expression is manifested in other aspects of the phenotype, remain to be elucidated. Our findings are relevant for understanding the consequences of outbreeding and hybridization and open up several questions for future studies.

The transcriptome of the avian malaria parasite Plasmodium ashfordi displays host-specific gene expression.

Malaria parasites (Plasmodium spp.) include some of the world's most widespread and virulent pathogens. Our knowledge of the molecular mechanisms these parasites use to invade and exploit their hosts other than in mice and primates is, however, extremely limited. It is therefore imperative to characterize transcriptome-wide gene expression from nonmodel malaria parasites and how this varies across individual hosts. Here, we used high-throughput Illumina RNA sequencing on blood from wild-caught Eurasian siskins experimentally infected with a clonal strain of the avian malaria parasite Plasmodium ashfordi (lineage GRW2). Using a bioinformatic multistep approach to filter out host transcripts, we successfully assembled the blood-stage transcriptome of P. ashfordi. A total of 11 954 expressed transcripts were identified, and 7860 were annotated with protein information. We quantified gene expression levels of all parasite transcripts across three hosts during two infection stages - peak and decreasing parasitemia. Interestingly, parasites from the same host displayed remarkably similar expression profiles during different infection stages, but showed large differences across hosts, indicating that P. ashfordi may adjust its gene expression to specific host individuals. We further show that the majority of transcripts are most similar to the human parasite Plasmodium falciparum, and a large number of red blood cell invasion genes were discovered, suggesting evolutionary conserved invasion strategies between mammalian and avian Plasmodium. The transcriptome of P. ashfordi and its host-specific gene expression advances our understanding of Plasmodium plasticity and is a valuable resource as it allows for further studies analysing gene evolution and comparisons of parasite gene expression.

The avian transcriptome response to malaria infection.

Malaria parasites are highly virulent pathogens which infect a wide range of vertebrates. Despite their importance, the way different hosts control and suppress malaria infections remains poorly understood. With recent developments in next-generation sequencing techniques, however, it is now possible to quantify the response of the entire transcriptome to infections. We experimentally infected Eurasian siskins (Carduelis spinus) with avian malaria parasites (Plasmodium ashfordi), and used high-throughput RNA-sequencing to measure the avian transcriptome in blood collected before infection (day 0), during peak parasitemia (day 21 postinfection), and when parasitemia was decreasing (day 31). We found considerable differences in the transcriptomes of infected and uninfected individuals, with a large number of genes differentially expressed during both peak and decreasing parasitemia stages. These genes were overrepresented among functions involved in the immune system, stress response, cell death regulation, metabolism, and telomerase activity. Comparative analyses of the differentially expressed genes in our study to those found in other hosts of malaria (human and mouse) revealed a set of genes that are potentially involved in highly conserved evolutionary responses to malaria infection. By using RNA-sequencing we gained a more complete view of the host response, and were able to pinpoint not only well-documented host genes but also unannotated genes with clear significance during infection, such as microRNAs. This study shows how the avian blood transcriptome shifts in response to malaria infection, and we believe that it will facilitate further research into the diversity of molecular mechanisms that hosts utilize to fight malaria infections.

Sex- and tissue-specific profiles of chemosensory gene expression in a herbivorous gall-inducing fly (Diptera: Cecidomyiidae).

BACKGROUND: The chemical senses of insects mediate behaviors that are closely linked to survival and reproduction. The order Diptera contains two model organisms, the vinegar fly Drosophila melanogaster and the mosquito Anopheles gambiae, whose chemosensory genes have been extensively studied. Representing a third dipteran lineage with an interesting phylogenetic position, and being ecologically distinct by feeding on plants, the Hessian fly (Mayetiola destructor Say, Diptera: Cecidomyiidae) genome sequence has recently become available. Among plant-feeding insects, the Hessian fly is unusual in 'reprogramming' the plant to create a superior food and in being the target of plant resistance genes, a feature shared by plant pathogens. Chemoreception is essential for reproductive success, including detection of sex pheromone and plant-produced chemicals by males and females, respectively. RESULTS: We identified genes encoding 122 odorant receptors (OR), 28 gustatory receptors (GR), 39 ionotropic receptors (IR), 32 odorant binding proteins, and 7 sensory neuron membrane proteins in the Hessian fly genome. We then mapped Illumina-sequenced transcriptome reads to the genome to explore gene expression in male and female antennae and terminal abdominal segments. Our results reveal that a large number of chemosensory genes have up-regulated expression in the antennae, and the expression is in many cases sex-specific. Sex-specific expression is particularly evident among the Or genes, consistent with the sex-divergent olfactory-mediated behaviors of the adults. In addition, the large number of Ors in the genome but the reduced set of Grs and divergent Irs suggest that the short-lived adults rely more on long-range olfaction than on short-range gustation. We also report up-regulated expression of some genes from all chemosensory gene families in the terminal segments of the abdomen, which play important roles in reproduction. CONCLUSIONS: We show that a large number of the chemosensory genes in the Hessian fly genome have sex- and tissue-specific expression profiles. Our findings provide the first insights into the molecular basis of chemoreception in plant-feeding flies, representing an important advance toward a more complete understanding of olfaction in Diptera and its links to ecological specialization.

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.

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.

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.

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.