Finding needles in a genomic haystack: targeted capture identifies clear signatures of selection in a nonmodel plant species

Teasing apart neutral and adaptive genomic processes and identifying loci that are targets of selection can be difficult, particularly for nonmodel species that lack a reference genome. However, identifying such loci and the factors driving selection have the potential to greatly assist conservation and restoration practices, especially for the management of species in the face of contemporary and future climate change. Here, we focus on assessing adaptive genomic variation within a nonmodel plant species, the narrow‐leaf hopbush (Dodonaea viscosa ssp. angustissima), commonly used for restoration in Australia. We used a hybrid‐capture target enrichment approach to selectively sequence 970 genes across 17 populations along a latitudinal gradient from 30°S to 36°S. We analysed 8462 single‐nucleotide polymorphisms (SNPs) for F_ST outliers as well as associations with environmental variables. Using three different methods, we found 55 SNPs with significant correlations to temperature and water availability, and 38 SNPs to elevation. Genes containing SNPs identified as under environmental selection were diverse, including aquaporin and abscisic acid genes, as well as genes with ontologies relating to responses to environmental stressors such as water deprivation and salt stress. Redundancy analysis demonstrated that only a small proportion of the total genetic variance was explained by environmental variables. We demonstrate that selection has led to clines in allele frequencies in a number of functional genes, including those linked to leaf shape and stomatal variation, which have been previously observed to vary along the sampled environmental cline. Using our approach, gene regions subject to environmental selection can be readily identified for nonmodel organisms.     

Whole‐genome resequencing uncovers molecular signatures of natural and sexual selection in wild bighorn sheep

The identification of genes influencing fitness is central to our understanding of the genetic basis of adaptation and how it shapes phenotypic variation in wild populations. Here, we used whole‐genome resequencing of wild Rocky Mountain bighorn sheep (Ovis canadensis) to >50‐fold coverage to identify 2.8 million single nucleotide polymorphisms (SNPs) and genomic regions bearing signatures of directional selection (i.e. selective sweeps). A comparison of SNP diversity between the X chromosome and the autosomes indicated that bighorn males had a dramatically reduced long‐term effective population size compared to females. This probably reflects a long history of intense sexual selection mediated by male–male competition for mates. Selective sweep scans based on heterozygosity and nucleotide diversity revealed evidence for a selective sweep shared across multiple populations at RXFP2, a gene that strongly affects horn size in domestic ungulates. The massive horns carried by bighorn rams appear to have evolved in part via strong positive selection at RXFP2. We identified evidence for selection within individual populations at genes affecting early body growth and cellular response to hypoxia; however, these must be interpreted more cautiously as genetic drift is strong within local populations and may have caused false positives. These results represent a rare example of strong genomic signatures of selection identified at genes with known function in wild populations of a nonmodel species. Our results also showcase the value of reference genome assemblies from agricultural or model species for studies of the genomic basis of adaptation in closely related wild taxa.     

Detecting selection along environmental gradients: analysis of eight methods and their effectiveness for outbreeding and selfing populations

Thanks to genome‐scale diversity data, present‐day studies can provide a detailed view of how natural and cultivated species adapt to their environment and particularly to environmental gradients. However, due to their sensitivity, up‐to‐date studies might be more sensitive to undocumented demographic effects such as the pattern of migration and the reproduction regime. In this study, we provide guidelines for the use of popular or recently developed statistical methods to detect footprints of selection. We simulated 100 populations along a selective gradient and explored different migration models, sampling schemes and rates of self‐fertilization. We investigated the power and robustness of eight methods to detect loci potentially under selection: three designed to detect genotype–environment correlations and five designed to detect adaptive differentiation (based on F_ST or similar measures). We show that genotype–environment correlation methods have substantially more power to detect selection than differentiation‐based methods but that they generally suffer from high rates of false positives. This effect is exacerbated whenever allele frequencies are correlated, either between populations or within populations. Our results suggest that, when the underlying genetic structure of the data is unknown, a number of robust methods are preferable. Moreover, in the simulated scenario we used, sampling many populations led to better results than sampling many individuals per population. Finally, care should be taken when using methods to identify genotype–environment correlations without correcting for allele frequency autocorrelation because of the risk of spurious signals due to allele frequency correlations between populations.     

Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds

