Genomic analysis of a migratory divide reveals candidate genes for migration and implicates selective sweeps in generating islands of differentiation

Differential gene flow, reductions in diversity following linked selection and/or features of the genome can structure patterns of genomic differentiation during the process of speciation. Possible sources of reproductive isolation are well studied between coastal and inland subspecies groups of Swainson's thrushes, with differences in seasonal migratory behaviour likely playing a key role in reducing hybrid fitness. We assembled and annotated a draft reference genome for this species and generated whole‐genome shotgun sequence data for populations adjacent to the hybrid zone between these groups. We documented substantial genomewide heterogeneity in relative estimates of genetic differentiation between the groups. Within population diversity was lower in areas of high relative differentiation, supporting a role for selective sweeps in generating this pattern. Absolute genetic differentiation was reduced in these areas, further suggesting that recurrent selective sweeps in the ancestral population and/or between divergent populations following secondary contact likely occurred. Relative genetic differentiation was also higher near centromeres and on the Z chromosome, suggesting that features of the genome also contribute to genomewide heterogeneity. Genes linked to migratory traits were concentrated in islands of differentiation, supporting previous suggestions that seasonal migration is under divergent selection between Swainson's thrushes. Differences in migratory behaviour likely play a central role in the speciation of many taxa; we developed the infrastructure here to permit future investigations into the role several candidate genes play in reducing gene flow between not only Swainson's thrushes but other species as well.     

Genetic changes in a novel breeding population of Brassica napus synthesized from hundreds of crosses between B. rapa and B. carinata

Introgression of genomic variation between and within related crop species is a significant evolutionary approach for population differentiation, genome reorganization and trait improvement. Using the Illumina Infinium Brassica 60K SNP array, we investigated genomic changes in a panel of advanced generation new‐type Brassica napus breeding lines developed from hundreds of interspecific crosses between 122 Brassica rapa and 74 Brassica carinata accessions, and compared them with representative accessions of their three parental species. The new‐type B. napus population presented rich genetic diversity and abundant novel genomic alterations, consisting of introgressions from B. rapa and B. carinata, novel allelic combinations, reconstructed linkage disequilibrium patterns and haplotype blocks, and frequent deletions and duplications (nonrandomly distributed), particularly in the C subgenome. After a much shorter, but very intensive, selection history compared to traditional B. napus, a total of 15 genomic regions with strong selective sweeps and 112 genomic regions with putative signals of selective sweeps were identified. Some of these regions were associated with important agronomic traits that were selected for during the breeding process, while others were potentially associated with restoration of genome stability and fertility after interspecific hybridization. Our results demonstrate how a novel method for population‐based crop genetic improvement can lead to rapid adaptation, restoration of genome stability and positive responses to artificial selection.     

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.     

Genetic hitchhiking,mitonuclear coadaptation,and the origins of mt DNA barcode gaps

DNA barcoding based on mitochondrial (mt) nucleotide sequences is an enigma. Neutral models of mt evolution predict DNA barcoding cannot work for recently diverged taxa, and yet, mt DNA barcoding accurately delimits species for many bilaterian animals. Meanwhile, mt DNA barcoding often fails for plants and fungi. I propose that because mt gene products must cofunction with nuclear gene products, the evolution of mt genomes is best understood with full consideration of the two environments that impose selective pressure on mt genes: the external environment and the internal genomic environment. Moreover, it is critical to fully consider the potential for adaptive evolution of not just protein products of mt genes but also of mt transfer RNAs and mt ribosomal RNAs. The tight linkage of genes on mt genomes that do not engage in recombination could facilitate selective sweeps whenever there is positive selection on any element in the mt genome, leading to the purging of mt genetic diversity within a population and to the rapid fixation of novel mt DNA sequences. Accordingly, the most important factor determining whether or not mt DNA sequences diagnose species boundaries may be the extent to which the mt chromosomes engage in recombination.     

Detecting selection signatures in three Iranian sheep breeds

