Inferring the chronological and biological age of individuals is fundamental to population ecology and our understanding of ageing itself, its evolution, and the biological processes that affect or even cause ageing. Epigenetic clocks based on DNA methylation (DNAm) at specific CpG sites show a strong correlation with chronological age in humans, and discrepancies between inferred and actual chronological age predict morbidity and mortality. Recently, a growing number of epigenetic clocks have been developed in non-model animals and we here review these studies. We also conduct a meta-analysis to assess the effects of different aspects of experimental protocol on the performance of epigenetic clocks for non-model animals. Two measures of performance are usually reported, the R2 of the association between the predicted and chronological age, and the mean/median absolute deviation (MAD) of estimated age from chronological age, and we argue that only the MAD reflects accuracy. R2 for epigenetic clocks based on the HorvathMammalMethylChip4 was higher and the MAD scaled to age range lower, compared with other DNAm quantification approaches. Scaled MAD tended to be lower among individuals in captive populations, and decreased with an increasing number of CpG sites. We conclude that epigenetic clocks can predict chronological age with relatively high accuracy, suggesting great potential in ecological epigenetics. We discuss general aspects of epigenetic clocks in the hope of stimulating further DNAm-based research on ageing, and perhaps more importantly, other key traits.

Africa is considered the "mother continent" from which the first hominins arose. The diversified wildlife and flora of Africa, which ranges from those in the scorching Sahara and Kalahari deserts to those in the vast Savannas and wet tropical forests, are also the most diverse of any continent. Although the continent's abundance and diversity of living resources have provided critical means of subsistence for its inhabitants, future utilization of this biodiversity will demand a fundamental understanding of genetic variation and its adaptive capabilities in the face of natural and man-made stressors. Molecular ecological insights have previously been gained from a variety of vertebrate species native to Africa, and some of these discoveries have larger evolutionary and conservation implications. Despite lagging in genomics research, African scientists are increasingly eager to use the increasingly accessible -omics technology to routinely sequence more animals and plants native to Africa. This overview, which focuses on Africa's vertebrate biodiversity, aims to provide a continental scale perspective on organismal and ecological adaptations discovered through prior genomics research, as well as what conceptually these findings suggest for future research.

Inbreeding depression refers to the reduced fitness of offspring produced by genetically-related individuals and is expected to be rare in large, outbred populations. When it occurs, marked fitness loss is possible as large populations can carry a substantial load of recessive harmful mutations which are normally sheltered at the heterozygous state. Using experimental cross data and genome-wide identity-by-descent (IBD) relationships from an outbred marine nine-spined stickleback (Pungitius pungitius) population, we documented a significant decrease in offspring survival probability with increasing parental IBD sharing associated with an average inbreeding load (B) of 10.5. Interestingly, we found that this relationship was also underlined by a positive effect of paternal inbreeding coefficient on offspring survival, suggesting that certain combinations of parental inbreeding and genetic relatedness among mates may promote offspring survival. Our results demonstrate the potential for substantial inbreeding load in an outbred population and emphasize the need to consider fine-scale genetic relatedness in future studies of inbreeding depression in the wild.

The rhizosphere is a vital soil compartment providing key plant-beneficial functions. However, little is known about the mechanisms driving viral diversity in the rhizosphere. Viruses can establish lytic or lysogenic interactions with their bacterial hosts. In the latter, they assume a dormant state integrated in the host genome and can be awakened by different perturbations that impact host cell physiology, triggering a viral bloom, which is potentially a fundamental mechanism driving soil viral diversity, as 22%-68% of soil bacteria are predicted to harbour dormant viruses. Here we assessed the viral bloom response in rhizospheric viromes by exposing them to three contrasting soil perturbation agents: earthworms, herbicide and antibiotic pollutant. The viromes were next screened for rhizosphere-relevant genes and also used as inoculant on microcosms incubations to test their impacts on pristine microbiomes. Our results show that while post-perturbation viromes diverged from control conditions, viral communities exposed to both herbicide and antibiotic pollutant were more similar to each other than those influenced by earthworms. The latter also favoured an increase in viral populations harbouring genes involved in plant-beneficial functions. Post-perturbation viromes inoculated on soil microcosms changed the diversity of pristine microbiomes, suggesting that viromes are important components of the soil ecological memory driving eco-evolutionary processes that determine future microbiome trajectories according to past events. Our findings demonstrate that viromes are active players in the rhizosphere and need to be considered in efforts to understand and control the microbial processes towards sustainable crop production.

Advancements in environmental DNA (eDNA) approaches have allowed for rapid and efficient species detections in diverse environments. Although most eDNA research is focused on leveraging genetic diversity to identify taxa, some recent studies have explored the potential for these approaches to detect within-species genetic variation, allowing for population genetic assessments and abundance estimates from environmental samples. However, we currently lack a framework outlining the key considerations specific to generating, analysing and applying eDNA data for these two purposes. Here, we discuss how various genetic markers differ with regard to genetic information and detectability in environmental samples and how analysis of eDNA samples differs from common tissue-based analyses. We then outline how it may be possible to obtain species absolute abundance estimates from eDNA by detecting intraspecific genetic variation in mixtures of DNA under multiple scenarios. We also identify the major causes contributing to allele detection and frequency errors in eDNA data, discuss their consequences for population-level analyses and outline bioinformatic approaches to detect and remove erroneous sequences. This review summarizes the key advances required to harness the full potential of eDNA-based intraspecific genetic variation to inform population-level questions in ecology, evolutionary biology and conservation management.

The evolution of Batesian mimicry - whereby harmless species avoid predation through their resemblance to harmful species - has long intrigued biologists. In rare cases, Batesian mimicry is linked to intraspecific colour variation, in which only some individuals within a population resemble a noxious 'model'. Here, we assess intraspecific colour variation within a widespread New Zealand stonefly, wherein highly melanized individuals of Zelandoperla closely resemble a chemically defended aposematic stonefly, Austroperla cyrene. We assess convergence in the colour pattern of these two species, compare their relative palatability to predators, and use genome-wide association mapping to assess the genetic basis of this resemblance. Our analysis reveals that melanized Zelandoperla overlap significantly with Austroperla in colour space but are significantly more palatable to predators, implying that they are indeed Batesian mimics. Analysis of 194,773 genome-wide SNPs reveals an outlier locus (ebony) strongly differentiating melanic versus non-melanic Zelandoperla. Genotyping of 338 specimens from a single Zelandoperla population indicates that ebony explains nearly 70% of the observed variance in melanism. As ebony has a well-documented role in insect melanin biosynthesis, our findings indicate this locus has a conserved function across deeply divergent hexapod lineages. Distributional records suggest a link between the occurrence of melanic Zelandoperla and the forested ecosystems where the model Austroperla is abundant, suggesting the potential for adaptive shifts in this system underpinned by environmental change.

It is a common observation that distinct lineages adapt to similar environments through similar "parallel" changes in phenotype. Less clear is whether parallel phenotypic change evolves through genetic variation in the same or distinct genes. This question of redundancy-the number of genetic answers to the same evolutionary problem-is an important yet poorly understood feature of adaptive traits, and is a major hurdle in our transition from understanding evolutionary responses to predicting them. In this issue of Molecular Ecology, Szukala et al. (2023) examine the redundancy of local adaptation in montane and alpine ecotypes that has occurred in the flowering plant Heliosperma pusillum. By comparing gene expression in each of four montane or alpine-adapted H. pusillum population pairs, the authors observe distinct genetic responses underlying each case of parallel local adaptation.

In epidemiology, endemicity characterizes sustained pathogen circulation in a geographical area, which involves a circulation that is not being maintained by external introductions. Because it could potentially shape the design of public health interventions, there is an interest in fully uncovering the endemic pattern of a disease. Here, we use a phylogeographic approach to investigate the endemic signature of rabies virus (RABV) circulation in Cambodia. Cambodia is located in one of the most affected regions by rabies in the world, but RABV circulation between and within Southeast Asian countries remains understudied. Our analyses are based on a new comprehensive data set of 199 RABV genomes collected between 2014 and 2017 as well as previously published Southeast Asian RABV sequences. We show that most Cambodian sequences belong to a distinct clade that has been circulating almost exclusively in Cambodia. Our results thus point towards rabies circulation in Cambodia that does not rely on external introductions. We further characterize within-Cambodia RABV circulation by estimating lineage dispersal metrics that appear to be similar to other settings, and by performing landscape phylogeographic analyses to investigate environmental factors impacting the dispersal dynamic of viral lineages. The latter analyses do not lead to the identification of environmental variables that would be associated with the heterogeneity of viral lineage dispersal velocities, which calls for a better understanding of local dog ecology and further investigations of the potential drivers of RABV spread in the region. Overall, our study illustrates how phylogeographic investigations can be performed to assess and characterize viral endemicity in a context of relatively limited data.

Untangling how factors such as environment, host, associations among bacterial species and dispersal predict microbial composition is a fundamental challenge. In this study, we use complementary machine-learning approaches to quantify the relative role of these factors in shaping microbiome variation of the blacklegged tick Ixodes scapularis. I. scapularis is the most important vector for Borrelia burgdorferi (the causative agent for Lyme disease) in the U.S. as well as a range of other important zoonotic pathogens. Yet the relative role of the interactions between pathogens and symbionts compared to other ecological forces is unknown. We found that positive associations between microbes where the occurrence of one microbe increases the probability of observing another, including between both pathogens and symbionts, was by far the most important factor shaping the tick microbiome. Microclimate and host factors played an important role for a subset of the tick microbiome including Borrelia (Borreliella) and Ralstonia, but for the majority of microbes, environmental and host variables were poor predictors at a regional scale. This study provides new hypotheses on how pathogens and symbionts might interact within tick species, as well as valuable predictions for how some taxa may respond to changing climate.

Rodent virus communities (viromes) are overrepresented with zoonotic viruses, and as such are a key host system for the study of zoonotic viruses. However, the extent of viral diversity beyond characterized zoonotic viruses, and the factors that modulate the viromes of rodents remain opaque. In this issue of Molecular Ecology, Raghwani et al. (2023) use rodents as a model to understand the role of seasonality in dictating virome abundance and composition-a factor known to play an important role in most animal one-host, one-pathogen systems. These data are not only highly relevant to rodents, but have broad applications across understanding and disentangling animal virome ecology.

The microbiome has been hypothesized as a driving force of phenotypic variation in host organisms that is capable of extending metabolic processes, altering development and in some cases, conferring novel functions that are critical for survival. Only a few studies have directly shown a causal role for the environmental microbiome in altering host phenotypic features. To assess the extent to which environmental microbes induce variation in host life-history traits and behaviour, we inoculated axenic Drosophila melanogaster with microbes isolated from drosophilid populations collected from two different field sites and generated two populations with distinct bacterial and fungal profiles. We show that microbes isolated from environmental sites with modest abiotic differences induce large variation in host reproduction, fatty acid levels, stress tolerance and sleep behaviour. Importantly, clearing microbes from each experimental population removed the phenotypic differences. The results support the causal role of environmental microbes as drivers of host phenotypic plasticity and potentially, rapid adaptation and evolution.

Continental East Asia has a mild Pleistocene climate and a complex recent geological history. Phylogeographic studies of animals over the last 30 years have produced several distinctive patterns. Glaciation refugia are numerous and are not restricted to any particular regions. Most of them are localized and species-specific, although several large refugia, for example the mountains of SW China, are shared by multiple species and have refugia-within-refugia. Furthermore, postglaciation range expansion events vary greatly in time, scale and direction. Large-scale south-to-north post-LGM expansions are few and mostly occurred in the northern regions. Additionally, several unique geographic features, including the three-step terrain of China and the northern arid belt, have significant impacts on many species histories. Overall, the impacts of Pleistocene glaciations, particularly the LGM, on species history vary drastically from nondetectable to significant. The impacts are the least for species from the southwestern region and are most dominant for species from the north. Geological events play a more significant role in shaping species history than Pleistocene climatic changes. Phylogeographic patterns among animals species are highly consistent with those of plants. Future phylogeographic endeavour in East Asia should be hypothesis-driven and seek processes that underlie common patterns. The wide use of genomic data allow accurate estimates of historical population processes and exploration of older history beyond the Pleistocene.

Polyploids recurrently emerge in angiosperms, but most polyploids are likely to go extinct before establishment due to minority cytotype exclusion, which may be specifically a constraint for dioecious plants. Here we test the hypothesis that a stable sex-determination system and spatial/ecological isolation facilitate the establishment of dioecious polyploids. We determined the ploidy levels of 351 individuals from 28 populations of the dioecious species Salix polyclona, and resequenced 190 individuals of S. polyclona and related taxa for genomic diversity analyses. The ploidy survey revealed a frequency 52% of tetraploids in S. polyclona, and genomic k-mer spectra analyses suggested an autopolyploid origin for them. Comparisons of diploid male and female genomes identified a female heterogametic sex-determining factor on chromosome 15, which probably also acts in the dioecious tetraploids. Phylogenetic analyses revealed two diploid clades and a separate clade/grade of tetraploids with a distinct geographic distribution confined to western and central China, where complex mountain systems create higher levels of environmental heterogeneity. Fossil-calibrated phylogenies showed that the polyploids emerged during 7.6-2.3 million years ago, and population demographic histories largely matched the geological and climatic history of the region. Our results suggest that inheritance of the sex-determining system from the diploid progenitor as intrinsic factor and spatial isolation as extrinsic factor may have facilitated the preservation and establishment of polyploid dioecious populations.

Fast-X evolution has been observed in a range of heteromorphic sex chromosomes. However, it remains unclear how early in the process of sex chromosome differentiation the Fast-X effect becomes detectible. Recently, we uncovered an extreme variation in sex chromosome heteromorphism across poeciliid fish species. The common guppy, Poecilia reticulata, Endler's guppy, P. wingei, swamp guppy, P. picta and para guppy, P. parae, appear to share the same XY system and exhibit a remarkable range of heteromorphism. Species outside this group lack this sex chromosome system. We combined analyses of sequence divergence and polymorphism data across poeciliids to investigate X chromosome evolution as a function of hemizygosity and reveal the causes for Fast-X effects. Consistent with the extent of Y degeneration in each species, we detect higher rates of divergence on the X relative to autosomes, a signal of Fast-X evolution, in P. picta and P. parae, species with high levels of X hemizygosity in males. In P. reticulata, which exhibits largely homomorphic sex chromosomes and little evidence of hemizygosity, we observe no change in the rate of evolution of X-linked relative to autosomal genes. In P. wingei, the species with intermediate sex chromosome differentiation, we see an increase in the rate of nonsynonymous substitutions on the older stratum of divergence only. We also use our comparative approach to test for the time of origin of the sex chromosomes in this clade. Taken together, our study reveals an important role of hemizygosity in Fast-X evolution.

The study of chromosomal inversion polymorphisms has received much recent attention, particularly in cases where inversions have drastic effects on phenotypes and fitness (e.g. lethality of homozygotes). Less attention has been paid to the question of the maintenance of inversion polymorphisms that show only weak effects. Here, we study the maintenance of such an inversion polymorphism that links 250 genes on chromosome Tgu11 in the zebra finch (Taeniopygia guttata). Based on data from over 6000 captive birds, we estimated the effects of this inversion on a wide range of fitness-related traits. We found that, compared with the ancestral allele A, the inverted allele D had small additive beneficial effects on male siring success and on female fecundity. These fitness-enhancing effects may explain the initial spread of the derived D allele (allele frequency 53%). However, individuals that were homozygous for D had a slightly lower survival rate, which may explain why the D allele has not spread to fixation. We used individual-based simulations to examine how an inversion polymorphism with such antagonistic fitness effects behaves over time. Our results indicate that polymorphisms become stabilized at an intermediate allele frequency if the inversion links an additively beneficial allele of small effect size to a recessive weakly deleterious mutation, overall resulting in weak net heterosis. Importantly, this conclusion remains valid over a wide range of selection coefficients against the homozygous DD (up to lethality), suggesting that the conditions needed to maintain the polymorphism may frequently be met. However, the simulations also suggest that in our zebra finch populations, the estimated recessive deleterious effect of the D allele (on survival in captivity) is not quite large enough to prevent fixation of the D allele in the long run. Estimates of fitness effects from free-living populations are needed to validate these results.

In the family Fagaceae, fertilization is delayed by several weeks to 1 year after pollination, leading to 1- or 2-year fruiting species depending on whether fruiting occurs in the same or the next year after flowering. To investigate physiological responses underlying the regulation of delayed fertilization, we monitored seasonal changes in genome-wide gene expression in tissues including leaves and buds over 2 years under natural conditions in one- (Quercus glauca) and 2-year fruiting species (Lithocarpus edulis). Genes associated with metabolic changes in response to winter cold, photosynthesis and cell proliferation, which are essential for survival and growth, showed highly conserved seasonal expression profiles between species. However, seasonal expression profiles diverged between species in genes associated with pollination, an important process contributing to the origin and maintenance of the reproductive barrier between plant species. By comparing seasonal progression of ovule development and gene expression in pistillate flowers, we revealed that ovules started developing after winter in the 2-year fruiting species, which could be linked to the activation of genes involved in fertilization and female gametophyte development after winter. These findings suggest that the 2-year fruiting species may have evolved a requirement of winter cold to prevent fertilization before winter and facilitate fertilization and embryo development in the following spring when temperature rises. This study offers new possibilities to explore the evolution of reproductive strategies in Fagaceae.

Anthropogenic changes have altered the historical distributions of many North American taxa. As environments shift, ecological and evolutionary processes can combine in complex ways to either stimulate or inhibit range expansion. Here, we examined the role of evolution in a rapid range expansion whose ecological context has been well-documented, Anna's Hummingbird (Calypte anna). Previous studies have suggested that the C. anna range expansion is the result of an ecological release facilitated by human-mediated environmental changes, where access to new food sources have allowed further filling of the abiotic niche. We examined the role of gene flow and adaptation during range expansion from their native California breeding range, north into Canada and east into New Mexico and Texas, USA. Using low coverage whole genome sequencing we found high genetic diversity, low divergence, and little evidence of selection on the northern and eastern expansion fronts. Additionally, there are no clear barriers to gene flow across the native and expanded range. The lack of selective signals between core and expanded ranges could reflect (i) an absence of novel selection pressure in the expanded range (supporting the ecological release hypothesis), (ii) swamping of adaptive variation due to high gene flow, or (iii) limitations of genome scans for detecting small shifts in allele frequencies across many loci. Nevertheless, our results provide an example where strong selection is not apparent during a rapid, contemporary range shift.

Major histocompatibility complex (MHC) genes are crucial for the adaptive immune response of jawed vertebrates. Their variation, reaching extreme levels, is driven mainly by an arms race between hosts and pathogens. One hypothesised mechanism contributing to MHC polymorphism is adaptive introgression, the exchange of genetic variants between hybridising species favoured by selection, yet its effect on MHC variation is poorly understood. Detection of adaptive MHC introgression, though challenging, may be facilitated by the analysis of species complexes forming multiple hybrid zones. Here, we investigated MHC introgression in six hybrid zones formed by seven species of Podarcis lizards inhabiting the Iberian Peninsula. To differentiate adaptive introgression from neutral introgression, we compared the patterns of gene exchange in MHC and genome-wide markers. We found elevated sharing of MHC alleles in the proximity of contact beyond the areas of detectable genome-wide admixture in most hybrid zones and, in half of them, asymmetric MHC exchange. In general, the elevated MHC allele sharing between species pairs with abutting ranges compared to geographically isolated species pairs also supports the prevalence of introgression. Collectively, our results demonstrate widespread MHC introgression in the Iberian Podarcis complex and suggest its adaptiveness. Contrary to previous results from Triturus newts, we did not observe differences in the rate of introgression between MHC classes. Our work adds support to the emerging view of adaptive introgression as a key mechanism shaping MHC diversity. It also raises questions about the effect of elevated MHC variation and factors leading to the asymmetry of adaptive introgression.

In nature, organisms have to cope with constantly changing environments. In certain conditions, it may be advantageous for the parents to pass on information about the environment, or resources to their offspring. Such transfers are known as parental effects, and they are well documented in plants and animals, but not in other eukaryotes, such as fungi. Many fungi disperse through spores, and fungal spores can potentially carry information or resources to the next generation. Understanding parental effects and their evolutionary consequences in fungi is of vital importance as they perform crucial ecosystem functions. In this study, we investigated whether parental effects are present in the filamentous fungus Neurospora crassa, how long do they last, whether the effects are adaptive, and what is their mechanism. We performed a fully factorial match/mismatch experiment for a good and a poor quality environment, in which we measured the initial growth of strains that experienced either a matched or mismatched environment in their previous generation. We found a strong silver-spoon effect in initial mycelium growth, which lasted for one generation, and increased fitness during competition experiments. By using deletion mutants that lacked key genes in epigenetic processes, we show that epigenetic mechanisms are not involved in this effect. Instead, we show that spore glycogen content, glucose availability and a radical transcription shift in spores are the main mechanisms behind this parental effect.

Evidence supporting non-DNA sequence-based inheritance in animals has increasingly been described in recent years, often under intergenerational inheritance or transgenerational epigenetic inheritance (TEI). Existence of the latter, a stronger indicator of germline transmission, has been demonstrated in invertebrates and mammals alike. The mechanisms and physiological implications of TEI, however, remain unclear. Here, in an unbiased approach, we compared existing transcriptomic data associated with so far available Drosophila models of inter- and trans-, and rodent models of inter-generational inheritance; observed phenotypic cross-species conservation and cross-generation directionality shift therein; and confirmed these observations experimentally in flies. Specifically, previous models of cold and diet-induced inheritance in both flies and mice were commonly associated with altered regulation of proteolysis genes. Besides, fly TEI models were in general characterized by opposite phenotypic regulation in transgenerational offsprings, compared to the ancestors. As insulin-producing cell (IPC) ablation was also associated with proteolysis gene dysregulation in one of the mouse models, we opted to use genetic ablation of IPCs in flies for the experimental validation. Remarkably, the ablation led to transcriptomic alterations across multiple generations, with dysregulated genes showing proteolysis enrichment. Similarly, phenotypic directionality changed in the opposite direction in transgenerational offsprings, in comparison of the ancestors. These results support evolutionary conservation, and both physiologically adaptive and maladaptive consequences of germline mediated epigenetic inheritance.

Fire is a major evolutionary and ecological driver that shapes biodiversity in forests. While above-ground community responses to fire have been well-documented, those below-ground are much less understood. However, below-ground communities, including fungi, play key roles in forests and facilitate the recovery of other organisms after fire. Here, we used internal transcribed spacer (ITS) meta-barcoding data from forests with three different times since fire [short (3 years), medium (13-19 years) and long (>26 years)] to characterize the temporal responses of soil fungal communities across functional groups, ectomycorrhizal exploration strategies and inter-guild associations. Our findings indicate that fire effects on fungal communities are strongest in the short to medium term, with clear distinctions between communities in forests with a short time (3 years) since fire, a medium time (13-19 years) and a long time (>26 years) since fire. Ectomycorrhizal fungi were disproportionately impacted by fire relative to saprotrophs, but the direction of the response varied depending on morphological structures and exploration strategies. For instance, short-distance ectomycorrhizal fungi increased with recent fire, while medium-distance (fringe) ectomycorrhizal fungi decreased. Further, we detected strong, negative inter-guild associations between ectomycorrhizal and saprotrophic fungi but only at medium and long times since fire. Given the functional significance of fungi, the temporal changes in fungal composition, inter-guild associations and functional groups after fire demonstrated in our study may have functional implications that require adaptive management to curtail.

