Balance between geographic,soil, and host tree parameters to shape soil microbiomes associated to clonal oak varies across soil zones along a European North–South transect

Tree root-associated microbiomes are shaped by geographic, soil physico-chemical, and host tree parameters. However, their respective impacts on microbiome variations in soils across larger spatial scales remain weakly studied. We out-planted saplings of oak clone DF159 (Quercus robur L.) as phytometer in four grassland field sites along a European North–South transect. After four years, we first compared the soil microbiomes of the tree root zone (RZ) and the tree root-free zone (RFZ). Then, we separately considered the total microbiomes of both zones, besides the microbiome with significant affinity to the RZ and compared their variability along the transect. Variations within the microbiome of the tree RFZ were shaped by geographic and soil physico-chemical changes, whereby bacteria responded more than fungi. Variations within both microbiomes of the tree RZ depended on the host tree and abiotic parameters. Based on perMANOVA and Mantel correlation tests, impacts of site specificities and geographic distance strongly decreased for the tree RZ affine microbiome. This pattern was more pronounced for fungi than bacteria. Shaping the microbiome of the soil zones in root proximity might be a mechanism mediating the acclimation of oaks to a wide range of environmental conditions across geographic regions.     

Host population genetics and biogeography structure the microbiome of the sponge Cliona delitrix

Sponges occur across diverse marine biomes and host internal microbial communities that can provide critical ecological functions. While strong patterns of host specificity have been observed consistently in sponge microbiomes, the precise ecological relationships between hosts and their symbiotic microbial communities remain to be fully delineated. In the current study, we investigate the relative roles of host population genetics and biogeography in structuring the microbial communities hosted by the excavating sponge Cliona delitrix. A total of 53 samples, previously used to demarcate the population genetic structure of C. delitrix, were selected from two locations in the Caribbean Sea and from eight locations across the reefs of Florida and the Bahamas. Microbial community diversity and composition were measured using Illumina‐based high‐throughput sequencing of the 16S rRNA V4 region and related to host population structure and geographic distribution. Most operational taxonomic units (OTUs) specific to Cliona delitrix microbiomes were rare, while other OTUs were shared with congeneric hosts. Across a large regional scale (>1,000 km), geographic distance was associated with considerable variability of the sponge microbiome, suggesting a distance–decay relationship, but little impact over smaller spatial scales (<300 km) was observed. Host population structure had a moderate effect on the structure of these microbial communities, regardless of geographic distance. These results support the interplay between geographic, environmental, and host factors as forces determining the community structure of microbiomes associated with C. delitrix. Moreover, these data suggest that the mechanisms of host regulation can be observed at the population genetic scale, prior to the onset of speciation.     

Landscape genomics of an obligate mutualism: Concordant and discordant population structures between the leafcutter ant Atta texana and its two main fungal symbiont types

To explore landscape genomics at the range limit of an obligate mutualism, we use genotyping‐by‐sequencing (ddRADseq) to quantify population structure and the effect of host–symbiont interactions between the northernmost fungus‐farming leafcutter ant Atta texana and its two main types of cultivated fungus. Genome‐wide differentiation between ants associated with either of the two fungal types is of the same order of magnitude as differentiation associated with temperature and precipitation across the ant's entire range, suggesting that specific ant–fungus genome–genome combinations may have been favoured by selection. For the ant hosts, we found a broad cline of genetic structure across the range, and a reduction of genetic diversity along the axis of range expansion towards the range margin. This population‐genetic structure was concordant between the ants and one cultivar type (M‐fungi, concordant clines) but discordant for the other cultivar type (T‐fungi). Discordance in population‐genetic structures between ant hosts and a fungal symbiont is surprising because the ant farmers codisperse with their vertically transmitted fungal symbionts. Discordance implies that (a) the fungi disperse also through between‐nest horizontal transfer or other unknown mechanisms, and (b) genetic drift and gene flow can differ in magnitude between each partner and between different ant–fungus combinations. Together, these findings imply that variation in the strength of drift and gene flow experienced by each mutualistic partner affects adaptation to environmental stress at the range margin, and genome–genome interactions between host and symbiont influence adaptive genetic differentiation of the host during range evolution in this obligate mutualism.     

