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.     

Genetic association of marbling score with intragenic nucleotide variants at selection signals of the bovine genome

Selection signals of Korean cattle might be attributed largely to artificial selection for meat quality. Rapidly increased intragenic markers of newly annotated genes in the bovine genome would help overcome limited findings of genetic markers associated with meat quality at the selection signals in a previous study. The present study examined genetic associations of marbling score (MS) with intragenic nucleotide variants at selection signals of Korean cattle. A total of 39 092 nucleotide variants of 407 Korean cattle were utilized in the association analysis. A total of 129 variants were selected within newly annotated genes in the bovine genome. Their genetic associations were analyzed using the mixed model with random polygenic effects based on identical-by-state genetic relationships among animals in order to control for spurious associations produced by population structure. Genetic associations of MS were found (P<3.88×10⁻⁴) with six intragenic nucleotide variants on bovine autosomes 3 (cache domain containing 1, CACHD1), 5 (like-glycosyltransferase, LARGE), 16 (cell division cycle 42 binding protein kinase alpha, CDC42BPA) and 21 (snurportin 1, SNUPN; protein tyrosine phosphatase, non-receptor type 9, PTPN9; chondroitin sulfate proteoglycan 4, CSPG4). In particular, the genetic associations with CDC42BPA and LARGE were confirmed using an independent data set of Korean cattle. The results implied that allele frequencies of functional variants and their proximity variants have been augmented by directional selection for greater MS and remain selection signals in the bovine genome. Further studies of fine mapping would be useful to incorporate favorable alleles in marker-assisted selection for MS of Korean cattle.     

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.     

Identification of positive selection signatures in pigs by comparing linkage disequilibrium variances

Selection affects the patterns of linkage disequilibrium (LD) around the site of a beneficial allele with an increase in LD among the hitchhiking alleles. Comparing the differences in regional LD between pig populations could help to identify putative genomic regions with potential adaptations for economic traits. In this study, using Illumina Porcine SNP60K BeadChip genotyping data from 207 Chinese indigenous, 117 South American village and 408 Large White pigs, we estimated the variation of genome‐wide LD between populations using the varld program. The top 0.1% standardized VarLD scores were used as a criterion for all comparisons, and compared with LD blocks, a total of four selection signatures on Sus scrofa chromosome (SSC) 7, 9, 13 and 14 were identified in all populations. These signatures overlapped with quantitative trait loci for linoleic acid content, age at puberty, number of muscle fibers per unit area, hip structure and body weight traits in pigs. Among them, one of the signatures (56.5–56.6 Mb on SSC7) in Large White pigs harbored the ADAMTSL3 gene, which is known to affect body length. The findings of this study seem to point toward recent selection in different pig populations. Further investigations are encouraged to confirm the selection signatures detected by varld in the present study.     

Assessing footprints of selection in commercial Atlantic salmon populations using microsatellite data

Relatively large rates of response to traits of economic importance have been observed in different selection experiments in salmon. Several QTL have been mapped in the salmon genome, explaining unprecedented levels of phenotypic variation. Owing to the relatively large selection intensity, individual loci may be indirectly selected, leaving molecular footprints of selection, together with increased inbreeding, as its likely relatives will share the selected loci. We used population differentiation and levels of linkage disequilibrium in chromosomes known to be harbouring QTL for body weight, infectious pancreatic necrosis resistance and infectious salmon anaemia resistance to assess the recent selection history at the genomic level in Atlantic salmon. The results clearly suggest that the marker SSA0343BSFU on chromosome 3 (body weight QTL) showed strong evidence of directional selection. It is intriguing that this marker is physically mapped to a region near the coding sequence of DVL2 , making it an ideal candidate gene to explain the rapid evolutionary response of this chromosome to selection for growth in Salmo salar. Weak evidence of diversifying selection was observed in the QTL associated with infectious pancreatic necrosis and infectious salmon anaemia resistance. Overall, this study showed that artificial selection has produced important changes in the Atlantic salmon genome, validating QTL in commercial salmon populations used for production purposes according to the recent selection history.     

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.     