Microbial rhodopsins are widely distributed in aquatic environments and may significantly contribute to phototrophy and energy budgets in global oceans. However, the study of freshwater rhodopsins has been largely limited. Here, we explored the diversity, ecological distribution, and expression of opsin genes that encode the apoproteins of type I rhodopsins in humic and clearwater lakes with contrasting physicochemical and optical characteristics. Using metagenomes and metagenome-assembled genomes, we recovered opsin genes from a wide range of taxa, mostly predicted to encode green light-absorbing proton pumps. Viral opsin and novel bacterial opsin clades were recovered. Opsin genes occurred more frequently in taxa from clearwater than from humic water, and opsins in some taxa have nontypical ion-pumping motifs that might be associated with physicochemical conditions of these two freshwater types. Analyses of the surface layer of 33 freshwater systems revealed an inverse correlation between opsin gene abundance and lake dissolved organic carbon (DOC). In humic water with high terrestrial DOC and light-absorbing humic substances, opsin gene abundance was low and dramatically declined within the first few meters, whereas the abundance remained relatively high along the bulk water column in clearwater lakes with low DOC, suggesting opsin gene distribution is influenced by lake optical properties and DOC. Gene expression analysis confirmed the significance of rhodopsin-based phototrophy in clearwater lakes and revealed different diel expressional patterns among major phyla. Overall, our analyses revealed freshwater opsin diversity, distribution and expression patterns, and suggested the significance of rhodopsin-based phototrophy in freshwater energy budgets, especially in clearwater lakes.

Populations on the edge of a species' distribution may represent an important source of adaptive diversity, yet these populations tend to be more fragmented and are more likely to be geographically isolated. Lack of genetic exchanges between such populations, due to barriers to animal movement, can not only compromise adaptive potential but also lead to the fixation of deleterious alleles. The south-eastern edge of chimpanzee distribution is particularly fragmented, and conflicting hypotheses have been proposed about population connectivity and viability. To address this uncertainty, we generated both mitochondrial and MiSeq-based microsatellite genotypes for 290 individuals ranging across western Tanzania. While shared mitochondrial haplotypes confirmed historical gene flow, our microsatellite analyses revealed two distinct clusters, suggesting two populations currently isolated from one another. However, we found evidence of high levels of gene flow maintained within each of these clusters, one of which covers an 18,000 km2 ecosystem. Landscape genetic analyses confirmed the presence of barriers to gene flow with rivers and bare habitats highly restricting chimpanzee movement. Our study demonstrates how advances in sequencing technologies, combined with the development of landscape genetics approaches, can resolve ambiguities in the genetic history of critical populations and better inform conservation efforts of endangered species.

Wild animals are under constant threat from a wide range of micro- and macroparasites in their environment. Animals make immune responses against parasites, and these are important in affecting the dynamics of parasite populations. Individual animals vary in their anti-parasite immune responses. Genetic polymorphism of immune-related loci contributes to inter-individual differences in immune responses, but most of what we know in this regard comes from studies of humans or laboratory animals; there are very few such studies of wild animals naturally infected with parasites. Here we have investigated the effect of single nucleotide polymorphisms (SNPs) in immune-related loci (the major histocompatibility complex [MHC], and loci coding for cytokines and Toll-like receptors) on a wide range of immune and infection phenotypes in UK wild house mice, Mus musculus domesticus. We found strong associations between SNPs in various MHC and cytokine-coding loci on both immune measures (antibody concentration and cytokine production) and on infection phenotypes (infection with mites, worms and viruses). Our study provides a comprehensive view of how polymorphism of immune-related loci affects immune and infection phenotypes in naturally infected wild rodent populations.

Gut microbial communities are complex and heterogeneous and play critical roles for animal hosts. Early-life disruptions to microbiome establishment can negatively impact host fitness and development. However, the consequences of such early-life disruptions remain unknown in wild birds. To help fill this gap, we investigated the effect of continuous early-life gut microbiome disruptions on the establishment and development of gut communities in wild Great tit (Parus major) and Blue tit (Cyanistes caeruleus) nestlings by applying antibiotics and probiotics. Treatment neither affected nestling growth nor their gut microbiome composition. Independent of treatment, nestling gut microbiomes of both species grouped by brood, which shared the highest numbers of bacterial taxa with both nest environment and their mother. Although fathers showed different gut communities than their nestlings and nests, they still contributed to structuring chick microbiomes. Lastly, we observed that the distance between nests increased inter-brood microbiome dissimilarity, but only in Great tits, indicating that species-specific foraging behaviour and/or microhabitat influence gut microbiomes. Overall, the strong maternal effect, driven by continuous recolonization from the nest environment and vertical transfer of microbes during feeding, appears to provide resilience towards early-life disruptions in nestling gut microbiomes.

Manipulation of host phenotypes by parasites is hypothesized to be an adaptive strategy enhancing parasite transmission across hosts and generations. Characterizing the molecular mechanisms of manipulation is important to advance our understanding of host-parasite coevolution. The trematode (Levinseniella byrdi) is known to alter the colour and behaviour of its amphipod host (Orchestia grillus) presumably increasing predation of amphipods which enhances trematode transmission through its life cycle. We sampled 24 infected and 24 uninfected amphipods from a salt marsh in Massachusetts to perform differential gene expression analysis. In addition, we constructed novel genomic tools for O. grillus including a de novo genome and transcriptome. We discovered that trematode infection results in upregulation of amphipod transcripts associated with pigmentation and detection of external stimuli, and downregulation of multiple amphipod transcripts implicated in invertebrate immune responses, such as vacuolar ATPase genes. We hypothesize that suppression of immune genes and the altered expression of genes associated with coloration and behaviour may allow the trematode to persist in the amphipod and engage in further biochemical manipulation that promotes transmission. The genomic tools and transcriptomic analyses reported provide new opportunities to discover how parasites alter diverse pathways underlying host phenotypic changes in natural populations.

Crassostrea gigas oysters represent a significant global food source with 4.7 million tons harvested per year. In 2001, the bacterium Vibrio aestuarianus subsp. francensis emerged as a pathogen that causes adult oyster mortality in France and Ireland. Its impact on oyster aquaculture has increased in Europe since its re-emergence in 2012. To better understand the evolutionary mechanisms leading to the emergence and persistence over time of this pathogen, we conducted a survey of mollusc diseases through national reference laboratories across Europe. We analysed 54 new genomes of Vibrio aestuarianus (Va) isolated from multiple environmental compartments since 2001, in areas with and without bivalve mortalities. We used a combination of comparative genomics and population genetics approaches and show that Va has a classical epidemic population structure from which the pathogenic Va francensis subspecies emerged and clonally expanded. Furthermore, we identified a specific cus-cop-containing island conferring copper resistance to Va francensis whose acquisition may have favoured the emergence of pathogenic lineages adapted and specialized to oysters.

Marine animals in the wild are often difficult to access, so they are studied in captivity. However, the implicit assumption that physiological processes of animals in artificial environments are not different from those in the wild has rarely been tested. Here, we investigate the extent to which an animal is impacted by captivity by comparing global gene expression in wild and captive crown-of-thorns starfish (COTS). In a preliminary analysis, we compared transcriptomes of three external tissues obtained from multiple wild COTS with a single captive COTS maintained in aquaria for at least 1 week. On average, an astonishingly large 24% of the coding sequences in the genome were differentially expressed. This led us to conduct a replicated experiment to test more comprehensively the impact of captivity on gene expression. Specifically, a comparison of 13 wild with 8 captive COTS coelomocyte transcriptomes revealed significant differences in the expression of 20% of coding sequences. Coelomocyte transcriptomes in captive COTS remain different from those in wild COTS for more than 30 days and show no indication of reverting back to a wild state (i.e. no evidence of acclimation). Genes upregulated in captivity include those involved in oxidative stress and energy metabolism, whereas genes downregulated are involved in cell signalling. These changes in gene expression indicate that being translocated and maintained in captivity has a marked impact on the physiology and health of these echinoderms. This study suggests that caution should be exercised when extrapolating results from captive aquatic invertebrates to their wild counterparts.

Speciation, the continuous process by which new species form, is often investigated by looking at the variation of nucleotide diversity and differentiation across the genome (hereafter genomic landscapes). A key challenge lies in how to determine the main evolutionary forces at play shaping these patterns. One promising strategy, albeit little used to date, is to comparatively investigate these genomic landscapes as progression through time by using a series of species pairs along a divergence gradient. Here, we resequenced 201 whole-genomes from eight closely related Populus species, with pairs of species at different stages along the divergence gradient to learn more about speciation processes. Using population structure and ancestry analyses, we document extensive introgression between some species pairs, especially those with parapatric distributions. We further investigate genomic landscapes, focusing on within-species (i.e. nucleotide diversity and recombination rate) and among-species (i.e. relative and absolute divergence) summary statistics of diversity and divergence. We observe relatively conserved patterns of genomic divergence across species pairs. Independent of the stage across the divergence gradient, we find support for signatures of linked selection (i.e. the interaction between natural selection and genetic linkage) in shaping these genomic landscapes, along with gene flow and standing genetic variation. We highlight the importance of investigating genomic patterns on multiple species across a divergence gradient and discuss prospects to better understand the evolutionary forces shaping the genomic landscapes of diversity and differentiation.

Climate is a fundamental abiotic factor that plays a key role in driving the evolution, distribution and population diversification of species. However, there have been few investigations of genomic signatures of adaptation to local climatic conditions in cladocerans. Here, we have provided the first high-quality chromosome-level genome assembly (~143 Mb, scaffold N50 12.6 Mb) of the waterflea, Daphnia galeata, and investigated genomic variation in 22 populations from Central Europe and Eastern China. Our ecological-niche models suggested that the historic distribution of D. galeata in Eurasia was significantly affected by Quaternary climate fluctuations. We detected pronounced genomic and morphometric divergences between European and Chinese D. galeata populations. Such divergences could be partly explained by genomic signatures of thermal adaptation to distinct climate regimes: a set of candidate single-nucleotide polymorphisms (SNPs) potentially associated with climate were detected. These SNPs were in genes significantly enriched in the Gene ontology terms "determination of adult lifespan" and "translation repressor activity", and especially, mthl5 and SOD1 involved in the IIS pathway, and EIF4EBP2 involved in the target of the rapamycin signalling pathway. Our study indicates that certain alleles might be associated with particular temperature regimes, playing a functional role in shaping the population structure of D. galeata at a large geographical scale. These results highlight the potential role of molecular variation in the response to climate variation, in the context of global climate change.

Telomeres are well known for their associations with lifespan and ageing across diverse taxa. Early-life telomere length can be influenced by developmental conditions and has been shown positively affect lifetime reproductive success in a limited number of studies. Whether these effects are caused by a change in lifespan, reproductive rate or perhaps most importantly reproductive senescence is unclear. Using long-term data on female breeding success from a threatened songbird (the hihi, Notiomystis cincta), we show that the early-life telomere length of individuals predicts the presence and rate of future senescence of key reproductive traits: clutch size and hatching success. In contrast, senescence of fledging success is not associated with early-life telomere length, which may be due to the added influence of biparental care at this stage. Early-life telomere length does not predict lifespan or lifetime reproductive success in this species. Females may therefore change their reproductive allocation strategy depending on their early developmental conditions, which we hypothesise are reflected in their early-life telomere length. Our results offer new insights on the role that telomeres play in reproductive senescence and individual fitness and suggest telomere length can be used as a predictor for future life history in threatened species.

Since Darwin put forward his opposing hypotheses to explain the successful establishment of species in areas outside their native ranges, the preadaptation and competition-relatedness hypotheses, known as Darwin's naturalization conundrum, numerous studies have sought to understand the relative importance of each. Here, we take advantage of well-characterized beetle communities across laurel forests of the Canary Islands for a first evaluation of the relative support for Darwin's two hypotheses within arthropods. We generated a mitogenome backbone tree comprising nearly half of the beetle genera recorded within the Canary Islands for the phylogenetic placement of native and introduced species sampled in laurel forests, using cytochrome c oxidase I (COI) sequences. For comparative purposes, we also assembled and phylogenetically placed a data set of COI sequences for introduced beetle species that were not sampled within laurel forests. Our results suggest a stronger effect of species preadaptation over resource competition, while also revealing an underappreciated shortfall in arthropod biodiversity data-knowledge of species as being native or introduced. We name this the Humboldtean shortfall and suggest that similar studies using arthropods should incorporate DNA barcode sequencing to mitigate this problem.

Many animals, including insects, exhibit plasticity of body colour in response to environmental changes. Varied expression of carotenoids, major cuticle pigments, significantly contributes to body colour flexibility. However, the molecular mechanisms by which environmental cues regulate carotenoid expression remain largely unknown. In this study, we used the ladybird Harmonia axyridis as a model to investigate the photoperiodic-responsive plasticity of elytra coloration and its endocrine regulation. It was found that H. axyridis females under long-day conditions develop elytra that are much redder than those under short-day conditions, resulting from the differential accumulation of carotenoids. Exogenous hormone application and RNAi-mediated gene knockdown indicate that carotenoid deposition was directed through the juvenile hormone (JH) receptor-mediated canonical pathway. Moreover, we characterized an SR-BI/CD36 (SCRB) gene SCRB10 as the carotenoid transporter responding to JH signalling and regulating the elytra coloration plasticity. Taken together, we propose that JH signalling transcriptionally regulates the carotenoid transporter gene for the photoperiodic coloration plasticity of elytra in the beetles, which reveals a novel role of the endocrine system in the regulation of carotenoid-associated animal body coloration under environmental stimuli.

Divergence in body shape is one of the most widespread and repeated patterns of morphological variation in fishes and is associated with habitat specification and swimming mechanics. Such ecological diversification is the first stage of the explosive adaptive radiation of cichlid fishes in the East African Rift Lakes. We use two hybrid crosses of cichlids (Metriaclima sp. × Aulonocara sp. and Labidochromis sp. × Labeotropheus sp., >975 animals total) to determine the genetic basis of body shape diversification that is similar to benthic-pelagic divergence across fishes. Using a series of both linear and geometric shape measurements, we identified 34 quantitative trait loci (QTL) that underlie various aspects of body shape variation. These QTL are spread throughout the genome, each explaining 3.2-8.6% of phenotypic variation, and are largely modular. Further, QTL are distinct both between these two crosses of Lake Malawi cichlids and compared to previously identified QTL for body shape in fishes such as sticklebacks. We find that body shape is controlled by many genes of small effect. In all, we find that convergent body shape phenotypes commonly observed across fish clades are most likely due to distinct genetic and molecular mechanisms.

Sexually dimorphic behaviour is pervasive across animals, with males and females exhibiting different mate selection, parental care, foraging, dispersal, and territorial strategies. However, the genetic underpinnings of sexually dimorphic behaviours are poorly understood. Here we investigate gene networks and expression patterns associated with sexually dimorphic imprinting-like learning in the butterfly Bicyclus anynana. In this species, both males and females learn visual preferences, but learn preferences for different traits and use different signals as salient, unconditioned cues. To identify genes and gene networks associated with this behaviour, we examined gene expression profiles of the brains and eyes of male and female butterflies immediately post training and compared them to the same tissues of naïve individuals. We found more differentially expressed genes and a greater number of associated gene networks in the eyes, indicating a role of the peripheral nervous system in visual imprinting-like learning. Females had higher chemoreceptor expression levels than males, supporting the hypothesized sexual dimorphic use of chemical cues during the learning process. In addition, genes that influence B. anynana wing patterns (sexual ornaments), such as invected, spalt, and apterous, were also differentially expressed in the brain and eye, suggesting that these genes may influence both sexual ornaments and the preferences for these ornaments. Our results indicate dynamic and sex-specific responses to social scenario in both the peripheral and central nervous systems and highlight the potential role of wing patterning genes in mate preference and learning across the Lepidoptera.

Deleterious variants are selected against but can linger in populations at low frequencies for long periods of time, decreasing fitness and contributing to disease burden in humans and other species. Deleterious variants occur at low frequency but distinguishing deleterious variants from low-frequency neutral variation is challenging based on population genomics data alone. As a result, we have little sense of the number and identity of deleterious variants in wild populations. For haplodiploid species, it has been hypothesised that deleterious alleles will be directly exposed to selection in haploid males, but selection can be masked in diploid females when deleterious variants are recessive, resulting in more efficient purging of deleterious mutations in males. Therefore, comparisons of the differences between haploid and diploid genomes from the same population may be a useful method for inferring rare deleterious variants. This study provides the first formal test of this hypothesis. Using wild populations of Northern paper wasps (Polistes fuscatus), we find that males have fewer missense and nonsense variants per generation than females from the same population. Allele frequency differences are especially pronounced for rare missense and nonsense variants and these differences lead to a lower mutational load in males than females. Based on these data we infer that many highly deleterious mutations are segregating in the paper wasp population. Stronger selection against deleterious alleles in haploid males may have implications for adaptation in other haplodiploid insects and provides evidence that wild populations harbour abundant deleterious variants.

Hybridisation and gene flow can have both deleterious and adaptive consequences for natural populations and species. To better understand the extent of hybridisation in nature and the balance between its beneficial and deleterious outcomes in a changing environment, information on naturally hybridising nonmodel organisms is needed. This requires the characterisation of the structure and extent of natural hybrid zones. Here, we study natural populations of five keystone mound-building wood ant species in the Formica rufa group across Finland. No genomic studies across the species group exist, and the extent of hybridisation and genomic differentiation in sympatry is unknown. Combining genome-wide and morphological data, we demonstrate more extensive hybridisation than was previously detected between all five species in Finland. Specifically, we reveal a mosaic hybrid zone between Formica aquilonia, F. rufa and F. polyctena, comprising further generation hybrid populations. Despite this, we find that F. rufa, F. aquilonia, F. lugubris and F. pratensis form distinct gene pools in Finland. We also find that hybrids occupy warmer microhabitats than the nonadmixed populations of cold-adapted F. aquilonia, and suggest that warm winters and springs, in particular, may benefit hybrids over F. aquilonia, the most abundant F. rufa group species in Finland. In summary, our results indicate that extensive hybridisation may create adaptive potential that could promote wood ant persistence in a changing climate. Additionally, they highlight the potentially significant ecological and evolutionary consequences of extensive mosaic hybrid zones, within which independent hybrid populations face an array of ecological and intrinsic selection pressures.

Individuals vary in their immune genotype, inbreeding coefficient f, immune responses, survival to adulthood, and adult longevity. However, whether immune genes predict survival or longevity, whether such relationships are mediated through immune responses, and how f affects immune genotype remain unclear. We use a wild song sparrow (Melospiza melodia) population in which survival to adulthood, adult longevity, and f were measured precisely, and in which immune responses have previously been assessed. We investigate four toll-like receptor (TLR) and the major histocompatibility complex (MHC) class IIB exon 2 genes. We test whether immune genes predict fitness (survival to adulthood or adult longevity); whether immune genes predict immune response; whether immune response predicts fitness and whether fitness, immune responses, or immune genotypes are correlated with f. We find that survival to adulthood is not associated with immune gene variation, but adult longevity is decreased by high MHC allele diversity (especially in birds that were relatively outbred), and by the presence of a specific MHC supertype. Immune responses were affected by specific immune genotypes. Survival to adulthood and adult longevity were not predicted by immune response, implying caution in the use of immune response as a predictor for fitness. We also found no relationship between f and immune genotype. This finding indicates that immune gene associations with longevity and immune response are not artefacts of f, and suggests that pathogen-mediated selection at functional loci can slow the loss of genetic variation arising from genetic drift and small population size.

Whole-genome duplication has shaped the evolution of angiosperms and other organisms, and is important for many crops. Structural reorganization of chromosomes and repatterning of gene expression are frequently observed in allopolyploids, with physiological and ecological consequences. Recurrent origins from different parental populations are widespread among polyploids, resulting in an array of lineages that provide excellent models to uncover mechanisms of adaptation to divergent environments in early phases of polyploid evolution. We integrate here transcriptomic and ecophysiological comparative studies to show that sibling allopolyploid marsh orchid species (Dactylorhiza, Orchidaceae) occur in different habitats (low nutrient fens vs. meadows with mesic soils) and are characterized by a complex suite of intertwined, pronounced ecophysiological differences between them. We uncover distinct features in leaf elemental chemistry, light-harvesting, photoprotection, nutrient transport and stomata activity of the two sibling allopolyploids, which appear to match their specific ecologies, in particular soil chemistry differences at their native sites. We argue that the phenotypic divergence between the sibling allopolyploids has a clear genetic basis, generating ecological barriers that maintain distinct, independent lineages, despite pervasive interspecific gene flow. This suggests that recurrent origins of polyploids bring about a long-term potential to trigger and maintain functional and ecological diversity in marsh orchids and other groups.

There is considerable evidence for local adaptation in nature, yet important questions remain regarding its genetic basis. How many loci are involved? What are their effect sizes? What is the relative importance of conditional neutrality versus genetic trade-offs? Here we address these questions in the self-pollinating, annual plant Arabidopsis thaliana. We used 400 recombinant inbred lines (RILs) derived from two locally adapted populations in Italy and Sweden, grew the RILs and parents at the parental locations, and mapped quantitative trait loci (QTL) for mean fitness (fruits/seedling planted). We previously published results from the first 3 years of the study, and here add five additional years, providing a unique opportunity to assess how temporal variation in selection might affect QTL detection and classification. We found 10 adaptive and one maladaptive QTL in Italy, and six adaptive and four maladaptive QTL in Sweden. The discovery of maladaptive QTL at both sites suggests that even locally adapted populations are not always at their genotypic optimum. Mean effect sizes for adaptive QTL, 0.97 and 0.55 fruits in Italy and Sweden, respectively, were large relative to the mean fitness of the RILs (approximately 8 fruits/seedling planted at both sites). Both genetic trade-offs (four cases) and conditional neutrality (seven cases) contribute to local adaptation in this system. The 8-year dataset provided greater power to detect QTL and to estimate their locations compared to our previous 3-year study, identifying one new genetic trade-off and resolving one genetic trade-off into two conditionally adaptive QTL.

As the global climate crisis continues, predictions concerning how wild populations will respond to changing climate conditions are informed by an understanding of how populations have responded and/or adapted to climate variables in the past. Changes in the local biotic and abiotic environment can drive differences in phenology, physiology, morphology and demography between populations leading to local adaptation, yet the molecular basis of adaptive evolution in wild non-model organisms is poorly understood. We leverage comparisons between two lineages of Calochortus venustus occurring along parallel transects that allow us to identify loci under selection and measure clinal variation in allele frequencies as evidence of population-specific responses to selection along climatic gradients. We identify targets of selection by distinguishing loci that are outliers to population structure and by using genotype-environment associations across transects to detect loci under selection from each of nine climatic variables. Despite gene flow between individuals of different floral phenotypes and between populations, we find evidence of ecological specialization at the molecular level, including genes associated with key plant functions linked to plant adaptation to California's Mediterranean climate. Single-nucleotide polymorphisms (SNPs) present in both transects show similar trends in allelic similarity across latitudes indicating parallel adaptation to northern climates. Comparisons between eastern and western populations across latitudes indicate divergent genetic evolution between transects, suggesting local adaptation to either coastal or inland habitats. Our study is among the first to show repeated allelic variation across climatic clines in a non-model organism.

Despite being the most important source of liquid freshwater on the planet, groundwater is severely threatened by climate change, agriculture, or industrial mining. It is thus extensively monitored for pollutants and declines in quantity. The organisms living in groundwater, however, are rarely the target of surveillance programmes and little is known about the fauna inhabiting underground habitats. The difficulties accessing groundwater, the lack of expertise, and the apparent scarcity of these organisms challenge sampling and prohibit adequate knowledge on groundwater fauna. Environmental DNA (eDNA) metabarcoding provides an approach to overcome these limitations but is largely unexplored. Here, we sampled water in 20 communal spring catchment boxes used for drinking water provisioning in Switzerland, with a high level of replication at both filtration and amplification steps. We sequenced a portion of the COI mitochondrial gene, which resulted in 4917 ASVs, yet only 3% of the reads could be assigned to a species, genus, or family with more than 90% identity. Careful evaluation of the unassigned reads corroborated that these sequences were true COI sequences belonging mostly to diverse eukaryotic groups, not present in the reference databases. Principal component analyses showed a strong correlation of the community composition with the surface land-use (agriculture vs. forest) and geology (fissured rock vs. unconsolidated sediment). While incomplete reference databases limit the assignment of taxa in groundwater eDNA metabarcoding, we showed that taxonomy-free approaches can reveal large hidden diversity and couple it with major land-use drivers, revealing their imprint on chemical and biological properties of groundwater.

