Genome scanning for detecting adaptive genes along environmental gradients in the Japanese conifer, Cryptomeria japonica

Climate-Associated Genetic Variation and Projected Genetic Offsets for Cryptomeria japonica D. Don Under Future Climate Scenarios

Revealing the spatial distribution of adaptive genetic variation is both a challenging and crucial task in evolutionary ecology, essential for understanding local adaptation within species, and in management, for predicting species responses to future climate change. This understanding is particularly important for long-lived tree species, which may not be able to migrate quickly enough to adapt to rapid climate changes and may need to rely on their standing genetic variation. In this study, we focused on Cryptomeria japonica, a major component of Japan's temperate forests and an important forestry species adapted to the humid environment of monsoon Asia. We extracted climate-associated genetic variation from the entire genome and evaluated its distribution and vulnerability under future climate scenarios using spatial modeling techniques. We analyzed 31,676 high-quality SNPs from 249 individuals across 22 natural populations of C. japonica, covering its entire distribution range. We identified 239 candidate climate-associated SNPs and found winter temperature, summer precipitation, and winter precipitation as the most significant factors explaining the genetic variation in these SNPs. The climate-associated genetic variation deviated from non-associated (neutral) genetic variation in the opposite (the Sea of Japan and Pacific Ocean) sides of Japanese archipelago, suggesting natural selection of different climate conditions in these regions. Difference in estimated allele frequency at the climate-associated loci (genetic offset) between the present and future (2090 in the SSP5-8.5 scenario) climate conditions was predicted to be larger in three areas (not only southwestern Japan but also coastal area on the Sea of Japan side and inland area on the Pacific Ocean side in northeastern Japan). This prediction implies the discrepancy between standing genetic variation at the present and that adaptive to the future climate in these areas, which underscores the necessity for proactive management to adjust the adaptive genetic variation.

Essential Oil Composition and Anti-Cholinesterase Properties of Cryptomeria japonica Foliage Harvested in São Miguel Island (Azores) in Two Different Seasons

The Azorean Cryptomeria japonica forest operations and wood industry generate considerable foliage biomass residues that are used for local essential oil (EO) production. However, research on seasonal variation of C. japonica EO remains scarce. In this study, the EOs from fresh Azorean C. japonica foliage (Az-CJF) collected in autumn (Aut) and spring (Spr) were obtained via hydrodistillation and investigated for their physical properties, yield, chemical composition, and bioactivities. Both EOs presented a strong odor, a yellowish color, a density around 0.9 g·mL-1, and similar yields (approximately 1% v/w, dry matter). Nevertheless, the GC-MS analyses showed a decrease in monoterpene hydrocarbons (MH) and an increase in oxygenated sesquiterpenes (OS) contents in Spr-EO compared with Aut-EO (16% vs. 35% for MH and 45% vs. 31% for OS, respectively). In addition, the predominant components were kaur-16-ene (23%) for Spr-EO and phyllocladene (19%) for Aut-EO, revealing that both EOs were rich in diterpene hydrocarbons (29% vs. 26%). Concerning its toxicity against brine shrimp, a low mortality (0-38%) was observed at a concentration range of 100-180 μg·mL-1. Regarding the anti-cholinesterase properties, both EOs were inactive against acetylcholinesterase but showed anti-butyrylcholinesterase activity superior to (-)-α-pinene, a major compound of Az-CJF EO (IC50 values: 84, 148, and 648 μg·mL-1 for Spr-EO, Aut-EO, and α-pinene, respectively). Overall, the results indicate the potential benefit of both seasonal EOs in Alzheimer's disease treatment. In conclusion, this study demonstrated that season strongly influences the Az-CJF EO quantitative composition and thus its bioactivity, aiding in the selection of the most high-quality raw materials for use in Azorean C. japonica EO aromatherapy industry.

Population structure of wild soybean (Glycine soja) based on SLAF-seq have implications for its conservation

Background

Glycine soja Sieb. & Zucc. is the wild ancestor from which the important crop plant soybean was bred. G. soja provides important germplasm resources for the breeding and improvement of cultivated soybean crops, however the species is threatened by habitat loss and fragmentation, and is experiencing population declines across its natural range. Understanding the patterns of genetic diversity in G. soja populations can help to inform conservation practices.

Methods

In this study, we analyzed the genetic diversity and differentiation of G. soja at different sites and investigated the gene flow within the species. We obtained 147 G. soja accessions collected from 16 locations across the natural range of the species from China, Korea and Japan. Samples were analyzed using SLAF-seq (Specific-Locus Amplified Fragment Sequencing).

Results

We obtained a total of 56,489 highly consistent SNPs. Our results suggested that G. soja harbors relatively high diversity and that populations of this species are highly differentiated. The populations harboring high genetic diversity, especially KR, should be considered first when devising conservation plans for the protection of G. soja, and in situ protection should be adopted in KR. G. soja populations from the Yangtze River, the Korean peninsula and northeastern China have a close relationship, although these areas are geographically disconnected. Other populations from north China clustered together. Analysis of gene flow suggested that historical migrations of G. soja may have occurred from the south northwards across the East-Asia land-bridge, but not across north China. All G. soja populations could be divided into one of two lineages, and these two lineages should be treated separately when formulating protection policies.

Epigenetic architecture of Pseudotaxus chienii: Revealing the synergistic effects of climate and soil variables

Whether conifers can withstand environmental changes especially temperature fluctuations has been controversial. Epigenetic analysis may provide new perspectives for solving the issue. Pseudotaxus chienii is an endangered gymnosperm species endemic to China. In this study, we have examined the genetic and epigenetic variations in its natural populations aiming to disentangle the synergistic effects of climate and soil on its population (epi)genetic differentiation by using amplified fragment length polymorphism (AFLP) and methylation-sensitive AFLP (MSAP) techniques. We identified 23 AFLP and 26, 7, and 5 MSAP outliers in P. chienii. Twenty-one of the putative adaptive AFLP loci were found associated with climate and/or soil variables including precipitation, temperature, K, Fe, Zn, and Cu, whereas 21, 7, and 4 MSAP outliers were significantly related to precipitation of wettest month (Bio13), precipitation driest of month (Bio14), percent tree cover (PTC), and soil Fe, Mn, and Cu compositions. Total precipitation and precipitation in the driest seasons were the most influential factors for genetic and epigenetic variation, respectively. In addition, a high full-methylation level and a strong correlation between genetic and epigenetic variation were detected in P. chienii. Climate is found of greater importance than soil in shaping adaptive (epi)genetic differentiation, and the synergistic effects of climate and climate-soil variables were also observed. The identified climate and soil variables should be considered when applying ex situ conservation.