The objective of genome mapping is to achieve valuable insight into the connection between gene variants (genotype) and observed traits (phenotype). Part of that objective is to understand the selective forces that have operated on a population. Finding links between genotype–phenotype changes makes it possible to identify selective sweeps by patterns of genetic variation and linkage disequilibrium. Based on Illumina 50KSNP chip data, two approaches, XP‐EHH (cross‐population extend haplotype homozygosity) and F_ST (fixation index), were carried out in this research to identify selective sweeps in the genome of three Iranian local sheep breeds: Baluchi (n = 86), Lori‐Bakhtiari (n = 45) and Zel (n = 45). Using both methods, 93 candidate genomic regions were identified as harboring putative selective sweeps. Bioinformatics analysis of the genomic regions showed that signatures of selection related to multiple candidate genes, such as HOXB9, HOXB13, ACAN, NPR2, TRIL, AOX1, CSF2, GHR, TNS2, SPAG8, HINT2, ALS2, AAAS, RARG, SYCP2, CAV1, PPP1R3D, PLA2G7, TTLL7 and C20orf10, that play a role in skeletal system and tail, sugar and energy metabolisms, growth, reproduction, immune and nervous system traits. Our findings indicated diverse genomic selection during the domestication of Iranian sheep breeds.     

Adaptation,ancestral variation and gene flow in a ‘Sky Island’ Drosophila species

Over time, populations of species can expand, contract, fragment and become isolated, creating subpopulations that must adapt to local conditions. Understanding how species maintain variation after divergence as well as adapt to these changes in the face of gene flow is of great interest, especially as the current climate crisis has caused range shifts and frequent migrations for many species. Here, we characterize how a mycophageous fly species, Drosophila innubila, came to inhabit and adapt to its current range which includes mountain forests in south‐western USA separated by large expanses of desert. Using population genomic data from more than 300 wild‐caught individuals, we examine four populations to determine their population history in these mountain forests, looking for signatures of local adaptation. In this first extensive study, establishing D. innubila as a key genomic "Sky Island" model, we find D. innubila spread northwards during the previous glaciation period (30–100 KYA) and have recently expanded even further (0.2–2 KYA). D. innubila shows little evidence of population structure, consistent with a recent establishment and genetic variation maintained since before geographic stratification. We also find some signatures of recent selective sweeps in chorion proteins and population differentiation in antifungal immune genes suggesting differences in the environments to which flies are adapting. However, we find little support for long‐term recurrent selection in these genes. In contrast, we find evidence of long‐term recurrent positive selection in immune pathways such as the Toll signalling system and the Toll‐regulated antimicrobial peptides.     

Whole genome resequencing reveals genomic regions associated with thermal adaptation in redband trout

Adaptation to local environments involves evolution of ecologically important traits and underlying physiological processes. Here, we used low coverage whole‐genome resequencing (lcWGR) on individuals to identify genome regions involved in thermal adaptation in wild redband trout Oncorhynchus mykiss gairdneri, a subspecies of rainbow trout that inhabits ecosystems ranging from cold montane forests to high elevation deserts. This study includes allele frequency‐based analyses for selective sweeps among populations, followed by multiple association tests for specific sets of phenotypes measured under thermal stress (acute and chronic survival/mortality; high or low cardiac performance groups). Depending on the groups in each set of analyses, sequencing reads covered 43%–75% of the genome at ≥15× and each analysis included millions of SNPs across the genome. In tests for selective sweeps among populations, a total of six chromosomal regions were significant. The further association tests for specific phenotypes revealed that the region on chromosome 4 was consistently the most significant and contains the cerk gene (ceramide kinase). This study provides insight into a potential genetic mechanism of local thermal adaptation and suggests cerk may be an important candidate gene. However, further validation of this cerk gene is necessary to determine if the association with cardiac performance results in a functional role to influence thermal performance when exposed to high water temperatures and hypoxic conditions.     

The adaptive genomic landscape of beak morphology in Darwin's finches