Little is known about when and how planktonic species arise and persist in the open ocean without apparent dispersal barriers. Pteropods are planktonic snails with thin shells susceptible to dissolution that are used as bio-indicators of ocean acidification. However, distinct evolutionary units respond to acidification differently, and defining species boundaries is therefore crucial for predicting the impact of changing ocean conditions. In this global population genomic study of the shelled pteropod Limacina bulimoides, we combined genetic (759,000 single nucleotide polymorphisms) and morphometric data from 161 individuals, revealing three major genetic lineages (FST  = 0.29-0.41): an "Atlantic lineage" sampled across the Atlantic, an "Indo-Pacific lineage" sampled in the North Pacific and Indian Ocean, and a "Pacific lineage" sampled in the North and South Pacific. A time-calibrated phylogeny suggests that the lineages diverged about 1 million years ago, with estimated effective population size remaining high (~10 million) throughout Pleistocene glacial cycles. We do not observe any signatures of recent hybridization, even in areas of sympatry in the North Pacific. While the lineages are reproductively isolated, they are morphologically cryptic, with overlapping shell shape and shell colour distributions. Despite showing that the circumglobal L. bulimoides consists of multiple species with smaller ranges than initially thought, we found that these pteropods still possess high levels of genetic variability. Our study adds to the growing evidence that speciation is often overlooked in the open ocean, and suggests the presence of distinct biological species within many other currently defined circumglobal planktonic species.

The plastic ability for a range of phenotypes to be exhibited by the same genotype allows organisms to respond to environmental variation and may modulate fitness in novel environments. Differing capacities for phenotypic plasticity within a population, apparent as genotype by environment interactions (GxE), can therefore have both ecological and evolutionary implications. Epigenetic gene regulation alters gene function in response to environmental cues without changes to the underlying genetic sequence and likely mediates phenotypic variation. DNA methylation is currently the most well described epigenetic mechanism and is related to transcriptional homeostasis in invertebrates. However, evidence quantitatively linking variation in DNA methylation with that of phenotype is lacking in some taxa, including reef-building corals. In this study, spatial and seasonal environmental variation in Bonaire, Caribbean Netherlands was utilized to assess relationships between physiology and DNA methylation profiles within genetic clones across different genotypes of Acropora cervicornis and A. palmata corals. The physiology of both species was highly influenced by environmental variation compared to the effect of genotype. GxE effects on phenotype were only apparent in A. cervicornis. DNA methylation in both species differed between genotypes and seasons and epigenetic variation was significantly related to coral physiological metrics. Furthermore, plastic shifts in physiology across seasons were significantly positively correlated with shifts in DNA methylation profiles in both species. These results highlight the dynamic influence of environmental conditions and genetic constraints on the physiology of two important Caribbean coral species. Additionally, this study provides quantitative support for the role of epigenetic DNA methylation in mediating phenotypic plasticity in invertebrates.

Most tick-borne pathogens (TBPs) are secondarily acquired by ticks during feeding on infected hosts, which imposes 'priority effect' constraints, as arrival order influences the establishment of new species in a microbial community. Here we tested whether once acquired, TBPs contribute to bacterial microbiota functioning by increasing community stability. For this, we used Hyalomma marginatum and Rhipicephalus bursa ticks collected from cattle in different locations of Corsica and combined 16S rRNA amplicon sequencing and co-occurrence network analysis, with high-throughput pathogen detection, and in silico removal of nodes to test for impact of rickettsial pathogens on network properties. Despite its low centrality, Rickettsia showed preferential connections in the networks, notably with a keystone taxon in H. marginatum, suggesting facilitation of Rickettsia colonisation by the keystone taxon. In addition, conserved patterns of community assembly in both tick species were affected by Rickettsia removal, suggesting that privileged connections of Rickettsia in the networks make this taxon a driver of community assembly. However, Rickettsia removal had minor impact on the conserved 'core bacterial microbiota' of H. marginatum and R. bursa. Interestingly, networks of the two tick species with Rickettsia have similar node centrality distribution, a property that is lost after Rickettsia removal, suggesting that this taxon drives specific hierarchical interactions between bacterial microbes in the microbiota. The study indicates that tick-borne Rickettsia play a significant role in the tick bacterial microbiota, despite their low centrality. These bacteria are influential and contribute to the conservation of the 'core bacterial microbiota' while also promoting community stability.

Maternal effect senescence, a decline in offspring viability with maternal age, has been documented across diverse animals, but its mechanisms remain largely unknown. Here, we test maternal effect senescence and explore its possible molecular mechanisms in a fish. We compared the levels of maternal mRNA transcripts of DNA repair genes and mtDNA copies in eggs and the levels of DNA damage in somatic and germline tissues between young and old female sticklebacks. We also tested, in an in vitro fertilization experiment, whether maternal age and sperm DNA damage level interactively influence the expression of DNA repair genes in early embryos. Old females transferred less mRNA transcripts of DNA repair genes into their eggs than did young females, but maternal age did not influence egg mtDNA density. Despite a higher level of oxidative DNA damage in the skeletal muscle, old females had a similar level of damage in the gonad to young females, suggesting the prioritization for germline maintenance during ageing. The embryos of both old and young mothers increased the expression of DNA repair genes in response to an increased level of oxidative DNA damage in sperm used for their fertilization. The offspring of old mothers showed higher rates of hatching, morphological deformity and post-hatching mortality and had smaller body size at maturity. These results suggest that maternal effect senescence may be mediated by reduced capacity of eggs to detect and repair DNA damages, especially prior to the embryonic genomic activation.

Rising mean and variance in temperatures elevates threats to endangered freshwater species such as lake sturgeon, Acipenser fulvescens. Previous research demonstrated that higher temperatures during development result in physiological consequences for lake sturgeon populations throughout Manitoba, Canada, with alteration of metabolic rate, thermal tolerance, transcriptional responses, growth and mortality. We acclimated lake sturgeon (30-60 days post fertilization, a period of high mortality) from northern and southern populations (56°02'46.5″N, 96°54'18.6″W and 50°17'52″N, 95°32'51″W, respectively, separated by approximately 650 km) within Manitoba to current (summer highs of 20-23°C) and future projected (+2-3°C) environmental temperatures of 16, 20 and 24°C for 30 days, and we measured gill transcriptional responses using RNAseq. Transcripts revealed SNPs consistent with genetically distinct populations and transcriptional responses altered by acclimation temperature. There were a higher number of differentially expressed transcripts observed in the southern, compared to the northern, population as temperatures increased, indicating enhanced transcriptional plasticity. Both lake sturgeon populations responded to elevated acclimation temperatures by downregulating the transcription of genes involved in protein synthesis and energy production. Furthermore, there were population-specific thresholds for the downregulation of processes promoting transcriptional plasticity as well as mitochondrial function as the northern population showed decreases at 20°C, while this capacity was not diminished until 24°C in the southern population. These transcriptional responses highlight the molecular impacts of increasing temperatures for divergent lake sturgeon populations during vulnerable developmental periods and the critical influence of transcriptome plasticity on acclimation capacity.

During the Late Pleistocene, major parts of North America were periodically covered by ice sheets. However, there are still questions about whether ice-free refugia were present in the Alexander Archipelago along the Southeast (SE) Alaska coast during the last glacial maximum (LGM). Numerous subfossils have been recovered from caves in SE Alaska, including American black (Ursus americanus) and brown (U. arctos) bears, which today are found in the Alexander Archipelago but are genetically distinct from mainland bear populations. Hence, these bear species offer an ideal system to investigate long-term occupation, potential refugial survival and lineage turnover. Here, we present genetic analyses based on 99 new complete mitochondrial genomes from ancient and modern brown and black bears spanning the last ~45,000 years. Black bears form two SE Alaskan subclades, one preglacial and another postglacial, that diverged >100,000 years ago. All postglacial ancient brown bears are closely related to modern brown bears in the archipelago, while a single preglacial brown bear is found in a distantly related clade. A hiatus in the bear subfossil record around the LGM and the deep split of their pre- and postglacial subclades fail to support a hypothesis of continuous occupancy in SE Alaska throughout the LGM for either species. Our results are consistent with an absence of refugia along the SE Alaska coast, but indicate that vegetation quickly expanded after deglaciation, allowing bears to recolonize the area after a short-lived LGM peak.

Suboptimal conditions during development can shorten telomeres, the protective DNA caps on the end of chromosomes. Shorter early-life telomere length (TL) can indicate reduced somatic maintenance, leading to lower survival and shorter lifespan. However, despite some clear evidence, not all studies show a relationship between early-life TL and survival or lifespan, which may be due to differences in biology or study design (e.g., survival period measured). In superb fairy-wrens (Malurus cyaneus), we assessed whether early-life TL predicts mortality across different life-history stages (fledgling, juvenile, adult). However, in contrast to a similar study on a congener, early-life TL did not predict mortality across any life stage in this species. We then performed a meta-analysis including 32 effect sizes from 23 studies (15 birds and three mammals) to quantify the effect of early-life TL on mortality whilst taking into consideration potential sources of biological and methodological variation. Overall, the effect of early-life TL on mortality was significant, corresponding to a 15% reduction in mortality risk with each standard deviation increase in TL. However, the effect became weaker when correcting for publication bias. Contrary to our predictions, there was no evidence that effects of early-life TL on mortality varied with species lifespan or the period over which survival was measured. However, negative effects of early-life TL on mortality risk were pervasive throughout life. These results imply that effects of early-life TL on mortality are more likely to be context-dependent than age-dependent, although substantial power and publication bias issues highlight the need for more research.

Male reproductive competition can select for condition-dependent, conspicuous traits that signal some aspect of fighting ability and facilitate assessment of potential rivals. However, the underlying mechanisms that link the signal to a male's current condition are difficult to investigate in wild populations, often requiring invasive experimental manipulation. Here, we use digital photographs and chest skin samples to investigate the mechanisms of a visual signal used in male competition in a wild primate, the red chest patch in geladas (Theropithecus gelada). We analysed photographs collected during natural (n = 144) and anaesthetized conditions (n = 38) to understand variability in male and female chest redness, and we used chest skin biopsies (n = 38) to explore sex differences in gene expression. Male and female geladas showed similar average redness, but males exhibited a wider within-individual range in redness under natural conditions. These sex differences were also reflected at the molecular level, with 10.5% of genes exhibiting significant sex differences in expression. Subadult males exhibited intermediate gene expression patterns between adult males and females, pointing to mechanisms underlying the development of the red chest patch. We found that genes more highly expressed in males were associated with blood vessel development and maintenance but not with androgen or oestrogen activity. Together, our results suggest male gelada redness variability is driven by increased blood vessel branching in the chest skin, providing a potential link between male chest redness and current condition as increased blood circulation to exposed skin could lead to heat loss in the cold, high-altitude environment of geladas.

Parasites with complex life cycles are known to induce phenotypic changes in their intermediate hosts to increase transmission to the final host. The magnitude of these changes could increase with the number of parasites, which would be beneficial to co-infecting parasites. Yet, adverse effects of high parasite load (i.e. many parasites in a single host) might stress both hosts and parasites (e.g. through an increased immune response). We investigated the consequences of parasite load on the transcriptional activity and morphology of the cestode Anomotaenia brevis and its intermediate host, the ant Temnothorax nylanderi. We demonstrated that many differentially expressed host genes shifted with parasite load, and their functions indicate a stronger immune response and fight against oxidative stress in heavily infected hosts. The expression of other host genes responded to infection in an all-or-nothing manner, as did the morphology of the host workers. However, the cestodes became smaller when they competed with other parasites for resources from a single host. Their expression profile further indicated shifts in host immune avoidance, starvation resistance and vesicle-mediated transport. In summary, our study reveals clear consequences of parasite load and highlights specific processes and traits affected by this.

Increasing antimicrobial resistance (AMR) poses a challenge for treatment of bacterial diseases. In real life, bacterial infections are typically embedded within complex multispecies communities and influenced by the environment, which can shape costs and benefits of AMR. However, knowledge of such interactions and their implications for AMR in vivo is limited. To address this knowledge gap, we investigated fitness-related traits of a pathogenic bacterium (Flavobacterium columnare) in its fish host, capturing the effects of bacterial antibiotic resistance, coinfections between bacterial strains and metazoan parasites (fluke Diplostomum pseudospathaceum) and antibiotic exposure. We quantified real-time replication and virulence of sensitive and resistant bacteria and demonstrate that both bacteria can benefit from coinfection in terms of persistence and replication, depending on the coinfecting partner and antibiotic presence. We also show that antibiotics can benefit resistant bacteria by increasing bacterial replication under coinfection with flukes. These results emphasize the importance of diverse, inter-kingdom coinfection interactions and antibiotic exposure in shaping costs and benefits of AMR, supporting their role as significant contributors to spread and long-term persistence of resistance.

In plants where seed dispersal is limited compared with pollen dispersal, hybridisation may enhance gene exchange and species dispersal. We provide genetic evidence of hybridisation contributing to the expansion of the rare Eucalyptus risdonii into the range of the widespread Eucalyptus amygdalina. These closely related tree species are morphologically distinct, and observations suggest that natural hybrids occur along their distribution boundaries and as isolated trees or in small patches within the range of E. amygdalina. Hybrid phenotypes occur outside the range of normal dispersal for E. risdonii seed, yet in some hybrid patches small individuals resembling E. risdonii occur and are hypothesised to be a result of backcrossing. Using 3362 genome-wide SNPs assessed from 97 individuals of E. risdonii and E. amygdalina and 171 hybrid trees, we show that (i) isolated hybrids match the genotypes expected of F1 /F2 hybrids, (ii) there is a continuum in the genetic composition among the isolated hybrid patches from patches dominated by F1 /F2 -like genotypes to those dominated by E. risdonii-backcross genotypes, and (iii) the E. risdonii-like phenotypes in the isolated hybrid patches are most-closely related to proximal larger hybrids. These results suggest that the E. risdonii phenotype has been resurrected in isolated hybrid patches established from pollen dispersal, providing the first steps in its invasion of suitable habitat by long-distance pollen dispersal and complete introgressive displacement of E. amygdalina. Such expansion accords with the population demographics, common garden performance data, and climate modelling which favours E. risdonii and highlights a role of interspecific hybridisation in climate change adaptation and species expansion.

Traditional agrosystems, where humans, crops and microbes have coevolved over long periods, can serve as models to understand the ecoevolutionary determinants of disease dynamics and help the engineering of durably resistant agrosystems. Here, we investigated the genetic and phenotypic relationship between rice (Oryza sativa) landraces and their rice blast pathogen (Pyricularia oryzae) in the traditional Yuanyang terraces of flooded rice paddies in China, where rice landraces have been grown and bred over centuries without significant disease outbreaks. Analyses of genetic subdivision revealed that indica rice plants clustered according to landrace names. Three new diverse lineages of rice blast specific to the Yuanyang terraces coexisted with lineages previously detected at the worldwide scale. Population subdivision in the pathogen population did not mirror pattern of population subdivision in the host. Measuring the pathogenicity of rice blast isolates on landraces revealed generalist life history traits. Our results suggest that the implementation of disease control strategies based on the emergence or maintenance of a generalist lifestyle in pathogens may sustainably reduce the burden of disease in crops.

Why species that in their core areas mainly reproduce sexually become enriched with clones in marginal populations ("geographic parthenogenesis") remains unclear. Earlier hypotheses have emphasized that selection might promote clonality because it protects locally adapted genotypes. On the other hand, it also hampers recombination and adaptation to changing conditions. The aim of the present study was to investigate the early stages of range expansion in a partially clonal species and what drives an increase in cloning during such expansion. We used genome-wide sequencing to investigate the origin and evolution of large clones formed in a macroalgal species (Fucus vesiculosus) during a recent expansion into the postglacial Baltic Sea. We found low but persistent clonality in core populations, while at range margins, large dominant clonal lineages had evolved repeatedly from different sexual populations. A range expansion model showed that even when asexual recruitment is less favourable than sexual recruitment in core populations, repeated bottlenecks at the expansion front can establish a genetically eroded clonal wave that spreads ahead of a sexual wave into the new area. Genetic variation decreases by drift following repeated bottlenecks at the expansion front. This results in the emerging clones having low expected heterozygosity, which corroborated our empirical observations. We conclude that Baker's Law (clones being favoured by uniparental reproductive assurance in new areas) can play an important role during range expansion in partially clonal species, resulting in a complex spatiotemporal mosaic of clonal and sexual lineages that might persist during thousands of generations.

Bottlenose dolphins (Tursiops spp.) are found in waters around Australia, with T. truncatus typically occupying deeper, more oceanic habitat, while T. aduncus occur in shallower, coastal waters. Little is known about the colonization history of T. aduncus along the Western Australian coastline; however, it has been hypothesized that extant populations are the result of an expansion along the coastline originating from a source in the north of Australia. To investigate the history of coastal T. aduncus populations in the area, we generated a genomic SNP dataset using a double-digest restriction-site-associated DNA (ddRAD) sequencing approach. The resulting dataset consisted of 103,201 biallelic SNPs for 112 individuals which were sampled from eleven coastal and two offshore sites between Shark Bay and Cygnet Bay, Western Australia. Our population genomic analyses showed a pattern consistent with the proposed source in the north with significant isolation by distance along the coastline, as well as a reduction in genomic diversity measures along the coastline with Shark Bay showing the most pronounced reduction. Our demographic analysis indicated that the expansion of T. aduncus along the coastline began around the last glacial maximum and progressed southwards with the Shark Bay population being founded only 13 kya. Our results are in line with coastal colonization histories inferred for Tursiops globally, highlighting the ability of delphinids to rapidly colonize novel coastal niches as habitat is released during glacial cycle-related global sea level and temperature changes.

Xylosandrus crassiusculus, a fungus-farming wood borer native to Southeastern Asia, is the most rapidly spreading invasive ambrosia species worldwide. Previous studies focusing on its genetic structure suggested the existence of cryptic genetic variation in this species. Yet, these studies used different genetic markers, focused on different geographical areas and did not include Europe. Our first goal was to determine the worldwide genetic structure of this species based on both mitochondrial and genomic markers. Our second goal was to study X. crassiusculus' invasion history on a global level and identify the origins of the invasion in Europe. We used a COI and RAD sequencing design to characterize 188 and 206 specimens worldwide, building the most comprehensive genetic data set for any ambrosia beetle to date. The results were largely consistent between markers. Two differentiated genetic clusters were invasive, albeit in different regions of the world. The markers were inconsistent only for a few specimens found exclusively in Japan. Mainland USA could have acted as a source for further expansion to Canada and Argentina through stepping stone expansion and bridgehead events. We showed that Europe was only colonized by Cluster II through a complex invasion history including several arrivals from multiple origins in the native area, and possibly including bridgehead from the United States. Our results also suggested that Spain was colonized directly from Italy through intracontinental dispersion. It is unclear whether the mutually exclusive allopatric distribution of the two clusters is due to neutral effects or due to different ecological requirements.

Individuals differ in the nature of the immune responses they produce, affecting disease susceptibility and ultimately health and fitness. These differences have been hypothesized to have an origin in events experienced early in life that then affect trajectories of immune development and responsiveness. Here, we investigate how early-life immune expression profiles influence life history outcomes in a natural population of field voles, Microtus agrestis, in which we are able to monitor variation between and within individuals through time by repeat sampling of individually marked animals. We analysed the co-expression of 20 immune genes in early life to create a correlation network consisting of three main clusters, one of which (containing Gata3, Il10 and Il17) was associated with later-life reproductive success and susceptibility to chronic bacterial (Bartonella) infection. More detailed analyses supported associations between early-life expression of Il17 and reproductive success later in life, and of Il10 expression early in life and later infection with Bartonella. We also found significant association between an Il17 genotype and the early-life expression of Il10. Our results demonstrate that immune expression profiles can be manifested during early life with effects that persist through adulthood and that shape the variability among individuals in susceptibility to infection and fitness widely seen in natural populations.

Although animal dispersal is known to play key roles in ecological and evolutionary processes such as colonization, population extinction and local adaptation, little is known about its genetic basis, particularly in vertebrates. Untapping the genetic basis of dispersal should deepen our understanding of how dispersal behaviour evolves, the molecular mechanisms that regulate it and link it to other phenotypic aspects in order to form the so-called dispersal syndromes. Here, we comprehensively combined quantitative genetics, genome-wide sequencing and transcriptome sequencing to investigate the genetic basis of natal dispersal in a known ecological and evolutionary model of vertebrate dispersal: the common lizard, Zootoca vivipara. Our study supports the heritability of dispersal in semi-natural populations, with less variation attributable to maternal and natal environment effects. In addition, we found an association between natal dispersal and both variation in the carbonic anhydrase (CA10) gene, and in the expression of several genes (TGFB2, SLC6A4, NOS1) involved in central nervous system functioning. These findings suggest that neurotransmitters (serotonin and nitric oxide) are involved in the regulation of dispersal and shaping dispersal syndromes. Several genes from the circadian clock (CRY2, KCTD21) were also differentially expressed between disperser and resident lizards, supporting that the circadian rhythm, known to be involved in long-distance migration in other taxa, might affect dispersal as well. Since neuronal and circadian pathways are relatively well conserved across vertebrates, our results are likely to be generalisable, and we therefore encourage future studies to further investigate the role of these pathways in shaping dispersal in vertebrates.

Understanding the immunogenetic basis of coronavirus (CoV) susceptibility in major pathogen reservoirs, such as bats, is central to inferring their zoonotic potential. Members of the cryptic Hipposideros bat species complex differ in CoV susceptibility, but the underlying mechanisms remain unclear. The genes of the major histocompatibility complex (MHC) are the best understood genetic basis of pathogen resistance, and differences in MHC diversity are one possible reason for asymmetrical infection patterns among closely related species. Here, we aimed to link asymmetries in observed CoV (CoV-229E, CoV-2B and CoV-2Bbasal) susceptibility to immunogenetic differences amongst four Hipposideros bat species. From the 2072 bats assigned to their respective species using the mtDNA cytochrome b gene, members of the most numerous and ubiquitous species, Hipposideros caffer D, were most infected with CoV-229E and SARS-related CoV-2B. Using a subset of 569 bats, we determined that much of the existent allelic and functional (i.e. supertype) MHC DRB class II diversity originated from common ancestry. One MHC supertype shared amongst all species, ST12, was consistently linked to susceptibility with CoV-229E, which is closely related to the common cold agent HCoV-229E, and infected bats and those carrying ST12 had a lower body condition. The same MHC supertype was connected to resistance to CoV-2B, and bats with ST12 were less likely be co-infected with CoV-229E and CoV-2B. Our work suggests a role of immunogenetics in determining CoV susceptibility in bats. We advocate for the preservation of functional genetic and species diversity in reservoirs as a means of mitigating the risk of disease spillover.

Parasite-mediated selection can rapidly drive up resistance levels in host populations, but fixation of resistance traits may be prevented by costs of resistance. Black bean aphids (Aphis fabae) benefit from increased resistance to parasitoids when carrying the defensive bacterial endosymbiont Hamiltonella defensa. However, due to fitness costs that come with symbiont infection, symbiont-conferred resistance may result in either a net benefit or a net cost to the aphid host, depending on parasitoid presence as well as on the general ecological context. Balancing selection may therefore explain why in natural aphid populations, H. defensa is often found at intermediate frequencies. Here we present a 2-year field study where we set out to look for signatures of balancing selection in natural aphid populations. We collected temporally well-resolved data on the prevalence of H. defensa in A. f. fabae and estimated the risk imposed by parasitoids using sentinel hosts. Despite a marked and consistent early-summer peak in parasitism risk, and significant changes in symbiont prevalence over time, we found just a weak correlation between parasitism risk and H. defensa frequency dynamics. H. defensa prevalence in the populations under study was, in fact, better explained by the number of heat days that previous aphid generations were exposed to. Our study grants an unprecedentedly well-resolved insight into the dynamics of endosymbiont and parasitoid communities of A. f. fabae populations, and it adds to a growing body of empirical evidence suggesting that not only parasitism risk, but rather multifarious selection is shaping H. defensa prevalence in the wild.

