Genome sequences and population genomics reveal climatic adaptation and genomic divergence between two closely related sweetgum species.

Understanding the genetic basis of population divergence and adaptation is an important goal in population genetics and evolutionary biology. However, the relative roles of demographic history, gene flow, and/or selective regime in driving genomic divergence, climatic adaptation, and speciation in non-model tree species are not yet fully understood. To address this issue, we generated whole-genome resequencing data of Liquidambar formosana and L. acalycina, which are broadly sympatric but altitudinally segregated in the Tertiary relict forests of subtropical China. We integrated genomic and environmental data to investigate the demographic history, genomic divergence, and climatic adaptation of these two sister species. We inferred a scenario of allopatric species divergence during the late Miocene, followed by secondary contact during the Holocene. We identified multiple genomic islands of elevated divergence that mainly evolved through divergence hitchhiking and recombination rate variation, likely fostered by long-term refugial isolation and recent differential introgression in low-recombination genomic regions. We also found some candidate genes with divergent selection signatures potentially involved in climatic adaptation and reproductive isolation. Our results contribute to a better understanding of how late Tertiary/Quaternary climatic change influenced speciation, genomic divergence, climatic adaptation, and introgressive hybridization in East Asia's Tertiary relict flora. In addition, they should facilitate future evolutionary, conservation genomics, and molecular breeding studies in Liquidambar, a genus of important medicinal and ornamental values.

Responses of selected plasma metabolites to a 2-day nutritional challenge of goats divergently selected for functional longevity.

Understanding the extent to which genetics × environment plays a role in shaping individual strategies to environmental challenges is of considerable interest for future selection of more resilient animals. The objective of this study was to evaluate the metabolic responses of goats, selected for functional longevity, to a nutritional challenge. This evaluation was based on plasma metabolites and was repeated across 2 experimental farms over 2 years. We carried out 6 different experimental trials from 2018 to 2022 (4 trials on site Bourges, 2018-2021 and 2 trials on site Grignon, 2021-2022) in which 267 first kidding goats, daughters of Alpine bucks divergently selected for functional longevity, longevity plus (n = 137), and longevity minus (n = 130), were exposed to a 2-d nutritional challenge in early lactation. The experiments consisted of a 5 or 7-d control period (prechallenge) on a standard lactation diet followed by a 2-d nutritional challenge with straw-only feeding, and then a 7- or 10-d recovery period on a standard lactation diet, for site Bourges and Grignon, respectively. During the challenge plasma metabolite composition was measured. Linear mixed models were used to analyze all traits, considering the individual as a random effect and the 2 × 2 treatments (i.e., genetic line and year nested in site) and litter size as fixed effects. The linear mixed-effects model using a piecewise arrangement was used to analyze the response/recovery profiles to the nutritional challenge. The estimates of the random effect for each individual, using the mixed models without the fixed effects of genetic line, were used in a sparse partial least square discriminant analysis (sPLS-DA) to compare the goat metabolism response to the challenge on a multivariate scale. The plasma glucose, BHB, nonesterified fatty acids (NEFA), and urea concentrations significantly responded to the 2-d nutritional challenge. Selection for functional longevity did not affect plasma glucose, NEFA, BHB, and urea response/recoveries to a 2-d nutritional challenge. However, site, trial, and litter size significantly affected these responses. Moreover, the plasma metabolites seem not to fully recover to prechallenge levels after the recovery phase. The sPLS-DA analysis did not discriminate between the 2 longevity lines. We observed meaningful between-individuals variability in plasma BHB, especially on the prechallenge and rate of response and rate of recovery from the 2-d nutritional challenge (CV = 26.2%, 36.1%, and 41.2%, repeatability = 0.749, 0.322, and 0.741, respectively). Plasma NEFA recovery from challenge also demonstrated high between-individuals variability (CV = 16.4%, repeatability = 0.323). Selection for functional longevity did not affect plasma metabolites responses to a 2-d nutritional challenge in dairy goats. Plasma NEFA and BHB response and recovery presented high between-individuals variability, indicating that individual adaptative characteristics to nutritional challenges are not related to the environmental conditions but to inherent individual characteristics.

