Common structuring principles of the Drosophila melanogaster microbiome on a continental scale and between host and substrate

Pervasive tissue-, genetic background-, and allele-specific gene expression effects in Drosophila melanogaster

The pervasiveness of gene expression variation and its contribution to phenotypic variation and evolution is well known. This gene expression variation is context dependent, with differences in regulatory architecture often associated with intrinsic and environmental factors, and is modulated by regulatory elements that can act in cis (linked) or in trans (unlinked) relative to the genes they affect. So far, little is known about how this genetic variation affects the evolution of regulatory architecture among closely related tissues during population divergence. To address this question, we analyzed gene expression in the midgut, hindgut, and Malpighian tubule as well as microbiome composition in the two gut tissues in four Drosophila melanogaster strains and their F1 hybrids from two divergent populations: one from the derived, European range and one from the ancestral, African range. In both the transcriptome and microbiome data, we detected extensive tissue- and genetic background-specific effects, including effects of genetic background on overall tissue specificity. Tissue-specific effects were typically stronger than genetic background-specific effects, although the two gut tissues were not more similar to each other than to the Malpighian tubules. An examination of allele specific expression revealed that, while both cis and trans effects were more tissue-specific in genes expressed differentially between populations than genes with conserved expression, trans effects were more tissue-specific than cis effects. Despite there being highly variable regulatory architecture, this observation was robust across tissues and genetic backgrounds, suggesting that the expression of trans variation can be spatially fine-tuned as well as or better than cis variation during population divergence and yielding new insights into cis and trans regulatory evolution.

Microbiota-mediated competition between Drosophila species

BACKGROUND

The influence of microbiota in ecological interactions, and in particular competition, is poorly known. We studied competition between two insect species, the invasive pest Drosophila suzukii and the model Drosophila melanogaster, whose larval ecological niches overlap in ripe, but not rotten, fruit.

RESULTS

We discovered D. suzukii females prevent costly interspecific larval competition by avoiding oviposition on substrates previously visited by D. melanogaster. More precisely, D. melanogaster association with gut bacteria of the genus Lactobacillus triggered D. suzukii avoidance. However, D. suzukii avoidance behavior is condition-dependent, and D. suzukii females that themselves carry D. melanogaster bacteria stop avoiding sites visited by D. melanogaster. The adaptive significance of avoiding cues from the competitor's microbiota was revealed by experimentally reproducing in-fruit larval competition: reduced survival of D. suzukii larvae only occurred if the competitor had its normal microbiota.

CONCLUSIONS

This study establishes microbiotas as potent mediators of interspecific competition and reveals a central role for context-dependent behaviors under bacterial influence. Video Abstract.

Effects of heat and pesticide stress on life history, physiology and the gut microbiome of two congeneric damselflies that differ in stressor tolerance

The combined impact of toxicants and warming on organisms is getting increased attention in ecotoxicology, but is still hard to predict, especially with regard to heat waves. Recent studies suggested that the gut microbiome may provide mechanistic insights into the single and combined stressor effects on their host. We therefore investigated effects of sequential exposure to a heat spike and a pesticide on both the phenotype (life history and physiology) and the gut microbiome composition of damselfly larvae. We compared the fast-paced Ischnura pumilio, which is more tolerant to both stressors, with the slow-paced I. elegans, to obtain mechanistic insights into species-specific stressor effects. The two species differed in gut microbiome composition, potentially contributing to their pace-of-life differences. Intriguingly, there was a general resemblance between the stressor response patterns in the phenotype and in the gut microbiome, whereby both species responded broadly similar to the single and combined stressors. The heat spike negatively affected the life history of both species (increased mortality, reduced growth rate), which could be explained not only by shared negative effects on physiology (inhibition of acetylcholinesterase, increase of malondialdehyde), but also by shared effects on gut bacterial species' abundances. The pesticide only had negative effects (reduced growth rate, reduced net energy budget) in I. elegans. The pesticide generated shifts in the bacterial community composition (e.g. increased abundance of Sphaerotilus and Enterobacteriaceae in the gut microbiome of I. pumilio), which potentially contributed to the relatively higher pesticide tolerance of I. pumilio. Moreover, in line with the response patterns in the host phenotype, the effects of the heat spike and the pesticide on the gut microbiome were mainly additive. By contrasting two species differing in stress tolerance, our results suggest that response patterns in the gut microbiome may improve our mechanistic understanding of single and combined stressor effects.