Understanding the drivers of evolutionary innovation provides a crucial perspective of how evolutionary processes unfold across taxa and ecological systems. It has been hypothesised that the Southern Ocean provided ecological opportunities for novelty in the past. However, the drivers of innovation are challenging to pinpoint as the evolutionary genetics of Southern Ocean fauna are influenced by Quaternary glacial-interglacial cycles, oceanic currents and species ecology. Here we examined the genome-wide single nucleotide polymorphisms of the Southern Ocean brittle stars Ophionotus victoriae (five arms, broadcaster) and O. hexactis (six arms, brooder). We found that O. victoriae and O. hexactis are closely-related species with interspecific gene flow. During the late Pleistocene, O. victoriae likely persisted in a connected deep water refugium and in situ refugia on the Antarctic continental shelf and around Antarctic islands; O. hexactis persisted exclusively within in situ island refugia. Within O. victoriae, contemporary gene flow linking to the Antarctic Circumpolar Current, regional gyres and other local oceanographic regimes was observed. Gene flow connecting West and East Antarctic islands near the Polar Front was also detected in O. hexactis. A strong association was detected between outlier loci and salinity in O. hexactis. Both O. victoriae and O. hexactis are associated with genome-wide increase in alleles at intermediate-frequencies; the alleles associated with this peak appear to be species specific, and these intermediate-frequency variants are far more excessive in O. hexactis. We hypothesise that the peak in alleles at intermediate frequencies could be related to adaptation in the recent past, linked to evolutionary innovations of increase in arm number and a switch to brooding from broadcasting, in O. hexactis.

A fundamental question in holobiont biology is the extent to which microbiomes are determined by host characteristics regulated by their genotype. Studies on the interactions of host genotype and microbiomes are emerging but disentangling the role that host genotype has in shaping microbiomes remains challenging in natural settings. Host genotypes tend to be segregated in space and affected by different environments. Here we overcome this challenge by studying an unusual situation where host asexual (5 clonal lineages) and sexual genotypes (15 non-clonal lineages) of the same species co-occur under the same environment. This allowed us to partition the influence of morphological traits and genotype in shaping host-associated bacterial communities. Lamina-associated bacteria of co-occurring kelp sexual non-clonal (Ecklonia radiata) and asexual clonal (E. brevipes) morphs were compared to test whether host genotype influences microbiomes beyond morphology. Similarity of bacterial composition and predicted functions were evaluated among individuals within a single clonal genotype or among non-clonal genotypes of each morph. Higher similarity in bacterial composition and inferred functions were found among identical clones of E. brevipes compared to other clonal genotypes or unique non-clonal E. radiata genotypes. Additionally, bacterial diversity and composition differed significantly between the two morphs and were related with one morphological trait in E. brevipes (haptera). Thus, factors regulated by the host genotype (e.g. secondary metabolite production) likely drive differences in microbial communities between morphs. The strong association of genotype and microbiome found here highlights the importance of genetic relatedness of hosts in determining variability in their bacterial symbionts.

Root-associated fungi could play a role in determining both the positive relationship between plant diversity and productivity in experimental grasslands, and its strengthening over time. This hypothesis assumes that specialized pathogenic and mutualistic fungal communities gradually assemble over time, enhancing plant growth more in species-rich than in species-poor plots. To test this hypothesis, we used high-throughput amplicon sequencing to characterize root-associated fungal communities in experimental grasslands of 1 and 15 years of age with varying levels of plant species richness. Specifically, we tested whether the relationship between fungal communities and plant richness and productivity becomes stronger with the age of the experimental plots. Our results showed that fungal diversity increased with plant diversity, but this relationship weakened rather than strengthened over the two time points. Contrastingly, fungal community composition showed increasing associations with plant diversity over time, suggesting a gradual build-up of specific fungal assemblages. Analyses of different fungal guilds showed that these changes were particularly marked in pathogenic fungi, whose shifts in relative abundance are consistent with the pathogen dilution hypothesis in diverse plant communities. Our results suggest that root-associated fungal pathogens play more specific roles in determining the diversity-productivity relationship than other root-associated plant symbionts.

Evolutionary genetics has long struggled with understanding how functional genes under selection remain polymorphic in natural populations. Taking as a starting point that natural selection is ultimately a manifestation of ecological processes, we spotlight an underemphasized and potentially ubiquitous ecological effect that may have fundamental effects on the maintenance of genetic variation. Negative frequency dependency is a well-established emergent property of density dependence in ecology, because the relative profitability of different modes of exploiting or utilizing limiting resources tends to be inversely proportional to their frequency in a population. We suggest that this may often generate negative frequency-dependent selection (NFDS) on major effect loci that affect rate-dependent physiological processes, such as metabolic rate, that are phenotypically manifested as polymorphism in pace-of-life syndromes. When such a locus under NFDS shows stable intermediate frequency polymorphism, this should generate epistatic selection potentially involving large numbers of loci with more minor effects on life-history (LH) traits. When alternative alleles at such loci show sign epistasis with a major effect locus, this associative NFDS will promote the maintenance of polygenic variation in LH genes. We provide examples of the kind of major effect loci that could be involved and suggest empirical avenues that may better inform us on the importance and reach of this process.

Biological introductions provide a natural ecological experiment unfolding in a recent historical timeframe to elucidate how evolutionary processes (such as founder effects, genetic diversity and adaptation) shape the genomic landscape of populations postintroduction. The Asian icefish, Protosalanx chinensis, is an economically important fishery resource, deliberately introduced into dozens of provinces across China for decades. However, while invading and disturbing the local ecosystem, many introduced populations declined, disappearing mysteriously in a very short time. The way in which various evolutionary forces integrate to result in invasion failure of an introduced population remains unknown. Here, we performed whole-genome sequencing of 10 species from the Salangidae family and 70 Asian icefish (Protosalanx chinensis) individuals from 7 geographic populations in China, aiming to characterize the evolutionary fate of introduced populations. Our results show that compared to other Salangidae species, P. chinensis has low genetic diversity, potentially due to the long-lasting decline in population size. In a recently introducted population, Lugu lake, severe sampling effects and a strong bottleneck further deteriorated the genomic landscape. Although the introduced population showed signs of reduced genetic load, the purging selection efficiency was low. Our selective sweep analysis revealed site frequency changes in candidate genes, including gata1a and hoxd4b, which could be associated with a decrease in dissolved oxygen in the deep-water plateau lake. These findings caution against the widespread introduction of P. chinensis in China and lay the groundwork for future use of this economically species.

Early lineage diversification is central to understand what mutational events drive species divergence. Particularly, gene misregulation in interspecific hybrids can inform about what genes and pathways underlie hybrid dysfunction. In Drosophila hybrids, how regulatory evolution impacts different reproductive tissues remains understudied. Here, we generate a new genome assembly and annotation in Drosophila willistoni and analyse the patterns of transcriptome divergence between two allopatrically evolved D. willistoni subspecies, their male sterile and female fertile hybrid progeny across testis, male accessory gland, and ovary. Patterns of transcriptome divergence and modes of regulatory evolution were tissue-specific. Despite no indication for cell-type differences in hybrid testis, this tissue exhibited the largest magnitude of expression differentiation between subspecies and between parentals and hybrids. No evidence for anomalous dosage compensation in hybrid male tissues was detected nor was a differential role for the neo- and the ancestral arms of the D. willistoni X chromosome. Compared to the autosomes, the X chromosome appeared enriched for transgressively expressed genes in testis despite being the least differentiated in expression between subspecies. Evidence for fine genome clustering of transgressively expressed genes suggests a role of chromatin structure on hybrid gene misregulation. Lastly, transgressively expressed genes in the testis of the sterile male progeny were enriched for GO terms not typically associated with sperm function, instead hinting at anomalous development of the reproductive tissue. Our thorough tissue-level portrait of transcriptome differentiation between recently diverged D. willistoni subspecies and their hybrids provides a more nuanced view of early regulatory changes during speciation.

The lowland Amazon rainforest houses some of the greatest tree diversity on Earth. While the vast majority of these species are rare, a small number are common and widespread and thus considered to play a disproportionate role in many of the global ecosystem services provided by the Amazon. However, the extent to which dominant Amazonian tree species actually include multiple clades, each on their own unique evolutionary trajectory, is unknown. Here we investigate the extent to which lineage divergence may be occurring within Protium subserratum (Burseraceae), a common and widespread tree species that is monophyletic with populations exhibiting genotypic and phenotypic differences associated with soil and geography. Utilizing a combination of phylogenomic and population genomic methods with sampling from across the range, we found that P. subserratum contains at least eight distinct clades. Specialization onto white-sand soils has evolved independently at least twice within the species; however, phenotype is not correlated with soil type. Finally, cryptic diversity at the base of the Andes is associated with elevational shifts. Together these results lend support to the hypothesis that common and widespread Amazon tree species may not represent evolutionary cohesive units. Instead, these dominant species may more commonly represent species complexes, undergoing evolutionary transitions on a trajectory to become multiple range-restricted, specialist species.

Epigenetic modifications, like DNA methylation, generate phenotypic diversity in fish and ultimately lead to adaptive evolutionary processes. Euryhaline marine species that migrate between salinity-contrasted habitats have received little attention regarding the role of salinity on whole-genome DNA methylation. Investigation of salinity-induced DNA methylation in fish will help to better understand the potential role of this process in salinity acclimation. Using whole-genome bisulfite sequencing, we compared DNA methylation patterns in European sea bass (Dicentrarchus labrax) juveniles in seawater and after freshwater transfer. We targeted the gill as a crucial organ involved in plastic responses to environmental changes. To investigate the function of DNA methylation in gills, we performed RNAseq and assessed DNA methylome-transcriptome correlations. We showed a negative correlation between gene expression levels and DNA methylation levels in promoters, first introns and first exons. A significant effect of salinity on DNA methylation dynamics with an overall DNA hypomethylation in freshwater-transferred fish compared to seawater controls was demonstrated. This suggests a role of DNA methylation changes in salinity acclimation. Genes involved in key functions as metabolism, ion transport and transepithelial permeability (junctional complexes) were differentially methylated and expressed between salinity conditions. Expression of genes involved in mitochondrial metabolism (tricarboxylic acid cycle) was increased, whereas the expression of DNA methyltransferases 3a was repressed. This study reveals novel links between DNA methylation, mainly in promoters and first exons/introns, and gene expression patterns following salinity change.

Numerous studies have investigated the spatiotemporal variability in water microbial communities, yet the effects of relic DNA on microbial community profiles, especially microeukaryotes, remain far from fully understood. Here, total and active bacterial and microeukaryotic community compositions were characterized using propidium monoazide (PMA) treatment coupled with high-throughput sequencing in a river-reservoir ecosystem. Beta diversity analysis showed a significant difference in community composition between both the PMA untreated and treated bacteria and microeukaryotes; however, the differentiating effect was much stronger for microeukaryotes. Relic DNA only resulted in underestimation of the relative abundances of Bacteroidota and Nitrospirota, while other bacterial taxa exhibited no significant changes. As for microeukaryotes, the relative abundances of some phytoplankton (e.g. Chlorophyta, Dinoflagellata and Ochrophyta) and fungi were greater after relic DNA removal, whereas Cercozoa and Ciliophora showed the opposite trend. Moreover, relic DNA removal weakened the size and complexity of cross-trophic microbial networks and significantly changed the relationships between environmental factors and microeukaryotic community composition. However, there was no significant difference in the rates of temporal community turnover between the PMA untreated and treated samples for either bacteria or microeukaryotes. Overall, our results imply that the presence of relic DNA in waters can give misleading information of the active microbial community composition, co-occurrence networks and their relationships with environmental conditions. More studies of the abundance, decay rate and functioning of nonviable DNA in freshwater ecosystems are highly recommended in the future.

Multisensory integration (MSI) combines information from more than one sensory modality to elicit behaviours distinct from unisensory behaviours. MSI is best understood in animals with complex brains and specialized centres for parsing different modes of sensory information, but dispersive larvae of sessile marine invertebrates utilize multimodal environmental sensory stimuli to base irreversible settlement decisions on, and most lack complex brains. Here, we examined the sensory determinants of settlement in actinula larvae of the hydrozoan Ectopleura crocea (Cnidaria), which possess a diffuse nerve net. A factorial settlement study revealed that photo-, chemo- and mechanosensory cues each influenced the settlement response in a complex and hierarchical manner that was dependent on specific combinations of cues, an indication of MSI. Additionally, sensory gene expression over development peaked with developmental competence to settle, which in actinulae, requires cnidocyte discharge. Transcriptome analyses also highlighted several deep homological links between cnidarian and bilaterian mechano-, chemo-, and photosensory pathways. Fluorescent in situ hybridization studies of candidate transcripts suggested cellular partitioning of sensory function among the few cell types that comprise the actinula nervous system, where ubiquitous polymodal sensory neurons expressing putative chemo- and photosensitivity interface with mechanoreceptive cnidocytes. We propose a simple multisensory processing circuit, involving polymodal chemo/photosensory neurons and mechanoreceptive cnidocytes, that is sufficient to explain MSI in actinulae settlement. Our study demonstrates that MSI is not exclusive to complex brains, but likely predated and contextualized their evolution.

Domesticated honeybees and wild bees are some of the most important beneficial insects for human and environmental health, but infectious diseases pose a serious risk to these pollinators, particularly following the emergence of the ectoparasitic mite Varroa destructor as a viral vector. The acquisition of this novel viral vector from the Asian honeybee Apis ceranae has fundamentally changed viral epidemiology in its new host, the western honeybee A. mellifera. While the recently discovered Lake Sinai Viruses (LSV) have been associated with weak honeybee colonies, they have not been associated with vector-borne transmission. By combining a large-scale multi-year survey of LSV in Chinese A. mellifera and A. cerana honeybee colonies with globally available LSV-sequence data, we investigate the global epidemiology of this virus. We find that globally distributed LSV is a highly diverse multi-strain virus, which is predominantly associated with the western honeybee A. mellifera. In contrast to the vector-borne deformed wing virus, LSV is not an emerging disease. Instead, demographic reconstruction and strong global and local population structure indicates that it is a highly variable multi-strain virus in a stable association with its main host, the western honeybee. Prevalence patterns in China suggest a potential role for migratory beekeeping in the spread of this pathogen, demonstrating the potential for disease transmission with the man-made transport of beneficial insects.

Differences in gut microbiome composition are linked with health, disease and ultimately host fitness; however, the molecular mechanisms underlying that relationship are not well characterized. Here, we modified the fish gut microbiota using antibiotic and probiotic feed treatments to address the effect of host microbiome on gene expression patterns. Chinook salmon (Oncorhynchus tshawytscha) gut gene expression was evaluated using whole transcriptome sequencing (RNA-Seq) on hindgut mucosa samples from individuals treated with antibiotic, probiotic and control diets to determine differentially expressed (DE) host genes. Fifty DE host genes were selected for further characterization using nanofluidic qPCR chips. We used 16S rRNA gene metabarcoding to characterize the rearing water and host gut microbiome (bacterial) communities. Daily administration of antibiotics and probiotics resulted in significant changes in fish gut and aquatic microbiota as well as more than 100 DE genes in the antibiotic and probiotic treatment fish, relative to healthy controls. Normal microbiota depletion by antibiotics mostly led to downregulation of different aspects of immunity and upregulation of apoptotic process. In the probiotic treatment, genes related to post-translation modification and inflammatory responses were up-regulated relative to controls. Our qPCR results revealed significant effects of treatment (antibiotic and probiotic) on rabep2, aifm3, manf, prmt3 gene transcription. Moreover, we found significant associations between members of Lactobacillaceae and Bifidobacteriaceae with host gene expression patterns. Overall, our analysis showed that the microbiota had significant impacts on many host signalling pathways, specifically targeting immune, developmental and metabolic processes. Our characterization of some of the molecular mechanisms involved in microbiome-host interactions will help develop new strategies for preventing/ treating microbiome disruption-related diseases.

Admixture between species is a cause for concern in wildlife management. Canids are particularly vulnerable to interspecific hybridisation, and genetic admixture has shaped their evolutionary history. Microsatellite DNA testing, relying on a small number of genetic markers and geographically restricted reference populations, has identified extensive domestic dog admixture in Australian dingoes and driven conservation management policy. But there exists a concern that geographic variation in dingo genotypes could confound ancestry analyses that use a small number of genetic markers. Here, we apply genome-wide single-nucleotide polymorphism (SNP) genotyping to a set of 402 wild and captive dingoes collected from across Australia and then carry out comparisons to domestic dogs. We then perform ancestry modelling and biogeographic analyses to characterise population structure in dingoes and investigate the extent of admixture between dingoes and dogs in different regions of the continent. We show that there are at least five distinct dingo populations across Australia. We observed limited evidence of dog admixture in wild dingoes. Our work challenges previous reports regarding the occurrence and extent of dog admixture in dingoes, as our ancestry analyses show that previous assessments severely overestimate the degree of domestic dog admixture in dingo populations, particularly in south-eastern Australia. These findings strongly support the use of genome-wide SNP genotyping as a refined method for wildlife managers and policymakers to assess and inform dingo management policy and legislation moving forwards.

Mitochondrial functions are intimately reliant on proteins and RNAs encoded in both the nuclear and mitochondrial genomes, leading to inter-genomic coevolution within taxa. Hybridization can break apart coevolved mitonuclear genotypes, resulting in decreased mitochondrial performance and reduced fitness. This hybrid breakdown is an important component of outbreeding depression and early-stage reproductive isolation. However, the mechanisms contributing to mitonuclear interactions remain poorly resolved. Here, we scored variation in developmental rate (a proxy for fitness) among reciprocal F2 interpopulation hybrids of the intertidal copepod Tigriopus californicus and used RNA sequencing to assess differences in gene expression between fast- and slow-developing hybrids. In total, differences in expression associated with developmental rate were detected for 2925 genes, whereas only 135 genes were differentially expressed as a result of differences in mitochondrial genotype. Upregulated expression in fast developers was enriched for genes involved in chitin-based cuticle development, oxidation-reduction processes, hydrogen peroxide catabolic processes and mitochondrial respiratory chain complex I. In contrast, upregulation in slow developers was enriched for DNA replication, cell division, DNA damage and DNA repair. Eighty-four nuclear-encoded mitochondrial genes were differentially expressed between fast- and slow-developing copepods, including 12 subunits of the electron transport system (ETS) which all had higher expression in fast developers than in slow developers. Nine of these genes were subunits of ETS complex I. Our results emphasize the major roles that mitonuclear interactions within the ETS, particularly in complex I, play in hybrid breakdown, and resolve strong candidate genes for involvement in mitonuclear interactions.

Seasonal environmental fluctuations provide formidable challenges for living organisms, especially small ectotherms such as butterflies. A common strategy to cope with harsh environments is to enter diapause, but some species avoid unsuitable conditions by migrating. Despite a growing understanding of migration in the life cycles of some butterfly species, it remains unknown how individuals register and store environmental cues to determine whether and where to migrate. Here, we explored how competition and host plant availability during larval development affect patterns of DNA methylation in the migratory painted lady (Vanessa cardui) butterfly. We identify a set of potentially functional methylome shifts associated with differences in the environment, indicating that DNA methylation is involved in the response to different conditions during larval development. By analysing the transcriptome for the same samples used for methylation profiling, we also uncovered a non-monotonic relationship between gene body methylation and gene expression. Our results provide a starting point for understanding the interplay between DNA methylation and gene expression in butterflies in general and how differences in environmental conditions during development can trigger unique epigenetic marks that might be important for behavioural decisions in the adult stage.

Individual animals in natural populations tend to host diverse parasite species concurrently over their lifetimes. In free-living ecological communities, organismal life histories shape interactions with their environment, which ultimately forms the basis of ecological succession. However, the structure and dynamics of mammalian parasite communities have not been contextualized in terms of primary ecological succession, in part because few datasets track occupancy and abundance of multiple parasites in wild hosts starting at birth. Here, we studied community dynamics of 12 subtypes of protozoan microparasites (Theileria spp.) in a herd of African buffalo. We show that Theileria communities followed predictable patterns of succession underpinned by four different parasite life history strategies. However, in contrast to many free-living communities, network complexity decreased with host age. Examining parasite communities through the lens of succession may better inform the effect of complex within host eco-evolutionary dynamics on infection outcomes, including parasite co-existence through the lifetime of the host.

After polyploid species are formed, interactions between diploid and polyploid lineages may generate additional diversity in novel cytotypes and phenotypes. In anurans, mate choice by acoustic communication is the primary method by which individuals identify their own species and assess suitable mates. As such, the evolution of acoustic signals is an important mechanism for contributing to reproductive isolation and diversification in this group. Here, we estimate the biogeographical history of the North American grey treefrog complex, consisting of the diploid Hyla chrysoscelis and the tetraploid Hyla versicolor, focusing specifically on the geographical origin of whole genome duplication and the expansion of lineages out of glacial refugia. We then test for lineage-specific differences in mating signals by applying comparative methods to a large acoustic data set collected over 52 years that includes >1500 individual frogs. Along with describing the overall biogeographical history and call diversity, we found evidence that the geographical origin of H. versicolor and the formation of the midwestern polyploid lineage are both associated with glacial limits, and that the southwestern polyploid lineage is associated with a shift in acoustic phenotype relative to the diploid lineage with which they share a mitochondrial lineage. In H. chrysoscelis, we see that acoustic signals are largely split by Eastern and Western lineages, but that northward expansion along either side of the Appalachian Mountains is associated with further acoustic diversification. Overall, results of this study provide substantial clarity on the evolution of grey treefrogs as it relates to their biogeography and acoustic communication.

The gut microbiota have important consequences for host biological processes and there is some evidence that they also affect fitness. However, the complex, interactive nature of ecological factors that influence the gut microbiota has scarcely been investigated in natural populations. We sampled the gut microbiota of wild great tits (Parus major) at different life stages allowing us to evaluate how microbiota varied with respect to a diverse range of key ecological factors of two broad types: (1) host state, namely age and sex, and the life history variables, timing of breeding, fecundity and reproductive success; and (2) the environment, including habitat type, the distance of the nest to the woodland edge, and the general nest and woodland site environments. The gut microbiota varied with life history and the environment in many ways that were largely dependent on age. Nestlings were far more sensitive to environmental variation than adults, pointing to a high degree of flexibility at an important time in development. As nestlings developed their microbiota from one to two weeks of life, they retained consistent (i.e., repeatable) among-individual differences. However these apparent individual differences were driven entirely by the effect of sharing the same nest. Our findings point to important early windows during development in which the gut microbiota are most sensitive to a variety of environmental drivers at multiple scales, and suggest reproductive timing, and hence potentially parental quality or food availability, are linked with the microbiota. Identifying and explicating the various ecological sources that shape an individual's gut bacteria is of vital importance for understanding the gut microbiota's role in animal fitness.

The 'good genes' hypothesis for the evolution of male secondary sexual traits poses that female preferences for such traits are driven by indirect genetic benefits. However, support for the hypothesis remains ambiguous, and, in particular, the genetic basis for the benefits has rarely been investigated. Here, we use seminatural populations of Trinidadian guppies to investigate whether sexually selected traits (orange, black and iridescent colouration, gonopodium length and body size) predict fitness measured as the number of grandoffspring, a metric that integrates across fitness components and sexes. Furthermore, we tested whether two potential sources of genetic benefits-major histocompatibility complex (MHC) genotypes and multilocus heterozygosity (MLH)-are significant predictors of fitness and of the size of sexually selected traits. We found a significant, nonlinear effect of the area of black pigmentation and male body size on the number of grandoffspring, suggesting stabilizing selection on black area, and nonlinear selection favouring small body size. MLH was heritable (h2  = 0.14) and significantly predicted the number of grandoffspring, indicating the potential for genetic benefits based on heterozygosity. We also found support for local heterozygosity effects, which may reflect a noneven distribution of genetic load across the genome. MHC genotype was not significantly associated with any tested fitness component, or with the load of Gyrodactylus parasites. Neither MHC nor MLH was significant predictor of sexually selected traits. Overall, our results highlight the role of heterozygosity in determining fitness, but do not provide support for male sexually selected traits being indicators of genetic quality.