Key triggers of adaptive genetic variability of sessile oak [Q. petraea (Matt.) Liebl.] from the Balkan refugia: outlier detection and association of SNP loci from ddRAD-seq data

Knowledge on the genetic composition of Quercus petraea in south-eastern Europe is limited despite the species' significant role in the re-colonisation of Europe during the Holocene, and the diverse climate and physical geography of the region. Therefore, it is imperative to conduct research on adaptation in sessile oak to better understand its ecological significance in the region. While large sets of SNPs have been developed for the species, there is a continued need for smaller sets of SNPs that are highly informative about the possible adaptation to this varied landscape. By using double digest restriction site associated DNA sequencing data from our previous study, we mapped RAD-seq loci to the Quercus robur reference genome and identified a set of SNPs putatively related to drought stress-response. A total of 179 individuals from eighteen natural populations at sites covering heterogeneous climatic conditions in the southeastern natural distribution range of Q. petraea were genotyped. The detected highly polymorphic variant sites revealed three genetic clusters with a generally low level of genetic differentiation and balanced diversity among them but showed a north-southeast gradient. Selection tests showed nine outlier SNPs positioned in different functional regions. Genotype-environment association analysis of these markers yielded a total of 53 significant associations, explaining 2.4-16.6% of the total genetic variation. Our work exemplifies that adaptation to drought may be under natural selection in the examined Q. petraea populations.

Changes in annual transcriptome dynamics of a clone of Japanese cedar (Cryptomeria japonica D. Don) planted under different climate conditions

Environmental responses are critical for plant growth and survival under different climate conditions. To elucidate the underlying biological mechanisms of environmental responses in Japanese cedar (Cryptomeria japonica D. Don), the annual transcriptome dynamics of common clonal trees (Godai1) planted at three different climate sites (Yamagata, Ibaraki, and Kumamoto Prefectures) were analyzed using microarrays. Both principal component analysis (PCA) and hierarchical clustering of the microarray data indicated the transition to dormant transcriptome status occurred earlier and the transition to active growth status later in the colder region. Interestingly, PCA also indicated that the transcriptomes of trees grown under three different conditions were similar during the growth period (June to September), whereas the transcriptomes differed between sites during the dormant period (January to March). In between-site comparisons, analyses of the annual expression profiles of genes for sites 'Yamagata vs. Kumamoto', 'Yamagata vs. Ibaraki', and 'Ibaraki vs. Kumamoto' identified 1,473, 1,137, and 925 targets exhibiting significantly different expression patterns, respectively. The total of 2,505 targets that exhibited significantly different expression patterns in all three comparisons may play important roles in enabling cuttings to adapt to local environmental conditions. Partial least-squares regression analysis and Pearson correlation coefficient analysis revealed that air temperature and day length were the dominant factors controlling the expression levels of these targets. GO and Pfam enrichment analyses indicated that these targets include genes that may contribute to environmental adaptation, such as genes related to stress and abiotic stimulus responses. This study provided fundamental information regarding transcripts that may play an important role in adaptation to environmental conditions at different planting sites.

Reduced representation sequencing to understand the evolutionary history of Torrey pine (Pinus torreyana Parry) with implications for rare species conservation

Understanding the contribution of neutral and adaptive evolutionary processes to population differentiation is often necessary for better informed management and conservation of rare species. In this study, we focused on Pinus torreyana Parry (Torrey pine), one of the world's rarest pines, endemic to one island and one mainland population in California. Small population size, low genetic diversity, and susceptibility to abiotic and biotic stresses suggest Torrey pine may benefit from inter-population genetic rescue to preserve the species' evolutionary potential. We leveraged reduced representation sequencing to tease apart the respective contributions of stochastic and deterministic evolutionary processes to population differentiation. We applied these data to model spatial and temporal demographic changes in effective population sizes and genetic connectivity, to identify loci possibly under selection, and evaluate genetic rescue as a potential conservation strategy. Overall, we observed exceedingly low standing variation within both Torrey pine populations, reflecting consistently low effective population sizes across time, and limited genetic differentiation, suggesting maintenance of gene flow between populations following divergence. However, genome scans identified more than 2000 candidate SNPs potentially under divergent selection. Combined with previous observations indicating population phenotypic differentiation, this indicates natural selection has likely contributed to the evolution of population genetic differences. Thus, while reduced genetic diversity, small effective population size, and genetic connectivity between populations suggest genetic rescue could mitigate the adverse effects of rarity, evidence for adaptive differentiation suggests genetic mixing could disrupt adaptation. Further work evaluating the fitness consequences of inter-population admixture is necessary to empirically evaluate the trade-offs associated with genetic rescue in Torrey pine.

Genome-environment association methods comparison supports omnigenic adaptation to ecological niche in malaria vector mosquitoes

The concept of a fundamental ecological niche is central to questions of geographic distribution, population demography, species conservation, and evolutionary potential. However, robust inference of genomic regions associated with evolutionary adaptation to particular environmental conditions remains difficult due to the myriad of potential confounding processes that can generate heterogeneous patterns of variation across the genome. Here, we interrogate the potential role of genome environment association (GEA) testing as an initial step in building an understanding of the genetic basis of ecological niche. We leverage publicly available genomic data from the Anopheles gambiae 1000 Genomes (Ag1000g) Consortium to test the ability of multiple analytically unique GEA methods to handle confounding patterns of genetic variation, control false positive rates, and discern associations with broadly relevant climate variables from random allele frequency patterns throughout the genome. We found evidence supporting the ability of commonly implemented GEA methods to account for confounding patterns of spatial and genetic variation, and control false positive rates. However, we fail to find evidence supporting the ability of GEA tests to reject signals of adaptation to randomly simulated environmental variables, indicating that discerning between true signals of genome environment adaptation and genome environment correlations resulting from alternative evolutionary processes, remains challenging. Because signals of environmental adaptation are so diffuse and confounded throughout the genome, we argue that genomic adaptation to ecological niche is likely best understood under an omnigenic model wherein highly interconnected, genome-wide gene regulatory networks shape genomic adaptation to key environmental conditions.