Metacommunity theory for transmission of heritable symbionts within insect communities

Microbial organisms are ubiquitous in nature and often form communities closely associated with their host, referred to as the microbiome. The microbiome has strong influence on species interactions, but microbiome studies rarely take interactions between hosts into account, and network interaction studies rarely consider microbiomes. Here, we propose to use metacommunity theory as a framework to unify research on microbiomes and host communities by considering host insects and their microbes as discretely defined “communities of communities” linked by dispersal (transmission) through biotic interactions. We provide an overview of the effects of heritable symbiotic bacteria on their insect hosts and how those effects subsequently influence host interactions, thereby altering the host community. We suggest multiple scenarios for integrating the microbiome into metacommunity ecology and demonstrate ways in which to employ and parameterize models of symbiont transmission to quantitatively assess metacommunity processes in host‐associated microbial systems. Successfully incorporating microbiota into community‐level studies is a crucial step for understanding the importance of the microbiome to host species and their interactions.     

Fungal disease and temperature alter skin microbiome structure in an experimental salamander system

Pathogens compete with host microbiomes for space and resources. Their shared environment impacts pathogen–microbiome–host interactions, which can lead to variation in disease outcome. The skin microbiome of red‐backed salamanders (Plethodon cinereus) can reduce infection by the pathogen Batrachochytrium dendrobatidis (Bd) at moderate infection loads, with high species richness and high abundance of competitors as putative mechanisms. However, it is unclear if the skin microbiome can reduce epizootic Bd loads across temperatures. We conducted a laboratory experiment to quantify skin microbiome and host responses (P. cinereus: n = 87) to Bd at mimicked epizootic loads across temperatures (13, 17 and 21°C). We quantified skin microbiomes using 16S rRNA gene metabarcoding and identified operational taxonomic units (OTUs) taxonomically similar to culturable bacteria known to kill Bd (anti‐Bd OTUs). Prior to pathogen exposure, temperature changed the microbiome (OTU richness decreased by 12% and the abundance of anti‐Bd OTUs increased by 18% per degree increase in temperature), but these changes were not predictive of disease outcome. After exposure, Bd changed the microbiome (OTU richness decreased by 0.1% and the abundance of anti‐Bd OTUs increased by 0.2% per 1% increase in Bd load) and caused high host mortality across temperatures (35/45: 78%). Temperature indirectly impacted microbiome change and mortality through its direct effect on pathogen load. We did not find support for the microbiome impacting Bd load or host survival. Our research reveals complex host, pathogen, microbiome and environmental interactions to demonstrate that during epizootic events the microbiome will be unlikely to reduce pathogen invasion, even for putatively Bd‐resistant species.     

Genome assembly of Scorias spongiosa and comparative genomics provide insights into ecological adaptation of honeydew-dependent sooty mould fungi

Sooty moulds are fungi of economic importance and with unique lifestyle mainly growing on insect honeydew. However, the lack of genomic data hinders investigation of genetic mechanisms underlying their ecological adaptation. With long-read sequencing technology, we generated the genome of Scorias spongiosa, an extraordinary sooty mould fungus associated with honeydew of colony aphids and producing large fruiting bodies. A 24.21 Mb high-quality genome assembly with a N50 length of 3.37 Mb was obtained. The genome contained 7758 protein coding genes, 97.13% of which were homologous to known genes, and approximately 0.29 Mb repeat sequences. Comparative genomics showed S. spongiosa lost relatively more gene families and contained fewer species-specific genes and gene families, with many CAZyme families and sugar transporters reduced or absent. This study not only promotes understanding of the ecological adaptation of sooty moulds, but also provides valuable genomic data resource for future comparative genomic and genetic studies.     

Host genetic variation impacts microbiome composition across human body sites