Differential patterns in runs of homozygosity in two mice lines under divergent selection for environmental variability for birth weight.

Runs of homozygosity (ROH) are defined as long continuous homozygous stretches in the genome which are assumed to originate from a common ancestor. It has been demonstrated that divergent selection for variability in mice is possible and that low variability in birth weight is associated with robustness. To analyse ROH patterns and ROH-based genomic inbreeding, two mouse lines that were divergently selected for birth weight variability for 26 generations were used, with: 752 individuals for the high variability line (H-Line), 766 individuals for the low variability line (L-Line) and 74 individuals as a reference population. Individuals were genotyped using the high density Affymetrix Mouse Diversity Genotyping Array. ROH were identified using both the sliding windows (SW) and the consecutive runs (CR) methods. Inbreeding coefficients were calculated based on pedigree (FPED ) information, on ROH identified using the SW method (FROHSW ) and on ROH identified using the CR method (FROHCR ). Differences in genomic inbreeding were not consistent across generations and these parameters did not show clear differences between lines. Correlations between FPED and FROH were high, particularly for FROHSW . Moreover, correlations between FROHSW and FPED were even higher when ROH were identified with no restrictions in the number of heterozygotes per ROH. The comparison of FROH estimates between either of the selected lines were based on significant differences at the chromosome level, mainly in chromosomes 3, 4, 6, 8, 11, 15 and 19. ROH-based inbreeding estimates that were computed using longer homozygous segments had a higher relationship with FPED . Differences in robustness between lines were not attributable to a higher homozygosis in the L-Line, but maybe to the different distribution of ROH at the chromosome level between lines. The analysis identified a set of genomic regions for future research to establish the genomic basis of robustness.

Species Persistence with Hybridization in Toad-Headed Lizards Driven by Divergent Selection and Low Recombination.

Speciation plays a central role in evolutionary studies, and particularly how reproductive isolation (RI) evolves. The origins and persistence of RI are distinct processes that require separate evaluations. Treating them separately clarifies the drivers of speciation and then it is possible to link the processes to understand large-scale patterns of diversity. Recent genomic studies have focused predominantly on how species or RI originate. However, we know little about how species persist in face of gene flow. Here, we evaluate a contact zone of two closely related toad-headed lizards (Phrynocephalus) using a chromosome-level genome assembly and population genomics. To some extent, recent asymmetric introgression from Phrynocephalus putjatai to P. vlangalii reduces their genomic differences. However, their highly divergent regions (HDRs) have heterogeneous distributions across the genomes. Functional gene annotation indicates that many genes within HDRs are involved in reproduction and RI. Compared with allopatric populations, contact areas exhibit recent divergent selection on the HDRs and a lower population recombination rate. Taken together, this implies that divergent selection and low genetic recombination help maintain RI. This study provides insights into the genomic mechanisms that drive RI and two species persistence in the face of gene flow during the late stage of speciation.

Genomic data and common garden experiments reveal climate-driven selection on ecophysiological traits in two Mediterranean oaks.

Improving our knowledge of how past climate-driven selection has acted on present-day trait population divergence is essential to understand local adaptation processes and improve our predictions of evolutionary trajectories in the face of altered selection pressures resulting from climate change. In this study, we investigated signals of selection on traits related to drought tolerance and growth rates in two Mediterranean oak species (Quercus faginea and Q. lusitanica) with contrasting distribution ranges and climatic niches. We genotyped 182 individuals from 24 natural populations of the two species using restriction-site-associated DNA sequencing and conducted a thorough functional characterization in 1602 seedlings from 21 populations cultivated in common garden experiments under contrasting watering treatments. Our genomic data revealed that both Q. faginea and Q. lusitanica have very weak population genetic structure, probably as a result of high rates of pollen-mediated gene flow among populations and large effective population sizes. In contrast, common garden experiments showed evidence of climate-driven divergent selection among populations on traits related to leaf morphology, physiology and growth in both species. Overall, our study suggests that climate is an important selective factor for Mediterranean oaks and that ecophysiological traits have evolved in drought-prone environments even in a context of very high rates of gene flow among populations.