A dataset consisting of 787 animals with high‐density SNP chip genotypes (346 774 SNPs) and 939 animals with medium‐density SNP chip genotypes (33 828 SNPs) from eight indigenous Swiss sheep breeds was analyzed to characterize population structure, quantify genomic inbreeding based on runs of homozygosity and identify selection signatures. In concordance with the recent known history of these breeds, the highest genetic diversity was observed in Engadine Red sheep and the lowest in Valais Blacknose sheep. Correlation between F_PED and F_ROH was around 0.50 and thereby lower than that found in similar studies in cattle. Mean F_ROH estimates from medium‐density data and HD data were highly correlated (0.95). Signatures of selection and candidate gene analysis revealed that the most prominent signatures of selection were found in the proximity of genes associated with body size (NCAPG, LCORL, LAP3, SPP1, PLAG1, ALOX12, TP53), litter size (SPP1), milk production (ABCG2, SPP1), coat color (KIT, ASIP, TBX3) and horn status (RXFP2). For the Valais Blacknose sheep, the private signatures in proximity of genes/QTL influencing body size, coat color and fatty acid composition were confirmed based on runs of homozygosity analysis. These private signatures underline the genetic uniqueness of the Valais Blacknose sheep breed. In conclusion, we identified differences in the genetic make‐up of Swiss sheep breeds, and we present relevant candidate genes responsible for breed differentiation in locally adapted breeds.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

Exome resequencing reveals signatures of demographic and adaptive processes across the genome and range of black cottonwood (Populus trichocarpa)

Extant variation in temperate and boreal plant species has been influenced by both demographic histories associated with Pleistocene glacial cycles and adaptation to local climate. We used sequence capture to investigate the role of these neutral and adaptive processes in shaping diversity in black cottonwood (Populus trichocarpa). Nucleotide diversity and Tajima's D were lowest at replacement sites and highest at intergenic sites, while LD showed the opposite pattern. With samples grouped into three populations arrayed latitudinally, effective population size was highest in the north, followed by south and centre, and LD was highest in the south followed by the north and centre, suggesting a possible northern glacial refuge. F_ST outlier analysis revealed that promoter, 5′‐UTR and intronic sites were enriched for outliers compared with coding regions, while no outliers were found among intergenic sites. Codon usage bias was evident, and genes with synonymous outliers had 30% higher average expression compared with genes containing replacement outliers. These results suggest divergent selection related to regulation of gene expression is important to local adaptation in P. trichocarpa. Finally, within‐population selective sweeps were much more pronounced in the central population than in putative northern and southern refugia, which may reflect the different demographic histories of the populations and concomitant effects on signatures of genetic hitchhiking from standing variation.     

A clever solution to a vexing problem

F_ST (as well as related measures such as G_ST) has long been used both as a measure of the relative amount of genetic variation between populations and as an indicator of the amount of gene flow among populations. Unfortunately, F_ST and its clones are also sensitive to mutation, particularly when the mutation rate per locus is greater than the migration rate among populations. Relatively high mutation rates cause estimates of F_ST and G_ST to be much lower than researchers sometimes expect, when migration rates are low in the studied species. Several recent suggestions for dealing with this problem have been unsatisfactory for one reason or another, and no general solution exists (if we are not to abandon these otherwise useful measures of differentiation). In an important article in this issue, Jinliang Wang (2015) shows that it is possible to identify whether the genetic markers in a given study are likely to give estimates of F_ST that are strongly affected by mutation. The proposed test is simple and elegant, and with it, molecular ecologists can determine whether the F_ST from their makers can be depended on for further inference about their species’ genome and the demographic forces which shaped its patterns.     

QST–FST comparisons with unbalanced half‐sib designs