Detection of signatures of selective sweeps in the Blonde d'Aquitaine cattle breed

Identifying recent positive selection signatures in domesticated animals could provide information on genome response to strong directional selection from domestication and artificial selection and therefore could help in identifying mutations responsible for improved traits. We used genotyping data generated using Illumina's BovineSNP50 Genotyping BeadChips to identify selection signatures in the Blonde d'Aquitaine breed, a well‐muscled French beef breed. For this purpose, we employed a hidden Markov model‐based test, which detects selection by studying local variations in the allele frequency spectrum along the genome, within a single population. Three regions containing selective sweeps were identified. Annotation of genes located within these regions revealed interesting candidate genes. For example, myostatin (also known as GDF8), a known muscle growth factor inhibitor, is located within the selection signature region found on chromosome 2. In addition, we have identified chromosomal regions that show some evidence of selection within QTL regions for economically important traits. The results of this study could help to better understand the mechanisms related to the selection of the Blonde d'Aquitaine breed.     

Identification of genetic associations of ECHS1 gene with milk fatty acid traits in dairy cattle

Our previous genome‐wide association study identified 83 genome‐wide significant SNPs and 20 novel promising candidate genes for milk fatty acids in Chinese Holstein. Among them, the enoyl‐CoA hydratase, short chain 1 (ECHS1) and enoyl‐CoA hydratase and 3‐hydroxyacyl CoA dehydrogenase (EHHADH) genes were located near two SNPs and one SNP respectively, and they play important roles in fatty acid metabolism pathways. We herein validated whether the two genes have genetic effects on milk fatty acid traits in dairy cattle. By re‐sequencing the full‐length coding region, partially adjacent introns and 3000 bp up/downstream flanking sequences, we identified 12 SNPs in ECHS1: two in exons, four in the 3′ flanking region and six in introns. The g.25858322C>T SNP results in an amino acid replacement from leucine to phenylalanine and changes the secondary structure of the ECHS1 protein, and single‐locus association analysis showed that it was significantly associated with three milk fatty acids (P = 0.0002–0.0013). The remaining 11 SNPs were found to be significantly associated with at least one milk fatty acid (P = <0.0001–0.0040). Also, we found that two haplotype blocks, consisting of nine and two SNPs respectively, were significantly associated with eight milk fatty acids (P = <0.0001–0.0125). However, none of polymorphisms was observed in the EHHADH gene. In conclusion, our findings are the first to indicate that the ECHS1 gene has a significant genetic impact on long‐chain unsaturated and medium‐chain saturated fatty acid traits in dairy cattle, although the biological mechanism is still undetermined and requires further in‐depth validation.     

No selection for change in polyandry under experimental evolution

What drives mating system variation is a major question in evolutionary biology. Female multiple mating (polyandry) has diverse evolutionary consequences, and there are many potential benefits and costs of polyandry. However, our understanding of its evolution is biased towards studies enforcing monandry in polyandrous species. What drives and maintains variation in polyandry between individuals, genotypes, populations and species remains poorly understood. Genetic variation in polyandry may be actively maintained by selection, or arise by chance if polyandry is selectively neutral. In Drosophila pseudoobscura, there is genetic variation in polyandry between and within populations. We used isofemale lines to found replicate populations with high or low initial levels of polyandry and tracked polyandry under experimental evolution over seven generations. Polyandry remained relatively stable, reflecting the starting frequencies of the experimental populations. There were no clear fitness differences between high versus low polyandry genotypes, and there was no signature of balancing selection. We confirmed these patterns in direct comparisons between evolved and ancestral females and found no consequences of polyandry for female fecundity. The absence of differential selection even when initiating populations with major differences in polyandry casts some doubt on the importance of polyandry for female fitness.     

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.     

Microevolution in time and space: SNP analysis of historical DNA reveals dynamic signatures of selection in Atlantic cod