Divergent selection and genetic introgression shape the genome landscape of heterosis in hybrid rice.

The successful application of heterosis in hybrid rice has dramatically improved rice productivity, but the genetic mechanism for heterosis in the hybrid rice remains unclear. In this study, we generated two populations of rice F1 hybrids with present-day commercial hybrid parents, genotyped the parents with 50k SNP chip and genome resequencing, and recorded the phenotype of ∼2,000 hybrids at three field trials. By integrating these data with the collected genotypes of ∼4,200 rice landraces and improved varieties that were reported previously, we found that the male and female parents have different levels of genome introgressions from other rice subpopulations, including indica, aus, and japonica, therefore shaping heterotic loci in the hybrids. Among the introgressed exogenous genome, we found that heterotic loci, including Ghd8/DTH8, Gn1a, and IPA1 existed in wild rice, but were significantly divergently selected among the rice subpopulations, suggesting these loci were subject to environmental adaptation. During modern rice hybrid breeding, heterotic loci were further selected by removing loci with negative effect and fixing loci with positive effect and pyramid breeding. Our results provide insight into the genetic basis underlying the heterosis of elite hybrid rice varieties, which could facilitate a better understanding of heterosis and rice hybrid breeding.

HDAC1: An Essential and Conserved Member of the Diverse Zn2+-Dependent HDAC Family Driven by Divergent Selection Pressure.

Zn2+-dependent histone deacetylases (HDACs) are enzymes that regulate gene expression by removing acetyl groups from histone proteins. These enzymes are essential in all living systems, playing key roles in cancer treatment and as potential pesticide targets. Previous phylogenetic analyses of HDAC in certain species have been published. However, their classification and evolutionary origins across biological kingdoms remain unclear, which limits our understanding of them. In this study, we collected the HDAC sequences from 1451 organisms and performed analyses. The HDACs are found to diverge into three classes and seven subclasses under divergent selection pressure. Most subclasses show species specificity, indicating that HDACs have evolved with high plasticity and diversification to adapt to different environmental conditions in different species. In contrast, HDAC1 and HDAC3, belonging to the oldest class, are conserved and crucial in major kingdoms of life, especially HDAC1. These findings lay the groundwork for the future application of HDACs.

Increased Floral Rewards due to Local Adaptation Drives Plant Ecological Speciation via Learned Preferences of Pollinators.

AbstractIn animal-pollinated plants, the growth environment and pollination environment are two important agents of natural selection. However, their simultaneous effects on plant speciation remain underexplored. Here, we report a theoretical finding that if plants' local adaptation to the growth environment increases their floral rewards for pollinators, it can strongly facilitate ecological speciation in plants. We consider two evolving plant traits, vegetative and floral signal traits, in a population genetic model for two plant populations under divergent selection from different growth environments. The vegetative trait determines plants' local adaptation. Locally adapted plants reward pollinators better than maladapted plants. By associative learning, pollinators acquire learned preferences for floral signal traits expressed by better-rewarding plants. If pollinators' learned preferences become divergent between populations, floral signal divergence occurs and plants develop genetic associations between vegetative and floral signal traits, leading to ecological speciation via a two-allele mechanism. Interestingly, speciation is contingent on whether novel floral signal variants arise before or after plant populations become locally adapted to the growth environment. Our results suggest that simultaneous selection from growth and pollination environments might be important for the ecological speciation of animal-pollinated plants.

Sympatric genetic divergence between early- and late-season weedy rice populations.