Imprints of selection in peripheral and ecologically marginal central-eastern European Scots pine populations

Knowledge of the molecular mechanisms underlying the stress response in plants is essential to understand evolutionary processes that result in long-term persistence of populations. Populations inhabiting marginal ecological conditions at the distribution range periphery may have preserved imprints of natural selection that have shaped functional genetic variation of the species. Our aim was to evaluate the extent of selection processes in the extremely fragmented, peripheral and isolated populations of Scots pine in central-eastern Europe. Autochthonous populations of the Carpathian Mts. and the Pannonian Basin were sampled and drought stress-related candidate genes were re-sequenced. Neutrality tests and outlier detection approaches were applied to infer the effect and direction of selection. Populations retained high genetic diversity by preserving a high number of alleles and haplotypes, many of them being population specific. Neutrality tests and outlier detection highlighted nucleotide positions that are under divergent selection and may be involved in local adaptation. The detected genetic pattern confirms that natural selection has played an important role in shaping modern-day genetic variation in marginal Scots pine populations, allowing for the long-term persistence of populations. Selection detected at functional regions possibly acts to maintain diversity and counteract the effect of genetic erosion.

Genomic insights on the contribution of balancing selection and local adaptation to the long-term survival of a widespread living fossil tree, Cercidiphyllum japonicum

'Living fossils' are testimonies of long-term sustained ecological success, but how demographic history and natural selection contributed to their survival, resilience, and persistence in the face of Quaternary climate fluctuations remains unclear. To better understand the interplay between demographic history and selection in shaping genomic diversity and evolution of such organisms, we assembled the whole genome of Cercidiphyllum japonicum, a widespread East Asian Tertiary relict tree, and resequenced 99 individuals of C. japonicum and its sister species, Cercidiphyllum magnificum (Central Japan). We dated this speciation event to the mid-Miocene, and the intraspecific lineage divergence of C. japonicum (China vs Japan) to the Early Pliocene. Throughout climatic upheavals of the late Tertiary/Quaternary, population bottlenecks greatly reduced the genetic diversity of C. japonicum. However, this polymorphism loss was likely counteracted by, first, long-term balancing selection at multiple chromosomal and heterozygous gene regions, potentially reflecting overdominance, and, second, selective sweeps at stress response and growth-related genes likely involved in local adaptation. Our findings contribute to a better understanding of how living fossils have survived climatic upheaval and maintained an extensive geographic range; that is, both types of selection could be major factors contributing to the species' survival, resilience, and persistence.

Ten Years of Provenance Trials and Application of Multivariate Random Forests Predicted the Most Preferable Seed Source for Silviculture of Abies sachalinensis in Hokkaido, Japan

Research highlights: Using 10-year tree height data obtained after planting from the range-wide provenance trials of Abies sachalinensis, we constructed multivariate random forests (MRF), a machine learning algorithm, with climatic variables. The constructed MRF enabled prediction of the optimum seed source to achieve good performance in terms of height growth at every planting site on a fine scale. Background and objectives: Because forest tree species are adapted to the local environment, local seeds are empirically considered as the best sources for planting. However, in some cases, local seed sources show lower performance in height growth than that showed by non-local seed sources. Tree improvement programs aim to identify seed sources for obtaining high-quality timber products by performing provenance trials. Materials and methods: Range-wide provenance trials for one of the most important silvicultural species, Abies sachalinensis, were established in 1980 at nine transplanting experimental sites. We constructed an MRF to estimate the responses of tree height at 10 years after planting at eight climatic variables at 1 km × 1 km resolution. The model was applied for prediction of tree height throughout Hokkaido Island. Results: Our model showed that four environmental variables were major factors affecting height growth—winter solar radiation, warmth index, maximum snow depth, and spring solar radiation. A tree height prediction map revealed that local seeds showed the best performance except in the southernmost region and several parts of northern regions. Moreover, the map of optimum seed provenance suggested that deployment of distant seed sources can outperform local sources in the southernmost and northern regions. Conclusions: We predicted that local seeds showed optimum growth, whereas non-local seeds had the potential to outperform local seeds in some regions. Several deployment options were proposed to improve tree growth.

Genetic structure of Australian glass shrimp, Paratya australiensis, in relation to altitude

Paratya australiensis Kemp (Decapoda: Atyidae) is a widely distributed freshwater shrimp in eastern Australia. The species has been considered as an important stream organism for studying genetics, dispersal, biology, behaviour and evolution in atyids and is a major food source for stream dwelling fishes. Paratya australiensis is a cryptic species complex consisting of nine highly divergent mitochondrial DNA lineages. Previous studies in southeast Queensland showed that “lineage 4” favours upstream sites at higher altitudes, with cooler water temperatures. This study aims to identify putative selection and population structure between high elevation and low elevation populations of this lineage at relatively small spatial scales. Sample localities were selected from three streams: Booloumba Creek, Broken Bridge Creek and Obi Obi Creek in the Conondale Range, southeast Queensland. Six sample localities, consisting of 142 individuals in total were sequenced using double digest Restriction Site Associated DNA-sequencing (ddRAD-seq) technique. Among the 142 individuals, 131 individuals shared 213 loci. Outlier analysis on 213 loci showed that 27 loci were putatively under selection between high elevation and low elevation populations. Outlier analysis on individual streams was also done to test for parallel patterns of adaptation, but there was no evidence of a parallel pattern. Population structure was observed using both the 27 outliers and 186 neutral loci and revealed similar population structure in both cases. Therefore, we cannot differentiate between selection and drift here. The highest genetic differentiation was observed between high elevation and low elevation populations of Booloumba Creek, with small levels of differentiation in the other two streams.