BackgroundThe composition of bacteria in and on the human body varies widely across human individuals, and has been associated with multiple health conditions. While microbial communities are influenced by environmental factors, some degree of genetic influence of the host on the microbiome is also expected. This study is part of an expanding effort to comprehensively profile the interactions between human genetic variation and the composition of this microbial ecosystem on a genome- and microbiome-wide scale.ResultsHere, we jointly analyze the composition of the human microbiome and host genetic variation. By mining the shotgun metagenomic data from the Human Microbiome Project for host DNA reads, we gathered information on host genetic variation for 93 individuals for whom bacterial abundance data are also available. Using this dataset, we identify significant associations between host genetic variation and microbiome composition in 10 of the 15 body sites tested. These associations are driven by host genetic variation in immunity-related pathways, and are especially enriched in host genes that have been previously associated with microbiome-related complex diseases, such as inflammatory bowel disease and obesity-related disorders. Lastly, we show that host genomic regions associated with the microbiome have high levels of genetic differentiation among human populations, possibly indicating host genomic adaptation to environment-specific microbiomes.ConclusionsOur results highlight the role of host genetic variation in shaping the composition of the human microbiome, and provide a starting point toward understanding the complex interaction between human genetics and the microbiome in the context of human evolution and disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0759-1) contains supplementary material, which is available to authorized users.     

Micro-coevolution of host genetics with gut microbiome in three Chinese ethnic groups

Understanding the micro-coevolution of the human gut microbiome with host genetics is challenging but essential in both evolutionary and medical studies. To gain insight into the interactions between host genetic variation and the gut microbiome, we analyzed both the human genome and gut microbiome collected from a cohort of 190 students in the same boarding college and representing 3 ethnic groups, Uyghur, Kazakh, and Han Chinese. We found that differences in gut microbiome were greater between genetically distinct ethnic groups than those genetically closely related ones in taxonomic composition, functional composition, enterotype stratification, and microbiome genetic differentiation. We also observed considerable correlations between host genetic variants and the abundance of a subset of gut microbial species. Notably, interactions between gut microbiome species and host genetic variants might have coordinated effects on specific human phenotypes. Bacteroides ovatus, previously reported to modulate intestinal immunity, is significantly correlated with the host genetic variant rs12899811 (meta-P = 5.55 × 10⁻⁵), which regulates the VPS33B expression in the colon, acting as a tumor suppressor of colorectal cancer. These results advance our understanding of the micro-coevolution of the human gut microbiome and their interactive effects with host genetic variation on phenotypic diversity.     

Genome evolution and host‐microbiome shifts correspond with intraspecific niche divergence within harmful algal bloom‐forming Microcystis aeruginosa

Intraspecific niche divergence is an important driver of species range, population abundance and impacts on ecosystem functions. Genetic changes are the primary focus when studying intraspecific divergence; however, the role of ecological interactions, particularly host‐microbiome symbioses, is receiving increased attention. The relative importance of these evolutionary and ecological mechanisms has seen only limited evaluation. To address this question, we used Microcystis aeruginosa, the globally distributed cyanobacterium that dominates freshwater harmful algal blooms. These blooms have been increasing in occurrence and intensity worldwide, causing major economic and ecological damages. We evaluated 46 isolates of M. aeruginosa and their microbiomes, collected from 14 lakes in Michigan, USA, that vary over 20‐fold in phosphorus levels, the primary limiting nutrient in freshwater systems. Genomes of M. aeruginosa diverged along this phosphorus gradient in genomic architecture and protein functions. Fitness in low‐phosphorus lakes corresponded with additional shifts within M. aeruginosa including genome‐wide reductions in nitrogen use, an expansion of phosphorus assimilation genes and an alternative life history strategy of nonclonal colony formation. In addition to host shifts, despite culturing in common‐garden conditions, host‐microbiomes diverged along the gradient in taxonomy, but converged in function with evidence of metabolic interdependence between the host and its microbiome. Divergence corresponded with a physiological trade‐off between fitness in low‐phosphorus environments and growth rate in phosphorus‐rich conditions. Co‐occurrence of genotypes adapted to different nutrient environments in phosphorus‐rich lakes may have critical implications for understanding how M. aeruginosa blooms persist after initial nutrient depletion. Ultimately, we demonstrate that the intertwined effects of genome evolution, host life history strategy and ecological interactions between a host and its microbiome correspond with an intraspecific niche shift with important implications for whole ecosystem function.     