Adaptive phenotypes are shaped by a combination of genetic and environmental forces, but how they interact remains poorly understood. Here, we utilize the cichlid oral jaw apparatus to better understand these gene-by-environment effects. First, we employed RNA-seq in bony and ligamentous tissues important for jaw opening to identify differentially expressed genes between species and across foraging environments. We used two Lake Malawi species adapted to different foraging habitats along the pelagic-benthic ecomorphological axis. Our foraging treatments were designed to force animals to employ either suction or biting/scraping, which broadly mimic pelagic or benthic modes of feeding. We found a large number of differentially expressed genes between species, and while we identified relatively few differences between environments, species differences were far more pronounced when they were challenged with a pelagic versus benthic foraging mode. Expression data carried the signature of genetic assimilation, and implicated cell cycle regulation in shaping the jaw across species and environments. Next, we repeated the foraging experiment and performed ATAC-seq procedures on nuclei harvested from the same tissues. Cross-referencing results from both analyses revealed subsets of genes that were both differentially expressed and differentially accessible. This reduced dataset implicated notable candidate genes including the Hedgehog effector, KIAA0586 and the ETS transcription factor, etv4, which connects environmental stress and craniofacial morphogenesis. Taken together, these data provide novel insights into the epigenetic, genetic and cellular bases of species- and environment-specific bone shapes.

Optimized nutrient utilization is crucial for the progression of microorganisms in competing communities. Here we investigate how different budding yeast species and ecological isolates have established divergent preferences for two alternative sugar substrates: Glucose, which is fermented preferentially by yeast, and galactose, which is alternatively used upon induction of the relevant GAL metabolic genes. We quantified the dose-dependent induction of the GAL1 gene encoding the central galactokinase enzyme and found that a very large diversification exists between different yeast ecotypes and species. The sensitivity of GAL1 induction correlates with the growth performance of the respective yeasts with the alternative sugar. We further define some of the mechanisms, which have established different glucose/galactose consumption strategies in representative yeast strains by modulating the activity of the Gal3 inducer. (1) Optimal galactose consumers, such as Saccharomyces uvarum, contain a hyperactive GAL3 promoter, sustaining highly sensitive GAL1 expression, which is not further improved upon repetitive galactose encounters. (2) Desensitized galactose consumers, such as S. cerevisiae Y12, contain a less sensitive Gal3 sensor, causing a shift of the galactose response towards higher sugar concentrations even in galactose experienced cells. (3) Galactose insensitive sugar consumers, such as S. cerevisiae DBVPG6044, contain an interrupted GAL3 gene, causing extremely reluctant galactose consumption, which is, however, improved upon repeated galactose availability. In summary, different yeast strains and natural isolates have evolved galactose utilization strategies, which cover the whole range of possible sensitivities by modulating the expression and/or activity of the inducible galactose sensor Gal3.

Insects have adapted to a multitude of environmental conditions, including the presence of xenobiotic noxious substances. Environmental microorganisms, particularly rich on ephemeral resources, employ these noxious chemicals in a chemical warfare against predators and competitors, driving co-evolutionary adaptations. In order to analyse how environmental microbes may be driving such evolutionary adaptations, we experimentally evolved Drosophila melanogaster populations by exposing larvae to the toxin-producing mould Aspergillus nidulans that infests the flies' breeding substrate. To disentangle the effects of the mycotoxin Sterigmatocystin from other substrate modifications inflicted by the mould, we used the following four selection regimes: (i) control without fungus, (ii) A. nidulans wild type, (iii) a mutant of A. nidulans ΔlaeA with impaired toxin production, (iv) synthetic Sterigmatocystin. Experimental evolution was carried out in five independent D. melanogaster populations each, for a total of 11 generations. We further combined our evolution experiment with transcriptome analysis to identify evolutionary shifts in gene expression due to the selection regimes and mould confrontation. Populations that evolved in presence of the toxin-producing mould or the pure mycotoxin rapidly adapted to the respective conditions and showed higher viability in subsequent confrontations. Yet, mycotoxin-selected populations had no advantage in A. nidulans wild type confrontation. Moreover, distinctive changes in gene expression related to the selection-regime contrast were only associated with the toxin-producing-fungus regime and comprised a narrow set of genes. Thus, it needs the specific conditions of the selection agent to enable adaptation to the fungus.

Corals are important models for understanding invertebrate host-microbe interactions; however, to fully discern mechanisms involved in these relationships, experimental approaches for manipulating coral-bacteria associations are needed. Coral-associated bacteria affect holobiont health via nutrient cycling, metabolic exchanges and pathogen exclusion, yet it is not fully understood how bacterial community shifts affect holobiont health and physiology. In this study, a combination of antibiotics (ampicillin, streptomycin and ciprofloxacin) was used to disrupt the bacterial communities of 14 colonies of the reef framework-building corals Pocillopora meandrina and P. verrucosa, originally collected from Panama and hosting diverse algal symbionts (family Symbiodiniaceae). Symbiodiniaceae photochemical efficiencies and holobiont oxygen consumption (as proxies for coral health) were measured throughout a 5-day exposure. Antibiotics altered bacterial community composition and reduced alpha and beta diversity, however, several bacteria persisted, leading to the hypothesis that these bacteria are either antibiotics resistant or occupy internal niches that are shielded from antibiotics. While antibiotics did not affect Symbiodiniaceae photochemical efficiency, antibiotics-treated corals had lower oxygen consumption rates. RNAseq revealed that antibiotics increased expression of Pocillopora immunity and stress response genes at the expense of cellular maintenance and metabolism functions. Together, these results reveal that antibiotic disruption of corals' native bacteria negatively impacts holobiont health by decreasing oxygen consumption and activating host immunity without directly impairing Symbiodiniaceae photosynthesis, underscoring the critical role of coral-associated bacteria in holobiont health. They also provide a baseline for future experiments that manipulate Pocillopora corals' symbioses by first reducing the diversity and complexity of coral-associated bacteria.

Environmental variation is increasingly recognized as an important driver of diversity in marine species despite the lack of physical barriers to dispersal and the presence of pelagic stages in many taxa. A robust understanding of the genomic and ecological processes involved in structuring populations is lacking for most marine species, often hindering management and conservation action. Cunner (Tautogolabrus adspersus) is a temperate reef fish with both pelagic early life-history stages and strong site-associated homing as adults; the species is also of interest for use as a cleaner fish in salmonid aquaculture in Atlantic Canada. We aimed to characterize genomic and geographic differentiation of cunner in the Northwest Atlantic. To achieve this, a chromosome-level genome assembly for cunner was produced and used to characterize spatial population structure throughout Atlantic Canada using whole-genome sequencing. The genome assembly spanned 0.72 Gbp and 24 chromosomes; whole-genome sequencing of 803 individuals from 20 locations from Newfoundland to New Jersey identified approximately 11 million genetic variants. Principal component analysis revealed four regional Atlantic Canadian groups. Pairwise FST and selection scans revealed signals of differentiation and selection at discrete genomic regions, including adjacent peaks on chromosome 10 across multiple pairwise comparisons (i.e. FST 0.5-0.75). Redundancy analysis suggested association of environmental variables related to benthic temperature and oxygen range with genomic structure. Results suggest regional scale diversity in this temperate reef fish and can directly inform the collection and translocation of cunner for aquaculture applications and the conservation of wild populations throughout the Northwest Atlantic.

While shaping of plant microbiome composition through 'host filtering' is well documented in legume-rhizobium symbioses, it is less clear to what extent different varieties and genotypes of the same plant species differentially influence symbiont community diversity and composition. Here, we compared how clover host varieties and genotypes affect the structure of Rhizobium populations in root nodules under conventional field and controlled greenhouse conditions. We first grew four Trifolium repens (white clover) F2 crosses and one variety in a conventional field trial and compared differences in root nodule Rhizobium leguminosarum symbiovar trifolii (Rlt) genotype diversity using high-throughput amplicon sequencing of chromosomal housekeeping (rpoB and recA) genes and auxiliary plasmid-borne symbiosis genes (nodA and nodD). We found that Rlt nodule diversities significantly differed between clover crosses, potentially due to host filtering. However, variance in Rlt diversity largely overlapped between crosses and was also explained by the spatial distribution of plants in the field, indicative of the role of local environmental conditions for nodule diversity. To test the effect of host filtering, we conducted a controlled greenhouse trial with a diverse Rlt inoculum and several host genotypes. We found that different clover varieties and genotypes of the same variety selected for significantly different Rlt nodule communities and that the strength of host filtering (deviation from the initial Rhizobium inoculant composition) was positively correlated with the efficiency of symbiosis (rate of plant greenness colouration). Together, our results suggest that selection by host genotype and local growth conditions jointly influence white clover Rlt nodule diversity and community composition.

Understanding the molecular basis of repeated evolution improves our ability to predict evolution across the tree of life. Only since the last decade has high-throughput sequencing enabled comparative genome scans to thoroughly examine the repeatability of genetic changes driving repeated phenotypic evolution. The Asian corn borer (ACB), Ostrinia furnacalis (Guenée), and the European corn borer (ECB), Ostrinia nubilalis (Hübner), are two closely related moths displaying repeatable phenological adaptation to a wide range of climates on two separate continents, largely manifesting as changes in the timing of diapause induction and termination across latitude. Candidate genes underlying diapause variation in North American ECB have been previously identified. Here, we sampled seven ACB populations across 23 degrees of latitude in China to elucidate the genetic basis of diapause variation and evolutionary mechanisms driving parallel clinal responses in the two species. Using pooled whole-genome sequencing (Pool-seq) data, population genomic analyses revealed hundreds of single nucleotide polymorphisms (SNP) whose allele frequencies covaried with mean diapause phenotypes along the cline. Genes involved in circadian rhythm were over-represented among candidate genes with strong signatures of spatially varying selection. Only one of two circadian clock genes associated with diapause evolution in ECB showed evidence of reuse in ACB (period [per]), but per alleles were not shared between species nor with their outgroup, implicating independent mutational paths. Nonetheless, evidence of adaptive introgression was discovered at putative diapause loci located elsewhere in the genome, suggesting that de novo mutations and introgression might both underlie the repeated phenological evolution.

Genes of the major histocompatibility complex (MHC) play important roles in vertebrate immunocompetence. MHC genes thus offer females indirect benefits to mate choice through the production of offspring of an optimal MHC genotype. Females may choose males with specific MHC haplotypes, dissimilar MHC genotypes, MHC heterozygous males or MHC-diverse males. We tested these four alternatives for both female social and paternal choice in wild golden snub-nosed monkeys (Rhinopithecus roxellana) by examining overall genetic variability (via microsatellites) and four MHC-genes (DRB1, DRB2, DQA1 and DQB1). Monte Carlo randomization tests showed that MHC dissimilarity was favoured for social choice (males to which females were socially affiliated) and intermediate MHC dissimilarity was favoured in paternal choice (fathers of offspring). No evidence of inbreeding avoidance was found for either social or paternal mates. We found that MHC heterozygotes, higher microsatellite multilocus heterozygosity and higher microsatellites diversity were favoured for social mates, and higher microsatellite diversity was favoured for paternal mates. Independent of male age, we found that the formation of male-female social pairings is significantly predicted by compatibility based on the sharing of MHC haplotypes. However, we found no evidence of independent genetic effects on the duration of male-female social pairings, male social status (achieving OMU leader male status or not), the number of females with which individual leader males paired, the likelihood of potential male-female pairings producing offspring, or whether males fathered offspring or not. Overall, our findings suggest different genetic factors are involved in social and paternal choice in R. roxellana.

Phenotypic plasticity can play an important role in the ability of animals to tolerate environmental stress, but the nature and magnitude of plastic responses are often specific to the developmental timing of exposure. Here, we examine changes in gene expression in the diaphragm of highland deer mice (Peromyscus maniculatus) in response to hypoxia exposure at different stages of development. In highland deer mice, developmental plasticity in diaphragm function may mediate changes in several respiratory traits that influence aerobic metabolism and performance under hypoxia. We generated RNAseq data from diaphragm tissue of adult deer mice exposed to (1) life-long hypoxia (before conception to adulthood), (2) post-natal hypoxia (birth to adulthood), (3) adult hypoxia (6-8 weeks only during adulthood) or (4) normoxia. We found five suites of co-regulated genes that are differentially expressed in response to hypoxia, but the patterns of differential expression depend on the developmental timing of exposure. We also identified four transcriptional modules that are associated with important respiratory traits. Many of the genes in these transcriptional modules bear signatures of altitude-related selection, providing an indirect line of evidence that observed changes in gene expression may be adaptive in hypoxic environments. Our results demonstrate the importance of developmental stage in determining the phenotypic response to environmental stressors.

Infectious diseases of wildlife continue to pose a threat to biodiversity worldwide, yet pathogens are far from uniform in virulence or host disease outcome. Within the same pathogen species, virulence can vary considerably depending on strain or lineage, in turn eliciting variable host responses. One pathogen that has caused extensive biodiversity loss is the amphibian-killing fungus, Batrachochytrium dendrobatidis (Bd), which is comprised of a globally widespread hypervirulent lineage (Bd-GPL), and multiple geographically restricted, enzootic lineages. Whereas host immunogenomic responses to Bd-GPL have been characterized in a number of amphibian species, immunogenomic responses to geographically restricted, enzootic Bd lineages are less clear. To examine lineage-specific host immune responses to Bd, we exposed a species of pumpkin toadlet, Brachycephalus pitanga, which is endemic to Brazil's Southern Atlantic Forest, to either the Bd-GPL or the enzootic Bd-Asia-2/Brazil (hereafter Bd-Brazil) lineage. Using temporal samples from early, mid, and late infection stages, we quantified functional immunogenomic responses over the course of infection using differential gene expression tests and coexpression network analyses. Host immune responses varied significantly with Bd lineage. Relative to controls, toadlet responses to Bd-Brazil were weak at early infection (25 genes significantly differentially expressed), peaked by mid-stage infection (414 genes), and were nearly fully resolved by late-stage infection (nine genes). In contrast, responses to Bd-GPL were magnified and delayed; toadlets significantly differentially expressed 111 genes early, 87 genes at mid-stage infection, and 726 genes by late-stage infection relative to controls. Given that infection intensity did not vary between mid- and late-stage disease in either Bd-Brazil or Bd-GPL treatments, this suggests that pumpkin toadlets may be at least partially tolerant to the enzootic Bd-Brazil lineage. In contrast, late-stage immune activation against Bd-GPL was consistent with immune dysregulation previously observed in other species. Our results demonstrate that both the timing of immune response and the particular immune pathways activated are specific to Bd lineage. Within regions where multiple Bd lineages co-occur, and given continued global Bd movement, these differential host responses may influence not only individual disease outcome, but transmission dynamics at the population and community levels.

Seasonal migration of Nearctic-Neotropical passerine birds may have profound effects on the diversity and abundance of their host-associated microbiota. Migratory birds experience seasonal change in environments and diets throughout the course of the annual cycle that, along with recurrent biological events such as reproduction, may significantly impact their microbiota. In this study, we characterize the intestinal microbiota of four closely related species of migratory Catharus thrushes at three time points of their migratory cycle: during spring migration, on the summer breeding territories and during fall migration. Using observations replicated over 3 years, we determined that microbial community diversity of Catharus thrushes was significantly different across distinct time periods of the annual cycle, whereas community composition was more similar within than across years. Elevated alpha diversity in the summer birds compared to either migratory period indicated that birds may harbour a reduced microbiota during active migration. We also found that community composition of the microbiota did not substantially differ between host species. Finally, we recovered two phyla, Cyanobacteria and Planctomycetota, which are not commonly described from birds, that were in relatively high abundance in specific years. This study contributes to our growing understanding of how microbiota in wild birds vary throughout disparate ecological conditions and reveals potential axes across which an animal's microbial flexibility adapts to variable environments and recurrent biological conditions throughout the annual cycle.

The outcome of natural enemy attack in insects is commonly impacted by the presence of defensive microbial symbionts residing within the host. The thermal environment is a factor known to affect symbiont-mediated traits in insects. Lower temperatures, for instance, have been shown to reduce Spiroplasma-mediated protection in Drosophila. Our understanding of protective symbiosis requires a deeper understanding of environment-symbiont-protection links. Here, we dissect the effect of the thermal environment on Spiroplasma-mediated protection against Leptopilina boulardi in Drosophila melanogaster by examining the effect of temperature before, during and after wasp attack on fly survival and wasp success. We observed that the developmental temperature of the mothers of attacked larvae, but not the temperature of the attacked larvae themselves during or after wasp attack, strongly determines the protective influence of Spiroplasma. Cooler maternal environments were associated with weaker Spiroplasma protection of their progeny. The effect of developmental temperature on Spiroplasma-mediated protection is probably mediated by a reduction in Spiroplasma titre. These results indicate that historical thermal environment is a stronger determinant of protection than current environment. Furthermore, protection is a character with transgenerational nongenetic variation probably to produce complex short-term responses to selection. In addition, the cool sensitivity of the Spiroplasma-Drosophila symbioses contrasts with the more common failure of symbioses at elevated temperatures, indicating a need to understand the mechanistic basis of low temperature sensitivity on this symbiosis.

The extremely high species diversity of soil bacterial community has fascinated and puzzled community ecologists. Although theory predicts that fluctuations in environments can facilitate diversity maintenance, the effects of fluctuating temperature on species diversity have rarely been investigated in species-rich microbial communities. Here, we examined whether fluctuating temperature had positive effects on species diversity relative to constant temperatures in soil bacterial communities, and investigated the effects of fluctuating temperature on bacterial performances (changes in relative abundance). We performed a temperature manipulation experiment with soils collected from temperate and subtropical zones, where the soils were subjected to constant high, low or fluctuating temperatures. We found that fluctuating temperatures showed significant positive effects on species diversity. The time-averaged effect of fluctuating temperatures (i.e., averaging out the differences between species in their environment-dependent performances) appeared to delay species loss in both the temperate and the subtropical communities. In addition, we found that the performances of temperature-responsive species at fluctuating temperatures significantly deviated from their time-weighted average performances at constant high and low temperatures, which was defined as fluctuation-dependent effects in our study. Intriguingly, fluctuation-dependent effects beyond time-averaged effect led to an opposite trend: differences in temperature-responsive species' performances decreased in the temperate communities, but increased in the subtropical communities. Our findings provide new insights into diversity maintenance in soil bacterial communities, and imply that the effects of fluctuating temperature on species diversity in soil bacterial community might vary across latitude.

Recent studies have highlighted associations between diseases and host microbiota. It remains extremely challenging - especially under natural conditions - to clarify whether host microbiota promote future infections, or whether changes in host microbiota result from infections. Nonetheless, deciphering between these two processes is essential for highlighting the role of microbes in disease progression. We longitudinally surveyed, in the wild, the microbiota of individual fish hosts (Leuciscus burdigalensis) both before and after infection by a crustacean ectoparasite (Tracheliastes polycolpus). We found a striking association between parasite infection and the host microbiota composition restricted to the fins the parasite anchored. We clearly demonstrated that infections by the parasite induced a shift in (and did not result from) the host fin microbiota. Furthermore during infection, the microbiota of infected fins got similar to the microbiota of the adult stage, and the free-living infective stage of the parasite with a predominance of the Burkholderiaceae bacteria family. This suggests that some Burkholderiaceae bacteria are involved in a coinfection process and possibly facilitate T. polycolpus infection. In this study, we reveal novel mechanistic insights for understanding the role of the microbiota in host-parasite interactions, which has implications for predicting the progression of diseases in natural host populations.

Skin microbiota play an important role in protecting bat hosts from the fungal pathogen Pseudogymnoascus destructans, which has caused dramatic bat population declines and extinctions. Recent studies have provided insights into the bacterial communities of bat skin, but variation in skin bacterial community structure in the context of the seasonal dynamics of fungal invasion, as well as the processes that drive such variation, remain largely unexplored. In this study, we characterized bat skin microbiota over the course of the bat hibernation and active season stages and used a neutral model of community ecology to determine the relative roles of neutral and selective processes in driving microbial community variation. Our results showed significant seasonal shifts in skin community structure, as well as less diverse microbiota in hibernation than in the active season. Skin microbiota were influenced by the environmental bacterial reservoir. During both the hibernation and active season stages, more than 78% of ASVs in bat skin microbiota were consistent with neutral distribution, implying that neutral processes, that is, dispersal or ecological drift contributing the most to shifts in skin microbiota. In addition, the neutral model showed that some ASVs were actively selected by the bats from the environmental bacterial reservoir, accounting for approximately 20% and 31% of the total community during hibernation and active season stages, respectively. Overall, this research provides insights into the assemblage of bat-associated bacterial communities and will aid in the development of conservation strategies against fungal disease.

Birds are among the most colourful terrestrial vertebrates, with various plumage colours and patterns. We conducted a genome-wide association study (GWAS) on an intercross F2 population of Pekin ducks and mallards (n = 722) and identified a 1.57-Mb genetic region (Chr11: 20,176,480-21,750,101 bp) related to duck melanism. Fine mapping by linkage disequilibrium (LD) and FST analysis narrowed the final candidate region to a region of 22,500 bp (Chr11: 20,677,500-20,700,000 bp) including three coding genes, TCF25, MC1R and TUBB3. Combined with transcriptome and qRT-PCR analysis, MC1R was identified as the unique genetic locus responsible for black plumage in ducks, and it was significantly more highly expressed in the feather bulbs of black ducks. We also identified 52G > A (Chr11: 20,696,354G > A) and 376G > A (Chr11: 20,696,678G > A) mutations in the MC1R coding region that have been widely studied in ducks. In addition, structural variations (SVs) were screened by nanopore sequencing, and no significant SV was found to be associated with the duck black plumage trait. However, we identified four novel single nucleotide polymorphisms in the MC1R regulator region (Chr11: 20,678,412G > A, Chr11: 20,679,236G > A, Chr11: 20,692,496 A > G and Chr11: 20,692,791 A > G) that had a strong association with the black plumage phenotype of ducks and combined with potential changes in transcription binding affinities. The luciferase reporter gene assay demonstrated that Chr11: 20,678,412G > A and Chr11: 20,679,236G > A led to significant promoter activity changes. Our research emphasizes the importance of MC1R regulatory region mutation in determining the duck black plumage phenotype, and these results expand our understanding of the genetic mechanism underlying duck plumage colour.

Bacterioplankton underpin biogeochemical cycles and an improved understanding of the patterns and drivers of variability in their distribution is needed to determine their wider functioning and importance. Sharp environmental gradients and dispersal barriers associated with ocean fronts are emerging as key determinants of bacterioplankton biodiversity patterns. We examined how the development of the Celtic Sea Front (CF), a tidal mixing front on the Northwest European Shelf affects bacterioplankton communities. We performed 16S-rRNA metabarcoding on 60 seawater samples collected from three depths (surface, 20 m and seafloor), across two research cruises (May and September 2018), encompassing the intra-annual range of the CF intensity. Communities above the thermocline of stratified frontal waters were clearly differentiated and less diverse than those below the thermocline and communities in the well-mixed waters of the Irish Sea. This effect was much more pronounced in September, when the CF was at its peak intensity. The stratified zone likely represents a stressful environment for bacterioplankton due to a combination of high temperatures and low nutrients, which fewer taxa can tolerate. Much of the observed variation was driven by Synechococcus spp. (cyanobacteria), which were more abundant within the stratified zone and are known to thrive in warm oligotrophic waters. Synechococcus spp. are key contributors to global primary productivity and carbon cycling and, as such, variability driven by the CF is likely to influence regional biogeochemical processes. However, further studies are required to explicitly link shifts in community structure to function and quantify their wider importance to pelagic ecosystems.