Beak shape in Darwin's ground finches (Geospiza) is emblematic of natural selection and adaptive radiation, yet our understanding of the genetic basis of beak shape variation, and thus the genetic target of natural selection, is still evolving. Here we reveal the genomic architecture of beak shape variation using genomewide comparisons of four closely related and hybridizing species across 13 islands subject to parallel natural selection. Pairwise contrasts among species were used to identify a large number of genomic loci that are consistently related to species differences across a complex landscape. These loci are associated with hundreds of genes that have enriched GO categories significantly associated with development. One genomic region of particular interest is a section of Chromosome 1A with many candidate genes and increased linkage. The distinct, pointed beak shape of the cactus finch is linked to an excess of intermediate frequency alleles and increased heterozygosity in significant SNPs, but not across the rest of the genome. Alleles associated with pointier beaks among species were associated with pointier‐beaked populations within each species, thus establishing a common basis for natural selection, species divergence and adaptive radiation. The adaptive genomic landscape for Darwin's finches mirrors theoretical expectations based on morphological variation. The implication that a large number of genes are actively maintained to facilitate beak variation across parallel populations with documented interspecies admixture challenges our understanding of evolutionary processes in the wild.     

A draft genome of sweet cherry (Prunus avium L.) reveals genome‐wide and local effects of domestication

Sweet cherry (Prunus avium L.) trees are both economically important fruit crops but also important components of natural forest ecosystems in Europe, Asia and Africa. Wild and domesticated trees currently coexist in the same geographic areas with important questions arising on their historical relationships. Little is known about the effects of the domestication process on the evolution of the sweet cherry genome. We assembled and annotated the genome of the cultivated variety “Big Star*” and assessed the genetic diversity among 97 sweet cherry accessions representing three different stages in the domestication and breeding process (wild trees, landraces and modern varieties). The genetic diversity analysis revealed significant genome‐wide losses of variation among the three stages and supports a clear distinction between wild and domesticated trees, with only limited gene flow being detected between wild trees and domesticated landraces. We identified 11 domestication sweeps and five breeding sweeps covering, respectively, 11.0 and 2.4 Mb of the P. avium genome. A considerable fraction of the domestication sweeps overlaps with those detected in the related species, Prunus persica (peach), indicating that artificial selection during domestication may have acted independently on the same regions and genes in the two species. We detected 104 candidate genes in sweep regions involved in different processes, such as the determination of fruit texture, the regulation of flowering and fruit ripening and the resistance to pathogens. The signatures of selection identified will enable future evolutionary studies and provide a valuable resource for genetic improvement and conservation programs in sweet cherry.     

阿拉善荒漠区一年生植物层片物种多样性及其分布特征

一年生植物层片是阿拉善荒漠区荒漠与绿洲生态系统的恒有植物层片,对荒漠生态系统的稳定和当地畜牧业的生产都有重要作用．对该层片的物种多样性及其在群落中的分布特征分析表明,荒漠与绿洲生态系统中,一年生植物层片植物种有12科35属61种,可归为一年生小禾草、一年生猪毛菜、一年生蒿类和一年生杂类草4个类群,并可划分成12种区系地理成分,其中有10个阿拉善或近阿拉善特有种．一年生植物层片广泛分布于土质、沙砾质、沙质、砾质、石质、盐土荒漠及各种类型的绿洲群落中,常一种或几种组成层片,均匀或斑块状分布,或沿水线分布．此外,还分析和评价了一年生植物层片在群落中的分布差异、生态适应模式与功能．     

Recurrent selection explains parallel evolution of genomic regions of high relative but low absolute differentiation in a ring species

Recent technological developments allow investigation of the repeatability of evolution at the genomic level. Such investigation is particularly powerful when applied to a ring species, in which spatial variation represents changes during the evolution of two species from one. We examined genomic variation among three subspecies of the greenish warbler ring species, using genotypes at 13 013 950 nucleotide sites along a new greenish warbler consensus genome assembly. Genomic regions of low within‐group variation are remarkably consistent between the three populations. These regions show high relative differentiation but low absolute differentiation between populations. Comparisons with outgroup species show the locations of these peaks of relative differentiation are not well explained by phylogenetically conserved variation in recombination rates or selection. These patterns are consistent with a model in which selection in an ancestral form has reduced variation at some parts of the genome, and those same regions experience recurrent selection that subsequently reduces variation within each subspecies. The degree of heterogeneity in nucleotide diversity is greater than explained by models of background selection, but is consistent with selective sweeps. Given the evidence that greenish warblers have had both population differentiation for a long period of time and periods of gene flow between those populations, we propose that some genomic regions underwent selective sweeps over a broad geographic area followed by within‐population selection‐induced reductions in variation. An important implication of this ‘sweep‐before‐differentiation’ model is that genomic regions of high relative differentiation may have moved among populations more recently than other genomic regions.     

