Phenotypic and Genomic Local Adaptation across Latitude and Altitude in Populus trichocarpa

Response of wild aquatic insect communities to thermal variation through comparative landscape transcriptomics

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.

Clinal variation in natural populations of Drosophila melanogaster: An old debate about natural selection and neutral processes

Distinguishing between environmental adaptations and neutral processes poses a challenge in population genetics and evolutionary studies, particularly when phenomena can be explained by both processes. Clines are genotypic or phenotypic characters correlated with environmental variables, because of that correlation, they are used as examples of spatially varying selection. At the same time, many genotypic clines can be explained by demographic history, like isolation by distance or secondary contact zones. Clines have been extensively studied in Drosophila melanogaster, especially in North America and Australia, where they are attributed to both differential selection and various demographic processes. This review explores existing literature supporting this conclusion and suggests new approaches to better understand the influence of these processes on clines. These innovative approaches aim to shed light on the longstanding debate regarding the importance of natural selection versus neutral processes in maintaining variation in natural populations.

Clinal variation in quantitative traits but not in evolutionary potential along elevational and latitudinal gradients in the widespread Anthyllis vulneraria

PREMISE

Strong elevational and latitudinal gradients allow the study of genetic differentiation in response to similar environmental changes. However, it is uncertain whether the environmental changes along the two types of gradients result in similar genetically based changes in quantitative traits. Peripheral arctic and alpine populations are thought to have less evolutionary potential than more central populations do.

METHODS

We studied quantitative traits of the widespread Anthyllis vulneraria in a common garden. Plants originated from 20 populations along a 2000-m elevational gradient from the lowlands to the elevational limit of the species in the Alps, and from 20 populations along a 2400-km latitudinal gradient from the center of the distribution of the species in Central Europe to its northern distributional margin.

RESULTS

Most traits showed similar clinal variations with elevation and latitude of origin, and the magnitude of all measured traits in relation to mean annual temperature was similar. Higher QST values than FST values in several traits indicated diversifying selection, but for others QST was smaller than FST. Genetic diversity of quantitative traits and neutral molecular markers was not correlated. Plasticity in response to favorable conditions declined with elevation and less strongly with latitude of origin, but the evolvability of traits did not.

CONCLUSIONS

The clinal variation suggests adaptive differentiation of quantitative traits along the two gradients. The evolutionary potential of peripheral populations is not necessarily reduced, but lower plasticity may threaten their survival under rapidly changing climatic conditions.

Populus cathayana genome and population resequencing provide insights into its evolution and adaptation

Populus cathayana Rehder, an indigenous poplar species of ecological and economic importance, is widely distributed in a high-elevation range from southwest to northeast China. Further development of this species as a sustainable poplar resource has been hindered by a lack of genome information the at the population level. Here, we produced a chromosome-level genome assembly of P. cathayana, covering 406.55 Mb (scaffold N50 = 20.86 Mb) and consisting of 19 chromosomes, with 35 977 protein-coding genes. Subsequently, we made a genomic variation atlas of 438 wild individuals covering 36 representative geographic areas of P. cathayana, which were divided into four geographic groups. It was inferred that the Northwest China regions served as the genetic diversity centers and a population bottleneck happened during the history of P. cathayana. By genotype-environment association analysis, 947 environment-association loci were significantly associated with temperature, solar radiation, precipitation, and altitude variables. We identified local adaptation genes involved in DNA repair and UV radiation response, among which UVR8, HY5, and CUL4 had key roles in high-altitude adaptation of P. cathayana. Predictions of adaptive potential under future climate conditions showed that P. cathayana populations in areas with drastic climate change were anticipated to have greater maladaptation risk. These results provide comprehensive insights for understanding wild poplar evolution and optimizing adaptive potential in molecular breeding.

Discerning structure versus speciation in phylogeographic analysis of Seepage Salamanders (Desmognathus aeneus) using demography, environment, geography, and phenotype