The environment experienced during early life is a crucial factor in the life of many organisms. This early life environment has been shown to have profound effects on morphology, physiology and fitness. However, the molecular mechanisms that mediate these effects are largely unknown, even though they are essential for our understanding of the processes that induce phenotypic variation in natural populations. DNA methylation is an epigenetic mechanism that has been suggested to explain such environmentally induced phenotypic changes early in life. To investigate whether DNA methylation changes are associated with experimentally induced early developmental effects, we cross-fostered great tit (Parus major) nestlings and manipulated their brood sizes in a natural study population. We assessed experimental brood size effects on pre-fledging biometry and behaviour. We linked this to genome-wide DNA methylation levels of CpG sites in erythrocyte DNA, using 122 individuals and an improved epiGBS2 laboratory protocol. Brood enlargement caused developmental stress and negatively affected nestling condition, predominantly during the second half of the breeding season, when conditions are harsher. Brood enlargement, however, affected nestling DNA methylation in only one CpG site and only if the hatch date was taken into account. In conclusion, this study shows that nutritional stress in enlarged broods does not associate with direct effects on genome-wide DNA methylation. Future studies should assess whether genome-wide DNA methylation variation may arise later in life as a consequence of phenotypic changes during early development.

Fungal communities are diverse and abundant in coastal waters, yet, their ecological roles and adaptations remain largely unknown. To address these gaps, ITS2 metabarcoding and metatranscriptomic analyses were used to capture the whole suite of fungal diversity and their metabolic potential in water column and sediments in the Yellow Sea during August and October 2019. ITS2 metabarcoding described successfully the abundance of Dikarya during August and October at the different examined habitats, but strongly underrepresented or failed to identify other fungal taxa, including zoosporic and early-diverging lineages, that were abundant in the mycobiome as uncovered by metatranscriptomes. Metatranscriptomics also revealed enriched expression of genes annotated to zoosporic fungi (e.g., chytrids) mainly in the surface water column in October. This enriched expression was correlated with the two-fold increase in chlorophyll-a intensity attributed to phytoplanktonic species which are known to be parasitized by chytrids. The concurrent high expression of genes related to calcium signalling and GTPase activity suggested that these metabolic traits facilitate the parasitic lifestyle of chytrids. Similarly, elevated expression of phagosome genes annotated to Rozellomycota, an early-diverging fungal phylum not fully detected with ITS2 metabarcoding, suggested that this taxon utilizes a suite of feeding modes, including phagotrophy in this coastal setting. Our data highlight the necessity of using combined approaches to accurately describe the community structure of coastal mycobiome. We also provide in-depth insights into the fungal ecological roles in coastal waters, and report potential metabolic mechanisms utilized by fungi to cope with environmental stresses that occur during distinct seasonal months in coastal ecosystems.

Speciation is a central topic in evolutionary biology. However, how genomic divergence originates and accumulates in the face of gene flow during ecological adaptation remains poorly understood. Closely related species that have adapted to distinct environments but inhabit some overlapping ranges provide an ideal system to evaluate this issue. Here, we combine population genomics and species distribution models (SDMs) to examine genomic divergences between two sister plant species, Medicago ruthenica and M. archiducis-nicolai, that occur in northern China and the northeast Qinghai-Tibet Plateau, respectively, with overlapping distributions in the border of the two regions. M. ruthenica and M. archiducis-nicolai are well-delimited based on population genomic data, although hybrids exist in sympatric sampling locations. Coalescent simulations and SDMs suggest that the two species diverged from each other in the Quaternary but have been in continuous contact with gene flow occurring between the two species since then. We also discovered positive selection signatures associated with genes both outside and within genomic islands in both species that are probably involved in adaptations to arid and high-altitude environments. Our findings provide insights into how natural selection and climatic changes in the Quaternary initiated and maintained interspecific divergence of these two sister species.

Within the Balkan Peninsula, topographic and climatic agents have promoted biodiversity and shaped the speciation history of many ectotherms. Here, we targeted an iconic European reptile, the nose-horned viper species-complex (Vipera ammodytes), and explored its spatial and temporal evolution. We (i) utilized genome-wide single nucleotide polymorphisms to infer genetic structure and build a time-calibrated species-tree, and (ii) applied species distribution modelling with niche-divergence tests among major phylogenomic clades. Geographically structured genetic diversity was found. Cycles of recurrent isolation and expansion during glacial-interglacial periods led to allopatric speciation and to secondary contacts and formation of multiple hybrid zones throughout the Balkan Peninsula. Deep divergence is still detected among populations separated by old and imminent biogeographical barriers (Pindos Mountain Range, the Cyclades islands, etc.), but in most cases speciation is incomplete. At the other end of the speciation continuum, we recognize two well-differentiated lineages, currently lacking any evidence of gene flow; one is distributed in the Northwestern Balkans and the other in the Southeastern Balkans, further expanding into Asia. Despite their split 5 million years ago, there is no evidence of ecological divergence, as speciation probably occurred in niche-pockets of analogous environments. These two lineages probably represent different species, while V. transcaucasiana does not merit species status. By comparing the genomic phylogenies to an updated mitochondrial one, we propose an evolutionary scenario that resolves all mitonuclear conflicts, according to which the history of the V. ammodytes species-complex was shaped by complex processes, including a major event of introgressive hybridization with asymmetric mitochondrial capture.

In oceanic ecosystems, the nature of barriers to gene flow and the processes by which populations may become isolated are different from the terrestrial environment, and less well understood. In this study we investigate a highly mobile species (the sperm whale, Physeter macrocephalus) that is genetically differentiated between an open North Atlantic population and the populations in the Mediterranean Sea. We apply high-resolution single nucleotide polymorphism (SNP) analysis to study the nature of barriers to gene flow in this system, assessing the putative boundary into the Mediterranean (Strait of Gibraltar and Alboran Sea region), and including novel analyses on structuring among sperm whale populations within the Mediterranean basin. Our data support a recent founding of the Mediterranean population, around the time of the last glacial maximum, and show concerted historical demographic profiles in both the Atlantic and the Mediterranean. In each region there is evidence for a population decline around the time of the founder event. The largest decline was seen within the Mediterranean Sea where effective population size is substantially lower (especially in the eastern basin). While differentiation is strongest at the Atlantic/Mediterranean boundary, there is also weaker but significant differentiation between the eastern and western basins of the Mediterranean Sea. We propose, however, that the mechanisms are different. While post-founding gene flow was reduced between the Mediterranean and Atlantic populations, within the Mediterranean an important factor differentiating the basins is probably a greater degree of admixture between the western basin and the North Atlantic and some level of isolation between the western and eastern Mediterranean basins. Subdivision within the Mediterranean Sea exacerbates conservation concerns and will require consideration of what distinct impacts may affect populations in the two basins.

Phylogeographic patterns in large mammals result from natural environmental factors and anthropogenic effects, which in some cases include domestication. The grey wolf was once widely distributed across the Holarctic, but experienced phylogeographic shifts and demographic declines during the Holocene. In the 19th-20th centuries, the species became extirpated from large parts of Europe due to direct extermination and habitat loss. We reconstructed the evolutionary history of the extinct Western European wolves based on the mitogenomic composition of 78 samples from France (Neolithic-20th century) in the context of other populations of wolves and dogs worldwide. We found a close genetic similarity of French wolves from ancient, medieval and recent populations, which suggests the long-term continuity of maternal lineages. MtDNA haplotypes of the French wolves showed large diversity and fell into two main haplogroups of modern Holarctic wolves. Our worldwide phylogeographic analysis indicated that haplogroup W1, which includes wolves from Eurasia and North America, originated in Northern Siberia. Haplogroup W2, which includes only European wolves, originated in Europe ~35 kya and its frequency was reduced during the Holocene due to an expansion of haplogroup W1 from the east. Moreover, we found that dog haplogroup D, currently restricted to Europe and the Middle East, was nested within the wolf haplogroup W2. This suggests European origin of haplogroup D, probably as a result of an ancient introgression from European wolves. Our results highlight the dynamic evolutionary history of European wolves during the Holocene, with a partial lineage replacement and introgressive hybridization with local dog populations.

The blacklegged tick (Ixodes scapularis (Journal of the Academy of Natural Sciences of Philadelphia, 1821, 2, 59)) is a vector of Borrelia burgdorferi sensu stricto (s.s.) (International Journal of Systematic Bacteriology, 1984, 34, 496), the causative bacterial agent of Lyme disease, part of a slow-moving epidemic of Lyme borreliosis spreading across the northern hemisphere. Well-known geographical differences in the vectorial capacity of these ticks are associated with genetic variation. Despite the need for detailed genetic information in this disease system, previous phylogeographical studies of these ticks have been restricted to relatively few populations or few genetic loci. Here we present the most comprehensive phylogeographical study of genome-wide markers in I. scapularis, conducted by using 3RAD (triple-enzyme restriction-site associated sequencing) and surveying 353 ticks from 33 counties throughout the species' range. We found limited genetic variation among populations from the Northeast and Upper Midwest, where Lyme disease is most common, and higher genetic variation among populations from the South. We identify five spatially associated genetic clusters of I. scapularis. In regions where Lyme disease is increasing in frequency, the I. scapularis populations genetically group with ticks from historically highly Lyme-endemic regions. Finally, we identify 10 variable DNA sites that contribute the most to population differentiation. These variable sites cluster on one of the chromosome-scale scaffolds for I. scapularis and are within identified genes. Our findings illuminate the need for additional research to identify loci causing variation in the vectorial capacity of I. scapularis and where additional tick sampling would be most valuable to further understand disease trends caused by pathogens transmitted by I. scapularis.

The key role of symbiotic skin bacteria communities in amphibian resistance to emerging pathogens is well recognized, but factors leading to their dysbiosis are not fully understood. In particular, the potential effects of population translocations on the composition and diversity of hosts' skin microbiota have received little attention, although such transfers are widely carried out as a strategy for amphibian conservation. To characterize the potential reorganization of the microbiota over such a sudden environmental change, we conducted a common-garden experiment simulating reciprocal translocations of yellow-spotted salamander larvae across three lakes. We sequenced skin microbiota samples collected before and 15 days after the transfer. Using a database of antifungal isolates, we identified symbionts with known function against the pathogen Batrachochytrium dendrobatidis, a major driver of amphibian declines. Our results indicate an important reorganization of bacterial assemblages throughout ontogeny, with strong changes in composition, diversity and structure of the skin microbiota in both control and translocated individuals over the 15 days of monitoring. Unexpectedly, the diversity and community structure of the microbiota were not significantly affected by the translocation event, thus suggesting a strong resilience of skin bacterial communities to environmental change-at least across the time-window studied here. A few phylotypes were more abundant in the microbiota of translocated larvae, but no differences were found among pathogen-inhibiting symbionts. Taken together, our results support amphibian translocations as a promising strategy for this endangered animal class, with limited impact on their skin microbiota.

Nature has evolved a wealth of sex determination (SD) mechanisms, driven by both genetic and environmental factors. Recent studies of SD in fishes have shown that not all taxa fit the classic paradigm of sex chromosome evolution and diverse SD methods can be found even among closely related species. Here, we apply a suite of genomic approaches to investigate sex-biased genomic variation in eight species of Sebastes rockfish found in the northeast Pacific Ocean. Using recently assembled chromosome-level rockfish genomes, we leverage published sequence data to identify disparate sex chromosomes and sex-biased loci in five species. We identify two putative male sex chromosomes in S. diaconus, a single putative sex chromosome in the sibling species S. carnatus and S. chrysomelas, and an unplaced sex determining contig in the sibling species S. miniatus and S. crocotulus. Our study provides evidence for disparate means of sex determination within a recently diverged set of species and sheds light on the diverse origins of sex determination mechanisms present in the animal kingdom.

Polymorphic species are useful models for investigating the evolutionary processes driving diversification. Such processes include colonization history as well as contemporary selection, gene flow, and genetic drift, which can vary between intraspecific morphs as a function of their distinct life histories. The interactive and relative influence of such evolutionary processes on morph differentiation critically informs morph-specific management decisions and our understanding of incipient speciation. We therefore investigated how geographic distance, environmental conditions, and colonization history interacted with morph migratory capacity in the highly polymorphic fish species, Arctic Charr (Salvelinus alpinus). Using an 87 k SNP chip we genetically characterized recently evolved anadromous, resident, and landlocked charr collected from 45 locations across a secondary contact zone of three charr glacial lineages in eastern Canada. A strong pattern of isolation by distance across all populations suggested geographic distance principally shaped genetic structure. Landlocked populations had lower genetic diversities and higher genetic differentiation than anadromous populations. However, effective population size was generally temporally stable in landlocked populations in comparison to anadromous populations. Genetic diversity positively correlated with latitude, potentially indicating southern anadromous populations' vulnerability to climate change and greater introgression between the Arctic and Atlantic glacial lineages in northern Labrador. Local adaptation was suggested by the observation of several environmental variables strongly associating with functionally relevant outlier genes including a region on chromosome AC21 potentially associated with anadromy. Our results demonstrate that gene flow, colonization history, and local adaptation uniquely interact to influence the genetic variation and evolutionary trajectory of populations.

Industrialization-including urbanization, participation in the global food chain and consumption of heavily processed foods-is thought to drive substantial shifts in the human microbiome. While diet strongly influences stool microbiome composition, the influence of diet on the oral microbiome is largely speculative. Multiple ecologically distinct surfaces in the mouth, each harbouring a unique microbial community, pose a challenge to assessing changes in the oral microbiome in the context of industrialization, as the results depend on the oral site under study. Here, we investigated whether microbial communities of dental plaque, the dense biofilm on non-shedding tooth surfaces, are distinctly different across populations with dissimilar subsistence strategies and degree of industrialized market integration. Using a metagenomic approach, we compared the dental plaque microbiomes of Baka foragers and Nzime subsistence agriculturalists in Cameroon (n = 46) with the dental plaque and calculus microbiomes of highly industrialized populations in North America and Europe (n = 38). We found that differences in microbial taxonomic composition between populations were minimal, with high conservation of abundant microbial taxa and no significant differences in microbial diversity related to dietary practices. Instead, we find that the major source of variation in dental plaque microbial species composition is related to tooth location and oxygen availability, which may be influenced by toothbrushing or other dental hygiene measures. Our results support that dental plaque, in contrast to the stool microbiome, maintains an inherent stability against ecological perturbations in the oral environment.

Transition from sexual reproduction to parthenogenesis constitutes a major life-history change with deep evolutionary consequences for sex-related traits, which are expected to decay. The pea aphid Acyrthosiphon pisum shows intraspecific reproductive polymorphism, with cold-resistant cyclically parthenogenetic (CP) lineages that alternate sexual and asexual generations and cold-sensitive obligately parthenogenetic (OP) lineages that produce only asexual females but still males. Here, the genotyping of 219 pea aphid lineages collected in cold-winter and mild-winter regions revealed contrasting population structures. Samples from cold-winter regions consisted mostly of distinct multilocus genotypes (MLGs) usually represented by a single sample (101 different MLGs for 111 samples) and were all phenotyped as CP. In contrast, fewer MLGs were found in mild-winter regions (28 MLGs for 108 samples), all but one being OP. Since the males produced by OP lineages are unlikely to pass on their genes (sexual females being rare in mild-winter regions), we tested the hypothesis that their traits could degenerate due to lack of selection by comparing male production and male reproductive success between OP and CP lineages. Male production was indeed reduced in OP lineages, but a less clear pattern was observed for male reproductive success: females mated with OP males laid fewer eggs (fertilized or not) but OP and CP males fertilized the same proportion of eggs. These differences may stem from the type of selective forces: male production may be counter-selected whereas male performances may evolve under the slower process of relaxed selection. The overall effective reproductive capacity of OP males could result from recent sex loss in OP lineages or underestimated reproductive opportunities.

Caraway (Carum carvi L.) is a crop species that is gaining in importance in Europe, especially as a condiment and medicinal plant. Here, we present the plant-pollinator network of caraway in a central European agricultural landscape, focusing on two diverse potential pollinator taxa, Diptera: Brachycera (= true flies) and Hymenoptera (sawflies, bees, and wasps). We specifically studied qualitative differences in interactions between the two insect taxa as well as the intraday and intraseasonal variability of the network. Insect and pollen plant species determination was done via morphological identification and DNA (meta)barcoding. In total, 121 species representing 33 families of Hymenoptera and Brachycera were found to carry caraway pollen. These taxa included many nonhoneybee and nonhoverfly species, showing a wide taxonomic breadth of potential pollinators and a higher network complexity than previously anticipated. There are distinct qualitative differences between Brachycera and Hymenoptera networks, suggesting complementary roles of both taxa in the pollination of native and crop plants. Strong intraday differences in potential pollinator diversity make it necessary to collect insects and pollen at different times of the day to compile complete plant-pollinator networks. Intraseasonal analyses of the plant-pollinator network of caraway show the potential of caraway as an important food source for insect species with an activity peak in late summer.

Understanding how microbial communities are shaped across spatial dimensions is of fundamental importance in microbial ecology. However, most studies on soil biogeography have focused on the topsoil microbiome, while the factors driving the subsoil microbiome distribution are largely unknown. Here we used 16S rRNA amplicon sequencing to analyse the factors underlying the bacterial β-diversity along vertical (0-240 cm of soil depth) and horizontal spatial dimensions (~500,000 km2 ) in the U.S. Corn Belt. With these data we tested whether the horizontal or vertical spatial variation had stronger impacts on the taxonomic (Bray-Curtis) and phylogenetic (weighted Unifrac) β-diversity. Additionally, we assessed whether the distance-decay (horizontal dimension) was greater in the topsoil (0-30 cm) or subsoil (in each 30 cm layer from 30-240 cm) using Mantel tests. The influence of geographic distance versus edaphic variables on the bacterial communities from the different soil layers was also compared. Results indicated that the phylogenetic β-diversity was impacted more by soil depth, while the taxonomic β-diversity changed more between geographic locations. The distance-decay was lower in the topsoil than in all subsoil layers analysed. Moreover, some subsoil layers were influenced more by geographic distance than any edaphic variable, including pH. Although different factors affected the topsoil and subsoil biogeography, niche-based models explained the community assembly of all soil layers. This comprehensive study contributed to elucidating important aspects of soil bacterial biogeography including the major impact of soil depth on the phylogenetic β-diversity, and the greater influence of geographic distance on subsoil than on topsoil bacterial communities in agroecosystems.

Global habitat degradation heightens the need to better understand patterns of genetic connectivity and diversity of marine biota across geographical ranges to guide conservation efforts. Corals across the Red Sea are subject to pronounced environmental differences, but studies so far suggest that animal populations are largely connected, excepting evidence for a genetic break between the northern-central and southern regions. Here, we investigated population structure and holobiont assemblage of two common pocilloporid corals, Pocillopora verrucosa and Stylophora pistillata, across the Red Sea. We found little evidence for population differentiation in P. verrucosa, except for the southernmost site. Conversely, S. pistillata exhibited a complex population structure with evidence for within-reef and regional genetic differentiation, in line with differences in their reproductive mode (P. verrucosa is a broadcast spawner and S. pistillata is a brooder). Analysis for genomic loci under positive selection identified 85 sites (18 of which were in coding sequences) that distinguished the southern P. verrucosa population from the remainder of the Red Sea population. By comparison, we found 128 loci (24 of which were residing in coding sequences) in S. pistillata with evidence for local adaptation at various sites. Functional annotation of the underlying proteins revealed putative roles in the response to stress, lipid metabolism, transport, cytoskeletal rearrangement, and ciliary function (among others). Microbial assemblages of both coral species showed pervasive association with microalgal symbionts from the genus Symbiodinium (former clade A) and bacteria from the genus Endozoicomonas that exhibited significant differences according to host genotype and environment. The disparity of population genetic and holobiont assemblage patterns even between closely related species (family Pocilloporidae) highlights the need for multispecies investigations to better understand the role of the environment in shaping evolutionary trajectories. It further emphasizes the importance of networks of reef reserves to achieve conservation of genetic variants critical to the future survival of coral ecosystems.

Low-vagility species may hold strong genetic signatures of past biogeographical processes but are also vulnerable to habitat loss. Flightless grasshoppers of the morabine group were once widespread in southeastern Australia, including Tasmania, but are becoming restricted to remnant patches of vegetation, with local ranges impacted by agriculture and development as well as management. Habitat fragmentation can generate genetically differentiated "island" populations with low genetic variation. However, following revegetation, populations could be re-established, and gene flow increased. Here we characterize single nucleotide polymorphism-based genetic variation in a widespread chromosomal race of the morabine Vandiemenella viatica (race 19) to investigate the genetic health of remnant populations and to provide guidelines for restoration efforts. We update the distribution of this race to new sites in Victoria and Tasmania, and show that V. viatica populations from northern Tasmania and eastern Victoria have reduced genetic variation compared to other mainland populations. In contrast, there was no effect of habitat fragment size on genetic variation. Tasmanian V. viatica populations fell into two groups, one connected genetically to eastern Victoria and the other connected to southwestern Victoria. Mainland populations showed isolation by distance. These patterns are consistent with expectations from past biogeographical processes rather than local recent population fragmentation and emphasize the importance of small local reserves in preserving genetic variation. The study highlights how genomic analyses can combine information on genetic variability and population structure to identify biogeographical patterns within a species, which in turn can inform decisions on potential source populations for translocations.

Network analysis is an effective tool to describe and quantify the ecological interactions between plants and root-associated fungi. Mycoheterotrophic plants, such as orchids, critically rely on mycorrhizal fungi for nutrients to survive, so investigating the structure of those intimate interactions brings new insights into the plant community assembly and coexistence. So far, there is little consensus on the structure of those interactions, described either as nested (generalist interactions), modular (highly specific interactions) or of both topologies. Biotic factors (e.g., mycorrhizal specificity) were shown to influence the network structure, while there is less evidence of abiotic factor effects. By using next-generation sequencing of the orchid mycorrhizal fungal (OMF) community associated to with plant individuals belonging to 17 orchid species, we assessed the structure of four orchid-OMF networks in two European regions under contrasting climatic conditions (Mediterranean vs. Continental). Each network contained four to 12 co-occurring orchid species, including six species shared among the regions. All four networks were both nested and modular, and fungal communities were different between co-occurring orchid species, despite multiple sharing of fungi across some orchids. Co-occurring orchid species growing in Mediterranean climate were associated with more dissimilar fungal communities, consistent with a more modular network structure compared to the Continental ones. OMF diversity was comparable among orchid species since most orchids were associated with multiple rarer fungi and with only a few highly dominant ones in the roots. Our results provide useful highlights into potential factors involved in structuring plant-mycorrhizal fungus interactions in different climatic conditions.

Deforestation threatens the integrity of the Amazon biome and the ecosystem services it provides, including greenhouse gas mitigation. Forest-to-pasture conversion has been shown to alter the flux of methane gas (CH4 ) in Amazonian soils, driving a switch from acting as a sink to a source of atmospheric CH4 . This study aimed to better understand this phenomenon by investigating soil microbial metagenomes, focusing on the taxonomic and functional structure of methane-cycling communities. Metagenomic data from forest and pasture soils were combined with measurements of in situ CH4 fluxes and soil edaphic factors and analysed using multivariate statistical approaches. We found a significantly higher abundance and diversity of methanogens in pasture soils. As inferred by co-occurrence networks, these microorganisms seem to be less interconnected within the soil microbiota in pasture soils. Metabolic traits were also different between land uses, with increased hydrogenotrophic and methylotrophic pathways of methanogenesis in pasture soils. Land-use change also induced shifts in taxonomic and functional traits of methanotrophs, with bacteria harbouring genes encoding the soluble form of methane monooxygenase enzyme (sMMO) depleted in pasture soils. Redundancy analysis and multimodel inference revealed that the shift in methane-cycling communities was associated with high pH, organic matter, soil porosity and micronutrients in pasture soils. These results comprehensively characterize the effect of forest-to-pasture conversion on the microbial communities driving the methane-cycling microorganisms in the Amazon rainforest, which will contribute to the efforts to preserve this important biome.