Drosophila melanogaster microbiome is shaped by strict filtering and neutrality along a latitudinal cline

Microbiomes affect many aspects of host biology, but the eco-evolutionary forces that shape their diversity in natural populations remain poorly understood. Geographical gradients, such as latitudinal clines, generate predictable patterns in biodiversity at macroecological scales, but whether these macroscale processes apply to host-microbiome interactions is an open question. To address this question, we sampled the microbiomes of 13 natural populations of Drosophila melanogaster along a latitudinal cline in the eastern United States. The microbiomes were surprisingly consistent across the cline, as latitude did not predict either alpha or beta diversity. Only a narrow taxonomic range of bacteria were present in all microbiomes, indicating that strict taxonomic filtering by the host and neutral ecological dynamics are the primary factors shaping the fly microbiome. Our findings reveal the complexity of eco-evolutionary interactions shaping microbial variation in D. melanogaster and highlight the need for additional sampling of the microbiomes in natural populations along environmental gradients.

Patterns of Microbiome Composition Vary Across Spatial Scales in a Specialist Insect

Microbial communities associated with animals vary based on both intrinsic and extrinsic factors. Of many possible determinants affecting microbiome composition, host phylogeny, host diet, and local environment are the most important. How these factors interact across spatial scales is not well understood. Here, we seek to identify the main influences on microbiome composition in a specialist insect, the western corn rootworm (WCR; Diabrotica virgifera virgifera), by analyzing the bacterial communities of adults collected from their obligate host plant, corn (Zea mays), across several geographic locations and comparing the patterns in communities to its congeneric species, the northern corn rootworm (NCR; Diabrotica barberi). We found that bacterial communities of WCR and NCR shared a portion of their bacterial communities even when collected from disparate locations. However, within each species, the location of collection significantly influenced the composition of their microbiome. Correlations of geographic distance between sites with WCR bacterial community composition revealed different patterns at different spatial scales. Community similarity decreased with increased geographic distance at smaller spatial scales (~25 km between the nearest sites). At broad spatial scales (>200 km), community composition was not correlated with distances between sites, but instead reflected the historical invasion path of WCR across the United States. These results suggest bacterial communities are structured directly by dispersal dynamics at small, regional spatial scales, while landscape-level genetic or environmental differences may drive community composition across broad spatial scales in this specialist insect.

Bowel Movement: Integrating Host Mobility and Microbial Transmission Across Host Taxa

The gut microbiota of animals displays a high degree of plasticity with respect to environmental or dietary adaptations and is shaped by factors like social interactions, diet diversity or the local environment. But the contribution of these drivers varies across host taxa and our ability to explain microbiome variability within wild populations remains limited. Terrestrial animals have divergent mobility ranges and can either crawl, walk or fly, from a couple of centimeters toward thousands of kilometers. Animal movement has been little regarded in host microbiota frameworks, though it can directly influence major drivers of the host microbiota: (1) Aggregation movement can enhance social transmissions, (2) foraging movement can extend range of diet diversity, and (3) dispersal movement determines the local environment of a host. Here, I would like to outline how movement behaviors of different host taxa matter for microbial acquisition across mammals, birds as well as insects. Host movement can have contrasting effects and either reduce or enlarge spatial scale. Increased dispersal movement could dissolve local effects of sampling location, while aggregation could enhance inter-host transmissions and uniformity among social groups. Host movement can also extend the boundaries of microbial dispersal limitations and connect habitat patches across plant-pollinator networks, while the microbiota of wild populations could converge toward a uniform pattern when mobility is interrupted in captivity or laboratory settings. Hence, the implementation of host movement would be a valuable addition to the metacommunity concept, to comprehend microbial dispersal within and across trophic levels.

Long-term gut microbiome dynamics in Drosophila melanogaster reveal environment-specific associations between bacterial taxa at the family level