Sympatric genetic divergence is the most appealing and controversial pattern in the theory of ecological speciation. Examples that support sympatric genetic divergence in plant species are extremely rare. Solid evidence of sympatric genetic divergence will provide deep insights for revealing the underlying mechanisms of ecological speciation. We analysed the total genomic DNA sequences of 120 weedy rice (WR; Oryza sativa f. spontanea) plants, representing three WR population pairs separately from three early- and late-season rice fields, in comparison with those of the co-occurring rice cultivars and other rice materials. We detected substantial genetic divergence within the pairs of the sympatric early- and late-season WR populations, although genetic divergence was unevenly distributed across the genomes. Restricted gene flow was determined between the sympatric WR populations, resulting in their distinct genetic structures. We also detected relatively low genetic diversity that was likely to be associated with stronger selection in early-season WR populations. Our findings provide strong evidence for sympatric genetic divergence between the WR populations in the same fields but in different seasons. We conclude that temporal isolation plays an important role in creating genetic divergence between sympatric populations/species in plants.

Evidence of climate-driven selection on tree traits and trait plasticity across the climatic range of a riparian foundation species.

Selection on quantitative traits by heterogeneous climatic conditions can lead to substantial trait variation across a species range. In the context of rapidly changing environments, however, it is equally important to understand selection on trait plasticity. To evaluate the role of selection in driving divergences in traits and their associated plasticities within a widespread species, we compared molecular and quantitative trait variation in Populus fremontii (Fremont cottonwood), a foundation riparian distributed throughout Arizona. Using SNP data and genotypes from 16 populations reciprocally planted in three common gardens, we first performed QST -FST analyses to detect selection on traits and trait plasticity. We then explored the environmental drivers of selection using trait-climate and plasticity-climate regressions. Three major findings emerged: (1) There was significant genetic variation in traits expressed in each of the common gardens and in the phenotypic plasticity of traits across gardens, both of which were heritable. (2) Based on QST -FST comparisons, there was evidence of selection in all traits measured; however, this result varied from no effect in one garden to highly significant in another, indicating that detection of past selection is environmentally dependent. We also found strong evidence of divergent selection on plasticity across environments for two traits. (3) Traits and/or their plasticity were often correlated with population source climate (R2 up to .77 and .66, respectively). These results suggest that steep climate gradients across the Southwest have played a major role in shaping the evolution of divergent phenotypic responses in populations and genotypes now experiencing climate change.

Habitat-linked genetic variation supports microgeographic adaptive divergence in an island-endemic bird species.

We investigated the potential mechanisms driving habitat-linked genetic divergence within a bird species endemic to a single 250-km2 island. The island scrub-jay (Aphelocoma insularis) exhibits microgeographic divergence in bill morphology across pine-oak ecotones on Santa Cruz Island, California (USA), similar to adaptive differences described in mainland congeners over much larger geographic scales. To test whether individuals exhibit genetic differentiation related to habitat type and divergence in bill length, we genotyped over 3000 single nucleotide polymorphisms in 123 adult island scrub-jay males from across Santa Cruz Island using restriction site-associated DNA sequencing. Neutral landscape genomic analyses revealed that genome-wide genetic differentiation was primarily related to geographic distance and differences in habitat composition. We also found 168 putatively adaptive loci associated with habitat type using multivariate redundancy analysis while controlling for spatial effects. Finally, two genome-wide association analyses revealed a polygenic basis to variation in bill length with multiple loci detected in or near genes known to affect bill morphology in other birds. Our findings support the hypothesis that divergent selection at microgeographic scales can cause adaptive divergence in the presence of ongoing gene flow.

Evolution of physical linkage between loci controlling ecological traits and mating preferences.

Coupling of multiple barriers to gene-flow, such as divergent local adaptation and reproductive isolation, facilitates speciation. However, alleles at loci that contribute to barrier effects can be dissociated by recombination. Models of linkage between diverging alleles often consider elements that reduce recombination, such as chromosomal inversions and alleles that modify recombination rate between existing loci. In contrast, here, we consider the evolution of linkage due to the close proximity of loci on the same chromosome. Examples of such physical linkage exist in several species, but in other cases, strong associations are maintained without physical linkage. We use an individual-based model to study the conditions under which the physical linkage between loci controlling ecological traits and mating preferences might be expected to evolve. We modelled a single locus controlling an ecological trait that acts also as a mating cue. Mating preferences are controlled by multiple loci, formed by mutations that are randomly placed in the "genome", within varying distances from the ecological trait locus, allowing us to examine which genomic architectures spread across the population. Our model reveals that stronger physical linkage is favoured when mating preferences and selection are weaker. Under such conditions mating among divergent phenotypes is more frequent, and matching ecological trait and mating preference alleles are more likely to become dissociated by recombination, favouring the evolution of genetic linkage. While most theoretical studies on clustering of divergent loci focus on how physical linkage influences speciation, we show how physical linkage itself can arise, establishing conditions that can favour speciation.