Q_ST, a measure of quantitative genetic differentiation among populations, is an index that can suggest local adaptation if Q_ST for a trait is sufficiently larger than the mean F_ST of neutral genetic markers. A previous method by Whitlock and Guillaume derived a simulation resampling approach to statistically test for a difference between Q_ST and F_ST, but that method is limited to balanced data sets with offspring related as half‐sibs through shared fathers. We extend this approach (i) to allow for a model more suitable for some plant populations or breeding designs in which offspring are related through mothers (assuming independent fathers for each offspring; half‐sibs by dam); and (ii) by explicitly allowing for unbalanced data sets. The resulting approach is made available through the R package QstFstComp.     

Local adaptation through genetic differentiation in highly fragmented Tilia cordata populations

We assessed the level of geographic differentiation of Tilia cordata in Denmark based on tests of 91 trees selected from 12 isolated populations. We used quantitative analysis of spring phenology and population genetic analysis based on SSR markers to infer the likely historical genetic processes within and among populations. High genetic variation within and among populations was observed in spring phenology, which correlated with spring temperatures at the origin of the tested T. cordata trees. The population genetic analysis revealed significant differentiation among the populations, but with no clear sign of isolation by distance. We infer the findings as indications of ongoing fine scale selection in favor of local growth conditions made possible by limited gene flow among the small and fragmented populations. This hypothesis fits well with reports of limited fruiting in the investigated Danish T. cordata populations, while the species is known for its ability to propagate vegetatively by root suckers. Our results suggest that both divergent selection and genetic drift may have played important roles in forming the genetic patterns of T. cordata at its northern distribution limit. However, we also speculate that epigenetic mechanism arising from the original population environment could have created similar patterns in regulating the spring phenology.     

Genetic diversity and genomic signatures of selection among cattle breeds from Siberia,eastern and northern Europe

Domestication in the near eastern region had a major impact on the gene pool of humpless taurine cattle (Bos taurus). As a result of subsequent natural and artificial selection, hundreds of different breeds have evolved, displaying a broad range of phenotypic traits. Here, 10 Eurasian B. taurus breeds from different biogeographic and production conditions, which exhibit different demographic histories and have been under artificial selection at various intensities, were investigated using the Illumina BovineSNP50 panel to understand their genetic diversity and population structure. In addition, we scanned genomes from eight breeds for signatures of diversifying selection. Our population structure analysis indicated six distinct breed groups, the most divergent being the Yakutian cattle from Siberia. Selection signals were shared (experimental P‐value < 0.01) with more than four breeds on chromosomes 6, 7, 13, 16 and 22. The strongest selection signals in the Yakutian cattle were found on chromosomes 7 and 21, where a miRNA gene and genes related to immune system processes are respectively located. In general, genomic regions indicating selection overlapped with known QTL associated with milk production (e.g. on chromosome 19), reproduction (e.g. on chromosome 24) and meat quality (e.g. on chromosome 7). The selection map created in this study shows that native cattle breeds and their genetic resources represent unique material for future breeding.     

A population genomic scan in Chorthippus grasshoppers unveils previously unknown phenotypic divergence

Understanding the genetics of speciation and the processes that drive it is a central goal of evolutionary biology. Grasshoppers of the Chorthippus species group differ strongly in calling song (and corresponding female preferences) but are exceedingly similar in other characteristics such as morphology. Here, we performed a population genomic scan on three Chorthippus species (Chorthippus biguttulus, C. mollis and C. brunneus) to gain insight into the genes and processes involved in divergence and speciation in this group. Using an RNA‐seq approach, we examined functional variation between the species by calling SNPs for each of the three species pairs and using F_ST‐based approaches to identify outliers. We found approximately 1% of SNPs in each comparison to be outliers. Between 37% and 40% of these outliers were nonsynonymous SNPs (as opposed to a global level of 17%) indicating that we recovered loci under selection. Among the outliers were several genes that may be involved in song production and hearing as well as genes involved in other traits such as food preferences and metabolism. Differences in food preferences between species were confirmed with a behavioural experiment. This indicates that multiple phenotypic differences implicating multiple evolutionary processes (sexual selection and natural selection) are present between the species.     