The influence of the microbiome on its host is well-documented, but the interplay of its members is not yet well-understood. Even for simple microbiomes, the interaction among members of the microbiome is difficult to study. Longitudinal studies provide a promising approach to studying such interactions through the temporal covariation of different taxonomic units. By contrast to most longitudinal studies, which span only a single host generation, we here present a post hoc analysis of a whole-genome dataset of 81 samples that follows microbiome composition for up to 180 host generations, which cover nearly 10 years. The microbiome diversity remained rather stable in replicated Drosophila melanogaster populations exposed to two different temperature regimes. The composition changed, however, systematically across replicates of the two temperature regimes. Significant associations between families, mostly specific to one temperature regime, indicate functional interdependence of different microbiome components. These associations also involve moderately abundant families, which emphasizes their functional importance, and highlights the importance of looking beyond the common constituents of the Drosophila microbiome.

Horizontal gene transfer-mediated bacterial strain variation affects host fitness in Drosophila

Background

How microbes affect host fitness and environmental adaptation has become a fundamental research question in evolutionary biology. To better understand the role of microbial genomic variation for host fitness, we tested for associations of bacterial genomic variation and Drosophila melanogaster offspring number in a microbial Genome Wide Association Study (GWAS).

Results

We performed a microbial GWAS, leveraging strain variation in the genus Gluconobacter, a genus of bacteria that are commonly associated with Drosophila under natural conditions. We pinpoint the thiamine biosynthesis pathway (TBP) as contributing to differences in fitness conferred to the fly host. While an effect of thiamine on fly development has been described, we show that strain variation in TBP between bacterial isolates from wild-caught D. melanogaster contributes to variation in offspring production by the host. By tracing the evolutionary history of TBP genes in Gluconobacter, we find that TBP genes were most likely lost and reacquired by horizontal gene transfer (HGT).

Conclusion

Our study emphasizes the importance of strain variation and highlights that HGT can add to microbiome flexibility and potentially to host adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01124-y.

Genome-Inferred Correspondence between Phylogeny and Metabolic Traits in the Wild Drosophila Gut Microbiome

Annotated genome sequences provide valuable insight into the functional capabilities of members of microbial communities. Nevertheless, most studies on the microbiome in animal guts use metagenomic data, hampering the assignment of genes to specific microbial taxa. Here, we make use of the readily culturable bacterial communities in the gut of the fruit fly Drosophila melanogaster to obtain draft genome sequences for 96 isolates from wild flies. These include 81 new de novo assembled genomes, assigned to three orders (Enterobacterales, Lactobacillales, and Rhodospirillales) with 80% of strains identified to species level using average nucleotide identity and phylogenomic reconstruction. Based on annotations by the RAST pipeline, among-isolate variation in metabolic function partitioned strongly by bacterial order, particularly by amino acid metabolism (Rhodospirillales), fermentation, and nucleotide metabolism (Lactobacillales) and arginine, urea, and polyamine metabolism (Enterobacterales). Seven bacterial species, comprising 2-3 species in each order, were well-represented among the isolates and included ≥5 strains, permitting analysis of metabolic functions in the accessory genome (i.e., genes not present in every strain). Overall, the metabolic function in the accessory genome partitioned by bacterial order. Two species, Gluconobacter cerinus (Rhodospirillales) and Lactiplantibacillus plantarum (Lactobacillales) had large accessory genomes, and metabolic functions were dominated by amino acid metabolism (G. cerinus) and carbohydrate metabolism (La. plantarum). The patterns of variation in metabolic capabilities at multiple phylogenetic scales provide the basis for future studies of the ecological and evolutionary processes shaping the diversity of microorganisms associated with natural populations of Drosophila.

Global population genetic structure and demographic trajectories of the black soldier fly, Hermetia illucens

Background

The black soldier fly (Hermetia illucens) is the most promising insect candidate for nutrient-recycling through bioconversion of organic waste into biomass, thereby improving sustainability of protein supplies for animal feed and facilitating transition to a circular economy. Contrary to conventional livestock, genetic resources of farmed insects remain poorly characterised. We present the first comprehensive population genetic characterisation of H. illucens. Based on 15 novel microsatellite markers, we genotyped and analysed 2862 individuals from 150 wild and captive populations originating from 57 countries on seven subcontinents.

Results