The effect of seed-dispersal timing on seedling recruitment is modulated by environmental conditions that vary across altitude in a threatened palm.

BACKGROUND AND AIMS: The timing of seed dispersal determines the environmental conditions that plants face during early life stages. In seasonal environments, selection is expected to favour dispersal timing that is matched to environmental conditions suitable for successful recruitment. Our aim here was to test whether the timing of seed dispersal influences seedling establishment success in two populations of Euterpe edulis that are located at contrasting altitudes, have different seed-dispersal phenologies and are subjected to distinct climatic conditions. METHODS: We sowed E. edulis seeds in contrasting altitudes on different dates, and monitored seed germination, emergence and seedling establishment at each altitude over 4 years. At the high-altitude site, five seed-dispersal cohorts were established during the natural dispersal period. At the low-altitude site, three seed-dispersal cohorts were established during natural dispersal, and two were established either before or after natural dispersal. KEY RESULTS: At the high-altitude site, seed-dispersal timing did not affect seed germination, seedling emergence or seedling establishment success. In contrast, at the low-altitude site, late seed dispersal near the end of the wet season resulted in a lower probability of seedling establishment, possibly due to the exposure of seeds, germinants and seedlings to unfavourable drought conditions. In addition, at the low-altitude site, the natural seed-dispersal period was poorly matched to favourable environmental conditions for seedling establishment. CONCLUSIONS: The greater effect of seed-dispersal timing on seedling establishment at the low-altitude site is probably related to a more seasonal and drought-prone environment that favours a restricted period of seed dispersal. The magnitude of the effect of dispersal timing on seedling establishment success was modulated by environmental conditions that vary across altitude. Furthermore, reproductive phenology appears to be subject to more intense selection at the lower limit of the altitudinal range, due to a more restrictive window of opportunity for successful seedling establishment.

Impact of selection for birth weight variability on reproductive longevity: A mice model.

Uniformity, understood as a similar performance in relevant livestock traits, such as birth weight within the litter, is being included as one of the selection objectives in breeding programmes, especially for polytocous livestock species. A divergent selection experiment for birth weight within-litter variability in mice during 23 generations showed that homogeneous animals were better for litter size, survival and feed efficiency but less heavy than heterogeneous animals. The aim of this study was to compare the reproductive longevity in both divergent lines as time to the end of the reproductive period. Two generations from both lines with an initial number of 43 females and 43 males were mated one to one and stayed together to have consecutive parturitions until the end of the reproductive life. Females were discarded when the time elapsed from the last parturition was longer than 63 days. The time to the end of the reproductive period between both lines was compared by fitting a Cox proportional hazard regression model adjusting for line, generation and its interaction. The rate of parturitions in both lines was also compared using a Prentice-Williams-Peterson model adjusted for the same effects. The low variability line was associated with a higher parturition rate, e.g., adjusted hazard ratio was 2.93 (95% CI 2.17-3.94). The Cox model showed that the low variability females also presented benefits of time to the end of the reproductive period, with an adjusted hazard ratio of 0.26 (95% CI 0.16-0.41). The median of reproductive days was 55.50 in the high variability line whilst the median was 252.50 days for the low variability females. The homogeneous line presented important reproductive advantages suggesting higher robustness and animal welfare. Further research should confirm whether the findings presented here of a better performance in the low variability line could be properly applied to some livestock species.

Genetic Architecture of Maize Rind Strength Revealed by the Analysis of Divergently Selected Populations.