SLAF-seq Uncovers the Genetic Diversity and Adaptation of Chinese Elm (Ulmus parvifolia) in Eastern China

The Chinese elm is an important tree ecologically; however, little is known about its genetic diversity and adaptation mechanisms. In this study, a total of 107 individuals collected from seven natural populations in eastern China were investigated by specific locus amplified fragment sequencing (SLAF-seq). Based on the single nucleotide polymorphisms (SNPs) detected by SLAF-seq, genetic diversity and markers associated with climate variables were identified. All seven populations showed medium genetic diversity, with PIC values ranging from 0.2632 to 0.2761. AMOVA and Fst indicated that a low genetic differentiation existed among populations. Environmental association analyses with three climate variables (annual rainfall, annual average temperature, and altitude) resulted in, altogether, 43 and 30 putative adaptive loci by Bayenv2 and LFMM, respectively. Five adaptive genes were annotated, which were related to the functions of glycosylation, peroxisome synthesis, nucleic acid metabolism, energy metabolism, and signaling. This study was the first on the genetic diversity and local adaptation in Chinese elms, and the results will be helpful in future work on molecular breeding.

Gene Flow and Genetic Variation Explain Signatures of Selection across a Climate Gradient in Two Riparian Species

Many species occur across environmental gradients and it is expected that these species will exhibit some signals of adaptation as heterogeneous environments and localized gene flow may facilitate local adaptation. While riparian zones can cross climate gradients, many of which are being impacted by climate change, they also create microclimates for the vegetation, reducing environmental heterogeneity. Species with differing distributions in these environments provide an opportunity to investigate the importance of genetic connectivity in influencing signals of adaptation over relatively short geographical distance. Association analysis with genomic data was used to compare signals of selection to climate variables in two species that have differing distributions along a river traversing a climate gradient. Results demonstrate links between connectivity, standing genetic variation, and the development of signals of selection. In the restricted species, the combination of high gene flow in the middle and lower catchment and occurrence in a microclimate created along riverbanks likely mitigated the development of selection to most climatic variables. In contrast the more widely distributed species with low gene flow showed a stronger signal of selection. Together these results strengthen our knowledge of the drivers and scale of adaptation and reinforce the importance of connectivity across a landscape to maintain adaptive potential of plant species.

Inferring the demographic history of Japanese cedar, Cryptomeria japonica, using amplicon sequencing

The evolution of a species depends on multiple forces, such as demography and natural selection. To understand the trajectory and driving forces of evolution of a target species, it is first necessary to uncover that species' population history, such as past and present population sizes, subdivision and gene flow, by using appropriate genetic markers. Cryptomeria japonica is a long-lived monoecious conifer species that is distributed in Japan. There are two main lines (omote-sugi and ura-sugi), which are distinguished by apparent differences in morphological traits that may have contributed to their local adaptation. The evolution of these morphological traits seems to be related to past climatic changes in East Asia, but no precise estimate is available for the divergence time of these two lines and the subsequent population dynamics in this species. Here, we analyzed the nucleotide variations at 120 nuclear genes in 94 individuals by using amplicon sequencing in combination with high-throughput sequencing technologies. Our analysis indicated that the population on Yakushima Island, the southern distribution limit of C. japonica in Japan, diverged from the other populations 0.85 million years ago (MYA). The divergence time of the other populations on mainland Japan was estimated to be 0.32 MYA suggesting that the divergence of omote-sugi and ura-sugi might have occurred before the last glacial maximum. Although we found modest levels of gene flow between the present populations, the long-term isolation and environmental heterogeneity caused by climatic changes might have contributed to the differentiation of the lines and their local adaptation.

Genome-wide SNP analyses reveal high gene flow and signatures of local adaptation among the scalloped spiny lobster (Panulirus homarus) along the Omani coastline

BACKGROUND

The scalloped spiny lobster (Panulirus homarus) is a popular seafood commodity worldwide and an important export item from Oman. Annual catches in commercial fisheries are in serious decline, which has resulted in calls for the development of an integrated stock management approach. In Oman, the scalloped spiny lobster is currently treated as a single management unit (MU) or stock and there is an absence of information on the genetic population structure of the species that can inform management decisions, particularly at a fine-scale level. This work is the first to identify genome-wide single nucleotide polymorphisms (SNPs) for P. homarus using Diversity Arrays Technology sequencing (DArT-seq) and to elucidate any stock structure in the species.

RESULTS

After stringent filtering, 7988 high utility SNPs were discovered and used to assess the genetic diversity, connectivity and structure of P. homarus populations from Al Ashkharah, Masirah Island, Duqm, Ras Madrakah, Haitam, Ashuwaymiyah, Mirbat and Dhalkut landing sites. Pairwise FST estimates revealed low differentiation among populations (pairwise FST range = - 0.0008 - 0.0021). Analysis of genetic variation using putatively directional FST outliers (504 SNPs) revealed higher and significant pairwise differentiation (p < 0.01) for all locations, with Ashuwaymiyah being the most diverged population (Ashuwaymiyah pairwise FST range = 0.0288-0.0736). Analysis of population structure using Discriminant Analysis of Principal Components (DAPC) revealed a broad admixture among P. homarus, however, Ashuwaymiyah stock appeared to be potentially under local adaptive pressures. Fine scale analysis using Netview R provided further support for the general admixture of P. homarus.

CONCLUSIONS

Findings here suggested that stocks of P. homarus along the Omani coastline are admixed. Yet, fishery managers need to treat the lobster stock from Ashuwaymiyah with caution as it might be subject to local adaptive pressures. We emphasize further study with larger number of samples to confirm the genetic status of the Ashuwaymiyah stock. The approach utilised in this study has high transferability in conservation and management of other marine stocks with similar biological and ecological attributes.

Inferring spatial patterns and drivers of population divergence of Neolitsea sericea (Lauraceae), based on molecular phylogeography and landscape genomics