Little is known about how quickly natural populations adapt to changes in their environment and how temporal and spatial variation in selection pressures interact to shape patterns of genetic diversity. We here address these issues with a series of genome scans in four overfished populations of Atlantic cod (Gadus morhua) studied over an 80‐year period. Screening of >1000 gene‐associated single‐nucleotide polymorphisms (SNPs) identified 77 loci that showed highly elevated levels of differentiation, likely as an effect of directional selection, in either time, space or both. Exploratory analysis suggested that temporal allele frequency shifts at certain loci may correlate with local temperature variation and with life history changes suggested to be fisheries induced. Interestingly, however, largely nonoverlapping sets of loci were temporal outliers in the different populations and outliers from the 1928 to 1960 period showed almost complete stability during later decades. The contrasting microevolutionary trajectories among populations resulted in sequential shifts in spatial outliers, with no locus maintaining elevated spatial differentiation throughout the study period. Simulations of migration coupled with observations of temporally stable spatial structure at neutral loci suggest that population replacement or gene flow alone could not explain all the observed allele frequency variation. Thus, the genetic changes are likely to at least partly be driven by highly dynamic temporally and spatially varying selection. These findings have important implications for our understanding of local adaptation and evolutionary potential in high gene flow organisms and underscore the need to carefully consider all dimensions of biocomplexity for evolutionarily sustainable management.     

Cunningham's skinks show low genetic connectivity and signatures of divergent selection across its distribution

Establishing corridors of connecting habitat has become a mainstay conservation strategy to maintain gene flow and facilitate climate‐driven range shifts. Yet, little attention has been given to ascertaining the extent to which corridors will benefit philopatric species, which might exhibit localized adaptation. Measures of genetic connectivity and adaptive genetic variation across species’ ranges can help fill this knowledge gap. Here, we characterized the spatial genetic structure of Cunningham's skink (Egernia cunninghami), a philopatric species distributed along Australia's Great Dividing Range, and assessed evidence of localized adaptation. Analysis of 4,274 SNPs from 94 individuals sampled at four localities spanning 500 km and 4° of latitude revealed strong genetic structuring at neutral loci (mean F_ST ± SD = 0.603 ± 0.237) among the localities. Putatively neutral SNPs and those under divergent selection yielded contrasting spatial patterns, with the latter identifying two genetically distinct clusters. Given low genetic connectivity of the four localities, we suggest that the natural movement rate of this species is insufficient to keep pace with spatial shifts to its climate envelope, irrespective of habitat availability. In addition, our finding of localized adaptation highlights the risk of outbreeding depression should the translocation of individuals be adopted as a conservation management strategy.     

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.     

Genetic signatures of natural selection in response to air pollution in red spruce (Picea rubens,Pinaceae)

One of the most important drivers of local adaptation for forest trees is climate. Coupled to these patterns, however, are human‐induced disturbances through habitat modification and pollution. The confounded effects of climate and disturbance have rarely been investigated with regard to selective pressure on forest trees. Here, we have developed and used a population genetic approach to search for signals of selection within a set of 36 candidate genes chosen for their putative effects on adaptation to climate and human‐induced air pollution within five populations of red spruce (Picea rubens Sarg.), distributed across its natural range and air pollution gradient in eastern North America. Specifically, we used F_ST outlier and environmental correlation analyses to highlight a set of seven single nucleotide polymorphisms (SNPs) that were overly correlated with climate and levels of sulphate pollution after correcting for the confounding effects of population history. Use of three age cohorts within each population allowed the effects of climate and pollution to be separated temporally, as climate‐related SNPs (n = 7) showed the strongest signals in the oldest cohort, while pollution‐related SNPs (n = 3) showed the strongest signals in the youngest cohorts. These results highlight the usefulness of population genetic scans for the identification of putatively nonneutral evolution within genomes of nonmodel forest tree species, but also highlight the need for the development and application of robust methodologies to deal with the inherent multivariate nature of the genetic and ecological data used in these types of analyses.     

