Fluctuations in temperature are recognized as a potent driver of selection pressure, fostering genomic variations that are crucial for the adaptation and survival of organisms under selection. Notably, water temperature is a pivotal factor influencing aquatic organism persistence. By comprehending how aquatic organisms respond to shifts in water temperature, we can understand their potential physiological adaptations to environmental change in one or multiple species. This, in turn, contributes to the formulation of biologically relevant guidelines for the landscape scale transcriptome profile of organisms in lotic systems. Here, we investigated the distinct responses of seven stream stonefly species, collected from four geographical regions across Japan, to variations in temperature, including atmospheric and water temperatures. We achieved this by assessing the differences in gene expression through RNA-sequencing within individual species and exploring the patterns of community-genes among different species. We identified 735 genes that exhibited differential expressions across the temperature gradient. Remarkably, the community displayed expression levels differences of respiration and metabolic genes. Additionally, the diversity in molecular functions appeared to be linked to spatial variation, with water temperature differences potentially contributing to the overall functional diversity of genes. We found 22 community-genes with consistent expression patterns among species in response to water temperature variations. These genes related to respiration, metabolism and development exhibited a clear gradient providing robust evidence of divergent adaptive responses to water temperature. Our findings underscore the differential adaptation of stonefly species to local environmental conditions, suggesting that shared responses in gene expression may occur across multiple species under similar environmental conditions. This study emphasizes the significance of considering various species when assessing the impacts of environmental changes on aquatic insect communities and understanding potential mechanisms to cope with such changes.

OBJECTIVE: Intraspecific variation in drought resistance traits, such as drought escape, appear to be frequent within wild, ruderal forb species. Understanding how these traits are arrayed across the landscape, particularly in association with climate, is critical to developing forbs for wildland restoration programs. Use of forbs is requisite for maintaining biological diversity and ecological services.
METHODS: Using 6074 greenhouse-grown Chaenactis douglasii seedlings from 95 wild, seed-sourced populations across the western United States, we recorded bolting phenology and estimated genome size using flow cytometry. Mixed-effects regression models were used to assess whether climate of seed origin was predictive for bolting phenology and genome size.
RESULTS: Variation in bolting, reflecting an annual vs. perennial lifespan in this species, was observed in 8.7% of the plants, with bolting plants disproportionately occurring in locations with warm, arid climates. Populations with increasing heat and aridity were positively correlated with observed bolting (r = 0.61, p < 0.0001). About one-third (22%) of the total (61%) lifespan variation was attributed to seed source climate and annual heat moisture index, a measure of aridity. Genome size had no significant effect on bolting. Projected climate modeling for mid-century (2041-2070) supports an increasing occurrence of annual lifespan.
CONCLUSIONS: Our analyses support a drought escape, bet-hedging strategy in C. douglasii. Populations exposed to greater aridity exhibited a higher proportion of individuals with an annual lifespan. Drought escape leading to an annual lifespan can affect how seeds are propagated and deployed for climate-informed restoration.

BACKGROUND: White clover (Trifolium repens) is a globally important perennial forage legume. This species also serves as an eco-evolutionary model system for studying within-species chemical defense variation; it features a well-studied polymorphism for cyanogenesis (HCN release following tissue damage), with higher frequencies of cyanogenic plants favored in warmer locations worldwide. Using a newly generated haplotype-resolved genome and two other long-read assemblies, we tested the hypothesis that copy number variants (CNVs) at cyanogenesis genes play a role in the ability of white clover to rapidly adapt to local environments. We also examined questions on subgenome evolution in this recently evolved allotetraploid species and on chromosomal rearrangements in the broader IRLC legume clade.
RESULTS: Integration of PacBio HiFi, Omni-C, Illumina, and linkage map data yielded a completely de novo genome assembly for white clover (created without a priori sequence assignment to subgenomes). We find that white clover has undergone extensive transposon diversification since its origin but otherwise shows highly conserved genome organization and composition with its diploid progenitors. Unlike some other clover species, its chromosomal structure is conserved with other IRLC legumes. We further find extensive evidence of CNVs at the major cyanogenesis loci; these contribute to quantitative variation in the cyanogenic phenotype and to local adaptation across wild North American populations.
CONCLUSIONS: This work provides a case study documenting the role of CNVs in local adaptation in a plant species, and it highlights the value of pan-genome data for identifying contributions of structural variants to adaptation in nature.