The relative roles of geography, climate and ecology in driving population divergence and (incipient) speciation has so far been largely neglected in studies addressing the evolution of East Asia's island flora. Here, we employed chloroplast and ribosomal DNA sequences and restriction site-associated DNA sequencing (RADseq) loci to investigate the phylogeography and drivers of population divergence of Neolitsea sericea. These data sets support the subdivision of N. sericea populations into the Southern and Northern lineages across the 'Tokara gap'. Two distinct sublineages were further identified for the Northern lineage of N. sericea from the RADseq data. RADseq was also used along with approximate Bayesian computation to show that the current distribution and differentiation of N. sericea populations resulted from a combination of relatively ancient migration and successive vicariant events that likely occurred during the mid to late Pleistocene. Landscape genomic analyses showed that, apart from geographic barriers, barrier, potentially local adaptation to different climatic conditions appears to be one of the major drivers for lineage diversification of N. sericea.

Inferences of population structure and demographic history for Taxodium distichum, a coniferous tree in North America, based on amplicon sequencing analysis

PREMISE OF THE STUDY

Studies of natural genetic variation can elucidate the genetic basis of phenotypic variation and the past population structure of species. Our study species, Taxodium distichum, is a unique conifer that inhabits the flood plains and swamps of North America. Morphological and ecological differences in two varieties, T. distichum var. distichum (bald cypress) and T. distichum var. imbricarium (pond cypress), are well known, but little is known about the level of genetic differentiation between the varieties and the demographic history of local populations.

METHODS

We analyzed nucleotide polymorphisms at 47 nuclear loci from 96 individuals collected from the Mississippi River Alluvial Valley (MRAV), and Gulf Coastal populations in Texas, Louisiana, and Florida using high-throughput DNA sequencing. Standard population genetic statistics were calculated, and demographic parameters were estimated using a composite-likelihood approach.

KEY RESULTS

Taxodium distichum in North America can be divided into at least three genetic groups, bald cypress in the MRAV and Texas, bald cypress in Florida, and pond cypress in Florida. The levels of genetic differentiation among the groups were low but significant. Several loci showed the signatures of positive selection, which might be responsible for local adaptation or varietal differentiation.

CONCLUSIONS

Bald cypress was genetically differentiated into two geographical groups, and the boundary was located between the MRAV and Florida. This differentiation could be explained by population expansion from east to west. Despite the overlap of the two varieties' ranges, they were genetically differentiated in Florida. The estimated demographic parameters suggested that pond cypress split from bald cypress during the late Miocene.

Local adaptation to temperature and precipitation in naturally fragmented populations of Cephalotaxus oliveri, an endangered conifer endemic to China

Cephalotaxus oliveri is an endangered tertiary relict conifer endemic to China. The species survives in a wide range from west to east with heterogeneous climatic conditions. Precipitation and temperature are main restrictive factors for distribution of C. oliveri. In order to comprehend the mechanism of adaptive evolution to climate variation, we employed ISSR markers to detect adaptive evolution loci, to identify the association between variation in temperature and precipitation and adaptive loci, and to investigate the genetic structure for 22 C. oliveri natural populations. In total, 14 outlier loci were identified, of which five were associated with temperature and precipitation. Among outlier loci, linkage disequilibrium (LD) was high (42.86%), which also provided strong evidence for selection. In addition, C. oliveri possessed high genetic variation (93.31%) and population differentiation, which may provide raw material to evolution and accelerate local adaptation, respectively. Ecological niche modeling showed that global warming will cause a shift for populations of C. oliveri from south to north with a shrinkage of southern areas. Our results contribute to understand the potential response of conifers to climatic changes, and provide new insights for conifer resource management and conservation strategies.

Evolutionary Quantitative Genomics of Populus trichocarpa

Forest trees generally show high levels of local adaptation and efforts focusing on understanding adaptation to climate will be crucial for species survival and management. Here, we address fundamental questions regarding the molecular basis of adaptation in undomesticated forest tree populations to past climatic environments by employing an integrative quantitative genetics and landscape genomics approach. Using this comprehensive approach, we studied the molecular basis of climate adaptation in 433 Populus trichocarpa (black cottonwood) genotypes originating across western North America. Variation in 74 field-assessed traits (growth, ecophysiology, phenology, leaf stomata, wood, and disease resistance) was investigated for signatures of selection (comparing Q_ST -F_ST) using clustering of individuals by climate of origin (temperature and precipitation). 29,354 SNPs were investigated employing three different outlier detection methods and marker-inferred relatedness was estimated to obtain the narrow-sense estimate of population differentiation in wild populations. In addition, we compared our results with previously assessed selection of candidate SNPs using the 25 topographical units (drainages) across the P. trichocarpa sampling range as population groupings. Narrow-sense Q_ST for 53% of distinct field traits was significantly divergent from expectations of neutrality (indicating adaptive trait variation); 2,855 SNPs showed signals of diversifying selection and of these, 118 SNPs (within 81 genes) were associated with adaptive traits (based on significant Q_ST). Many SNPs were putatively pleiotropic for functionally uncorrelated adaptive traits, such as autumn phenology, height, and disease resistance. Evolutionary quantitative genomics in P. trichocarpa provides an enhanced understanding regarding the molecular basis of climate-driven selection in forest trees and we highlight that important loci underlying adaptive trait variation also show relationship to climate of origin. We consider our approach the most comprehensive, as it uncovers the molecular mechanisms of adaptation using multiple methods and tests. We also provide a detailed outline of the required analyses for studying adaptation to the environment in a population genomics context to better understand the species’ potential adaptive capacity to future climatic scenarios.

Population Structure, Genetic Variation, and Linkage Disequilibrium in Perennial Ryegrass Populations Divergently Selected for Freezing Tolerance