Genome‐wide detection of signatures of selection in Korean Hanwoo cattle

The Korean Hanwoo cattle have been intensively selected for production traits, especially high intramuscular fat content. It is believed that ancient crossings between different breeds contributed to forming the Hanwoo, but little is known about the genomic differences and similarities between other cattle breeds and the Hanwoo. In this work, cattle breeds were grouped by origin into four types and used for comparisons: the Europeans (represented by six breeds), zebu (Nelore), African taurine (N'Dama) and Hanwoo. All animals had genotypes for around 680 000 SNPs after quality control of genotypes. Average heterozygosity was lower in Nelore and N'Dama (0.22 and 0.21 respectively) than in Europeans (0.26–0.31, with Shorthorn as outlier at 0.24) and Hanwoo (0.29). Pairwise F_ST analyses demonstrated that Hanwoo are more related to European cattle than to Nelore, with N'Dama in an intermediate position. This finding was corroborated by principal components and unsupervised hierarchical clustering. Using genome‐wide smoothed F_ST, 55 genomic regions potentially under positive selection in Hanwoo were identified. Among these, 29 were regions also detected in previous studies. Twenty‐four regions were exclusive to Hanwoo, and a number of other regions were shared with one or two of the other groups. These regions overlap a number of genes that are related to immune, reproduction and fatty acid metabolism pathways. Further analyses are needed to better characterize the ancestry of the Hanwoo cattle and to define the genes responsible to the identified selection peaks.     

Bud phenology and growth are subject to divergent selection across a latitudinal gradient in Populus angustifolia and impact adaptation across the distributional range and associated arthropods

Temperate forest tree species that span large geographical areas and climatic gradients often have high levels of genetic variation. Such species are ideal for testing how neutral demographic factors and climate‐driven selection structure genetic variation within species, and how this genetic variation can affect ecological communities. Here, we quantified genetic variation in vegetative phenology and growth traits in narrowleaf cottonwood, Populus angustifolia, using three common gardens planted with genotypes originating from source populations spanning the species' range along the Rocky Mountains of North America (ca. 1700 km). We present three main findings. First, we found strong evidence of divergent selection (Q_ST > F_ST) on fall phenology (bud set) with adaptive consequences for frost avoidance. We also found evidence for selection on bud flush duration, tree height, and basal diameter, resulting in population differentiation. Second, we found strong associations with climate variables that were strongly correlated with latitude of origin. More strongly differentiated traits also showed stronger climate correlations, which emphasizes the role that climate has played in divergent selection throughout the range. We found population × garden interaction effects; for some traits, this accounted for more of the variance than either factor alone. Tree height was influenced by the difference in climate of the source and garden locations and declined with increasing transfer distance. Third, growth traits were correlated with dependent arthropod community diversity metrics. Synthesis. Overall, we conclude that climate has influenced genetic variation and structure in phenology and growth traits and leads to local adaptation in P. angustifolia, which can then impact dependent arthropod species. Importantly, relocation of genotypes far northward or southward often resulted in poor growth, likely due to a phenological mismatch with photoperiod, the proximate cue for fall growth cessation. Genotypes moved too far southward suffer from early growth cessation, whereas those moved too far northward are prone to fall frost and winter dieback. In the face of current and forecasted climate change, habitat restoration, forestry, and tree breeding efforts should utilize these findings to better match latitudinal and climatic source environments with management locations for optimal future outcomes.     

Detecting strong positive selection in the genome

New statistical tests have been developed in the past decade that enable us to infer evidence of recent strong positive selection from genome-wide data on single-nucleotide polymorphism and to localize the targets of selection in the genome. Based on these tests, past demographic events that led to distortions of the site-frequency spectrum of variation can be distinguished from selection, in particular if linkage disequilibrium is taken into account. These methods have been successfully applied to species from which complete sequence information and polymorphism data are available, including Drosophila melanogaster, humans, and several plant species. To make full use of the available data, however, the tests that were primarily designed for panmictic populations need to be extended to spatially structured populations.     

Reanalysis suggests that genomic islands of speciation are due to reduced diversity,not reduced gene flow

The metaphor of ‘genomic islands of speciation’ was first used to describe heterogeneous differentiation among loci between the genomes of closely related species. The biological model proposed to explain these differences was that the regions showing high levels of differentiation were resistant to gene flow between species, while the remainder of the genome was being homogenized by gene flow and consequently showed lower levels of differentiation. However, the conditions under which such differentiation can occur at multiple unlinked loci are restrictive; additionally, essentially, all previous analyses have been carried out using relative measures of divergence, which can be misleading when regions with different levels of recombination are compared. Here, we test the model of differential gene flow by asking whether absolute divergence is also higher in the previously identified ‘islands’. Using five species pairs for which full sequence data are available, we find that absolute measures of divergence are not higher in genomic islands. Instead, in all cases examined, we find reduced diversity in these regions, a consequence of which is that relative measures of divergence are abnormally high. These data therefore do not support a model of differential gene flow among loci, although islands of relative divergence may represent loci involved in local adaptation. Simulations using the program IMa2 further suggest that inferences of any gene flow may be incorrect in many comparisons. We instead present an alternative explanation for heterogeneous patterns of differentiation, one in which postspeciation selection generates patterns consistent with multiple aspects of the data.     

Large‐scale adaptive divergence in Boechera fecunda,an endangered wild relative of Arabidopsis

Many biological species are threatened with extinction because of a number of factors such as climate change and habitat loss, and their preservation depends on an accurate understanding of the extent of their genetic variability within and among populations. In this study, we assessed the genetic divergence of five quantitative traits in 10 populations of an endangered cruciferous species, Boechera fecunda, found in only several populations in each of two geographic regions (WEST and EAST) in southwestern Montana. We analyzed variation in quantitative traits, neutral molecular markers, and environmental factors and provided evidence that despite the restricted geographical distribution of this species, it exhibits a high level of genetic variation and regional adaptation. Conservation efforts therefore should be directed to the preservation of populations in each of these two regions without attempting transplantation between regions. Heritabilities and genetic coefficients of variation estimated from nested ANOVAs were generally high for leaf and rosette traits, although lower (and not significantly different from 0) for water‐use efficiency. Measures of quantitative genetic differentiation, Q_ST, were calculated for each trait from each pair of populations. For three of the five traits, these values were significantly higher between regions compared with those within regions (after adjustment for neutral genetic variation, F_ST). This suggested that natural selection has played an important role in producing regional divergence in this species. Our analysis also revealed that the B. fecunda populations appear to be locally adapted due, at least in part, to differences in environmental conditions in the EAST and WEST regions.     

Analysis of ROH patterns in the Noriker horse breed reveals signatures of selection for coat color and body size

Overlapping runs of homozygosity (ROH islands) shared by the majority of a population are hypothesized to be the result of selection around a target locus. In this study we investigated the impact of selection for coat color within the Noriker horse on autozygosity and ROH patterns. We analyzed overlapping homozygous regions (ROH islands) for gene content in fragments shared by more than 50% of horses. Long‐term assortative mating of chestnut horses and the small effective population size of leopard spotted and tobiano horses resulted in higher mean genome‐wide ROH coverage (S_ROH) within the range of 237.4–284.2 Mb, whereas for bay, black and roan horses, where rotation mating is commonly applied, lower autozygosity (S_ROH from 176.4–180.0 Mb) was determined. We identified seven common ROH islands considering all Noriker horses from our dataset. Specific islands were documented for chestnut, leopard spotted, roan and bay horses. The ROH islands contained, among others, genes associated with body size (ZFAT, LASP1 and LCORL/NCAPG), coat color (MC1R in chestnut and the factor PATN1 in leopard spotted horses) and morphogenesis (HOXB cluster in all color strains except leopard spotted horses). This study demonstrates that within a closed population sharing the same founders and ancestors, selection on a single phenotypic trait, in this case coat color, can result in genetic fragmentation affecting levels of autozygosity and distribution of ROH islands and enclosed gene content.     

Genetic and evolution analysis of extrafloral nectary in cotton

Extrafloral nectaries are a defence trait that plays important roles in plant–animal interactions. Gossypium species are characterized by cellular grooves in leaf midribs that secret large amounts of nectar. Here, with a panel of 215 G. arboreum accessions, we compared extrafloral nectaries to nectariless accessions to identify a region of Chr12 that showed strong differentiation and overlapped with signals from GWAS of nectaries. Fine mapping of an F₂ population identified GaNEC1, encoding a PB1 domain‐containing protein, as a positive regulator of nectary formation. An InDel, encoding a five amino acid deletion, together with a nonsynonymous substitution, was predicted to cause 3D structural changes in GaNEC1 protein that could confer the nectariless phenotype. mRNA‐Seq analysis showed that JA‐related genes are up‐regulated and cell wall‐related genes are down‐regulated in the nectary. Silencing of GaNEC1 led to a smaller size of foliar nectary phenotype. Metabolomics analysis identified more than 400 metabolites in nectar, including expected saccharides and amino acids. The identification of GaNEC1 helps establish the network regulating nectary formation and nectar secretion, and has implications for understanding the production of secondary metabolites in nectar. Our results will deepen our understanding of plant–mutualism co‐evolution and interactions, and will enable utilization of a plant defence trait in cotton breeding efforts.     

Evolutionary dynamics of the leaf phenological cycle in an oak metapopulation along an elevation gradient

It has been predicted that environmental changes will radically alter the selective pressures on phenological traits. Long‐lived species, such as trees, will be particularly affected, as they may need to undergo major adaptive change over only one or a few generations. The traits describing the annual life cycle of trees are generally highly evolvable, but nothing is known about the strength of their genetic correlations. Tight correlations can impose strong evolutionary constraints, potentially hampering the adaptation of multivariate phenological phenotypes. In this study, we investigated the evolutionary, genetic and environmental components of the timing of leaf unfolding and senescence within an oak metapopulation along an elevation gradient. Population divergence, estimated from in situ and common‐garden data, was compared to expectations under neutral evolution, based on microsatellite markers. This approach made it possible (1) to evaluate the influence of genetic correlation on multivariate local adaptation to elevation and (2) to identify traits probably exposed to past selective pressures due to the colder climate at high elevation. The genetic correlation was positive but very weak, indicating that genetic constraints did not shape the local adaptation pattern for leaf phenology. Both spring and fall (leaf unfolding and senescence, respectively) phenology timings were involved in local adaptation, but leaf unfolding was probably the trait most exposed to climate change‐induced selection. Our data indicated that genetic variation makes a much smaller contribution to adaptation than the considerable plastic variation displayed by a tree during its lifetime. The evolutionary potential of leaf phenology is, therefore, probably not the most critical aspect for short‐term population survival in a changing climate.     

Detection of selection signatures of population‐specific genomic regions selected during domestication process in Jinhua pigs

Chinese pigs have been undergoing both natural and artificial selection for thousands of years. Jinhua pigs are of great importance, as they can be a valuable model for exploring the genetic mechanisms linked to meat quality and other traits such as disease resistance, reproduction and production. The purpose of this study was to identify distinctive footprints of selection between Jinhua pigs and other breeds utilizing genome‐wide SNP data. Genotyping by genome reducing and sequencing was implemented in order to perform cross‐population extended haplotype homozygosity to reveal strong signatures of selection for those economically important traits. This work was performed at a 2% genome level, which comprised 152 006 SNPs genotyped in a total of 517 individuals. Population‐specific footprints of selective sweeps were searched for in the genome of Jinhua pigs using six native breeds and three European breeds as reference groups. Several candidate genes associated with meat quality, health and reproduction, such as GH1, CRHR2, TRAF4 and CCK, were found to be overlapping with the significantly positive outliers. Additionally, the results revealed that some genomic regions associated with meat quality, immune response and reproduction in Jinhua pigs have evolved directionally under domestication and subsequent selections. The identified genes and biological pathways in Jinhua pigs showed different selection patterns in comparison with the Chinese and European breeds.