Reconstructing the geographic origins of invasive species is critical for establishing effective management strategies. Frequently, molecular investigations are undertaken when the source population is not known, however; these analyses are constrained both by the amount of diversity present in the native region and by changes in the genetic background of the invading population following bottlenecks and/or hybridization events. Here we explore the geographical origins of the invasive winter moth (Operopthera brumata L.) that has caused widespread defoliation to forests, orchards, and crops in four discrete regions: Nova Scotia, British Columbia, Oregon, and the northeastern United States. It is not known whether these represent independent introductions to North America, or “stepping stone” spread among regions. Using a combination of Bayesian assignment and approximate Bayesian computation methods, we analyzed a population genetic dataset of 24 polymorphic microsatellite loci. We estimate that winter moth was introduced to North America on at least four occasions, with the Nova Scotian and British Columbian populations likely being introduced from France and Sweden, respectively; the Oregonian population likely being introduced from either the British Isles or northern Fennoscandia; and the population in the northeastern United States likely being introduced from somewhere in Central Europe. To our surprise, we found that hybridization has not played a large role in the establishment of winter moth populations even though previous reports have documented widespread hybridization between winter moth and a native congener. We discuss the impact of genetic bottlenecks on analyses meant to determine region of origin.

Natural hybrid zones provide opportunities for studies of the evolution of reproductive isolation in wild populations. Although recent investigations have found that the formation of neo-sex chromosomes is associated with reproductive isolation, the mechanisms remain unclear in most cases. Here, we assess the contemporary structure of gene flow in the contact zone between largely allopatric cytotypes of the dioecious plant Rumex hastatulus, a species with evidence of sex chromosome turn-over. Males to the west of the Mississippi river, USA, have an X and a single Y chromosome, whereas populations to the east of the river have undergone a chromosomal rearrangement giving rise to a larger X and two Y chromosomes. Using reduced-representation sequencing, we provide evidence that hybrids form readily and survive multiple backcross generations in the field, demonstrating the potential for ongoing gene flow between the cytotypes. Cline analysis of each chromosome separately captured no signals of difference in cline shape between chromosomes. However, principal component regression revealed a significant increase in the contribution of individual SNPs to inter-cytotype differentiation on the neo-X chromosome, but no correlation with recombination rate. Cline analysis revealed that the only SNPs with significantly shallower clines than the genome-average were located on the neo-X. Our data are consistent with a role for neo-sex chromosomes in reproductive isolation between R. hastatulus cytotypes. Our investigation highlights the importance of studying plant hybrid zones in species with sex chromosomes for understanding mechanisms of reproductive isolation and the role of gene flow in the spread of neo-X chromosomes.

Taxon-specific characteristics and extrinsic climatic and geological forces may both shape population differentiation and speciation. In geographically and taxonomically focused investigations, differentiation may occur synchronously as species respond to the same external conditions. Conversely, when evolution is investigated in taxa with largely varying traits, population differentiation and speciation is complex and shaped by interactions of Earth’s template and species-specific traits. As such, it is important to characterize evolutionary histories broadly across the tree of life, especially in geographic regions that are exceptionally diverse and under pressures from human activities such as in biodiversity hotspots. Here, using whole-genome sequencing data, we characterize genomic variation in populations of six Ethiopian Highlands forest bird species separated by a lowland biogeographic barrier, the Great Rift Valley (GRV). In all six species, populations on either side of the GRV exhibited significant but varying levels of genetic differentiation. Species’ dispersal ability was negatively correlated with levels of population differentiation. Isolation with migration models indicated varied patterns of population differentiation and connectivity among populations of the focal species. We found that demographic histories—estimated for each individual—varied by both species and population but were consistent between individuals of the same species and sampling region. We found that genomic diversity varied by half an order of magnitude across species, and that this variation could largely be explained by the harmonic mean of effective population size over the past 200,000 years. Overall, we found that even in highly dispersive species like birds, the GRV acts as a substantial biogeographic barrier.