Host plants and edaphic factors influence the distribution and diversity of ectomycorrhizal fungal fruiting bodies within rainforests from Tshopo,Democratic Republic of the Congo

Ectomycorrhizal fungi constitute an important component of forest ecosystems that enhances plant nutrition and resistance against stresses. Diversity of ectomycorrhizal (EcM) fungi is, however, affected by host plant diversity and soil heterogeneity. This study provides information about the influence of host plants and soil resources on the diversity of ectomycorrhizal fungal fruiting bodies from rainforests of the Democratic Republic of the Congo. Based on the presence of fungal fruiting bodies, significant differences in the number of ectomycorrhizal fungi species existed between forest stand types (p < 0.001). The most ectomycorrhizal species‐rich forest was the Gilbertiodendron dewevrei‐dominated forest (61 species). Of all 93 species of ectomycorrhizal fungi, 19 demonstrated a significant indicator value for particular forest stand types. Of all analysed edaphic factors, the percentage of silt particles was the most important parameter influencing EcM fungi host plant tree distribution. Both host trees and edaphic factors strongly affected the distribution and diversity of EcM fungi. EcM fungi may have developed differently their ability to successfully colonise root systems in relation to the availability of nutrients.     

How to assemble a beneficial microbiome in three easy steps

There is great interest in explaining how beneficial microbiomes are assembled. Antibiotic‐producing microbiomes are arguably the most abundant class of beneficial microbiome in nature, having been found on corals, arthropods, molluscs, vertebrates and plant rhizospheres. An exemplar is the attine ants, which cultivate a fungus for food and host a cuticular microbiome that releases antibiotics to defend the fungus from parasites. One explanation posits long‐term vertical transmission of Pseudonocardia bacteria, which (somehow) evolve new compounds in arms‐race fashion against parasites. Alternatively, attines (somehow) selectively recruit multiple, non‐coevolved actinobacterial genera from the soil, enabling a ‘multi‐drug’ strategy against parasites. We reconcile the models by showing that when hosts fuel interference competition by providing abundant resources, the interference competition favours the recruitment of antibiotic‐producing (and ‐resistant) bacteria. This partner‐choice mechanism is more effective when at least one actinobacterial symbiont is vertically transmitted or has a high immigration rate, as in disease‐suppressive soils.     

Geography is more important than host plant use for the population genetic structure of a generalist insect herbivore

Population divergence can occur due to mechanisms associated with geographic isolation and/or due to selection associated with different ecological niches. Much of the evidence for selection‐driven speciation has come from studies of specialist insect herbivores that use different host plant species; however, the influence of host plant use on population divergence of generalist herbivores remains poorly understood. We tested how diet breadth, host plant species and geographic distance influence population divergence of the fall webworm (Hyphantria cunea; FW). FW is a broadly distributed, extreme generalist herbivore consisting of two morphotypes that have been argued to represent two different species: black‐headed and red‐headed. We characterized the differentiation of FW populations at two geographic scales. We first analysed the influence of host plant and geographic distance on genetic divergence across a broad continental scale for both colour types. We further analysed the influence of host plant, diet breadth and geographic distance on divergence at a finer geographic scale focusing on red‐headed FW in Colorado. We found clear genetic and morphological distinction between red‐ and black‐headed FW, and Colorado FW formed a genetic cluster distinct from other locations. Although both geographic distance and host plant use were correlated with genetic distance, geographic distance accounted for up to 3× more variation in genetic distance than did host plant use. As a rare study investigating the genetic structure of a widespread generalist herbivore over a broad geographic range (up to 3,000 km), our study supports a strong role for geographic isolation in divergence in this system.     

Genome‐wide association studies on the phyllosphere microbiome: Embracing complexity in host–microbe interactions