Predator responses to warming can occur via phenotypic plasticity, evolutionary adaptation or a combination of both, changing their top-down effects on prey communities. However, we lack evidence of how warming-induced evolutionary changes in predators may influence natural food webs. Here, we ask whether wild fish subject to warming across multiple generations differ in their impacts on prey communities compared with their nearby conspecifics experiencing a natural thermal regime. We carried out a common garden mesocosm experiment with larval perch (Perca fluviatilis), originating from a heated or reference coastal environment, feeding on zooplankton communities under a gradient of experimental temperatures. Overall, in the presence of fish of heated origin, zooplankton abundance was higher and did not change with experimental warming, whereas in the presence of fish of unheated origin, it declined with experimental temperature. Responses in zooplankton taxonomic and size composition suggest that larvae of heated origin consume more large-sized taxa as the temperature increases. Our findings show that differences between fish populations, potentially representing adaptation to their long-term thermal environments, can affect the abundance, biomass, size and species composition of their prey communities. This suggests that rapid microevolution in predators to ongoing climate warming might have indirect cross-generational ecological consequences propagating through food webs.

Adaptive responses to climate change, based on heritable variation in stress tolerance, may be important for plant population persistence. It is unclear which populations will mount the strongest future adaptive responses. It may be fruitful to identify populations that have escaped trade-offs among performance traits, which can hinder adaptation. Barring strong genetic constraints, the extent of trade-offs may depend on spatial relationships among climate variables shaping different traits. Here, we test for climate-driven ecotypic variation and trade-offs among drought and freezing sensitivity, and growth, for Lemmon\'s willow (Salix lemmonii) in a common garden study of 90 genotypes from 38 sites in the Sierra Nevada, USA. Salix lemmonii exhibits ecotypic variation in leaf turgor loss point, a measure of drought sensitivity, from -0.95 to -0.74 MPa along a gradient of spring snowpack. We also find variation in spring freezing sensitivity with minimum May temperature. However, we find no trade-off, as the climatic gradients shaping these traits are spatially uncorrelated in our study region, despite being negatively correlated across the Sierra Nevada. Species may escape adaptive trade-offs in geographic regions where climate variables are spatially decoupled. These regions may represent valuable reservoirs of heritable adaptive phenotypic variation.

Human transport of species across oceans disrupts natural dispersal barriers and facilitates hybridization between previously allopatric species. The recent introduction of the North Pacific sea squirt, Ciona robusta, into the native range of the North Atlantic sea squirt, Ciona intestinalis, is a good example of this outcome. Recent studies have revealed an adaptive introgression in a single chromosomal region from the introduced into the native species. Here, we monitored this adaptive introgression over time, examining both the frequency of adaptive alleles at the core and the hitchhiking footprint in the shoulders of the introgression island by studying a thousand Ciona spp. individuals collected in 22 ports of the contact zone, 14 of which were sampled 20 generations apart. For that purpose, we developed a KASP multiplex genotyping approach, which proved effective in identifying native, nonindigenous and hybrid individuals and in detecting introgressed haplotypes. We found no early generation hybrids in the entire sample, and field observations suggest a decline in the introduced species. At the core region of the introgression sweep, where the frequency of C. robusta alleles is the highest and local adaptation genes must be, we observed stable frequencies of adaptive alien alleles in both space and time. In contrast, we observed erosion of C. robusta ancestry tracts in flanking chromosomal shoulders on the edges of the core, consistent with the second phase of a local sweep and a purge of hitchhiked incompatible mutations. We hypothesize that adaptive introgression may have modified the competition relationships between the native and invasive species in human-altered environments.

Contemporary symbioses in extreme environments can give an insight into mechanisms that stabilize species interactions during environmental change. The intertidal sea anemone, Anthopleura elegantissima, engages in a nutritional symbiosis with microalgae similar to tropical coral, but withstands more intense environmental fluctuations during tidal inundations. In this study, we compare baseline symbiotic traits and their sensitivity to thermal stress within and among anemone aggregations across the intertidal using a laboratory-based tank experiment to better understand how fixed genotypic and plastic environmental effects contribute to the successful maintenance of this symbiosis in extreme habitats. High intertidal anemones had lower baseline symbiont-to-host cell ratios under control conditions, but their symbionts had higher baseline photosynthetic efficiency compared to low intertidal anemone symbionts. Symbiont communities were identical across all samples, suggesting that shifts in symbiont density and photosynthetic performance could be an acclimatory mechanism to maintain symbiosis in different environments. Despite lower baseline symbiont-to-host cell ratios, high intertidal anemones maintained greater symbiont-to-host cell ratios under heat stress compared with low intertidal anemones, suggesting greater thermal tolerance of high intertidal holobionts. However, the thermal tolerance of clonal anemones acclimatized to different zones was not explained by tidal height alone, indicating additional environmental variables contribute to physiological differences. Host genotype significantly influenced anemone weight, but only explained a minor proportion of variation among symbiotic traits and their response to thermal stress, further implicating environmental history as the primary driver of holobiont tolerance. These results indicate that this symbiosis is highly plastic and may be able to acclimatize to climate change over ecological timescales, defying the convention that symbiotic organisms are more susceptible to environmental stress.