The strength of the stalk rind, measured as rind penetrometer resistance (RPR), is an important contributor to stalk lodging resistance. To enhance the genetic architecture of RPR, we combined selection mapping on populations developed by 15 cycles of divergent selection for high and low RPR with time-course transcriptomic and metabolic analyses of the stalks. Divergent selection significantly altered allele frequencies of 3,656 and 3,412 single- nucleotide polymorphisms (SNPs) in the high and low RPR populations, respectively. Surprisingly, only 110 (1.56%) SNPs under selection were common in both populations, while the majority (98.4%) were unique to each population. This result indicated that high and low RPR phenotypes are produced by biologically distinct mechanisms. Remarkably, regions harboring lignin and polysaccharide genes were preferentially selected in high and low RPR populations, respectively. The preferential selection was manifested as higher lignification and increased saccharification of the high and low RPR stalks, respectively. The evolution of distinct gene classes according to the direction of selection was unexpected in the context of parallel evolution and demonstrated that selection for a trait, albeit in different directions, does not necessarily act on the same genes. Tricin, a grass-specific monolignol that initiates the incorporation of lignin in the cell walls, emerged as a key determinant of RPR. Integration of selection mapping and transcriptomic analyses with published genetic studies of RPR identified several candidate genes including ZmMYB31, ZmNAC25, ZmMADS1, ZmEXPA2, ZmIAA41 and hk5. These findings provide a foundation for an enhanced understanding of RPR and the improvement of stalk lodging resistance.

Genomic signatures of spatially divergent selection at clownfish range margins.

Understanding how evolutionary forces interact to drive patterns of selection and distribute genetic variation across a species' range is of great interest in ecology and evolution, especially in an era of global change. While theory predicts how and when populations at range margins are likely to undergo local adaptation, empirical evidence testing these models remains sparse. Here, we address this knowledge gap by investigating the relationship between selection, gene flow and genetic drift in the yellowtail clownfish, Amphiprion clarkii, from the core to the northern periphery of the species range. Analyses reveal low genetic diversity at the range edge, gene flow from the core to the edge and genomic signatures of local adaptation at 56 single nucleotide polymorphisms in 25 candidate genes, most of which are significantly correlated with minimum annual sea surface temperature. Several of these candidate genes play a role in functions that are upregulated during cold stress, including protein turnover, metabolism and translation. Our results illustrate how spatially divergent selection spanning the range core to the periphery can occur despite the potential for strong genetic drift at the range edge and moderate gene flow from the core populations.

Population genomics reveals repeated signals of adaptive divergence in the Atlantic salmon of north-eastern Europe.

Our ability to examine genetic variation across entire genomes has enabled many studies searching for the genetic basis of local adaptation. These studies have identified numerous loci as candidates for differential local selection; however, relatively few have examined the overlap among candidate loci identified from independent studies of the same species in different geographic areas or evolutionary lineages. We used an allelotyping approach with a 220K SNP array to characterize the population genetic structure of Atlantic salmon in north-eastern Europe and ask whether the same genomic segments emerged as outliers among populations in different geographic regions. Genome-wide data recapitulated the phylogeographic structure previously inferred from mtDNA and microsatellite markers. Independent analyses of three genetically and geographically distinct groups of populations repeatedly inferred the same 17 haploblocks to contain loci under differential local selection. The most strongly supported of these replicated haploblocks had known strong associations with life-history variation or immune response in Atlantic salmon. Our results are consistent with these genomic segments harbouring large-effect loci which have a major role in Atlantic salmon diversification and are ideal targets for validation studies.

Parallel evolution of morphological traits and body shape in littoral and pelagic brook charr, Salvelinus fontinalis, along a gradient of interspecific competition.