Environmental sequencing shows that plants harbor complex communities of microbes that vary across environments. However, many approaches for mapping plant genetic variation to microbe‐related traits were developed in the relatively simple context of binary host–microbe interactions under controlled conditions. Recent advances in sequencing and statistics make genome‐wide association studies (GWAS) an increasingly promising approach for identifying the plant genetic variation associated with microbes in a community context. This review discusses early efforts on GWAS of the plant phyllosphere microbiome and the outlook for future studies based on human microbiome GWAS. A workflow for GWAS of the phyllosphere microbiome is then presented, with particular attention to how perspectives on the mechanisms, evolution and environmental dependence of plant–microbe interactions will influence the choice of traits to be mapped.     

Host specificity of microbiome assembly and its fitness effects in phytoplankton

Insights into symbiosis between eukaryotic hosts and their microbiomes have shifted paradigms on what determines host fitness, ecology, and behavior. Questions remain regarding the roles of host versus environment in shaping microbiomes, and how microbiome composition affects host fitness. Using a model system in ecology, phytoplankton, we tested whether microbiomes are host-specific, confer fitness benefits that are host-specific, and remain conserved in time in their composition and fitness effects. We used an experimental approach in which hosts were cleaned of bacteria and then exposed to bacterial communities from natural environments to permit recruitment of microbiomes. We found that phytoplankton microbiomes consisted of a subset of taxa recruited from these natural environments. Microbiome recruitment was host-specific, with host species explaining more variation in microbiome composition than environment. While microbiome composition shifted and then stabilized over time, host specificity remained for dozens of generations. Microbiomes increased host fitness, but these fitness effects were host-specific for only two of the five species. The shifts in microbiome composition over time amplified fitness benefits to the hosts. Overall, this work solidifies the importance of host factors in shaping microbiomes and elucidates the temporal dynamics of microbiome compositional and fitness effects.Subject terms: Microbial ecology, Freshwater ecology     

Finding flies in the mushroom soup: Host specificity of fungus‐associated communities revisited with a novel molecular method

Fruiting bodies of fungi constitute an important resource for thousands of other taxa. The structure of these diverse assemblages has traditionally been studied with labour‐intensive methods involving cultivation and morphology‐based species identification, to which molecular information might offer convenient complements. To overcome challenges in DNA extraction and PCR associated with the complex chemical properties of fruiting bodies, we developed a pipeline applicable for extracting amplifiable total DNA from soft fungal samples of any size. Our protocol purifies DNA in two sequential steps: (a) initial salt–isopropanol extraction of all nucleic acids in the sample is followed by (b) an extra clean‐up step using solid‐phase reversible immobilization (SPRI) magnetic beads. The protocol proved highly efficient, with practically all of our samples—regardless of biomass or other properties—being successfully PCR‐amplified using metabarcoding primers and subsequently sequenced. As a proof of concept, we apply our methods to address a topical ecological question: is host specificity a major characteristic of fungus‐associated communities, that is, do different fungus species harbour different communities of associated organisms? Based on an analysis of 312 fungal fruiting bodies representing 10 species in five genera from three orders, we show that molecular methods are suitable for studying this rich natural microcosm. Comparing to previous knowledge based on rearing and morphology‐based identifications, we find a species‐rich assemblage characterized by a low degree of host specialization. Our method opens up new horizons for molecular analyses of fungus‐associated interaction webs and communities. Fruiting bodies of fungi constitute an important resource for thousands of other taxa. The structure of these diverse assemblages has traditionally been studied with labour‐intensive methods involving cultivation and morphology‐based species identification, to which molecular information might offer convenient complements. To overcome challenges in DNA extraction and PCR associated with the complex chemical properties of fruiting bodies, we developed a pipeline applicable for extracting amplifiable total DNA from soft fungal samples of any size. Our protocol purifies DNA in two sequential steps: (a) initial salt–isopropanol extraction of all nucleic acids in the sample is followed by (b) an extra clean‐up step using solid‐phase reversible immobilization (SPRI) magnetic beads. The protocol proved highly efficient, with practically all of our samples—regardless of biomass or other properties—being successfully PCR‐amplified using metabarcoding primers and subsequently sequenced. As a proof of concept, we apply our methods to address a topical ecological question: is host specificity a major characteristic of fungus‐associated communities, that is, do different fungus species harbour different communities of associated organisms? Based on an analysis of 312 fungal fruiting bodies representing 10 species in five genera from three orders, we show that molecular methods are suitable for studying this rich natural microcosm. Comparing to previous knowledge based on rearing and morphology‐based identifications, we find a species‐rich assemblage characterized by a low degree of host specialization. Our method opens up new horizons for molecular analyses of fungus‐associated interaction webs and communities.     