When adaptive phenotypic variation or QTLs map within an inverted segment of a chromosome, researchers often despair because the suppression of crossing over will prevent the discovery of selective target genes that established the rearrangement. If an inversion polymorphism is old enough, then the accumulation of gene conversion tracts offers the promise that QTLs or selected loci within inversions can be mapped. The inversion polymorphism of Drosophila pseudoobscura is a model system to show that gene conversion analysis is a useful tool for mapping selected loci within inversions. D. pseudoobscura has over 30 different chromosomal arrangements on the third chromosome (Muller C) in natural populations and their frequencies vary with changes in environmental habitats. Statistical tests of five D. pseudoobscura gene arrangements identified outlier genes within inverted regions that had potentially heritable variation, either fixed amino acid differences or differential expression patterns. We use genome sequences of the inverted third chromosome (Muller C) to infer 98,443 gene conversion tracts for a total coverage of 142 Mb or 7.2 x coverage of the 19.7 Mb chromosome. We estimated gene conversion tract coverage in the 2,668 genes on Muller C and tested whether gene conversion coverage was similar among arrangements for outlier versus non-outlier loci. Outlier genes had lower gene conversion tract coverage among arrangements than the non-outlier genes suggesting that selection removes exchanged DNA in the outlier genes. These data support the hypothesis that the third chromosome in D. pseudoobscura captured locally adapted combinations of alleles prior to inversion mutation events.

Species that repeatedly evolve phenotypic clines across environmental gradients have been highlighted as ideal systems for characterizing the genomic basis of local environmental adaptation. However, few studies have assessed the importance of observed phenotypic clines for local adaptation: conspicuous traits that vary clinally may not necessarily be the most critical in determining local fitness. The present study was designed to fill this gap, using a plant species characterized by repeatedly evolved adaptive phenotypic clines. White clover is naturally polymorphic for its chemical defence cyanogenesis (HCN release with tissue damage); climate-associated cyanogenesis clines have evolved throughout its native and introduced range worldwide. We performed landscape genomic analyses on 415 wild genotypes from 43 locations spanning much of the North American species range to assess the relative importance of cyanogenesis loci vs. other genomic factors in local climatic adaptation. We find clear evidence of local adaptation, with temperature-related climatic variables best describing genome-wide differentiation between sampling locations. The same climatic variables are also strongly correlated with cyanogenesis frequencies and gene copy number variations (CNVs) at cyanogenesis loci. However, landscape genomic analyses indicate no significant contribution of cyanogenesis loci to local adaptation. Instead, several genomic regions containing promising candidate genes for plant response to seasonal cues are identified - some of which are shared with previously identified QTLs for locally adaptive fitness traits in North American white clover. Our findings suggest that local adaptation in white clover is likely determined primarily by genes controlling the timing of growth and flowering in response to local seasonal cues. More generally, this work suggests that caution is warranted when considering the importance of conspicuous phenotypic clines as primary determinants of local adaptation.

Identifying dormancy traits is important for predicting insect population success, particularly in a changing climate that could disrupt evolved traits. The mountain pine beetle (Dendroctonus ponderosae Hopkins) is native to North America, is responsible for millions of acres of tree mortality, and is expanding northward in Canada. Research has identified thermal traits important to epidemic-phase ecology that vary among populations. Genomic research identified 3 mountain pine beetle haplogroups representing Pleistocene glacial refugia. Significant variation in generation timing aligning with the haplogroups has been observed. The adult stage was previously identified as the likely cause of differences among populations, although the mechanism(s) remain unclear. We tested for an adult summer diapause that varies among populations from 2 haplogroups, southern Colorado (CO) (central haplogroup) and southern Idaho (ID) (eastern haplogroup) using respirometry and reproduction experiments. Warm temperatures (25 °C) resulted in reduced respiration rates of central haplogroup mountain pine beetle compared to a cool temperature treatment (15 °C), whereas respiration of the eastern haplogroup did not differ between the treatments. Mated pairs of central haplogroup mountain pine beetle reared/held at 15 °C were more likely to be classified with a higher reproductive success rating compared to pairs reared/held at 25 °C. These results support a facultative summer adult diapause in southern CO central haplogroup mountain pine beetle. Manifestation of this diapause was low/absent among adults from the northerly ID location. This diapause likely serves to maintain univoltinism shown to be important for mountain pine beetle epidemic-phase ecology. The variation occurring among haplogroups highlights the long-term, evolved processes driving local adaptations in mountain pine beetle.