We identified 16 well-distinguished genetic clusters indicating substantial global population structure. The data revealed genetic hotspots in central South America and successive northwards range expansions within the indigenous ranges of the Americas. Colonisations and naturalisations of largely unique genetic profiles occurred on all non-native continents, either preceded by demographically independent founder events from various single sources or involving admixture scenarios. A decisive primarily admixed Polynesian bridgehead population serially colonised the entire Australasian region and its secondarily admixed descendants successively mediated invasions into Africa and Europe. Conversely, captive populations from several continents traced back to a single North American origin and exhibit considerably reduced genetic diversity, although some farmed strains carry distinct genetic signatures. We highlight genetic footprints characteristic of progressing domestication due to increasing socio-economic importance of H. illucens, and ongoing introgression between domesticated strains globally traded for large-scale farming and wild populations in some regions.

Conclusions

We document the dynamic population genetic history of a cosmopolitan dipteran of South American origin shaped by striking geographic patterns. These reflect both ancient dispersal routes, and stochastic and heterogeneous anthropogenic introductions during the last century leading to pronounced diversification of worldwide structure of H. illucens. Upon the recent advent of its agronomic commercialisation, however, current human-mediated translocations of the black soldier fly largely involve genetically highly uniform domesticated strains, which meanwhile threaten the genetic integrity of differentiated unique local resources through introgression. Our in-depth reconstruction of the contemporary and historical demographic trajectories of H. illucens emphasises benchmarking potential for applied future research on this emerging model of the prospering insect-livestock sector.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01029-w.

Predicted Metabolic Function of the Gut Microbiota of Drosophila melanogaster

An important goal for many nutrition-based microbiome studies is to identify the metabolic function of microbes in complex microbial communities and their impact on host physiology. This research can be confounded by poorly understood effects of community composition and host diet on the metabolic traits of individual taxa. Here, we investigated these multiway interactions by constructing and analyzing metabolic models comprising every combination of five bacterial members of the Drosophila gut microbiome (from single taxa to the five-member community of Acetobacter and Lactobacillus species) under three nutrient regimes. We show that the metabolic function of Drosophila gut bacteria is dynamic, influenced by community composition, and responsive to dietary modulation. Furthermore, we show that ecological interactions such as competition and mutualism identified from the growth patterns of gut bacteria are underlain by a diversity of metabolic interactions, and show that the bacteria tend to compete for amino acids and B vitamins more frequently than for carbon sources. Our results reveal that, in addition to fermentation products such as acetate, intermediates of the tricarboxylic acid (TCA) cycle, including 2-oxoglutarate and succinate, are produced at high flux and cross-fed between bacterial taxa, suggesting important roles for TCA cycle intermediates in modulating Drosophila gut microbe interactions and the potential to influence host traits. These metabolic models provide specific predictions of the patterns of ecological and metabolic interactions among gut bacteria under different nutrient regimes, with potentially important consequences for overall community metabolic function and nutritional interactions with the host.IMPORTANCE Drosophila is an important model for microbiome research partly because of the low complexity of its mostly culturable gut microbiota. Our current understanding of how Drosophila interacts with its gut microbes and how these interactions influence host traits derives almost entirely from empirical studies that focus on individual microbial taxa or classes of metabolites. These studies have failed to capture fully the complexity of metabolic interactions that occur between host and microbe. To overcome this limitation, we reconstructed and analyzed 31 metabolic models for every combination of the five principal bacterial taxa in the gut microbiome of Drosophila This revealed that metabolic interactions between Drosophila gut bacterial taxa are highly dynamic and influenced by cooccurring bacteria and nutrient availability. Our results generate testable hypotheses about among-microbe ecological interactions in the Drosophila gut and the diversity of metabolites available to influence host traits.

Predicted Metabolic Function of the Gut Microbiota of Drosophila melanogaster