Numerous mechanisms can drive speciation, including isolation by adaptation, distance, and environment. These forces can promote genetic and phenotypic differentiation of local populations, the formation of phylogeographic lineages, and ultimately, completed speciation. However, conceptually similar mechanisms may also result in stabilizing rather than diversifying selection, leading to lineage integration and the long-term persistence of population structure within genetically cohesive species. Processes that drive the formation and maintenance of geographic genetic diversity while facilitating high rates of migration and limiting phenotypic differentiation may thereby result in population genetic structure that is not accompanied by reproductive isolation. We suggest that this framework can be applied more broadly to address the classic dilemma of "structure" versus "species" when evaluating phylogeographic diversity, unifying population genetics, species delimitation, and the underlying study of speciation. We demonstrate one such instance in the Seepage Salamander (Desmognathus aeneus) from the southeastern United States. Recent studies estimated up to 6.3% mitochondrial divergence and four phylogenomic lineages with broad admixture across geographic hybrid zones, which could potentially represent distinct species supported by our species-delimitation analyses. However, while limited dispersal promotes substantial isolation by distance, microhabitat specificity appears to yield stabilizing selection on a single, uniform, ecologically mediated phenotype. As a result, climatic cycles promote recurrent contact between lineages and repeated instances of high migration through time. Subsequent hybridization is apparently not counteracted by adaptive differentiation limiting introgression, leaving a single unified species with deeply divergent phylogeographic lineages that nonetheless do not appear to represent incipient species.

Divergent selection and climate adaptation fuel genomic differentiation between sister species of Sphagnum (peat moss)

BACKGROUND AND AIMS

New plant species can evolve through the reinforcement of reproductive isolation via local adaptation along habitat gradients. Peat mosses (Sphagnaceae) are an emerging model system for the study of evolutionary genomics and have well-documented niche differentiation among species. Recent molecular studies have demonstrated that the globally distributed species Sphagnum magellanicum is a complex of morphologically cryptic lineages that are phylogenetically and ecologically distinct. Here, we describe the architecture of genomic differentiation between two sister species in this complex known from eastern North America: the northern S. diabolicum and the largely southern S. magniae.

METHODS

We sampled plant populations from across a latitudinal gradient in eastern North America and performed whole genome and restriction-site associated DNA sequencing. These sequencing data were then analyzed computationally.

KEY RESULTS

Using sliding-window population genetic analyses we find that differentiation is concentrated within 'islands' of the genome spanning up to 400 kb that are characterized by elevated genetic divergence, suppressed recombination, reduced nucleotide diversity and increased rates of non-synonymous substitution. Sequence variants that are significantly associated with genetic structure and bioclimatic variables occur within genes that have functional enrichment for biological processes including abiotic stress response, photoperiodism and hormone-mediated signalling. Demographic modelling demonstrates that these two species diverged no more than 225 000 generations ago with secondary contact occurring where their ranges overlap.

CONCLUSIONS

We suggest that this heterogeneity of genomic differentiation is a result of linked selection and reflects the role of local adaptation to contrasting climatic zones in driving speciation. This research provides insight into the process of speciation in a group of ecologically important plants and strengthens our predictive understanding of how plant populations will respond as Earth's climate rapidly changes.

Pleistocene glaciations caused the latitudinal gradient of within-species genetic diversity

Intraspecific genetic diversity is a key aspect of biodiversity. Quaternary climatic change and glaciation influenced intraspecific genetic diversity by promoting range shifts and population size change. However, the extent to which glaciation affected genetic diversity on a global scale is not well established. Here we quantify nucleotide diversity, a common metric of intraspecific genetic diversity, in more than 38,000 plant and animal species using georeferenced DNA sequences from millions of samples. Results demonstrate that tropical species contain significantly more intraspecific genetic diversity than nontropical species. To explore potential evolutionary processes that may have contributed to this pattern, we calculated summary statistics that measure population demographic change and detected significant correlations between these statistics and latitude. We find that nontropical species are more likely to deviate from neutral expectations, indicating that they have historically experienced dramatic fluctuations in population size likely associated with Pleistocene glacial cycles. By analyzing the most comprehensive data set to date, our results imply that Quaternary climate perturbations may be more important as a process driving the latitudinal gradient in species richness than previously appreciated.

Hydraulic trade-offs underlie enhanced performance of polyploid trees under soil water deficit