Low temperature is one of the abiotic stresses seriously affecting the growth of perennial ryegrass (Lolium perenne L.), and freezing tolerance is a complex trait of major agronomical importance in northern and central Europe. Understanding the genetic control of freezing tolerance would aid in the development of cultivars of perennial ryegrass with improved adaptation to frost. The plant material investigated in this study was an experimental synthetic population derived from pair-crosses among five European perennial ryegrass genotypes, representing adaptations to a range of climatic conditions across Europe. A total number of 80 individuals (24 of High frost [HF]; 29 of Low frost [LF], and 27 of Unselected [US]) from the second generation of the two divergently selected populations and an unselected (US) control population were genotyped using 278 genome-wide SNPs derived from perennial ryegrass transcriptome sequences. Our studies investigated the genetic diversity among the three experimental populations by analysis of molecular variance and population structure, and determined that the HF and LF populations are very divergent after selection for freezing tolerance, whereas the HF and US populations are more similar. Linkage disequilibrium (LD) decay varied across the seven chromosomes and the conspicuous pattern of LD between the HF and LF population confirmed their divergence in freezing tolerance. Furthermore, two F st outlier methods; finite island model (fdist) by LOSITAN and hierarchical structure model using ARLEQUIN, both detected six loci under directional selection. These outlier loci are most probably linked to genes involved in freezing tolerance, cold adaptation, and abiotic stress. These six candidate loci under directional selection for freezing tolerance might be potential marker resources for breeding perennial ryegrass cultivars with improved freezing tolerance.

Neutral and Selective Processes Drive Population Differentiation for Iris hexagona

Gene flow among widespread populations can be reduced by geographical distance or by divergent selection resulting from local adaptation. In this study, we tested for the divergence of phenotypes and genotypes among 8 populations of Iris hexagona. Using a genotyping-by-sequencing approach, we generated a panel of 750 single nucleotide polymorphisms (SNPs) and used population genetic analyses to determine what may affect patterns of divergence across I. hexagona populations. Specifically, genetic differentiation was compared between populations at neutral and nonneutral SNPs and detected significant differences between the 2 types of markers. We then asked whether loci with the strongest degree of population genetic differentiation were also the loci with the strongest association to morphology or climate differences, allowing us to test if pollinators or climate drive population differentiation or some combination of both. We found 2 markers that were associated with morphology and 1 marker associated with 2 of the environmental variables, which were also identified in the outlier analysis. We then show that the SNPs putatively under selection were positively correlated with both geographic distance and phenotypic distance, albeit weakly to phenotypic distance. Moreover, neutral SNPs were only correlated with geographic distance and thus isolation-by-distance was observed for neutral SNPs. Our data suggest that both deterministic and neutral processes have contributed to the evolutionary trajectory of I. hexagona populations.

Evidence for cryptic northern refugia in the last glacial period in Cryptomeria japonica

BACKGROUND AND AIMS

Distribution shifts and natural selection during past climatic changes are important factors in determining the genetic structure of forest species. In particular, climatic fluctuations during the Quaternary appear to have caused changes in the distribution ranges of plants, and thus strongly affected their genetic structure. This study was undertaken to identify the responses of the conifer Cryptomeria japonica, endemic to the Japanese Archipelago, to past climatic changes using a combination of phylogeography and species distribution modelling (SDM) methods. Specifically, this study focused on the locations of refugia during the last glacial maximum (LGM).

METHODS

Genetic diversity and structure were examined using 20 microsatellite markers in 37 populations of C. japonica. The locations of glacial refugia were assessed using STRUCTURE analysis, and potential habitats under current and past climate conditions were predicted using SDM. The process of genetic divergence was also examined using the approximate Bayesian computation procedure (ABC) in DIY ABC to test the divergence time between the gene pools detected by the STRUCTURE analysis.

KEY RESULTS

STRUCTURE analysis identified four gene pools: northern Tohoku district; from Chubu to Chugoku district; from Tohoku to Shikoku district on the Pacific Ocean side of the Archipelago; and Yakushima Island. DIY ABC analysis indicated that the four gene pools diverged at the same time before the LGM. SDM also indicated potential northern cryptic refugia.

CONCLUSIONS

The combined evidence from microsatellites and SDM clearly indicates that climatic changes have shaped the genetic structure of C. japonica. The gene pool detected in northern Tohoku district is likely to have been established by cryptic northern refugia on the coast of the Japan Sea to the west of the Archipelago. The gene pool in Yakushima Island can probably be explained simply by long-term isolation from the other gene pools since the LGM. These results are supported by those of SDM and the predicted divergence time determined using ABC analysis.

On the prospect of identifying adaptive loci in recently bottlenecked populations

Identifying adaptively important loci in recently bottlenecked populations – be it natural selection acting on a population following the colonization of novel habitats in the wild, or artificial selection during the domestication of a breed – remains a major challenge. Here we report the results of a simulation study examining the performance of available population-genetic tools for identifying genomic regions under selection. To illustrate our findings, we examined the interplay between selection and demography in two species of Peromyscus mice, for which we have independent evidence of selection acting on phenotype as well as functional evidence identifying the underlying genotype. With this unusual information, we tested whether population-genetic-based approaches could have been utilized to identify the adaptive locus. Contrary to published claims, we conclude that the use of the background site frequency spectrum as a null model is largely ineffective in bottlenecked populations. Results are quantified both for site frequency spectrum and linkage disequilibrium-based predictions, and are found to hold true across a large parameter space that encompasses many species and populations currently under study. These results suggest that the genomic footprint left by selection on both new and standing variation in strongly bottlenecked populations will be difficult, if not impossible, to find using current approaches.

Genetic differentiation and evolutionary adaptation in Cryptomeria japonica

Local adaptation of plant species is a central issue for survival during global climate change, especially for long-lived forest trees, with their lengthy regeneration time and spatially limited gene flow. Identification of loci and/or genomic regions associated with local adaptation is necessary for knowledge of both evolution and molecular breeding for climate change. Cryptomeria japonica is an important species for forestry in Japan; it has a broad natural distribution and can survive in a range of different environments. The genetic structure of 14 natural populations of this species was investigated using 3930 SNP markers. Populations on the Pacific Ocean side of Japan are clearly different from those on the Japan Sea side, as discussed in previous studies. Structure analysis and population network trees show that peripheral populations, including the most northerly and southerly ones, have unique features. We found that the genetic differentiation coefficient is low, F_ST = 0.05, although it must account for the presence of important genes associated with adaptation to specific environments. In total, 208 outlier loci were detected, of which 43 were associated with environmental variables. Four clumped regions of outlier loci were detected in the genome by linkage analysis. Linkage disequilibrium (LD) was quite high in these clumps of outlier loci, which were found in linkage groups (LGs) 2, 7, 10, and 11, especially between populations of two varieties, and when interchromosomal LD was also detected. The LG7 region is characteristic of the Yakushima population, which is a large, isolated, peripheral population occupying a specific environment resulting from isolation combined with volcanic activity in the region. The detected LD may provide strong evidence for selection between varieties.