Acquisition of new genes often results in the emergence of novel functions and is a key step in lineage-specific adaptation. As a group of sessile crustaceans, barnacles establish permanent attachment through initial cement secretion at the larval phase followed by continuous cement secretion in juveniles and adults. However, the origins and evolution of barnacle larval and adult cement proteins remain poorly understood. By performing microdissection of larval cement glands, transcriptome and shotgun proteomics and immunohistochemistry validation, we identified 30 larval and 27 adult cement proteins of the epibiotic turtle barnacle Chelonibia testudinaria, of which the majority are stage- and barnacle-specific. While only two proteins, SIPC and CP100K, were expressed in both larvae and adults, detection of protease inhibitors and the cross-linking enzyme lysyl oxidase paralogs in larvae and adult cement. Other barnacle-specific cement proteins such as CP100k and CP52k likely share a common origin dating back at least to the divergence of Rhizocephala and Thoracica. Different CP52k paralogues could be detected in larval and adult cement, suggesting stage-specific cement proteins may arise from duplication followed by changes in expression timing of the duplicates. Interestingly, the biochemical properties of larval- and adult-specific CP52k paralogues exhibited remarkable differences. We conclude that barnacle larval and adult cement systems evolved independently, and both emerged from co-option of existing genes and de novo formation, duplication and functional divergence of lineage-specific cement protein genes. Our findings provide important insights into the evolutionary mechanisms of bioadhesives in sessile marine invertebrates.

A compelling interest in marine biology is to elucidate how species boundaries between sympatric free-spawning marine invertebrates such as bivalve molluscs are maintained in the face of potential hybridization. Hybrid zones provide the natural resources for us to study the underlying genetic mechanisms of reproductive isolation between hybridizing species. Against this backdrop, we examined the occurrence of introgressive hybridization (introgression) between two bivalves distributed in the western Pacific margin, Atrina japonica and Atrina lischkeana, based on single-nucleotide polymorphisms (SNPs) derived from restriction site-associated DNA sequencing. Using 1066 ancestry-informative SNP sites, we also investigated the extent of introgression within the genome to search for SNP sites with reduced interspecies gene flow. A series of our individual-level clustering analyses including the principal component analysis, Bayesian model-based clustering, and triangle plotting based on ancestry-heterozygosity relationships for an admixed population sample from the Seto Inland Sea (Japan) consistently suggested the presence of specimens with varying degrees of genomic admixture, thereby implying that the two species are not completely isolated. The Bayesian genomic cline analysis identified 10 SNP sites with reduced introgression, each of which was located within a genic region or an intergenic region physically close to a functional gene. No, or very few, heterozygotes were observed at these sites in the hybrid zone, suggesting that selection acts against heterozygotes. Accordingly, we raised the possibility that the SNP sites are within genomic regions that are incompatible between the two species. Our finding of restricted interspecies gene flow at certain genomic regions gives new insight into the maintenance of species boundaries in hybridizing broadcast-spawning molluscs.

Insect social parasites are characterized by exploiting the hosts' social behaviour. Why exactly hosts direct their caring behaviour towards these parasites and their offspring remains largely unstudied. One hypothesis is that hosts do not perceive their social environment as altered and accept the parasitic colony as their own. We used the ant Leptothorax acervorum, host of the dulotic, obligate social parasite Harpagoxenus sublaevis, to shed light on molecular mechanisms underlying behavioural exploitation by contrasting tissue-specific transcriptomes in young host workers. Host pupae were experimentally (re-)introduced into fragments of their original, another conspecific, heterospecific or parasitic colony. Brain and antennal mRNA was extracted and sequenced from adult ants after they had lived in the experimental colony for at least 50 days after eclosion. The resulting transcriptomes of L. acervorum revealed that ants were indeed affected by their social environment. Host brain transcriptomes were altered by the presence of social parasites, suggesting that the parasitic environment influences brain activity, which may be linked to behavioural changes. Transcriptional activity in the antennae changed most with the presence of unrelated individuals, regardless of whether they were conspecifics or parasites. This suggests early priming of odour perception, which was further supported by sensory perception of odour as an enriched function of differentially expressed genes. Furthermore, gene expression in the antennae, but not in the brain corresponded to ant worker behaviour before sampling. Our study demonstrated that the exploitation of social behaviours by brood parasites correlates with transcriptomic alterations in the central and peripheral nervous systems.

Infectious diseases can cause steep declines in wildlife populations, leading to changes in genetic diversity that may affect the susceptibility of individuals to infection and the overall resilience of populations to pathogen outbreaks. Here, we examine evidence for a genetic bottleneck in a population of American crows (Corvus brachyrhynchos) before and after the emergence of West Nile virus (WNV). More than 50% of marked birds in this population were lost over the 2-year period of the epizootic, representing a 10-fold increase in adult mortality. Using analyses of single-nucleotide polymorphisms (SNPs) and microsatellite markers, we tested for evidence of a genetic bottleneck and compared levels of inbreeding and immigration in the pre- and post-WNV populations. Counter to expectations, genetic diversity (allelic diversity and the number of new alleles) increased after WNV emergence. This was likely due to increases in immigration, as the estimated membership coefficients were lower in the post-WNV population. Simultaneously, however, the frequency of inbreeding appeared to increase: Mean inbreeding coefficients were higher among SNP markers, and heterozygosity-heterozygosity correlations were stronger among microsatellite markers, in the post-WNV population. These results indicate that loss of genetic diversity at the population level is not an inevitable consequence of a population decline, particularly in the presence of gene flow. The changes observed in post-WNV crows could have very different implications for their response to future pathogen risks, potentially making the population as a whole more resilient to a changing pathogen community, while increasing the frequency of inbred individuals with elevated susceptibility to disease.

Early events in the evolution of an ancestral lineage can shape the adaptive patterns of descendant species, but the evolutionary mechanisms driving initial adaptation from an ancestor remain largely unexplored. High-altitude adaptations have been extensively explored from the viewpoint of protein-coding genes; however, the contribution of noncoding regions remains relatively neglected. Here, we integrate genomic and transcriptomic data to investigate adaptive evolution in the ancestor of three high-altitude snowfinch species endemic to the Qinghai-Tibet Plateau. Our genome-wide scan for adaptation in the snowfinch ancestor identifies strong adaptation signals in functions of development and metabolism for the coding genes, but in functions of the nervous system development for noncoding regions. This pattern is exclusive to the snowfinch ancestor compared to a control ancestral lineage subject to weak selection. Changes in noncoding regions in the snowfinch ancestor, especially those nearest to coding genes, may be disproportionately associated with the differential expression of genes in the brain tissue compared to other tissues. Extensive gene expression in the brain tissue can be further altered via genetic regulatory networks of transcription factors harbouring potential accelerated regulatory regions (e.g., the development-related transcription factor YEATS4). Altogether, our study provides new evidence concerning how coding and noncoding sequences work through decoupled pathways in initial adaptation to the selective pressure of high-altitude environments. The analysis highlights the idea that noncoding sequences may be promising elements in facilitating the rapid evolution and adaptation to high altitudes.

The evergreen versus deciduous leaf habit is an important functional trait for adaptation of forest trees and has been hypothesized to be related to the evolutionary processes of the component species under paleoclimatic change, and potentially reflected in the dynamic history of evergreen broadleaved forests (EBLFs) in East Asia. However, knowledge about the shift of evergreen versus deciduous leaf with the impact of paleoclimatic change using genomic data remains rare. Here, we focus on the Litsea complex (Lauraceae), a key lineage with dominant species of EBLFs, to gain insights into how evergreen versus deciduous trait shifted, providing insights into the origin and historical dynamics of EBLFs in East Asia under Cenozoic climate change. We reconstructed a robust phylogeny of the Litsea complex using genome-wide single-nucleotide variants (SNVs) with eight clades resolved. Fossil-calibrated analyses, diversification rate shifts, ancestral habit, ecological niche modelling and climate niche reconstruction were employed to estimate its origin and diversification pattern. Taking into account studies on other plant lineages dominating EBLFs of East Asia, it was revealed that the prototype of EBLFs in East Asia probably emerged in the Early Eocene (55-50 million years ago [Ma]), facilitated by the greenhouse warming. As a response to the cooling and drying climate in the Middle to Late Eocene (48-38 Ma), deciduous habits were evolved in the dominant lineages of the EBLFs in East Asia. Up to the Early Miocene (23 Ma), the prevailing of East Asian monsoon increased the extreme seasonal precipitation and accelerated the emergence of evergreen habits of the dominant lineages, and ultimately shaped the vegetation resembling that of today.

The distribution of ecotypic variation in natural populations is influenced by neutral and adaptive evolutionary forces that are challenging to disentangle. This study provides a high-resolution portrait of genomic variation in Chinook salmon (Oncorhynchus tshawytscha) with emphasis on a region of major effect for ecotypic variation in migration timing. With a filtered data set of ~13 million single nucleotide polymorphisms (SNPs) from low-coverage whole genome resequencing of 53 populations (3566 barcoded individuals), we contrasted patterns of genomic structure within and among major lineages and examined the extent of a selective sweep at a major effect region underlying migration timing (GREB1L/ROCK1). Neutral variation provided support for fine-scale structure of populations, while allele frequency variation in GREB1L/ROCK1 was highly correlated with mean return timing for early and late migrating populations within each of the lineages (r2  = .58-.95; p < .001). However, the extent of selection within the genomic region controlling migration timing was much narrower in one lineage (interior stream-type) compared to the other two major lineages, which corresponded to the breadth of phenotypic variation in migration timing observed among lineages. Evidence of a duplicated block within GREB1L/ROCK1 may be responsible for reduced recombination in this portion of the genome and contributes to phenotypic variation within and across lineages. Lastly, SNP positions across GREB1L/ROCK1 were assessed for their utility in discriminating migration timing among lineages, and we recommend multiple markers nearest the duplication to provide highest accuracy in conservation applications such as those that aim to protect early migrating Chinook salmon. These results highlight the need to investigate variation throughout the genome and the effects of structural variants on ecologically relevant phenotypic variation in natural species.

The examination of genetic structure in the deep-ocean hadal zone has focused on divergence between tectonic trenches to understand how environment and geography may drive species divergence and promote endemism. There has been little attempt to examine localized genetic structure within trenches, partly because of logistical challenges associated with sampling at an appropriate scale, and the large effective population sizes of species that can be sampled adequately may mask underlying genetic structure. Here we examine genetic structure in the superabundant amphipod Hirondellea gigas in the Mariana Trench at depths of 8126-10,545 m. RAD sequencing was used to identify 3182 loci containing 43,408 single nucleotide polymorphisms (SNPs) across individuals after stringent pruning of loci to prevent paralogous multicopy genomic regions being erroneously merged. Principal components analysis of SNP genotypes resolved no genetic structure between sampling locations, consistent with a signature of panmixia. However, discriminant analysis of principal components identified divergence between all sites driven by 301 outlier SNPs in 169 loci and significantly associated with latitude and depth. Functional annotation of loci identified differences between singleton loci used in analysis and paralogous loci pruned from the data set and also between outlier and nonoutlier loci, all consistent with hypotheses explaining the role of transposable elements driving genome dynamics. This study challenges the traditional perspective that highly abundant amphipods within a trench form a single panmictic population. We discuss the findings in relation to eco-evolutionary and ontogenetic processes operating in the deep sea, and highlight key challenges associated with population genetic analysis in nonmodel systems with inherent large effective population sizes and genomes.

Phylogeographic studies uncover hidden pathways of divergence and inform conservation. Brown bears (Ursus arctos) have one of the broadest distributions of all land mammals, ranging from Eurasia to North America, and are an important model for evolutionary studies. Although several whole genomes were available for individuals from North America, Europe and Asia, limited whole-genome data were available from Central Asia, including the highly imperilled brown bears in the Gobi Desert. To fill this knowledge gap, we sequenced whole genomes from nine Asian brown bears from the Gobi Desert of Mongolia, Northern Mongolia and the Himalayas of Pakistan. We combined these data with published brown bear sequences from Europe, Asia and North America, as well as other bear species. Our goals were to determine the evolutionary relationships among brown bear populations worldwide, their genetic diversity and their historical demography. Our analyses revealed five major lineages of brown bears based on a filtered set of 684,081 single nucleotide polymorphisms. We found distinct evolutionary lineages of brown bears in the Gobi, Himalayas, northern Mongolia, Europe and North America. The lowest level of genetic diversity and the highest level of inbreeding were found in Pakistan, the Gobi Desert and Central Italy. Furthermore, the effective population size (Ne ) for all brown bears decreased over the last 70,000 years. Our results confirm the genetic distinctiveness and ancient lineage of brown bear subspecies in the Gobi Desert of Mongolia and the Himalayas of Pakistan and highlight their importance for conservation.

In organisms reproducing sexually, speciation occurs when increasing divergence results in pre- or post-zygotic reproductive isolation between lineages. Studies focusing on reproductive isolation origin in early stages of speciation are common and many rely on genomic scans to infer introgression providing limited information on the genomic architecture of reproductive isolation long-term maintenance. This study analyses a natural hybrid zone between two species in a late stage of speciation. We used ddRADseq genotyping in the contact between Podarcis bocagei and P. carbonelli to examine admixture extent, analyse hybrid zone stability and assess genome-wide variation in selection against introgression. We confirmed strong but incomplete reproductive isolation in a bimodal hybrid zone. New findings revealed population genetic structure within P. carbonelli in the contact zone; geographical and genomic clines analysis suggested strong selection against gene flow, but a relatively small proportion of the loci can introgress, mostly within the narrow contact zone. However, geographical clines revealed that a few introgressed loci show signs of potential positive selection, particularly into P. bocagei. Geographical clines also detected a signal of hybrid zone movement towards P. bocagei distribution. Genomic cline analysis revealed heterogeneous patterns of introgression among loci within the syntopy zone, but the majority maintain a strong association with the genomic background of origin. However, incongruences between both cline approaches were found, potentially driven by confounding effects on genomic clines. Last, an important role of the Z chromosome in reproductive isolation is suggested. Importantly, overall patterns of restricted introgression seem to result from numerous strong intrinsic barriers across the genome.

The distribution of genetic diversity is often heterogeneous in space, and it usually correlates with environmental transitions or historical processes that affect demography. The coast of Chile encompasses two biogeographic provinces and spans a broad environmental gradient together with oceanographic processes linked to coastal topography that can affect species' genetic diversity. Here, we evaluated the genetic connectivity and historical demography of four Scurria limpets, S. scurra, S. variabilis, S. ceciliana and S. araucana, between ca. 19° S and 53° S in the Chilean coast using genome-wide SNPs markers. Genetic structure varied among species which was evidenced by species-specific breaks together with two shared breaks. One of the shared breaks was located at 22-25° S and was observed in S. araucana and S. variabilis, while the second break around 31-34° S was shared by three Scurria species. Interestingly, the identified genetic breaks are also shared with other low-disperser invertebrates. Demographic histories show bottlenecks in S. scurra and S. araucana populations and recent population expansion in all species. The shared genetic breaks can be linked to oceanographic features acting as soft barriers to dispersal and also to historical climate, evidencing the utility of comparing multiple and sympatric species to understand the influence of a particular seascape on genetic diversity.

Many organisms possess multiple discrete genomes (i.e. nuclear and organellar), which are inherited separately and may have unique and even conflicting evolutionary histories. Phylogenetic reconstructions from these discrete genomes can yield different patterns of relatedness, a phenomenon known as cytonuclear discordance. In many animals, mitonuclear discordance (i.e. discordant evolutionary histories between the nuclear and mitochondrial genomes) has been widely documented, but its causes are often considered idiosyncratic and inscrutable. We show that a case of mitonuclear discordance in Todiramphus kingfishers can be explained by extensive genome-wide incomplete lineage sorting (ILS), likely a result of the explosive diversification history of this genus. For these kingfishers, quartet frequencies reveal that the nuclear genome is dominated by discordant topologies, with none of the internal branches in our consensus nuclear tree recovered in >50% of genome-wide gene trees. Meanwhile, a lack of inter-species shared ancestry, non-significant pairwise tests for gene flow, and little evidence for meaningful migration edges between species, leads to the conclusion that gene flow cannot explain the mitonuclear discordance we observe. This lack of evidence for gene flow combined with evidence for extensive genome-wide gene tree discordance, a hallmark of ILS, leads us to conclude that the mitonuclear discordance we observe likely results from ILS, specifically deep coalescence of the mitochondrial genome. Based on this case study, we hypothesize that similar demographic histories in other 'great speciator' taxa across the Indo-Pacific likely predispose these groups to high levels of ILS and high likelihoods of mitonuclear discordance.

Predicting how quickly populations expand their range and whether they will retain genetic diversity when they are introduced to new regions or track environmental conditions suited to their survival is an important applied and theoretical challenge. The literature suggests that long-distance dispersal, landscape heterogeneity and the evolution of dispersal influence populations' expansion rates and genetic diversity. We used individual-based spatially explicit simulations to examine these relationships for Tench (Tinca tinca), an invasive fish expanding its geographical range in eastern North America since the 1990s. Simulated populations varied greatly in expansion rates (1.1-28.6 patches year-1 ) and genetic diversity metrics, including changes in observed heterozygosity (-19 to +0.8%) and effective number of alleles (-0.32 to -0.01). Populations with greater dispersal distances expanded faster than those with smaller dispersal distances but exhibited considerable variation in expansion rate among local populations, implying less predictable expansions. However, they tended to retain genetic diversity as they expanded, suggesting more predictable evolutionary trajectories. In contrast, populations with smaller dispersal distances spread predictably more slowly but exhibited more variability among local populations in genetic diversity losses. Consistent with empirical data, populations spreading in a longer, narrower dispersal corridor lost more neutral genetic variation to the stochastic fixation of alleles. Given the unprecedented pace of anthropogenic environmental change and the increasing need to manage range-expanding populations, our results have conservation ramifications as they imply that the evolutionary trajectories of populations characterised by shorter dispersal distances spreading in narrower landscapes are more variable and, therefore, less predictable.

The processes governing soil bacteria biogeography are still not fully understood. It remains unknown how the importance of environmental filtering and dispersal differs between bacterial taxonomic and functional biogeography, and whether their importance is scale-dependent. We sampled soils across the Tibet plateau, with distances among plots ranging from 20 m to 1550 km. Taxonomic composition of bacterial community was characterized by 16S amplicon sequencing and functional community composition by qPCR targeting 9 functional groups involved in N dynamics. Factors representing climate, soil and plant community were measured to assess different facets of environmental dissimilarity. Both bacterial taxonomic and functional dissimilarities were more related to abiotic dissimilarity than biotic (vegetation) dissimilarity or distance. Taxonomic dissimilarity was mostly explained by differences in soil pH and mean annual temperature (MAT), while functional dissimilarity was linked to differences in soil N and P availabilities and N:P ratio. Soil pH and MAT remained the main determinants of taxonomic dissimilarity across spatial scales. In contrast, the explanatory variables of N-related functional dissimilarity varied across the scales, with soil moisture and organic matter having the highest role across short distances (<~330 km), and available P, N:P ratio and distance being important over long distances (>~660 km). Our results demonstrate how biodiversity dimension (taxonomic versus functional aspects) and spatial scale influence the factors driving soil bacterial biogeography.

Invasive species can impact native populations through competition, predation, habitat alteration, and disease transmission, but also genetically through hybridization. Potential outcomes of hybridization span the continuum from extinction to hybrid speciation and can be further complicated by anthropogenic habitat disturbance. Hybridization between the native green anole lizard (Anolis carolinensis) and a morphologically similar invader (A. porcatus) in south Florida provides an ideal opportunity to study interspecific admixture across a heterogeneous landscape. We used reduced-representation sequencing to describe introgression in this hybrid system and to test for a relationship between urbanization and non-native ancestry. Our findings indicate that hybridization between green anole lineages was probably a limited, historic event, producing a hybrid population characterized by a diverse continuum of ancestry proportions. Genomic cline analyses revealed rapid introgression and disproportionate representation of non-native alleles at many loci and no evidence for reproductive isolation between parental species. Three loci were associated with urban habitat characteristics; urbanization and non-native ancestry were positively correlated, although this relationship did not remain significant when accounting for spatial nonindependence. Ultimately, our study demonstrates the persistence of non-native genetic material even in the absence of ongoing immigration, indicating that selection favouring non-native alleles can override the demographic limitation of low propagule pressure. We also note that not all outcomes of admixture between native and non-native species should be considered intrinsically negative. Hybridization with ecologically robust invaders can lead to adaptive introgression, which may facilitate the long-term survival of native populations otherwise unable to adapt to anthropogenically mediated global change.

Speciation is often accompanied by frequent gene exchanges, which have been gradually recognized as a common phenomenon in nature. Although gene flow may influence different reproductive isolations, the specific mechanism of the process still requires more experimental evidence, especially in hybrid populations that have no significant differentiation and isolation. To address this challenge, this study aims to elucidate the underlying mechanisms of sympatry and parapatry in related species. Three sclerophyllous oaks (Quercus spinosa, Quercus aquifolioides and Quercus rehderiana), which are mainly distributed in the sympatry/parapatry of the East Himalaya-Hengduan Mountains and adjacent regions, were used to explore the population dynamics and evolutionary history. Based on 12,420 genome-wide single nucleotide polymorphism datasets, gene flow detection showed that the three species did not present any obvious genetic barriers. Evolutionary analysis indicated that the three species diverged during the Tertiary Period, and no migration events occurred in the early stages of species divergence. Combined with the data of 19 ecological factors, geological movements and climatic turbulence caused the rapid radiated differentiation of the three species during the Neocene, and similar selective pressures resulted in the same evolutionary pattern based on demographic history analysis. In addition, the predicted niche occupancy profiles and Generalized Dissimilarity Modelling revealed that the three species occupied distinct niches and had significant differences in ecological adaptation, which may explain the specific morphological characteristics of the different species. Therefore, we believe that the populations of the three related species underwent adaptive evolution in different habitats during the early stages of divergence. This study provides new experimental evidence of the formation patterns of parallel speciation.

Heat waves are predicted to be detrimental for organismal physiology with costs for survival that could be reflected in markers of biological state such as telomeres. Changes in early life telomere dynamics driven by thermal stress are of particular interest during the early post-natal stages of altricial birds because nestlings quickly shift from being ectothermic to endothermic after hatching. Telomeres of ectothermic and endothermic organisms respond differently to environmental temperature, but few investigations within species that transition from ectothermy to endothermy are available. Also, ambient temperature influences parental brooding behaviour, which will alter the temperature experienced by offspring and thereby, potentially, their telomeres. We exposed zebra finch nestlings to experimental heat waves and compared their telomere dynamics to that of a control group at 5, 12 and 80 days of age that encapsulate the transition from the ectothermic to the endothermic thermoregulatory stage; we also recorded parental brooding, offspring sex, mass, growth rates, brood size and hatch order. Nestling mass showed an inverse relationship with telomere length, and nestlings exposed to heat waves showed lower telomere attrition during their first 12 days of life (ectothermic stage) compared to controls. Additionally, parents of heated broods reduced the time they spent brooding offspring (at 5 days old) compared to controls. Our results indicate that the effect of heat waves on telomere dynamics likely varies depending on age and thermoregulatory stage of the offspring in combination with parental brooding behaviour during growth.