Introductions of invasive species to new environments often result in rapid rates of trait evolution. While in some cases these evolutionary transitions are adaptive and driven by natural selection, they can also result from patterns of genetic and phenotypic variation associated with the invasion history. Here, we examined the brown anole (Anolis sagrei), a widespread invasive lizard for which genetic data have helped trace the sources of non-native populations. We focused on the dewlap, a complex signaling trait known to be subject to multiple selective pressures. We measured dewlap reflectance, pattern, and size in 30 non-native populations across the southeastern United States. As well, we quantified environmental variables known to influence dewlap signal effectiveness, such as canopy openness. Further, we used genome-wide data to estimate genetic ancestry, perform association mapping, and test for signatures of selection. We found that among-population variation in dewlap characteristics was best explained by genetic ancestry. This result was supported by genome-wide association mapping, which identified several ancestry-specific loci associated with dewlap traits. Despite the strong imprint of this aspect of the invasion history on dewlap variation, we also detected significant relationships between dewlap traits and local environmental conditions. However, we found limited evidence that dewlap-associated genetic variants have been subject to selection. Our study emphasizes the importance of genetic ancestry and admixture in shaping phenotypes during biological invasion, while leaving the role of selection unresolved, likely due to the polygenic genetic architecture of dewlaps and selection acting on many genes of small effect.

Stony coral tissue loss disease (SCTLD) remains an unprecedented disease outbreak due to its high mortality rate and rapid spread throughout Florida’s Coral Reef and wider Caribbean. A collaborative effort is underway to evaluate disease intervention strategies that mitigate the spread of SCTLD across coral colonies and reefs. We conducted an in-situ experiment in Southeast Florida to assess molecular responses among SCTLD-affected Montastraea cavernosa pre- and post-application of the most widely-used intervention method, CoreRx Base 2B with amoxicillin. Through Tag-Seq gene expression profiling of apparently healthy, diseased, and treated corals, we identified modulation of metabolomic and immune pathways following antibiotic treatment. In a complementary ex-situ disease challenge experiment, we exposed nursery-cultured M. cavernosa and Orbicella faveolata fragments to SCTLD-affected donor corals to compare transcriptomic profiles among clonal individuals from unexposed controls, those exposed and displaying disease signs, and corals exposed and not displaying disease signs. Suppression of metabolic functional groups and activation of stress gene pathways as a result of SCTLD exposure were apparent in both species. Amoxicillin treatment led to a ‘reversal’ of the majority of gene pathways implicated in disease response, suggesting potential recovery of corals following antibiotic application. In addition to increasing our understanding of molecular responses to SCTLD, we provide resource managers with transcriptomic evidence that disease interventions with antibiotics appear to be successful and may help to modulate coral immune responses to SCTLD. These results contribute to feasibility assessments of intervention efforts following disease outbreaks and improved predictions of coral reef health in Southeast Florida.

Trade-offs between traits arise and reflect constraints imposed by the environment and physicochemical laws. Trade-off situations are expected to be highly relevant for sessile plants, which have to respond to changes in the environment to ensure survival. Despite increasing interest in determining the genetic and molecular basis of plant trade-offs, there are still gaps and differences with respect to how trade-offs are defined, how they are measured, and how their genetic architecture is dissected. The first step to fill these gaps is to establish what is meant by trade-offs. In this review we provide a classification of the existing definitions of trade-offs according to: (1) the measures used for their quantification, (2) the dependence of trade-offs on environment, and (3) whether data based on which they are inferred are from a single individual across different environments or a population of individuals in single or multiple environments. We then compare the approaches for quantification of trade-offs based on phenotypic, between-individual, and genetic correlations, and stress the need for developing further quantification indices particularly for trade-offs between multiple traits. Lastly, we highlight the genetic mechanisms underpinning trade-offs and experimental designs that facilitate their discovery in plants, with focus on usage of natural variability. This review also offers a perspective for future research aimed at identification of plant trade-offs, dissection of their genetic architecture, and development of strategies to overcome trade-offs, with applications in crop breeding.