Genomewide signatures of selection in Epichloë reveal candidate genes for host specialization

Host specialization is a key process in ecological divergence and speciation of plant‐associated fungi. The underlying determinants of host specialization are generally poorly understood, especially in endophytes, which constitute one of the most abundant components of the plant microbiome. We addressed the genetic basis of host specialization in two sympatric subspecies of grass‐endophytic fungi from the Epichloë typhina complex: subsp. typhina and clarkii. The life cycle of these fungi entails unrestricted dispersal of gametes and sexual reproduction before infection of a new host, implying that the host imposes a selective barrier on viability of the progeny. We aimed to detect genes under divergent selection between subspecies, experiencing restricted gene flow due to adaptation to different hosts. Using pooled whole‐genome sequencing data, we combined F_ST and D_XY population statistics in genome scans and detected 57 outlier genes showing strong differentiation between the two subspecies. Genomewide analyses of nucleotide diversity (π), Tajima's D and dN/dS ratios indicated that these genes have evolved under positive selection. Genes encoding secreted proteins were enriched among the genes showing evidence of positive selection, suggesting that molecular plant–fungus interactions are strong drivers of endophyte divergence. We focused on five genes encoding secreted proteins, which were further sequenced in 28 additional isolates collected across Europe to assess genetic variation in a larger sample size. Signature of positive selection in these isolates and putative identification of pathogenic function supports our findings that these genes represent strong candidates for host specialization determinants in Epichloë endophytes. Our results highlight the role of secreted proteins as key determinants of host specialization.     

Host‐derived population genomics data provides insights into bacterial and diatom composition of the killer whale skin

Recent exploration into the interactions and relationship between hosts and their microbiota has revealed a connection between many aspects of the host's biology, health and associated micro‐organisms. Whereas amplicon sequencing has traditionally been used to characterize the microbiome, the increasing number of published population genomics data sets offers an underexploited opportunity to study microbial profiles from the host shotgun sequencing data. Here, we use sequence data originally generated from killer whale Orcinus orca skin biopsies for population genomics, to characterize the skin microbiome and investigate how host social and geographical factors influence the microbial community composition. Having identified 845 microbial taxa from 2.4 million reads that did not map to the killer whale reference genome, we found that both ecotypic and geographical factors influence community composition of killer whale skin microbiomes. Furthermore, we uncovered key taxa that drive the microbiome community composition and showed that they are embedded in unique networks, one of which is tentatively linked to diatom presence and poor skin condition. Community composition differed between Antarctic killer whales with and without diatom coverage, suggesting that the previously reported episodic migrations of Antarctic killer whales to warmer waters associated with skin turnover may control the effects of potentially pathogenic bacteria such as Tenacibaculum dicentrarchi. Our work demonstrates the feasibility of microbiome studies from host shotgun sequencing data and highlights the importance of metagenomics in understanding the relationship between host and microbial ecology.     

Bacterial microbiomes of individual ectomycorrhizal Pinus sylvestris roots are shaped by soil horizon and differentially sensitive to nitrogen addition