The genetic architecture of freezing tolerance varies across the range of Arabidopsis thaliana

The capacity to tolerate freezing temperatures limits the geographical distribution of many plants, including several species of agricultural importance. However, the genes involved in freezing tolerance remain largely unknown. Here, we describe the variation in constitutive freezing tolerance that occurs among worldwide accessions of Arabidopsis thaliana. We found that although plants from high latitudes tend to be more freezing tolerant than plants from low latitudes, the environmental factors that shape cold adaptation differ across the species range. Consistent with this, we found that the genetic architecture of freezing tolerance also differs across its range. Conventional genome‐wide association studies helped identify a priori and other promising candidate genes. However, simultaneously modelling climate variables and freezing tolerance together pinpointed other excellent a priori candidate genes. This suggests that if the selective factor underlying phenotypic variation is known, multi‐trait mixed models may aid in identifying the genes that underlie adaptation.     

A roadmap of tandemly arrayed genes in the genomes of human, mouse, and rat

Tandemly arrayed genes (TAGs) play an important functional and physiological role in the genome. Most previous studies have focused on individual TAG families in a few species, yet a broad characterization of TAGs is not available. Here we identified all TAGs in the genomes of humans, mouse, and rat and performed a comprehensive analysis of TAG distribution, TAG sizes, TAG orientations and intergenic distances, and TAG functions. TAGs account for about 14-17% of all genes in the genome and nearly one-third of all duplicated genes, highlighting the predominant role that tandem duplication plays in gene duplication. For all species, TAG distribution is highly heterogeneous along chromosomes and some chromosomes are enriched with TAG forests, whereas others are enriched with TAG deserts. The majority of TAGs are of size 2 for all genomes, similar to the previous findings in Caenorhabditis elegans, Arabidopsis thaliana, and Oryza sativa, suggesting that it is a rather general phenomenon in eukaryotes. The comparison with the genome patterns shows that TAG members have a significantly higher proportion of parallel gene orientation in all species, corroborating Graham's claim that parallel orientation is the preferred form of orientation in TAGs. Moreover, TAG members with parallel orientation tend to be closer to each other than all neighboring genes in the genome with parallel orientation. The analyses of Gene Ontology function indicate that genes with receptor or binding activities are significantly overrepresented by TAGs. Computer simulation reveals that random gene rearrangements have little effect on the statistics of TAGs for all genomes. Finally, the average proportion of TAGs shows a trend of increase with the increase of family sizes, although the correlation between TAG proportions in individual families and family sizes is not significant.     

Selective sweeps reveal candidate genes for adaptation to drought and salt tolerance in common sunflower, Helianthus annuus

Here we report the results of an analysis of variation at 128 EST-based microsatellites in wild Helianthus annuus, using populations from the species' typical plains habitat in Kansas and Colorado, as well as two arid desert and two distinct brackish marsh areas in Utah. The test statistics lnRV and lnRH were used to find regions of the genome that were significantly less variable in one population relative to the others and thus are likely to contain genes under selection. A small but detectable percentage (1.5-6%) of genes showed evidence for selection from both statistics in any particular environment, and a total of 17 loci showed evidence of selection in at least one environment. Distance-based measures provided additional evidence of selection for 15 of the 17 loci. Global F(ST)-values were significantly higher for candidate loci, as expected under divergent selection. However, pairwise F(ST)-values were lower for populations that shared a selective sweep. Moreover, while spatially separated populations undergoing similar selective pressures showed evidence of divergence at some loci, they evolved in concert at other loci. Thus, this study illustrates how selective sweeps might contribute both to the integration of conspecific populations and to the differentiation of races or species.     

Localizing recent adaptive evolution in the human genome