Corals show spatial acclimatisation to local environment conditions. However, the various cellular mechanisms involved in local acclimatisation and variable bleaching patterns in corals remain to be thoroughly understood. In this study, the modulation of a protein implicated in cellular heat stress tolerance, the Heat shock protein 70, was compared at both gene (hsp70) and protein (Hsp70) expression level in bleaching tolerant near-coast Acropora muricata colonies and bleaching susceptible reef colonies, in the lagoon of Belle Mare (Mauritius). The relative Hsp70 levels varied significantly between colonies from the two different locations, colonies having different health conditions and the year of collection. Before the bleaching event of 2016, near-coast colonies had higher basal levels of both Hsp70 gene and protein compared to reef colonies. During the bleaching event, the near-coast colonies did not bleach and had significantly higher relative levels of both Hsp70 gene and protein compared to bleached reef colonies. No significant genetic differentiation between the two studied coral populations was observed and all the colonies analysed were associated with Symbiodiniaceae of the genus Symbiodinium (Clade A) irrespective of location and sampling period. These findings provide further evidence of the involvement of Hsp70 in conferring bleaching tolerance to corals. Moreover, the consistent expression differences of Hsp70 gene and protein between the near-coast and reef coral populations in a natural setting indicate that the modulation of this Hsp is involved in local acclimatisation of corals to their environments.

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 examine the role of evolution in a rapid range expansion whose ecological context has been well-documented, Anna’s Hummingbird (Calypte anna). Previous work suggests 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 examine the role of gene flow and adaptation during range expansion from their native California 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 few (if any) limits 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 extended 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.

Understanding how eco-evolutionary processes and environmental factors drive population differentiation and adaptation are key challenges in evolutionary biology of relevance for biodiversity protection. Differentiation requires at least partial reproductive separation, which may result from different modes of isolation such as geographic isolation (allopatry) or isolation by distance (IBD), resistance (IBR), and environment (IBE). Despite that multiple modes might jointly influence differentiation, studies that compare the relative contributions are scarce. Using RADseq, we analyse neutral and adaptive genetic diversity and structure in 11 pike (Esox lucius) populations along a latitudinal gradient (54.9 - 63.6°N), to investigate the relative effects of IBD, IBE and IBR, and to assess whether the effects differ between neutral and adaptive variation, or across structural levels. Patterns of neutral and adaptive variation differed, likely reflecting that they have been differently affected by stochastic and deterministic processes. The importance of the different modes of isolation differed between neutral and adaptive diversity, yet were consistent across structural levels. Neutral variation was influenced by interactions among all three modes of isolation, with IBR (seascape features) playing a central role, wheares adaptive variation was mainly influenced by IBE (environmental conditions). Taken together, this and previous studies suggest that it is common that multiple modes of isolation interactively shape patterns of genetic variation, and that their relative contributions differ among systems. To enable identification of general patterns and understand how various factors influence the relative contributions, it is important that several modes are simultaneously investigated in additional populations, species and environmental settings.

Pathogen-mediated selection (PMS) is thought to maintain the high level of allelic diversity observed in the major histocompatibility complex (MHC) class II genes. A comprehensive way to demonstrate contemporary selection is to examine associations between MHC variation and individual fitness. As individual fitness is hard to measure, many studies examine associations between MHC diversity and phenotypic traits which include direct or indirect measures of adaptive immunity thought to contribute to fitness. Here, we tested associations between MHC class II variation and five representative phenotypic traits measured in August: weight, strongyle faecal egg count, and plasma IgA, IgE and IgG immunoglobulin titres against the gastrointestinal nematode parasite Teladorsagia circumcincta in a free-living population of Soay sheep. We found no association between MHC class II variation and August weight or strongyle faecal egg count. We did however find associations between MHC class II variation and immunoglobulin levels which varied with age, isotype and sex. Our results suggest associations between MHC and phenotypic traits are more likely to be found for traits more closely associated with pathogen defence than integrative traits such as body weight and highlight a useful role of MHC-antibody associations in examining selection on MHC genes.