An important goal for many nutrition-based microbiome studies is to identify the metabolic function of microbes in complex microbial communities and their impact on host physiology. This research can be confounded by poorly understood effects of community composition and host diet on the metabolic traits of individual taxa. Here, we investigated these multiway interactions by constructing and analyzing metabolic models comprising every combination of five bacterial members of the Drosophila gut microbiome (from single taxa to the five-member community of Acetobacter and Lactobacillus species) under three nutrient regimes. We show that the metabolic function of Drosophila gut bacteria is dynamic, influenced by community composition, and responsive to dietary modulation. Furthermore, we show that ecological interactions such as competition and mutualism identified from the growth patterns of gut bacteria are underlain by a diversity of metabolic interactions, and show that the bacteria tend to compete for amino acids and B vitamins more frequently than for carbon sources. Our results reveal that, in addition to fermentation products such as acetate, intermediates of the tricarboxylic acid (TCA) cycle, including 2-oxoglutarate and succinate, are produced at high flux and cross-fed between bacterial taxa, suggesting important roles for TCA cycle intermediates in modulating Drosophila gut microbe interactions and the potential to influence host traits. These metabolic models provide specific predictions of the patterns of ecological and metabolic interactions among gut bacteria under different nutrient regimes, with potentially important consequences for overall community metabolic function and nutritional interactions with the host.IMPORTANCE Drosophila is an important model for microbiome research partly because of the low complexity of its mostly culturable gut microbiota. Our current understanding of how Drosophila interacts with its gut microbes and how these interactions influence host traits derives almost entirely from empirical studies that focus on individual microbial taxa or classes of metabolites. These studies have failed to capture fully the complexity of metabolic interactions that occur between host and microbe. To overcome this limitation, we reconstructed and analyzed 31 metabolic models for every combination of the five principal bacterial taxa in the gut microbiome of Drosophila This revealed that metabolic interactions between Drosophila gut bacterial taxa are highly dynamic and influenced by cooccurring bacteria and nutrient availability. Our results generate testable hypotheses about among-microbe ecological interactions in the Drosophila gut and the diversity of metabolites available to influence host traits.

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.

The discovery, distribution, and diversity of DNA viruses associated with Drosophila melanogaster in Europe

Drosophila melanogaster is an important model for antiviral immunity in arthropods, but very few DNA viruses have been described from the family Drosophilidae. This deficiency limits our opportunity to use natural host-pathogen combinations in experimental studies, and may bias our understanding of the Drosophila virome. Here, we report fourteen DNA viruses detected in a metagenomic analysis of 6668 pool-sequenced Drosophila, sampled from forty-seven European locations between 2014 and 2016. These include three new nudiviruses, a new and divergent entomopoxvirus, a virus related to Leptopilina boulardi filamentous virus, and a virus related to Musca domestica salivary gland hypertrophy virus. We also find an endogenous genomic copy of galbut virus, a double-stranded RNA partitivirus, segregating at very low frequency. Remarkably, we find that Drosophila Vesanto virus, a small DNA virus previously described as a bidnavirus, may be composed of up to twelve segments and thus represent a new lineage of segmented DNA viruses. Two of the DNA viruses, Drosophila Kallithea nudivirus and Drosophila Vesanto virus are relatively common, found in 2 per cent or more of wild flies. The others are rare, with many likely to be represented by a single infected fly. We find that virus prevalence in Europe reflects the prevalence seen in publicly available datasets, with Drosophila Kallithea nudivirus and Drosophila Vesanto virus the only ones commonly detectable in public data from wild-caught flies and large population cages, and the other viruses being rare or absent. These analyses suggest that DNA viruses are at lower prevalence than RNA viruses in D.melanogaster, and may be less likely to persist in laboratory cultures. Our findings go some way to redressing an earlier bias toward RNA virus studies in Drosophila, and lay the foundation needed to harness the power of Drosophila as a model system for the study of DNA viruses.

Demographic analyses of a new sample of haploid genomes from a Swedish population of Drosophila melanogaster

European and African natural populations of Drosophila melanogaster have been the focus of several studies aiming at inferring demographic and adaptive processes based on genetic variation data. However, in these analyses little attention has been given to gene flow between African and European samples. Here we present a dataset consisting of 14 fully sequenced haploid genomes sampled from a natural population from the northern species range (Umeå, Sweden). We co-analyzed this new data with an African population to compare the likelihood of several competing demographic scenarios for European and African populations and show that gene flow improves the fit of demographic models to data.

Genomic Analysis of European Drosophila melanogaster Populations Reveals Longitudinal Structure, Continent-Wide Selection, and Previously Unknown DNA Viruses

Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatiotemporal patterns of genetic variation in natural populations have been challenging, as comprehensive sampling is logistically difficult, and sequencing of entire populations costly. Here, we address these issues using a collaborative approach, sequencing 48 pooled population samples from 32 locations, and perform the first continent-wide genomic analysis of genetic variation in European Drosophila melanogaster. Our analyses uncover longitudinal population structure, provide evidence for continent-wide selective sweeps, identify candidate genes for local climate adaptation, and document clines in chromosomal inversion and transposable element frequencies. We also characterize variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.