The relationships between aerial organ morpho-anatomy of woody polyploid plants with their functional hydraulics under water stress remain largely understudied. We evaluated growth-associated traits, aerial organ xylem anatomy, and physiological parameters of diploid, triploid, and tetraploid genotypes of atemoyas (Annona cherimola x Annona squamosa), which belong to the woody perennial genus Annona (Annonaceae), testing their performance under long-term soil water reduction. The contrasting phenotypes of vigorous triploids and dwarf tetraploids consistently showed stomatal size-density trade-off. The vessel elements in aerial organs were ∼1.5 times wider in polyploids compared with diploids, and triploids displayed the lowest vessel density. Plant hydraulic conductance was higher in well-irrigated diploids while their tolerance to drought was lower. The phenotypic disparity of atemoya polyploids associated with contrasting leaf and stem xylem porosity traits that coordinate to regulate water balances between the trees and the belowground and aboveground environments. Polyploid trees displayed better performance under soil water scarcity, presenting as more sustainable agricultural and forestry genotypes to cope with water stress.

Seed type and origin-dependent seedling emergence patterns in Danthonia californica, a species commonly used in grassland restoration

Danthonia californica Bolander (Poaceae)is a native perennial bunchgrass commonly used in the restoration of prairie ecosystems in the western United States. Plants of this species simultaneously produce both chasmogamous (potentially outcrossed) and cleistogamous (obligately self-fertilized) seeds. Restoration practitioners almost exclusively use chasmogamous seeds for outplanting, which are predicted to perform better in novel environments due to their greater genetic diversity. Meanwhile, cleistogamous seeds may exhibit greater local adaptation to the conditions in which the maternal plant exists. We performed a common garden experiment at two sites in the Willamette Valley, Oregon, to assess the influence of seed type and source population (eight populations from a latitudinal gradient) on seedling emergence and found no evidence of local adaptation for either seed type. Cleistogamous seeds outperformed chasmogamous seeds, regardless of whether seeds were sourced directly from the common gardens (local seeds) or other populations (nonlocal seeds). Furthermore, average seed weight had a strong positive effect on seedling emergence, despite the fact that chasmogamous seeds had significantly greater mass than cleistogamous seeds. At one common garden, we observed that seeds of both types sourced from north of our planting site performed significantly better than local or southern-sourced seeds. We also found a significant seed type and distance-dependent interaction, with cleistogamous seedling emergence peaking approximately 125 km from the garden. These results suggest that cleistogamous seeds should be considered for greater use in D. californica restoration.

Genome-wide association study of 23 flowering phenology traits and 4 floral agronomic traits in tree peony (Paeonia section Moutan DC.) reveals five genes known to regulate flowering time

Tree peony is a unique traditional flower in China, with large, fragrant, and colorful flowers. However, a relatively short and concentrated flowering period limits the applications and production of tree peony. A genome-wide association study (GWAS) was conducted to accelerate molecular breeding for the improvement of flowering phenology traits and ornamental phenotypes in tree peony. A diverse panel of 451 tree peony accessions was phenotyped for 23 flowering phenology traits and 4 floral agronomic traits over 3 years. Genotyping by sequencing (GBS) was used to obtain a large number of genome-wide single-nucleotide polymorphisms (SNPs) (107 050) for the panel genotypes, and 1047 candidate genes were identified by association mapping. Eighty-two related genes were observed during at least 2 years for flowering, and seven SNPs repeatedly identified for multiple flowering phenology traits over multiple years were highly significantly associated with five genes known to regulate flowering time. We validated the temporal expression profiles of these candidate genes and highlighted their possible roles in the regulation of flower bud differentiation and flowering time in tree peony. This study shows that GWAS based on GBS can be used to identify the genetic determinants of complex traits in tree peony. The results expand our understanding of flowering time control in perennial woody plants. Identification of markers closely related to these flowering phenology traits can be used in tree peony breeding programs for important agronomic traits.

Modeling the effect of adaptation to future climate change on spring phenological trend of European beech (Fagus sylvatica L.)