Long term environmental variation often drives local adaptation and leads to trait differentiation across populations. Additionally, when traits change in an environment-dependent way through phenotypic plasticity, the underlying genetic variation will also be under selection, but only in the inducing environment. Both of these processes will create a landscape of differentiation across populations in trait means as well as their plasticity. However, studies uncovering environmental drivers of this variation are scarce. With this work, we studied drought responses in seedlings of a shrub species from the Cape Floristic Region, the common sugarbush (Protea repens). We measured morphological and physiological traits as well as whole transcriptomes in 8 populations that represent both the climatic and the geographic distribution of this species. We found that there is substantial variation in how populations respond to drought, but we also observe common patterns such as reduced leaf size and thickness and upregulation of stress- and down-regulation of growth-related gene groups. Both environmental heterogeneity and milder source site climates were associated with higher plasticity in various traits and co-expression gene networks. By uncovering associations between traits, trait plasticity, co-expression gene networks with source site climate, we showed that temperature plays a bigger role in shaping these patterns when compared to precipitation, in line with recent changes in the region due to climate change. We also found that traits respond to climatic variation in a context dependent manner: some associations between traits and climate were apparent only under certain growing conditions.

Molecular Ecology News and Views - Perspective on Jacobs and Elmer (2021) “Alternative splicing and gene expression play contrasting roles in the parallel phenotypic evolution of a salmonid fish”

Studies of natural hybrid zones can provide documentation of range shifts in response to climate change and identify loci important to reproductive isolation. Using a deep temporal (36-38 years) comparison of the black-capped (Poecile atricapillus) and Carolina (P. carolinensis) chickadee hybrid zone, we investigated movement of the western portion of the zone (western Missouri) and assessed whether loci and pathways underpinning reproductive isolation were similar to those in the eastern portion of the hybrid zone. Using 92 birds sampled along the hybrid zone transect in 2016 and 68 birds sampled between 1978 and 1980, we generated 11,669 SNPs via ddRADseq. These SNPs were used to assess movement of the hybrid zone through time and to evaluate variation in introgression among loci. We demonstrate that the interface has moved ~5 km to the northwest over the last 36-38 years, i.e., at only one-fifth the rate at which the eastern portion (e.g., Pennsylvania, Ohio) of the hybrid zone has moved. Temperature trends over the last 38 years reveal that eastern areas have warmed 50% more than western areas in terms of annual mean temperature, possibly providing an explanation for the slower movement of the hybrid zone in Missouri. Our results suggest hybrid zone movement in broadly distributed species, such as chickadees, will vary between areas in response to local differences in the impacts of climate change.

DNA metabarcoding has been widely used to access and monitor species. However, several challenges remain open for its mainstream application in ecological studies, particularly when dealing with a quantitative approach. In a from the Cover article in this issue of Molecular Ecology, Cédric et al. (2021) report species-level ichthyoplankton dynamics for 97 fish species from two Amazon river basins using a clever quantitative metabarcoding approach employing a probe capture method. They clearly show that most species spawned during the rainy season when the floods started, but interestingly, species from the same genus reproduced in distinct periods (i.e., inverse phenology). Opportunistically, Cédric et al. (2021) reported that during an intense hydrological anomaly, several species had a sharp reduction in spawning activity, demonstrating a quick response to environmental cues. This is an interesting result since the speed at which fish species can react to environmental changes, during the spawning period, is largely unknown. Thus, this study brings remarkable insights into basic life history information that is imperative for proposing strategies that could lead to a realistic framework for sustainable fisheries management practices and conservation, fundamental for an under-studied and threatened realm, such as the Amazon River basin.

Inducible defenses of prey are evolved under diverse and variable predation risks. However, during the co-evolution of prey and multiple predators, the responses of prey to antagonistic predation risks, which may put the prey into a dilemma of responding to predators, remain unclear. Based on antagonistic predation pressure from an invertebrate (Chaoborus larvae) and a vertebrate (Rhodeus ocellatus) predator, we studied the responses of multiple traits and transcriptomes of the freshwater crustacean Ceriodaphnia cornuta under multiple predation risks. Chaoborus predation risk altered the expression of genes encoding cuticle proteins and modulated the biosynthesis of steroid hormones, cutin, suberine, and wax, leading to the development of horns and increase in size at the late developmental stage. Meanwhile, fish predation risk primarily triggered genes encoding ribosomes and those involved in unsaturated fatty acid biosynthesis and cysteine and methionine metabolism, resulting in smaller individual size and earlier reproduction. Inducible responses of both transcriptome and individual traits revealed that predator-dependent unique responses were dominant and the dilemma of antagonistic responses was relatively limited. However, the unique individual traits in response to invertebrate predation could be significantly impaired by vertebrate predation risk, even though the unique responses to different predators were extremely weakly correlated and could be elicited simultaneously. These results indicate that diverse predator-dependent unique responses are favored by Ceriodaphnia during its co-evolution with multiple predators. Nonetheless, Ceriodaphnia is not a generalist that can fully adopt all predator-dependent unique responses simultaneously under multiple predation risks.