Mito‐nuclear discordance with evidence of shared ancestral polymorphism and selection in cactophilic species of Drosophila

The Drosophila serido haplogroup is a monophyletic group composed of the following four cryptic and cactophilic species that are endemic to eastern Brazil: D. borborema, D. gouveai, D. seriema and D. serido. Here, we investigate the mito‐nuclear discordance in these species found among the cytochrome c oxidase subunit I (COI) mitochondrial gene, the autosomal alpha‐Esterase‐5 (α‐Est5) and the X‐linked period gene (per). Our analysis indicates that shared polymorphisms in these three molecular markers may be explained by the maintenance of ancestral polymorphisms rather than introgressive hybridization. The primary structures of COI, per and α‐Est5 genes evolve primarily under purifying selection, but we detected some sites that evolved under positive selection in α‐Est5. Considering the high variability of cacti species in eastern Brazil and the role attributed to Drosophila esterases in digestion metabolism and/or the detoxification of several compounds found in cactus tissues, we conjecture about the role of natural selection triggered by host shifts as an important factor in the intraspecific diversification of the D. serido haplogroup.     

Spatial analysis of Listronotus maculicollis immature stages demonstrates strong associations with conspecifics and turfgrass damage but not with optimal hosts on golf course fairways

For insects that develop on few hosts and/or have immobile immature stages, optimal oviposition theory suggests that females should seek high‐quality hosts that maximize larval development and reduce competition from conspecifics. However, there is a growing amount of evidence that suggests female choice may often be at odds with their offspring's development. Listronotus maculicollis (Kirby) (Coleoptera: Curculionidae) is a serious pest of golf course turfgrass in eastern North America. The weevil develops on few hosts and demonstrates improved fitness traits when developing on Poa annua L. (Poaceae). However, previous population studies observed either weak or no correlations between the spatial dispersion of larval populations and P. annua in the field. In this study, populations on three golf course fairways were monitored over a 4‐year period (2009–2012) to determine whether the lack of spatial associations between preferred hosts and immatures was a result of spatial scale or the density and distribution of conspecifics. Spatial Analysis by Distance IndicEs (SADIE) was used to characterize the spatial dispersion of populations of individual stages (larvae and pupae), P. annua, and turfgrass damage. Life stages were aggregated in each observation, independent of population density or the spatial dispersion of hosts. The distribution of consecutive and non‐consecutive immature stages was found to be correlated in all years, suggesting that females do not avoid patches already occupied by conspecific eggs. Surprisingly, significant spatial associations were not found between larvae and P. annua when the host plant was relatively abundant. Hence, multiple mechanisms may drive L. maculicollis oviposition site‐selection behavior, and a flexible strategy may allow the weevil to persist in areas where P. annua is not the dominant species. Future studies are required to determine what other factors (e.g., natural enemy‐free space, egg or time limitations) influence oviposition behavior.     

ADP‐glucose pyrophosphorylase genes,associated with kernel weight,underwent selection during wheat domestication and breeding

ADP‐glucose pyrophosphorylase, comprising two small subunits and two large subunits, is considered a key enzyme in the endosperm starch synthesis pathway in wheat (Triticum aestivum L.). Two genes, TaAGP‐S1‐7A and TaAGP‐L‐1B, were investigated in this study. Haplotypes of these genes were associated with thousand kernel weight (TKW) in different populations. Mean TKWs of favoured haplotypes were significantly higher than those of nonfavoured ones. Two molecular markers developed to distinguish these haplotypes could be used in molecular breeding. Frequencies of favoured haplotypes were dramatically increased in cultivars released in China after the 1940s. These favoured haplotypes were also positively selected in six major wheat production regions globally. Selection of AGP‐S1 and AGP‐L‐1B in wheat mainly occurred during and after hexaploidization. Strong additive effects of the favoured haplotypes of with other genes for starch synthesis were also detected in different populations.