A fundamental goal of population genetic studies is to identify historical biogeographic patterns and understand the processes that generate them. However, localized demographic events can skew population genetic inference. Assessing populations with multiple types of genetic markers, each with unique mutation rates and responses to changes in population size, can help to identify potentially confounding population-specific demographic processes. Here, we compared population structure and connectivity inferred from microsatellites and restriction site-associated DNA loci among 17 populations of an arid-specialist lizard, the desert night lizard, Xantusia vigilis, in central California to test among historical processes structuring population genetic diversity. We found that both marker types yielded generally concordant insights into population genetic structure including a major phylogenetic break maintained between two populations separated by less than 10 km, suggesting that either marker type could be used to understand generalized demographic patterns across the region for management purposes. However, we also found that the effects of demography on marker discordance could be used to elucidate population histories and distinguish among competing biogeographic hypotheses. Our results suggest that comparisons of within-population diversity across marker types provide powerful opportunities for leveraging marker discordance, particularly for understanding the creation and maintenance of contact zones among clades.

Pesticide resistance represents a clear and trackable case of adaptive evolution with a strong societal impact. Understanding the factors associated with the evolution and spread of resistance is imperative to develop sustainable crop management strategies. The two-spotted spider mite Tetranychus urticae, a major crop pest with worldwide distribution and a polyphagous lifestyle, has evolved resistance to most classes of pesticides. Tetranychus urticae exists as either a green- or a red-coloured morph. However, the extent of genetic divergence and reproductive compatibility vary across populations of these colour morphs, complicating their taxonomic resolution at the species level. Here, we studied patterns of genetic differentiation and barriers to gene flow within and between morphs of T. urticae in order to understand the factors that influence the spread of resistance mutations across its populations. We derived multiple iso-female lines from Tetranychus populations collected from agricultural crops. We generated genomic and morphological data, characterized their bacterial communities and performed controlled crosses. Despite morphological similarities, we found large genomic differentiation between the morphs. This pattern was reflected in the incomplete, but strong postzygotic incompatibility in crosses between colour morphs, while crosses within morphs from different geographical locations were largely compatible. In addition, our results suggest recent/on-going gene flow between green-coloured T. urticae and T. turkestani. By screening the sequences of 10 resistance genes, we found evidence for multiple independent origins and for single evolutionary origins of target-site resistance mutations. Our results indicate that target-site mutations mostly evolve independently in populations on different geographical locations, and that these mutations can spread due to incomplete barriers to gene flow within and between populations.

The tropical Andes are one of the most important biodiversity hotspots on Earth, yet our understanding of how their biotas have responded to Quaternary climatic oscillations is extraordinarily limited and the alternative models proposed to explain their demographic dynamics have been seldom formally evaluated. Here, we test the hypothesis that the interplay between the spatial configuration of geographical barriers to dispersal and elevational displacements driven by Quaternary cooling-warming cycles has shaped the demographic trajectories of montane oak forests (Quercus humboldtii) from the Colombian Andes. Specifically, we integrate genomic data and environmental niche modelling at fine temporal resolution to test competing spatially explicit demographic and coalescent models, including scenarios considering (i) isotropic gene flow through the landscape, (ii) the hypothetical impact of contemporary barriers to dispersal (i.e., inter-Andean valleys), and (iii) distributional shifts of montane oak forests from the Last Glacial Maximum to the present. Although our data revealed a marked genetic fragmentation of montane oak forests, statistical support for isolation-with-migration models indicates that geographically separated populations from the different Andean Cordilleras regularly exchange gene flow. Accordingly, spatiotemporally explicit demographic analyses supported a model of flickering connectivity, with scenarios considering isotropic gene flow or currently unsuitable habitats as persistent barriers to dispersal providing a comparatively worse fit to empirical genomic data. Overall, these results emphasize the role of landscape heterogeneity on shaping spatial patterns of genomic variation in montane oak forests, rejecting the hypothesis of genetic continuity and supporting a significant impact of Quaternary climatic oscillations on their demographic trajectories.

Adaptation to local conditions is known to occur in seagrasses; however, knowledge of the genetic basis underlying this phenomenon remains scarce. Here, we analysed Posidonia oceanica from six sites within and around the Stagnone di Marsala, a semi-enclosed coastal lagoon where salinity and temperature exceed the generally described tolerance thresholds of the species. Sea surface temperatures (SSTs) were measured and plant samples were collected for the assessment of morphology, flowering rate and for screening genome-wide polymorphisms using double digest restriction-site-associated DNA sequencing. Results demonstrated more extreme SSTs and salinity levels inside the lagoon than the outer lagoon regions. Morphological results showed significantly fewer and shorter leaves and reduced rhizome growth of P. oceanica from the inner lagoon and past flowering events were recorded only for a meadow farthest away from the lagoon. Using an array of 51,329 single nucleotide polymorphisms, we revealed a clear genetic structure among the study sites and confirmed the genetic isolation and high clonality of the innermost site. In all, 14 outlier loci were identified and annotated with several proteins including those relate to plant stress response, protein transport and regulators of plant-specific developmental events. Especially, five outlier loci showed maximum allele frequency at the innermost site, likely reflecting adaptation to the extreme temperature and salinity regimes, possibly due to the selection of more resistant genotypes and the progressive restriction of gene flow. Overall, this study helps us to disentangle the genetic basis of seagrass adaptation to local environmental conditions and may support future works on assisted evolution in seagrasses.

Genomic structural variations (SVs) are widespread in plant and animal genomes and play important roles in phenotypic novelty and species adaptation. Frequent whole genome duplications followed by (re)diploidizations have resulted in high diversity of genome architecture among extant species. In this study, we identified abundant genomic SVs in the Panax genus that are hypothesized to have occurred through during the repeated polyploidizations/(re)diploidizations. Our genome-wide comparisons demonstrated that although these polyploidization-derived SVs have evolved at distinct evolutionary stages, a large number of SV-intersecting genes showed enrichment in functionally important pathways related to secondary metabolites, photosynthesis and basic cellular activities. In line with these observations, our metabolic analyses of these Panax species revealed high diversity of primary and secondary metabolites both at the tissue and interspecific levels. In particular, genomic SVs identified at ginsenoside biosynthesis genes, including copy number variation and large fragment deletion, appear to have played important roles in the evolution and diversification of ginsenosides. A further herbivore deterrence experiment demonstrated that, as major triterpenoidal saponins found exclusively in Panax, ginsenosides provide protection against insect herbivores. Our study provides new insights on how polyploidization-derived SVs have contributed to phenotypic novelty and plant adaptation.

Genetic differentiation between and within natural populations is the result of the joint effects of neutral and adaptative processes. In addition, the spatial arrangement of the landscape promotes connectivity or creates barriers to gene flow, directly affecting speciation processes. In this study, we carried out a landscape genomics analysis using NextRAD data from a montane forest specialist bird complex, the Mesoamerican Chestnut-capped/Green-striped Brushfinch of the genus Arremon. Specifically, we examined population genomic structure using different assignment methods and genomic differentiation and diversity, and we tested alternative genetic isolation hypotheses at the individual level (e.g., isolation by barrier, IBB; isolation by environment, IBE; isolation by resistance, IBR). We found well-delimited genomic structuring (K = 5) across Mesoamerican montane forests in the studied group. Individual-level genetic distances among major montane ranges were mainly explained by IBR hypotheses in this sedentary Neotropical taxon. Our results uncover genetic distances/differentiation and patterns of gene flow in allopatric species that support the role of tropical mountains as spatial landscape drivers of biodiversity. IBR clearly supports a pattern of conserved niche-tracking of suitable habitat conditions and topographic complexity throughout glacial-interglacial dynamics.

The repeated evolution of phenotypes provides clear evidence for the role of natural selection in driving evolutionary change. However, the evolutionary origin of repeated phenotypes can be difficult to disentangle as it can arise from a combination of factors such as gene flow, shared ancestral polymorphisms or mutation. Here, we investigate the presence of these evolutionary processes in the Hawaiian spiny-leg Tetragnatha adaptive radiation, which includes four microhabitat-specialists or ecomorphs, with different body pigmentation and size (Green, Large Brown, Maroon, and Small Brown). We investigated the evolutionary history of this radiation using 76 newly generated low-coverage, whole-genome resequenced samples, along with phylogenetic and population genomic tools. Considering the Green ecomorph as the ancestral state, our results suggest that the Green ecomorph likely re-evolved once, the Large Brown and Maroon ecomorphs evolved twice and the Small Brown evolved three times. We found that the evolution of the Maroon and Small Brown ecomorphs likely involved ancestral hybridization events, while the Green and Large Brown ecomorphs likely evolved through novel mutations, despite a high rate of incomplete lineage sorting in the dataset. Our findings demonstrate that the repeated evolution of ecomorphs in the Hawaiian spiny-leg Tetragnatha is influenced by multiple evolutionary processes.

The intra- and interspecific facets of biodiversity have traditionally been analysed separately, limiting our understanding of how evolution has shaped biodiversity, how biodiversity (as a whole) alters ecological dynamics and hence eco-evolutionary feedbacks at the community scale. Here, we propose using candidate genes phylogenetically-conserved across species and sustaining functional traits as an inclusive biodiversity unit transcending the intra- and interspecific boundaries. This framework merges knowledge from functional genomics and functional ecology, and we first provide guidelines and a concrete example for identifying phylogenetically-conserved candidate genes (PCCGs) within communities and for measuring biodiversity from PCCGs. We then explain how biodiversity measured at PCCGs can be linked to ecosystem functions, which unifies recent observations that both intra- and interspecific biodiversity are important for ecosystem functions. We then highlight the eco-evolutionary processes shaping PCCG diversity patterns and argue that their respective role can be inferred from concepts derived from population genetics. Finally, we explain how PCCGs may shift the field of eco-evolutionary dynamics from a focal-species approach to a more realistic focal-community approach. This framework provides a novel perspective to investigate the global ecosystem consequences of diversity loss across biological scales, and how these ecological changes further alter biodiversity evolution.

Adaptation to environmental conditions, and the mechanisms underlying these adaptations, can vary greatly among species. This variation can be attributed to a variety of factors including the strength of evolutionary processes like selection, gene flow, time since divergence, and/or genetic drift, as well as the interactions between these processes. A number of simulation and theoretical studies have helped elucidate the role of these processes on the genomic basis of adaptation (Schaal et al., 2022; Yeaman et al., 2016). However, complementary empirical studies to test these theoretical expectations for within-species adaptation have been limited due to the challenging nature of evaluating these processes in a comparative framework. To do this effectively, it is necessary to have systems where the range of environmental variation is similar between species, but where one or more of these evolutionary processes vary. In a From the Cover article in this issue of Molecular Ecology, Shi et al. (2022) provide an excellent example of a freshwater system where rates of gene flow differ between populations of six riverine species due to variation in spawning strategies (i.e., broadcast spawners = high gene flow, nest spawners = low gene flow), but all experience the same variation in environmental conditions across their distributions. The authors take a multivariate approach to evaluate the genomic basis of adaptation by using a combination of differentiation-based and genotype-environment association (GEA) methods. By comparing the amount of gene flow between species and the resulting genomic basis of local adaptation, they are able to infer how genomic architecture may be shaped by rates of gene flow. Their results identify a general pattern of increased genomic clustering in species with increasing levels of gene flow. However, two of six species did not follow this pattern, which could be due to additional factors not assessed. Additionally, they provide convincing evidence that the underlying evolutionary mechanisms that formed genomic clusters within each species vary. These deviations from a general pattern highlight how difficult evaluating these processes in natural populations are, particularly because species-specific responses can vary dramatically. Taken together, their comparative framework for assessing the genomic architecture of adaptation is unique, sheds important light on how evolutionary processes can impact adaptation, and provides robust empirical support of foundational theoretical and simulation studies.

The extent of parallel genomic responses to similar selective pressures depends on a complex array of environmental, demographic, and evolutionary forces. Laboratory experiments with replicated selective pressures yield mixed outcomes under controlled conditions and our understanding of genomic parallelism in the wild is limited to a few well-established systems. Here, we examine genomic signals of selection in the eelgrass Zostera marina across temperature gradients in adjacent embayments. Although we find many genomic regions with signals of selection within each bay there is very little overlap in signals of selection at the SNP level, despite most polymorphisms being shared across bays. We do find overlap at the gene level, potentially suggesting multiple mutational pathways to the same phenotype. Using polygenic models we find that some sets of candidate SNPs are able to predict temperature across both bays, suggesting that small but parallel shifts in allele frequencies may be missed by independent genome scans. Together, these results highlight the continuous rather than binary nature of parallel evolution in polygenic traits and the complexity of evolutionary predictability.

Hydrogen sulfide is a toxic gas that disrupts numerous biological processes, including energy production in the mitochondria, yet fish in the Poecilia mexicana species complex have independently evolved sulfide tolerance several times. Despite clear evidence for convergence at the phenotypic level in these fishes, it is unclear if the repeated evolution of hydrogen sulfide tolerance is the result of similar genomic changes. To address this gap, we used a targeted capture approach to sequence genes associated with sulfide processes and toxicity from five sulfidic and five nonsulfidic populations in the species complex. By comparing sequence variation in candidate genes to a reference set, we identified similar population structure and differentiation, suggesting that patterns of variation in most genes associated with sulfide processes and toxicity are due to demographic history and not selection. But the presence of tree discordance for a subset of genes suggests that several loci are evolving divergently between ecotypes. We identified two differentiation outlier genes that are associated with sulfide detoxification in the mitochondria that have signatures of selection in all five sulfidic populations. Further investigation into these regions identified long, shared haplotypes among sulfidic populations. Together, these results reveal that selection on standing genetic variation in putatively adaptive genes may be driving phenotypic convergence in this species complex.

Mountains are renowned for their bountiful biodiversity. Explanations on the origin of such abundant life are usually regarded to their orogenic history. However, ancient mountain systems with geological stability also exhibit astounding levels of number of species and endemism, as illustrated by the Brazilian Quartzitic Mountains (BQM) in Eastern South America. Thus, cycles of climatic changes over the last couple million years are usually assumed to play an important role in the origin of mountainous biota. These climatic oscillations potentially isolated and reconnected adjacent populations, a phenomenon known as flickering connectivity, accelerating speciation events due to range fragmentation, dispersion, secondary contact, and hybridization. To evaluate the role of the climatic fluctuations on the diversification of the BQM biota, we estimated the ancient demography of distinct endemic species of animals and plants using hierarchical approximate Bayesian computation analysis and Ecological Niche Modelling. Additionally, we evaluated if climatic oscillations have driven a genetic spatial congruence in the genetic structure of codistributed species from the Espinhaço Range, one of the main BQM areas. Our results show that the majority of plant lineages underwent a synchronous expansion over the Last Glacial Maximum (LGM, c. 21 thousand years ago), although we could not obtain a clear demographic pattern for the animal lineages. We also obtained a signal of a congruent phylogeographic break between lineages endemic to the Espinhaço Range, suggesting how ancient climatic oscillations might have driven the evolutionary history of the Espinhaço's biota.

Across the animal kingdom, newly independent juveniles form social associations that influence later fitness, mate choice and gene flow, but little is known about the ontogeny of social environments, particularly in wild populations. Here we test whether associations among young animals form randomly or are influenced by environmental or genetic conditions established by parents. Parents' decisions determine natal birth sites, which could affect who independent young initially encounter; secondly, mate choice determines genetic condition (e.g. inbreeding) of young and the parental care they receive, which can affect sociability. However, genetic and environmental factors are confounded unless related offspring experience different natal environments. Therefore, we used a long-term genetic pedigree, breeding records and social network data from three cohorts of a songbird with high extra-pair paternity (hihi, Notiomystis cincta) to disentangle (1) how nest location and relatedness contribute to association structure once juveniles disperse away from birth sites, and (2) if juvenile and/or parental inbreeding predicts individual sociability. We detected positive spatial autocorrelation: hihi that fledged closer by were more likely to associate even after dispersing, irrespective of genetic relatedness. Juvenile inbreeding did not predict sociability, but those raised by more inbred fathers formed more, stronger, associations, which did not depend on whether that male was the genetic parent or not. These results suggest that the natal environment created by parents, rather than focal genetic condition, establishes the foundation for social associations. Overall, we highlight how social inheritance may play an important role in population dynamics and evolutionary potential in wild animals.

Fire has shaped global ecosystems for millennia by directly killing organisms and indirectly altering habitats and resources. All terrestrial ecosystems, including fire-prone ecosystems, rely on soil-inhabiting fungi, where they play vital roles in ecological processes. Yet our understanding of how fire regimes influence soil fungi remains limited and our knowledge of these interactions in semiarid landscapes is virtually absent. We collected soil samples and vegetation measurements from sites across a gradient in time-since-fire ages (0-75 years-since-fire) and fire frequency (burnt 0-5 times during the recent 29-year period) in a semiarid heathland of south-eastern Australia. We characterized fungal communities using ITS amplicon-sequencing and assigned fungi taxonomically to trophic guilds. We used structural equation models to examine direct, indirect and total effects of time-since-fire and fire frequency on total fungal, ectomycorrhizal, saprotrophic and pathogenic richness. We used multivariate analyses to investigate how total fungal, ectomycorrhizal, saprotrophic and pathogenic species composition differed between post-fire successional stages and fire frequency classes. Time-since-fire was an important driver of saprotrophic richness; directly, saprotrophic richness increased with time-since-fire, and indirectly, saprotrophic richness declined with time-since-fire (resulting in a positive total effect), mediated through the impact of fire on substrates. Frequently burnt sites had lower numbers of saprotrophic and pathogenic species. Post-fire successional stages and fire frequency classes were characterized by distinct fungal communities, with large differences in ectomycorrhizal species composition. Understanding the complex responses of fungal communities to fire can be improved by exploring how the effects of fire flow through ecosystems. Diverse fire histories may be important for maintaining the functional diversity of fungi in semiarid regions.

The impact of post-divergence gene flow in speciation has been documented across a range of taxa in recent years, and may have been especially widespread in highly mobile, wide-ranging marine species, such as cetaceans. Here, we studied individual genomes from nine species across the three families of the toothed whale superfamily Delphinoidea (Delphinidae, Phocoenidae and Monodontidae). To investigate the role of post-divergence gene flow in the speciation process, we used a multifaceted approach, including (i) phylogenomics, (ii) the distribution of shared derived alleles and (iii) demographic inference. We found the divergence of lineages within Delphinoidea did not follow a process of pure bifurcation, but was much more complex. Sliding-window phylogenomics reveal a high prevalence of discordant topologies within the superfamily, with further analyses indicating these discordances arose due to both incomplete lineage sorting and gene flow. D-statistics and f-branch analyses supported gene flow between members of Delphinoidea, with the vast majority of gene flow occurring as ancient interfamilial events. Demographic analyses provided evidence that introgressive gene flow has likely ceased between all species pairs tested, despite reports of contemporary interspecific hybrids. Our study provides the first steps towards resolving the large complexity of speciation within Delphinoidea; we reveal the prevalence of ancient interfamilial gene flow events prior to the diversification of each family, and suggest that contemporary hybridisation events may be disadvantageous, as hybrid individuals do not appear to contribute to the parental species' gene pools.

The current advances of environmental DNA (eDNA) bring profound changes to ecological monitoring and provide unique insights on the biological diversity of ecosystems. The very nature of eDNA data is challenging yet also revolutionizing how biological monitoring information is analysed. In particular, new metrics and approaches should take full advantage of the extent and detail of molecular data produced by genetic methods. In this perspective, machine learning algorithms are particularly promising as they can capture complex relationships between the multiple environmental pressures and the diversity of biological communities. We investigated the potential of a new generation of biomonitoring tools that implement machine-learning techniques to fully exploit eDNA datasets. We trained a machine learning model to discriminate between reference and impacted communities of freshwater macroinvertebrates and assessed its performances using a large eDNA dataset collected at 64 standard federal monitoring sites across Switzerland. We show that a model trained on eDNA is significantly better than a naive model and performs similarly to a model trained on traditional data. Our proof-of-concept shows that such a combination of eDNA and machine learning approaches has the potential to complement or even replace traditional environmental monitoring, and could be scaled along temporal or spatial dimensions.

The broad scale distribution of population-specific genetic diversity (GDP ) across taxa remains understudied relative to species diversity gradients, despite its relevance for systematic conservation planning. We used nuclear DNA data collected from 3678 vertebrate populations across the Americas to assess the role of environmental and spatial variables in structuring the distribution of GDP , a key component of adaptive potential in the face of environmental change. We specifically assessed non-linear trends for a metric of GDP, expected heterozygosity (HE ), and found more evidence for spatial hotspots and cold spots in HE rather than a strict pattern with latitude. We also detected inconsistent relationships between HE and environmental variables, where only 11 of 30 environmental comparisons among taxa groups were statistically significant at the .05 level, and the shape of significant trends differed substantially across vertebrate groups. Only one of six taxonomic groups, freshwater fishes, consistently showed significant relationships between HE and most (four of five) environmental variables. The remaining groups had statistically significant relationships for either two (amphibians, reptiles), one (birds, mammals), or no variables (anadromous fishes). Our study highlights gaps in the theoretical foundation upon which macrogenetic predictions have been made thus far in the literature, as well as the nuances for assessing broad patterns in GDP among vertebrate groups. Overall, our results suggest a disconnect between patterns of species and genetic diversity, and underscores that large-scale factors affecting genetic diversity may not be the same factors as those shaping taxonomic diversity. Thus, careful spatial and taxonomic-specific considerations are needed for applying macrogenetics to conservation planning.

With continued global change, recovery of species listed under the Endangered Species Act is increasingly challenging. One rare success was the recovery and delisting of the Channel Island fox (Urocyon littoralis) after 90%-99% population declines in the 1990s. While their demographic recovery was marked, less is known about their genetic recovery. To address genetic changes, we conducted the first multi-individual and population-level direct genetic comparison of samples collected before and after the recent bottlenecks. Using whole-exome sequencing, we found that already genetically depauperate populations were further degraded by the 1990s declines and remain low, particularly on San Miguel and Santa Rosa Islands, which underwent the most severe bottlenecks. The two other islands that experienced recent bottlenecks (Santa Cruz, and Santa Catalina islands) showed mixed results based on multiple metrics of genetic diversity. Previous island fox genomics studies showed low genetic diversity before the declines and no change after the demographic recovery, thus this is the first study to show a decrease in genetic diversity over time in U. littoralis. Additionally, we found that divergence between populations consistently increased over time, complicating prospects for using inter-island translocation as a conservation tool. The Santa Catalina subspecies is now federally listed as threatened, yet other de-listed subspecies are still recovering genetic variation which may limit their ability to adapt to changing environmental conditions. This study further demonstrates that species conservation is more complex than population size and that some island fox populations are not yet 'out of the woods'.

Understanding how interactions among microevolutionary forces generate genetic population structure of exploited species is vital to the implementation of management policies that facilitate persistence. Philopatry displayed by many coastal shark species can impact gene flow and facilitate selection, and has direct implications for the spatial scales of management. Here, genetic structure of the blacktip shark (Carcharhinus limbatus) was examined using a mixed-marker approach employing mitochondrial control region sequences and 4339 SNP-containing loci generated using ddRAD-Seq. Genetic variation was assessed among young-of-the-year sampled in 11 sites in waters of the United States in the western North Atlantic Ocean, including the Gulf of Mexico. Spatial and environmental analyses detected 68 nuclear loci putatively under selection, enabling separate assessments of neutral and adaptive genetic structure. Both mitochondrial and neutral SNP data indicated three genetically distinct units-the Atlantic, eastern Gulf, and western Gulf-that align with regional stocks and suggest regional philopatry by males and females. Heterogeneity at loci putatively under selection, associated with temperature and salinity, was observed among sites within Gulf units, suggesting local adaptation. Furthermore, five pairs of siblings were identified in the same site across timescales corresponding with female reproductive cycles. This indicates that females re-used a site for parturition, which has the potential to facilitate the sorting of adaptive variation among neighbouring sites. The results demonstrate differential impacts of microevolutionary forces at varying spatial scales and highlight the importance of conserving essential habitats to maintain sources of adaptive variation that may buffer species against environmental change.