Temperate trees could cope with climate change through phenotypic plasticity of phenological key events or adaptation in situ via selection on genetic variation. However, the relative contribution of local adaptation and phenotypic plasticity to phenological change is unclear for many ecologically important tree species. Here, we analyzed the leaf-out data of European beech (Fagus sylvatica L.) from 50 provenances planted in 7 trial sites. We first constructed a function between chilling accumulation (CA) and photoperiod-associated heat requirement (PHR) of leaf-out date for each provenance and quantified the relationship between parameters of the CA-PHR function and climatic variables at provenance origins by using the random forest model. Furthermore, we used the provenance-specific CA-PHR function to simulate future leaf-out dates under two climate change scenarios (RCP 4.5 and 8.5) and two assumptions (no adaptation and adaptation). The results showed that both CA, provenance, and their interactions affected the PHR of leaf-out. The provenances from southeastern Europe exhibited a stronger response of PHR to CA and thus flushed earlier than northwestern provenances. The parameters of the CA-PHR function were connected with climatic variables (e.g., mean diurnal temperature range, temperature seasonality) at the originating sites of each provenance. If only considering the phenotypic plasticity, the leaf-out date of European beech in 2070-2099 will advance by 6.8 and 9.0 days on average relative to 1951-2020 under RCP 4.5 and RCP 8.5, respectively. However, if F. sylvatica adapts to future climate change by adopting the current strategy, the advance of the leaf-out date will weaken by 1.4 and 3.4 days under RCP 4.5 and RCP 8.5, respectively. Our results suggest that the European beech could slow down its spring phenological advances and reduce its spring frost risk if it adopts the current strategy to adapt to future climate change.

Assisted migration is plausible for a boreal tree species under climate change: A quantitative and population genetics study of trembling aspen (Populus tremuloides Michx.) in western Canada

A novel method was tested for improving tree breeding strategies that integrate quantitative and population genetics based on range‐wide reciprocal transplant experiments. Five reciprocal common garden tests of Populus tremuloides were investigated including 6450 trees across western Canada focusing on adaptation traits and growth. Both genetic parameters and home‐site transplant models were evaluated. We found a genetic trade‐off between growth and early spring leaf flush and late fall senescence. Coefficients of phenotypic variation (CVp) of cell lysis (CL), a measure of freezing injury, shrank from 0.28 to 0.10 during acclimation in the fall, and the CVp slope versus the freezing temperature was significantly different from zero (R² = 0.33, p = .02). There was more between‐population genetic variation in fall phenology than in spring leaf phenology. We suggest that P. tremuloides demonstrated a discrepancy between the ecological optimum and the physiological optimum minimum winter temperature. The sub‐optimal growing condition of P. tremuloides is potentially caused by the warmer ecological optimum than the physiological optimum. Assisted migration and breeding of fast growers to reforest cooler plantation sites can improve productivity. Transferring the study populations to less than 4°C of extreme minimum temperature appears safe for reforestation aligning with the historical recolonization direction of the species. This is equivalent to a 5–10° latitudinal northward movement. Fall frost hardiness is an effective criterion for family selection in the range tested in this study.

Genome-Wide Identification, Evolution, and Expression Analysis of GASA Gene Family in Prunus mume

The Gibberellic Acid Stimulated Arabidopsis/Gibberellin Stimulated Transcript (GASA/GAST) gene family is a group of plant-specific genes encoding cysteine-rich peptides essential to plant growth, development, and stress responses. Although GASA family genes have been identified in various plant species, their functional roles in Prunus mume are still unknown. In this study, a total of 16 PmGASA genes were identified via a genome-wide scan in Prunus mume and were grouped into three major gene clades based on the phylogenetic tree. All PmGASA proteins possessed the conserved GASA domain, consisting of 12-cysteine residues, but varied slightly in protein physiochemical properties and motif composition. With evolutionary analysis, we observed that duplications and purifying selection are major forces driving PmGASA family gene evolution. By analyzing PmGASA promoters, we detected a number of hormonal-response related cis-elements and constructed a putative transcriptional regulatory network for PmGASAs. To further understand the functional role of PmGASA genes, we analyzed the expression patterns of PmGASAs across different organs and during various biological processes. The expression analysis revealed the functional implication of PmGASA gene members in gibberellic acid-, abscisic acid-, and auxin-signaling, and during the progression of floral bud break in P. mume. To summarize, these findings provide a comprehensive understanding of GASA family genes in P. mume and offer a theoretical basis for future research on the functional characterization of GASA genes in other woody perennials.

Integrated transcriptome and small RNA sequencing in revealing miRNA-mediated regulatory network of floral bud break in Prunus mume