Temporally isolated lineages of pink salmon reveal unique signatures of selection on distinct pools of standing genetic variation

A species' genetic diversity bears the marks of evolutionary processes that have occurred throughout its history. However, robust detection of selection in wild populations is difficult and often impeded by lack of replicate tests. Here, we investigate selection in pink salmon (Oncorhynchus gorbuscha) using genome scans coupled with inference from a haploid-assisted linkage map. Pink salmon have a strict 2-year semelparous life history which has resulted in temporally isolated (allochronic) lineages that remain sympatric through sharing of spawning habitats in alternate years. The lineages differ in a range of adaptive traits, suggesting different genetic backgrounds. We used genotyping by sequencing of haploids to generate a high-density linkage map with 7035 loci and screened an existing panel of 8036 loci for signatures of selection. The linkage map enabled identification of novel genomic regions displaying signatures of parallel selection shared between lineages. Furthermore, 24 loci demonstrated divergent selection and differences in genetic diversity between lineages, suggesting that adaptation in the 2 lineages has arisen from different pools of standing genetic variation. Findings have implications for understanding asynchronous population abundances as well as predicting future ecosystem impacts from lineage-specific responses to climate change.

Cross-species outlier detection reveals different evolutionary pressures between sister species

Lodgepole pine (Pinus contorta var. latifolia) and jack pine (Pinus banksiana) hybridize in western Canada, an area of recent mountain pine beetle range expansion. Given the heterogeneity of the environment, and indications of local adaptation, there are many unknowns regarding the response of these forests to future outbreaks. To better understand this we aim to identify genetic regions that have adaptive potential. We used data collected on 472 single nucleotide polymorphism (SNP) loci from 576 tree samples collected across 13 lodgepole pine-dominated sites and four jack pine-dominated sites. We looked at the relationship of genetic diversity with the environment, and we identified candidate loci using both frequency-based (arlequin and bayescan) and correlation-based (matsam and bayenv) methods. We found contrasting relationships between environmental variation and genetic diversity for the species. While we identified a number of candidate outliers (34 in lodgepole pine, 25 in jack pine, and 43 interspecific loci), we did not find any loci in common between lodgepole and jack pine. Many of the outlier loci identified were correlated with environmental variation. Using rigorous criteria we have been able to identify potential outlier SNPs. We have also found evidence of contrasting environmental adaptations between lodgepole and jack pine which could have implications for beetle spread risk.

Detecting SNPs underlying domestication-related traits in soybean

BACKGROUND

Cultivated soybean (Glycine max) experienced a severe genetic bottleneck during its domestication and a further loss in diversity during its subsequent selection. Here, a panel of 65 wild (G. soja) and 353 cultivated accessions was genotyped at 552 single-nucleotide polymorphism loci to search for signals of selection during and after domestication.

RESULTS

The wild and cultivated populations were well differentiated from one another. Application of the Fst outlier test revealed 64 loci showing evidence for selection. Of these, 35 related to selection during domestication, while the other 29 likely gradually became monomorphic as a result of prolonged selection during post domestication. Two of the SNP locus outliers were associated with testa color.

CONCLUSIONS

Identifying genes controlling domestication-related traits is important for maintaining the diversity of crops. SNP locus outliers detected by a combined forward genetics and population genetics approach can provide markers with utility for the conservation of wild accessions and for trait improvement in the cultivated genepool.

A dense linkage map for Chinook salmon (Oncorhynchus tshawytscha) reveals variable chromosomal divergence after an ancestral whole genome duplication event

Comparisons between the genomes of salmon species reveal that they underwent extensive chromosomal rearrangements following whole genome duplication that occurred in their lineage 58−63 million years ago. Extant salmonids are diploid, but occasional pairing between homeologous chromosomes exists in males. The consequences of re-diploidization can be characterized by mapping the position of duplicated loci in such species. Linkage maps are also a valuable tool for genome-wide applications such as genome-wide association studies, quantitative trait loci mapping or genome scans. Here, we investigated chromosomal evolution in Chinook salmon (Oncorhynchus tshawytscha) after genome duplication by mapping 7146 restriction-site associated DNA loci in gynogenetic haploid, gynogenetic diploid, and diploid crosses. In the process, we developed a reference database of restriction-site associated DNA loci for Chinook salmon comprising 48528 non-duplicated loci and 6409 known duplicated loci, which will facilitate locus identification and data sharing. We created a very dense linkage map anchored to all 34 chromosomes for the species, and all arms were identified through centromere mapping. The map positions of 799 duplicated loci revealed that homeologous pairs have diverged at different rates following whole genome duplication, and that degree of differentiation along arms was variable. Many of the homeologous pairs with high numbers of duplicated markers appear conserved with other salmon species, suggesting that retention of conserved homeologous pairing in some arms preceded species divergence. As chromosome arms are highly conserved across species, the major resources developed for Chinook salmon in this study are also relevant for other related species.

Genome-wide distribution of genetic diversity and linkage disequilibrium in a mass-selected population of maritime pine

Background

The accessibility of high-throughput genotyping technologies has contributed greatly to the development of genomic resources in non-model organisms. High-density genotyping arrays have only recently been developed for some economically important species such as conifers. The potential for using genomic technologies in association mapping and breeding depends largely on the genome wide patterns of diversity and linkage disequilibrium in current breeding populations. This study aims to deepen our knowledge regarding these issues in maritime pine, the first species used for reforestation in south western Europe.

Results