The parallel evolution of similar ecotypes in response to comparable environmental conditions is believed to reveal the importance of divergent selection in phenotypic diversifying processes. Systems characterized by the presence of multiple replicate populations expressing resource polymorphism thus provide an ideal opportunity to address the occurrence and factors affecting the parallel evolution of ecotypes. Previous studies have shown that brook charr (Salvelinus fontinalis) exhibit resource polymorphism in some Canadian Shield lakes, where a littoral ecotype feeds mainly on zoobenthos and a pelagic ecotype feeds mostly on zooplankton. Using morphological traits and geometric morphometric analyses on 18 native brook charr populations, we explicitly tested (i) whether brook charr ecotypes show parallel evolution across populations (i.e. the same morphological traits discriminate ecotypes among lakes) and (ii) whether interspecific competition decreases the amplitude of morphological differentiation between ecotypes, if any, because brook charr experience some level of competitive exclusion from the littoral habitat in the presence of creek chub or white sucker. We observed a low level of parallel evolution, where the littoral ecotype was overall stouter with longer fins and smaller eyes than the pelagic ecotype. Interspecific competition had no clear impacts on the amplitude of morphological differentiation. We also observed that inter-lake morphological differences are greater than between ecotypes within lakes, suggesting an important effect of local environmental factors on population morphology. Early-stage of diversification as well as phenotypic plasticity and morphological integration could explain why resource polymorphism is still subtle in brook charr populations.

Combining population genomics and fitness QTLs to identify the genetics of local adaptation in Arabidopsis thaliana.

Evidence for adaptation to different climates in the model species Arabidopsis thaliana is seen in reciprocal transplant experiments, but the genetic basis of this adaptation remains poorly understood. Field-based quantitative trait locus (QTL) studies provide direct but low-resolution evidence for the genetic basis of local adaptation. Using high-resolution population genomic approaches, we examine local adaptation along previously identified genetic trade-off (GT) and conditionally neutral (CN) QTLs for fitness between locally adapted Italian and Swedish A. thaliana populations [Ågren J, et al. (2013) Proc Natl Acad Sci USA 110:21077-21082]. We find that genomic regions enriched in high FST SNPs colocalize with GT QTL peaks. Many of these high FST regions also colocalize with regions enriched for SNPs significantly correlated to climate in Eurasia and evidence of recent selective sweeps in Sweden. Examining unfolded site frequency spectra across genes containing high FST SNPs suggests GTs may be due to more recent adaptation in Sweden than Italy. Finally, we collapse a list of thousands of genes spanning GT QTLs to 42 genes that likely underlie the observed GTs and explore potential biological processes driving these trade-offs, from protein phosphorylation, to seed dormancy and longevity. Our analyses link population genomic analyses and field-based QTL studies of local adaptation, and emphasize that GTs play an important role in the process of local adaptation.

Positive selection alone is sufficient for whole genome differentiation at the early stage of speciation process in the fall armyworm.

BACKGROUND: The process of speciation involves differentiation of whole genome sequences between a pair of diverging taxa. In the absence of a geographic barrier and in the presence of gene flow, genomic differentiation may occur when the homogenizing effect of recombination is overcome across the whole genome. The fall armyworm is observed as two sympatric strains with different host-plant preferences across the entire habitat. These two strains exhibit a very low level of genetic differentiation across the whole genome, suggesting that genomic differentiation occurred at an early stage of speciation. In this study, we aim at identifying critical evolutionary forces responsible for genomic differentiation in the fall armyworm. RESULTS: These two strains exhibit a low level of genomic differentiation (FST = 0.0174), while 99.2% of 200 kb windows have genetically differentiated sequences (FST > 0). We found that the combined effect of mild positive selection and genetic linkage to selectively targeted loci are responsible for the genomic differentiation. However, a single event of very strong positive selection appears not to be responsible for genomic differentiation. The contribution of chromosomal inversions or tight genetic linkage among positively selected loci causing reproductive barriers is not supported by our data. Phylogenetic analysis shows that the genomic differentiation occurred by sub-setting of genetic variants in one strain from the other. CONCLUSIONS: From these results, we concluded that genomic differentiation may occur at the early stage of a speciation process in the fall armyworm and that mild positive selection targeting many loci alone is sufficient evolutionary force for generating the pattern of genomic differentiation. This genomic differentiation may provide a condition for accelerated genomic differentiation by synergistic effects among linkage disequilibrium generated by following events of positive selection. Our study highlights genomic differentiation as a key evolutionary factor connecting positive selection to divergent selection.