MicroRNAs is one class of small non-coding RNAs that play important roles in plant growth and development. Though miRNAs and their target genes have been widely studied in many plant species, their functional roles in floral bud break and dormancy release in woody perennials is still unclear. In this study, we applied transcriptome and small RNA sequencing together to systematically explore the transcriptional and post-transcriptional regulation of floral bud break in P. mume. Through expression profiling, we identified a few candidate genes and miRNAs during different developmental stage transitions. In total, we characterized 1,553 DEGs associated with endodormancy release and 2,084 DEGs associated with bud flush. Additionally, we identified 48 known miRNAs and 53 novel miRNAs targeting genes enriched in biological processes such as floral organ morphogenesis and hormone signaling transudation. We further validated the regulatory relationship between differentially expressed miRNAs and their target genes combining computational prediction, degradome sequencing, and expression pattern analysis. Finally, we integrated weighted gene co-expression analysis and constructed miRNA-mRNA regulatory networks mediating floral bud flushing competency. In general, our study revealed the miRNA-mediated networks in modulating floral bud break in P. mume. The findings will contribute to the comprehensive understanding of miRNA-mediated regulatory mechanism governing floral bud break and dormancy cycling in wood perennials.

Transcriptomic responses and physiological changes to cold stress among natural populations provide insights into local adaptation of weeping forsythia

Genetic mechanisms of species local adaptation are an emerging topic of great interest in evolutionary biology and molecular ecology. In this study, we compared the changes of physiological and phenotypic indexes and gene expression of four weeping forsythia populations under cold stress through a common garden experiment. Physiological and phenotypic results showed that there were differences in cold tolerance among populations. cold tolerance of high the latitude population (HBWZ) was the strongest, followed by the middle latitude population (SXWL), while the low latitude populations (SXHM) and (SXLJ) expressed the weakest cold tolerance. We identified significant differences in gene expression of cold tolerance related pathways and ontologies, including genes of oxylipin and isoquinoline alkaloid biosynthetic process, galactose, tyrosine and unsaturated fatty acids metabolism, among these populations under the same experimental temperature treatments. Even under the same degree of stress, there were notable differences in gene expression among natural populations. In this study, we present a working model of weeping forsythia populations which evolved in the context of different intensities of cold stress. Our study provides new insights for comprehending the genetic mechanisms of local adaptation for non-model species.

Local Adaptation in Marine Foundation Species at Microgeographic Scales

AbstractNearshore foundation species in coastal and estuarine systems (e.g., salt marsh grasses, mangroves, seagrasses, corals) drive the ecological functions of ecosystems and entire biomes by creating physical structure that alters local abiotic conditions and influences species interactions and composition. The resilience of foundation species and the ecosystem functions they provide depends on their phenotypic and genetic responses to spatial and temporal shifts in environmental conditions. In this review, we explore what is known about the causes and consequences of adaptive genetic differentiation in marine foundation species over spatial scales shorter than dispersal capabilities (i.e., microgeographic scales). We describe the strength of coupling field and laboratory experiments with population genetic techniques to illuminate patterns of local adaptation, and we illustrate this approach by using several foundation species. Among the major themes that emerge from our review include (1) adaptive differentiation of marine foundation species repeatedly evolves along vertical (i.e., elevation or depth) gradients, and (2) mating system and phenology may facilitate this differentiation. Microgeographic adaptation is an understudied mechanism potentially underpinning the resilience of many sessile marine species, and this evolutionary mechanism likely has particularly important consequences for the ecosystem functions provided by foundation species.

Genomic signatures of natural selection at phenology-related genes in a widely distributed tree species Fagus sylvatica L

BACKGROUND

Diversity among phenology-related genes is predicted to be a contributing factor in local adaptations seen in widely distributed plant species that grow in climatically variable geographic areas, such as forest trees. European beech (Fagus sylvatica L.) is widespread, and is one of the most important broadleaved tree species in Europe; however, its potential for adaptation to climate change is a matter of uncertainty, and little is known about the molecular basis of climate change-relevant traits like bud burst.

RESULTS