Monarch butterflies are known for their spectacular annual migration in eastern North America, with millions of monarchs flying up to 4,500 kilometers to overwintering sites in central Mexico. Monarchs also live west of the Rocky Mountains, where they travel shorter distances to overwinter along the Pacific Coast. It is often assumed that eastern and western monarchs form distinct evolutionary units, but genomic studies to support this notion are lacking. We used a tethered flight mill to show that migratory eastern monarchs have greater flight performance than western monarchs, consistent with their greater migratory distances. However, analyzing more than 20 million SNPs in 43 monarch genomes, we found no evidence for genomic differentiation between eastern and western monarchs. Genomic analysis also showed identical and low levels of genetic diversity, and demographic analyses indicated similar effective population sizes and ongoing gene flow between eastern and western monarchs. Gene expression analysis of a subset of candidate genes during active flight revealed differential gene expression related to non-muscular motor activity. Our results demonstrate that eastern and western monarchs maintain migratory differences despite ongoing gene flow, and suggest that migratory differences between eastern and western monarchs are not driven by select major-effects alleles. Instead, variation in migratory distance and destination may be driven by environmentally induced differential gene expression, or by many alleles of small effect.

Host switching allows parasites to expand their niches. However, successful switching may require suites of adaptations and also may decrease performance on the old host. As a result, reductions in gene flow accompany many host switches, driving speciation. Because host switches tend to be rapid, it is difficult to study their demographic parameters in real-time. Fundamental factors that control subsequent parasite evolution, such as the size of the switching population or the extent of immigration from the original host, remain largely unknown. To shed light on the host switching process, we explored the history of independent switches by two ectoparasitic honey bee mites (Varroa destructor and V. jacobsoni). Both switched to the western honey bee (Apis mellifera) after it was brought into contact with their ancestral host (Apis cerana), ~70 and ~12 years ago, respectively. Varroa destructor subsequently caused worldwide collapses of honey bee populations. Using whole-genome sequencing on 63 mites collected in their native ranges from both the ancestral and novel hosts, we were able to reconstruct the known temporal dynamics of the switch. We further found multiple previously undiscovered mitochondrial lineages on the novel host, along with the genetic equivalent of tens of individuals that were involved in the initial host switch. Despite being greatly reduced, some gene flow remains between mites adapted to different hosts. Our findings suggest that while reproductive isolation may facilitate the fixation of traits beneficial for exploitation of the new host, ongoing genetic exchange may allow genetic amelioration of inbreeding effects.

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. Many rely on indirect measures of introgression providing limited information on the genomic architecture of reproductive isolation maintenance in the long term. This study focuses on direct measures of introgression between two species in a late stage of speciation. We used ddRADseq genotyping in a natural hybrid zone between Podarcis bocagei and P. carbonelli to examine admixture extent, estimate effective selection, analyse hybrid zone stability, and assess variation in selection against introgression across the genome. Hybridization was confirmed the narrow and bimodal hybrid zone demonstrating the existence of strong mechanisms of reproductive isolation. Simulations suggested that simple premating barriers were not enough to explain the observed distribution of admixture classes, pointing out the role of post-mating isolation. A geographic cline approach confirmed strong reproductive isolation and high effective selection preventing extensive introgression outside of the contact zone; and 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 most of the loci do not introgress more or less than the genomic average maintaining a strong association with the genomic background of origin. However, genomic clines can be driven by confounding effects resulting in incongruences between both cline approaches. Importantly, overall patterns of introgression seem to result from strong intrinsic barriers across the genome, without clear “islands of differentiation”. Last, an important role of the Z chromosome in reproductive isolation is suggested.