Plant roots select non‐random communities of fungi and bacteria from the surrounding soil that have effects on their health and growth, but we know little about the factors influencing their composition. We profiled bacterial microbiomes associated with individual ectomycorrhizal Pinus sylvestris roots colonized by different fungi and analyzed differences in microbiome structure related to soils from distinct podzol horizons and effects of short‐term additions of N, a growth‐limiting nutrient commonly applied as a fertilizer, but known to influence patterns of carbon allocation to roots. Ectomycorrhizal roots growing in soil from different horizons harboured distinct bacterial communities. The fungi colonizing individual roots had a strong effect on the associated bacterial communities. Even closely related species within the same ectomycorrhizal genus had distinct bacterial microbiomes in unfertilized soil, but fertilization removed this specificity. Effects of N were rapid and context dependent, being influenced by both soil type and the particular ectomycorrhizal fungi involved. Fungal community composition changed in soil from all horizons, but bacteria only responded strongly to N in soil from the B horizon where community structure was different and bacterial diversity was significantly reduced, possibly reflecting changed carbon allocation patterns.     

Standing geographic variation in eclosion time and the genomics of host race formation in Rhagoletis pomonella fruit flies

Taxa harboring high levels of standing variation may be more likely to adapt to rapid environmental shifts and experience ecological speciation. Here, we characterize geographic and host‐related differentiation for 10,241 single nucleotide polymorphisms in Rhagoletis pomonella fruit flies to infer whether standing genetic variation in adult eclosion time in the ancestral hawthorn (Crataegus spp.)‐infesting host race, as opposed to new mutations, contributed substantially to its recent shift to earlier fruiting apple (Malus domestica). Allele frequency differences associated with early vs. late eclosion time within each host race were significantly related to geographic genetic variation and host race differentiation across four sites, arrayed from north to south along a 430‐km transect, where the host races co‐occur in sympatry in the Midwest United States. Host fruiting phenology is clinal, with both apple and hawthorn trees fruiting earlier in the North and later in the South. Thus, we expected alleles associated with earlier eclosion to be at higher frequencies in northern populations. This pattern was observed in the hawthorn race across all four populations; however, allele frequency patterns in the apple race were more complex. Despite the generally earlier eclosion timing of apple flies and corresponding apple fruiting phenology, alleles on chromosomes 2 and 3 associated with earlier emergence were paradoxically at lower frequency in the apple than hawthorn host race across all four sympatric sites. However, loci on chromosome 1 did show higher frequencies of early eclosion‐associated alleles in the apple than hawthorn host race at the two southern sites, potentially accounting for their earlier eclosion phenotype. Thus, although extensive clinal genetic variation in the ancestral hawthorn race exists and contributed to the host shift to apple, further study is needed to resolve details of how this standing variation was selected to generate earlier eclosing apple fly populations in the North.     

Population bottlenecks constrain host microbiome diversity and genetic variation impeding fitness

It is becoming increasingly clear that microbial symbionts influence key aspects of their host’s fitness, and vice versa. This may fundamentally change our thinking about how microbes and hosts interact in influencing fitness and adaptation to changing environments. Here we explore how reductions in population size commonly experienced by threatened species influence microbiome diversity. Consequences of such reductions are normally interpreted in terms of a loss of genetic variation, increased inbreeding and associated inbreeding depression. However, fitness effects of population bottlenecks might also be mediated through microbiome diversity, such as through loss of functionally important microbes. Here we utilise 50 Drosophila melanogaster lines with different histories of population bottlenecks to explore these questions. The lines were phenotyped for egg-to-adult viability and their genomes sequenced to estimate genetic variation. The bacterial 16S rRNA gene was amplified in these lines to investigate microbial diversity. We found that 1) host population bottlenecks constrained microbiome richness and diversity, 2) core microbiomes of hosts with low genetic variation were constituted from subsets of microbiomes found in flies with higher genetic variation, 3) both microbiome diversity and host genetic variation contributed to host population fitness, 4) connectivity and robustness of bacterial networks was low in the inbred lines regardless of host genetic variation, 5) reduced microbial diversity was associated with weaker evolutionary responses of hosts in stressful environments, and 6) these effects were unrelated to Wolbachia density. These findings suggest that population bottlenecks reduce hologenomic variation (combined host and microbial genetic variation). Thus, while the current biodiversity crisis focuses on population sizes and genetic variation of eukaryotes, an additional focal point should be the microbial diversity carried by the eukaryotes, which in turn may influence host fitness and adaptability with consequences for the persistence of populations.