Using a new map merging algorithm, we first established a 1,712 cM composite linkage map (comprising 1,838 SNP markers in 12 linkage groups) by bringing together three already available genetic maps. Using rigorous statistical testing based on kernel density estimation and resampling we identified cold and hot spots of recombination. In parallel, 186 unrelated trees of a mass-selected population were genotyped using a 12k-SNP array. A total of 2,600 informative SNPs allowed to describe historical recombination, genetic diversity and genetic structure of this recently domesticated breeding pool that forms the basis of much of the current and future breeding of this species. We observe very low levels of population genetic structure and find no evidence that artificial selection has caused a reduction in genetic diversity. By combining these two pieces of information, we provided the map position of 1,671 SNPs corresponding to 1,192 different loci. This made it possible to analyze the spatial pattern of genetic diversity (H_e) and long distance linkage disequilibrium (LD) along the chromosomes. We found no particular pattern in the empirical variogram of H_e across the 12 linkage groups and, as expected for an outcrossing species with large effective population size, we observed an almost complete lack of long distance LD.

Conclusions

These results are a stepping stone for the development of strategies for studies in population genomics, association mapping and genomic prediction in this economical and ecologically important forest tree species.

Environmental versus geographical determinants of genetic structure in two subalpine conifers

Alpine ecosystems are facing rapid human-induced environmental changes, and so more knowledge about tree adaptive potential is needed. This study investigated the relative role of isolation by distance (IBD) versus isolation by adaptation (IBA) in explaining population genetic structure in Abies alba and Larix decidua, based on 231 and 233 single nucleotide polymorphisms (SNPs) sampled across 36 and 22 natural populations, respectively, in the Alps and Apennines. Genetic structure was investigated for both geographical and environmental groups, using analysis of molecular variance (AMOVA). For each species, nine environmental groups were defined using climate variables selected from a multiple factor analysis. Complementary methods were applied to identify outliers based on these groups, and to test for IBD versus IBA. AMOVA showed weak but significant genetic structure for both species, with higher values in L. decidua. Among the potential outliers detected, up to two loci were found for geographical groups and up to seven for environmental groups. A stronger effect of IBD than IBA was found in both species; nevertheless, once spatial effects had been removed, temperature and soil in A. alba, and precipitation in both species, were relevant factors explaining genetic structure. Based on our findings, in the Alpine region, genetic structure seems to be affected by both geographical isolation and environmental gradients, creating opportunities for local adaptation.

Demonstration of genome-wide association studies for identifying markers for wood property and male strobili traits in Cryptomeria japonica

Genome-wide association studies (GWAS) are an alternative to bi-parental QTL mapping in long-lived perennials. In the present study, we examined the potential of GWAS in conifers using 367 unrelated plus trees of Cryptomeria japonica D. Don, which is the most widely planted and commercially important tree species in Japan, and tried to detect significant associations between wood property traits and quantity of male strobili on the one hand, and 1,032 single nucleotide polymorphisms (SNPs) assigned to 1,032 genes on the other. Association analysis was performed with the mixed linear model taking into account kinship relationships and subpopulation structure. In total, 6 SNPs were found to have significant associations with the variations in phenotype. These SNPs were not associated with the positions of known genes and QTLs that have been reported to date, thus they may identify novel QTLs. These 6 SNPs were all found in sequences showing similarities with known genes, although further analysis is required to dissect the ways in which they affect wood property traits and abundance of male strobili. These presumptive QTL loci provide opportunities for improvement of C. japonica, based on a marker approach. The results suggest that GWAS has potential for use in future breeding programs in C. japonica.

Signatures of selection in the Iberian honey bee (Apis mellifera iberiensis) revealed by a genome scan analysis of single nucleotide polymorphisms

Understanding the genetic mechanisms of adaptive population divergence is one of the most fundamental endeavours in evolutionary biology and is becoming increasingly important as it will allow predictions about how organisms will respond to global environmental crisis. This is particularly important for the honey bee, a species of unquestionable ecological and economical importance that has been exposed to increasing human-mediated selection pressures. Here, we conducted a single nucleotide polymorphism (SNP)-based genome scan in honey bees collected across an environmental gradient in Iberia and used four FST -based outlier tests to identify genomic regions exhibiting signatures of selection. Additionally, we analysed associations between genetic and environmental data for the identification of factors that might be correlated or act as selective pressures. With these approaches, 4.4% (17 of 383) of outlier loci were cross-validated by four FST -based methods, and 8.9% (34 of 383) were cross-validated by at least three methods. Of the 34 outliers, 15 were found to be strongly associated with one or more environmental variables. Further support for selection, provided by functional genomic information, was particularly compelling for SNP outliers mapped to different genes putatively involved in the same function such as vision, xenobiotic detoxification and innate immune response. This study enabled a more rigorous consideration of selection as the underlying cause of diversity patterns in Iberian honey bees, representing an important first step towards the identification of polymorphisms implicated in local adaptation and possibly in response to recent human-mediated environmental changes.

Imprints of natural selection along environmental gradients in phenology-related genes of Quercus petraea

We explored single nucleotide polymorphism (SNP) variation in candidate genes for bud burst from Quercus petraea populations sampled along gradients of latitude and altitude in Western Europe. SNP diversity was monitored for 106 candidate genes, in 758 individuals from 32 natural populations. We investigated whether SNP variation reflected the clinal pattern of bud burst observed in common garden experiments. We used different methods to detect imprints of natural selection (FST outlier, clinal variation at allelic frequencies, association tests) and compared the results obtained for the two gradients. FST outlier SNPs were found in 15 genes, 5 of which were common to both gradients. The type of selection differed between the two gradients (directional or balancing) for 3 of these 5. Clinal variations were observed for six SNPs, and one cline was conserved across both gradients. Association tests between the phenotypic or breeding values of trees and SNP genotypes identified 14 significant associations, involving 12 genes. The results of outlier detection on the basis of population differentiation or clinal variation were not very consistent with the results of association tests. The discrepancies between these approaches may reflect the different hierarchical levels of selection considered (inter- and intrapopulation selection). Finally, we obtained evidence for convergent selection (similar for gradients) and clinal variation for a few genes, suggesting that comparisons between parallel gradients could be used to screen for major candidate genes responding to natural selection in trees.