We explored single nucleotide polymorphisms (SNP) at candidate genes related to bud burst in beech individuals sampled across 47 populations from Europe. SNP diversity was monitored for 380 candidate genes using a sequence capture approach, providing 2909 unlinked SNP loci. We used two complementary analytical methods to find loci significantly associated with geographic variables, climatic variables (expressed as principal components), or phenotypic variables (spring and autumn phenology, height, survival). Redundancy analysis (RDA) was used to detect candidate markers across two spatial scales (entire study area and within subregions). We revealed 201 candidate SNPs at the broadest scale, 53.2% of which were associated with phenotypic variables. Additive polygenic scores, which provide a measure of the cumulative signal across significant candidate SNPs, were correlated with a climate variable (first principal component, PC1) related to temperature and precipitation availability, and spring phenology. However, different genotype-environment associations were identified within Southeastern Europe as compared to the entire geographic range of European beech.

CONCLUSIONS

Environmental conditions play important roles as drivers of genetic diversity of phenology-related genes that could influence local adaptation in European beech. Selection in beech favors genotypes with earlier bud burst under warmer and wetter habitats within its range; however, selection pressures may differ across spatial scales.

The genetic basis of adaptation in phenology in an introduced population of Black Cottonwood (Populus trichocarpa, Torr. & Gray)

BACKGROUND

Entering and exiting winter dormancy present important trade-offs between growth and survival at northern latitudes. Many forest trees display local adaptation across latitude in traits associated with these phenology transitions. Transfers of a species outside its native range introduce the species to novel combinations of environmental conditions potentially requiring different combinations of alleles to optimize growth and survival. In this study, we performed genome wide association analyses and a selection scan in a P. trichocarpa mapping population derived from crossings between clones collected across the native range and introduced into Sweden. GWAS analyses were performed using phenotypic data collected across two field seasons and in a controlled phytotron experiment.

RESULTS

We uncovered 584 putative candidate genes associated with spring and autumn phenology traits as well as with growth. Many regions harboring variation significantly associated with the initiation of leaf shed and leaf autumn coloring appeared to have been evolving under positive selection in the native environments of P. trichocarpa. A comparison between the candidate genes identified with results from earlier GWAS analyses performed in the native environment found a smaller overlap for spring phenology traits than for autumn phenology traits, aligning well with earlier observations that spring phenology transitions have a more complex genetic basis than autumn phenology transitions.

CONCLUSIONS

In a small and structured introduced population of P. trichocarpa, we find complex genetic architectures underlying all phenology and growth traits, and identify multiple putative candidate genes despite the limitations of the study population.

The genomics of rapid climatic adaptation and parallel evolution in North American house mice

Parallel changes in genotype and phenotype in response to similar selection pressures in different populations provide compelling evidence of adaptation. House mice (Mus musculus domesticus) have recently colonized North America and are found in a wide range of environments. Here we measure phenotypic and genotypic differentiation among house mice from five populations sampled across 21° of latitude in western North America, and we compare our results to a parallel latitudinal cline in eastern North America. First, we show that mice are genetically differentiated between transects, indicating that they have independently colonized similar environments in eastern and western North America. Next, we find genetically-based differences in body weight and nest building behavior between mice from the ends of the western transect which mirror differences seen in the eastern transect, demonstrating parallel phenotypic change. We then conduct genome-wide scans for selection and a genome-wide association study to identify targets of selection and candidate genes for body weight. We find some genomic signatures that are unique to each transect, indicating population-specific responses to selection. However, there is significant overlap between genes under selection in eastern and western house mouse transects, providing evidence of parallel genetic evolution in response to similar selection pressures across North America.

Natural population re-sequencing detects the genetic basis of local adaptation to low temperature in a woody plant

Total of 14 SNPs associated with overwintering-related traits and 75 selective regions were detected. Important candidate genes were identified and a possible network of cold-stress responses in woody plants was proposed. Local adaptation to low temperature is essential for woody plants to against changeable climate and safely survive the winter. To uncover the specific molecular mechanism of low temperature adaptation in woody plants, we sequenced 134 core individuals selected from 494 paper mulberry (Broussonetia papyrifera), which naturally distributed in different climate zones and latitudes. The population structure analysis, PCA analysis and neighbor-joining tree analysis indicated that the individuals were classified into three clusters, which showed forceful geographic distribution patterns because of the adaptation to local climate. Using two overwintering phenotypic data collected at high latitudes of 40°N and one bioclimatic variable, genome-phenotype and genome-environment associations, and genome-wide scans were performed. We detected 75 selective regions which possibly undergone temperature selection and identified 14 trait-associated SNPs that corresponded to 16 candidate genes (including LRR-RLK, PP2A, BCS1, etc.). Meanwhile, low temperature adaptation was also supported by other three trait-associated SNPs which exhibiting significant differences in overwintering traits between alleles within three geographic groups. To sum up, a possible network of cold signal perception and responses in woody plants were proposed, including important genes that have been confirmed in previous studies while others could be key potential candidates of woody plants. Overall, our results highlighted the specific and complex molecular mechanism of low temperature adaptation and overwintering of woody plants.