Identifying genomic locations that have experienced selective sweeps is an important first step toward understanding the molecular basis of adaptive evolution. Using statistical methods that account for the confounding effects of population demography, recombination rate variation, and single-nucleotide polymorphism ascertainment, while also providing fine-scale estimates of the position of the selected site, we analyzed a genomic dataset of 1.2 million human single-nucleotide polymorphisms genotyped in African-American, European-American, and Chinese samples. We identify 101 regions of the human genome with very strong evidence (p < 10⁻⁵) of a recent selective sweep and where our estimate of the position of the selective sweep falls within 100 kb of a known gene. Within these regions, genes of biological interest include genes in pigmentation pathways, components of the dystrophin protein complex, clusters of olfactory receptors, genes involved in nervous system development and function, immune system genes, and heat shock genes. We also observe consistent evidence of selective sweeps in centromeric regions. In general, we find that recent adaptation is strikingly pervasive in the human genome, with as much as 10% of the genome affected by linkage to a selective sweep.     

A Chromosome-level Genome Assembly of Wild Castor Provides New Insights into its Adaptive Evolution in Tropical Desert

Wild castor grows in the high-altitude tropical desert of the African Plateau,a region known for high ultraviolet radiation,strong light,and extremely dry condition.To investigate the potential genetic basis of adaptation to both highland and tropical deserts,we generated a chromosome-level genome sequence assembly of the wild castor accession WT05,with a genome size of 316 Mb,a scaffold N50 of 31.93 Mb,and a contig N50 of 8.96 Mb,respectively.Compared with cultivated castor and other Euphorbiac...     

Selection signatures for fiber production in commercial species: A review

Natural fibers derived from diverse animal species have gained increased attention in recent years due to their favorable environmental effects, long-term sustainability benefits, and remarkable physical and mechanical properties that make them valuable raw materials used for textile and non-textile production. Domestication and selective breeding for the economically significant fiber traits play an imperative role in shaping the genomes and, thus, positively impact the overall productivity of the various fiber-producing species. These selection pressures leave unique footprints on the genome due to alteration in the allelic frequencies at specific loci, characterizing selective sweeps. Recent advances in genomics have enabled the discovery of selection signatures across the genome using a variety of methods. The increased demand for ‘green products’ manufactured from natural fibers necessitates a detailed investigation of the genomes of the various fiber-producing plant and animal species to identify the candidate genes associated with important fiber attributes such as fiber diameter/fineness, color, length, and strength, among others. The objective of this review is to present a comprehensive overview of the concept of selection signature and selective sweeps, discuss the main methods used for its detection, and address the selection signature studies conducted so far in the diverse fiber-producing animal species.     

Divergent selection and drift shape the genomes of two avian sister species spanning a saline–freshwater ecotone

The role of species divergence due to ecologically based divergent selection—or ecological speciation—in generating and maintaining biodiversity is a central question in evolutionary biology. Comparison of the genomes of phylogenetically related taxa spanning a selective habitat gradient enables discovery of divergent signatures of selection and thereby provides valuable insight into the role of divergent ecological selection in speciation. Tidal marsh ecosystems provide tractable opportunities for studying organisms' adaptations to selective pressures that underlie ecological divergence. Sharp environmental gradients across the saline–freshwater ecotone within tidal marshes present extreme adaptive challenges to terrestrial vertebrates. Here, we sequence 20 whole genomes of two avian sister species endemic to tidal marshes—the saltmarsh sparrow (Ammospiza caudacutus) and Nelson's sparrow (A. nelsoni)—to evaluate the influence of selective and demographic processes in shaping genome‐wide patterns of divergence. Genome‐wide divergence between these two recently diverged sister species was notably high (genome‐wide F_ST = 0.32). Against a background of high genome‐wide divergence, regions of elevated divergence were widespread throughout the genome, as opposed to focused within islands of differentiation. These patterns may be the result of genetic drift resulting from past tidal march colonization events in conjunction with divergent selection to different environments. We identified several candidate genes that exhibited elevated divergence between saltmarsh and Nelson's sparrows, including genes linked to osmotic regulation, circadian rhythm, and plumage melanism—all putative candidates linked to adaptation to tidal marsh environments. These findings provide new insights into the roles of divergent selection and genetic drift in generating and maintaining biodiversity.