Genomic analyses provide insights into peach local adaptation and responses to climate change

The environment has constantly shaped plant genomes, but the genetic bases underlying how plants adapt to environmental influences remain largely unknown. We constructed a high-density genomic variation map of 263 geographically representative peach landraces and wild relatives. A combination of whole-genome selection scans and genome-wide environmental association studies (GWEAS) was performed to reveal the genomic bases of peach adaptation to diverse climates. A total of 2092 selective sweeps that underlie local adaptation to both mild and extreme climates were identified, including 339 sweeps conferring genomic pattern of adaptation to high altitudes. Using genome-wide environmental association studies (GWEAS), a total of 2755 genomic loci strongly associated with 51 specific environmental variables were detected. The molecular mechanism underlying adaptive evolution of high drought, strong UVB, cold hardiness, sugar content, flesh color, and bloom date were revealed. Finally, based on 30 yr of observation, a candidate gene associated with bloom date advance, representing peach responses to global warming, was identified. Collectively, our study provides insights into molecular bases of how environments have shaped peach genomes by natural selection and adds candidate genes for future studies on evolutionary genetics, adaptation to climate changes, and breeding.

Locally adapted gut microbiomes mediate host stress tolerance

While evidence for the role of the microbiome in shaping host stress tolerance is becoming well-established, to what extent this depends on the interaction between the host and its local microbiome is less clear. Therefore, we investigated whether locally adapted gut microbiomes affect host stress tolerance. In the water flea Daphnia magna, we studied if the host performs better when receiving a microbiome from their source region than from another region when facing a stressful condition, more in particular exposure to the toxic cyanobacteria Microcystis aeruginosa. Therefore, a reciprocal transplant experiment was performed in which recipient, germ-free D. magna, isolated from different ponds, received a donor microbiome from sympatric or allopatric D. magna that were pre-exposed to toxic cyanobacteria or not. We tested for effects on host life history traits and gut microbiome composition. Our data indicate that Daphnia interact with particular microbial strains mediating local adaptation in host stress tolerance. Most recipient D. magna individuals performed better when inoculated with sympatric than with allopatric microbiomes. This effect was most pronounced when the donors were pre-exposed to the toxic cyanobacteria, but this effect was also pond and genotype dependent. We discuss how this host fitness benefit is associated with microbiome diversity patterns.

Clair JB, Crepeau MW, Salzberg SL, Langley CH, Allen B, Neale DB. Dissecting the Polygenic Basis of Cold Adaptation Using Genome-Wide Association of Traits and Environmental Data in Douglas-fir

Understanding the genomic and environmental basis of cold adaptation is key to understand how plants survive and adapt to different environmental conditions across their natural range. Univariate and multivariate genome-wide association (GWAS) and genotype-environment association (GEA) analyses were used to test associations among genome-wide SNPs obtained from whole-genome resequencing, measures of growth, phenology, emergence, cold hardiness, and range-wide environmental variation in coastal Douglas-fir (Pseudotsuga menziesii). Results suggest a complex genomic architecture of cold adaptation, in which traits are either highly polygenic or controlled by both large and small effect genes. Newly discovered associations for cold adaptation in Douglas-fir included 130 genes involved in many important biological functions such as primary and secondary metabolism, growth and reproductive development, transcription regulation, stress and signaling, and DNA processes. These genes were related to growth, phenology and cold hardiness and strongly depend on variation in environmental variables such degree days below 0c, precipitation, elevation and distance from the coast. This study is a step forward in our understanding of the complex interconnection between environment and genomics and their role in cold-associated trait variation in boreal tree species, providing a baseline for the species' predictions under climate change.

Integrating Genome-Wide Association Analysis With Transcriptome Sequencing to Identify Candidate Genes Related to Blooming Time in Prunus mume

Prunus mume is one of the most important woody perennials for edible and ornamental use. Despite a substantial variation in the flowering phenology among the P. mume germplasm resources, the genetic control for flowering time remains to be elucidated. In this study, we examined five blooming time-related traits of 235 P. mume landraces for 2 years. Based on the phenotypic data, we performed genome-wide association studies, which included a combination of marker- and gene-based association tests, and identified 1,445 candidate genes that are consistently linked with flowering time across multiple years. Furthermore, we assessed the global transcriptome change of floral buds from the two P. mume cultivars exhibiting contrasting bloom dates and detected 617 associated genes that were differentially expressed during the flowering process. By integrating a co-expression network analysis, we screened out 191 gene candidates of conserved transcriptional pattern during blooming across cultivars. Finally, we validated the temporal expression profiles of these candidates and highlighted their putative roles in regulating floral bud break and blooming time in P. mume. Our findings are important to expand the understanding of flowering time control in woody perennials and will boost the molecular breeding of novel varieties in P. mume.

Ecological Factors Generally Not Altitude Related Played Main Roles in Driving Potential Adaptive Evolution at Elevational Range Margin Populations of Taiwan Incense Cedar (Calocedrus formosana)

Population diversification can be shaped by a combination of environmental factors as well as geographic isolation interacting with gene flow. We surveyed genetic variation of 243 samples from 12 populations of Calocedrus formosana using amplified fragment length polymorphism (AFLP) and scored a total of 437 AFLP fragments using 11 selective amplification primer pairs. The AFLP variation was used to assess the role of gene flow on the pattern of genetic diversity and to test environments in driving population adaptive evolution. This study found the relatively lower level of genetic diversity and the higher level of population differentiation in C. formosana compared with those estimated in previous studies of conifers including Cunninghamia konishii, Keteleeria davidiana var. formosana, and Taiwania cryptomerioides occurring in Taiwan. BAYESCAN detected 26 F_ST outlier loci that were found to be associated strongly with various environmental variables using multiple univariate logistic regression, latent factor mixed model, and Bayesian logistic regression. We found several environmentally dependent adaptive loci with high frequencies in low- or high-elevation populations, suggesting their involvement in local adaptation. Ecological factors, including relative humidity and sunshine hours, that are generally not altitude related could have been the most important selective drivers for population divergent evolution in C. formosana. The present study provides fundamental information in relation to adaptive evolution and can be useful for assisted migration program of C. formosana in the future conservation of this species.

Latitudinal clines in bud flush phenology reflect genetic variation in chilling requirements in balsam poplar, Populus balsamifera

PREMISE

Boreal and northern temperate forest trees possess finely tuned mechanisms of dormancy, which match bud phenology with local seasonality. After winter dormancy, the accumulation of chilling degree days (CDD) required for rest completion before the accumulation of growing degree days (GDD) during quiescence is an important step in the transition to spring bud flush. While bud flush timing is known to be genetically variable within species, few studies have investigated variation among genotypes from different climates in response to variable chilling duration.

METHODS

We performed a controlled environment study using dormant cuttings from 10 genotypes of Populus balsamifera, representing a broad latitudinal gradient (43-58°N). We exposed cuttings to varying amounts of chilling (0-10 weeks) and monitored subsequent GDD to bud flush at a constant forcing temperature.

RESULTS

Chilling duration strongly accelerated bud flush timing, with increasing CDD resulting in fewer GDD to flush. Genotypic variation for bud flush was significant and stratified by latitude, with southern genotypes requiring more GDD to flush than northern genotypes. The latitudinal cline was pronounced under minimal chilling, whereas genotypic variation in GDD to bud flush converged as CDD increased.

CONCLUSIONS

We demonstrate that increased chilling lessens GDD to bud flush in a genotype-specific manner. Our results emphasize that latitudinal clines in bud flush reflect a critical genotype-by-environment interaction, whereby differences in bud flush between southern vs. northern genotypes depend on chilling. Our results suggest selection has shaped chilling requirements and depth of rest as an adaptive strategy to avoid precocious flush in climates with midwinter warming.