Contributions of flowering time genes to sunflower domestication and improvement

Determining the identity and distribution of molecular changes leading to the evolution of modern crop species provides major insights into the timing and nature of historical forces involved in rapid phenotypic evolution. In this study, we employed an integrated candidate gene strategy to identify loci involved in the evolution of flowering time during early domestication and modern improvement of the sunflower (Helianthus annuus). Sunflower homologs of many genes with known functions in flowering time were isolated and cataloged. Then, colocalization with previously mapped quantitative trait loci (QTLs), expression, or protein sequence differences between wild and domesticated sunflower, and molecular evolutionary signatures of selective sweeps were applied as step-wise criteria for narrowing down an original pool of 30 candidates. This process led to the discovery that five paralogs in the flowering locus T/terminal flower 1 gene family experienced selective sweeps during the evolution of cultivated sunflower and may be the causal loci underlying flowering time QTLs. Our findings suggest that gene duplication fosters evolutionary innovation and that natural variation in both coding and regulatory sequences of these paralogs responded to a complex history of artificial selection on flowering time during the evolution of cultivated sunflower.     

Association Mapping and the Genomic Consequences of Selection in Sunflower

The combination of large-scale population genomic analyses and trait-based mapping approaches has the potential to provide novel insights into the evolutionary history and genome organization of crop plants. Here, we describe the detailed genotypic and phenotypic analysis of a sunflower (Helianthus annuus L.) association mapping population that captures nearly 90% of the allelic diversity present within the cultivated sunflower germplasm collection. We used these data to characterize overall patterns of genomic diversity and to perform association analyses on plant architecture (i.e., branching) and flowering time, successfully identifying numerous associations underlying these agronomically and evolutionarily important traits. Overall, we found variable levels of linkage disequilibrium (LD) across the genome. In general, islands of elevated LD correspond to genomic regions underlying traits that are known to have been targeted by selection during the evolution of cultivated sunflower. In many cases, these regions also showed significantly elevated levels of differentiation between the two major sunflower breeding groups, consistent with the occurrence of divergence due to strong selection. One of these regions, which harbors a major branching locus, spans a surprisingly long genetic interval (ca. 25 cM), indicating the occurrence of an extended selective sweep in an otherwise recombinogenic interval.     

Selection under domestication: evidence for a sweep in the rice waxy genomic region

Rice (Oryza sativa) was cultivated by Asian Neolithic farmers >11,000 years ago, and different cultures have selected for divergent starch qualities in the rice grain during and after the domestication process. An intron 1 splice donor site mutation of the Waxy gene is responsible for the absence of amylose in glutinous rice varieties. This mutation appears to have also played an important role in the origin of low amylose, nonglutinous temperate japonica rice varieties, which form a primary component of Northeast Asian cuisines. Waxy DNA sequence analyses indicate that the splice donor mutation is prevalent in temperate japonica rice varieties, but rare or absent in tropical japonica, indica, aus, and aromatic varieties. Sequence analysis across a 500-kb genomic region centered on Waxy reveals patterns consistent with a selective sweep in the temperate japonicas associated with the mutation. The size of the selective sweep (>250 kb) indicates very strong selection in this region, with an inferred selection coefficient that is higher than similar estimates from maize domestication genes or wild species. These findings demonstrate that selection pressures associated with crop domestication regimes can exceed by one to two orders of magnitude those observed for genes under even strong selection in natural systems.     

Comparative population genomics identified genomic regions and candidate genes associated with fruit domestication traits in peach

Crop evolution is a long‐term process involving selection by natural evolutionary forces and anthropogenic influences; however, the genetic mechanisms underlying the domestication and improvement of fruit crops have not been well studied to date. Here, we performed a population structure analysis in peach (Prunus persica) based on the genome‐wide resequencing of 418 accessions and confirmed the presence of an obvious domestication event during evolution. We identified 132 and 106 selective sweeps associated with domestication and improvement, respectively. Analysis of their tissue‐specific expression patterns indicated that the up‐regulation of selection genes during domestication occurred mostly in fruit and seeds as opposed to other organs. However, during the improvement stage, more up‐regulated selection genes were identified in leaves and seeds than in the other organs. Genome‐wide association studies (GWAS) using 4.24 million single nucleotide polymorphisms (SNPs) revealed 171 loci associated with 26 fruit domestication traits. Among these loci, three candidate genes were highly associated with fruit weight and the sorbitol and catechin content in fruit. We demonstrated that as the allele frequency of the SNPs associated with high polyphenol composition decreased during peach evolution, alleles associated with high sugar content increased significantly. This indicates that there is genetic potential for the breeding of more nutritious fruit with enhanced bioactive polyphenols without disturbing a harmonious sugar and acid balance by crossing with wild species. This study also describes the development of the genomic resources necessary for evolutionary research in peach and provides the large‐scale characterization of key agronomic traits in this crop species.     

Genomics of invasion: diversity and selection in introduced populations of monkeyflowers (Mimulus guttatus)

Global trade and travel is irreversibly changing the distribution of species around the world. Because introduced species experience drastic demographic events during colonization and often face novel environmental challenges from their native range, introduced populations may undergo rapid evolutionary change. Genomic studies provide the opportunity to investigate the extent to which demographic, historical and selective processes shape the genomic structure of introduced populations by analysing the signature that these processes leave on genomic variation. Here, we use next‐generation sequencing to compare genome‐wide relationships and patterns of diversity in native and introduced populations of the yellow monkeyflower (Mimulus guttatus). Genome resequencing data from 10 introduced populations from the United Kingdom (UK) and 12 native M. guttatus populations in North America (NA) demonstrated reduced neutral genetic diversity in the introduced range and showed that UK populations are derived from a geographic region around the North Pacific. A selective‐sweep analysis revealed site frequency changes consistent with selection on five of 14 chromosomes, with genes in these regions showing reduced silent site diversity. While the target of selection is unknown, genes associated with flowering time and biotic and abiotic stresses were located within the swept regions. The future identification of the specific source of origin of introduced UK populations will help determining whether the observed selective sweeps can be traced to unsampled native populations or occurred since dispersal across the Atlantic. Our study demonstrates the general potential of genome‐wide analyses to uncover a range of evolutionary processes affecting invasive populations.     

A Genomic Scan for Selection Reveals Candidates for Genes Involved in the Evolution of Cultivated Sunflower (Helianthus annuus)

Genomic scans for selection are a useful tool for identifying genes underlying phenotypic transitions. In this article, we describe the results of a genome scan designed to identify candidates for genes targeted by selection during the evolution of cultivated sunflower. This work involved screening 492 loci derived from ESTs on a large panel of wild, primitive (i.e., landrace), and improved sunflower (Helianthus annuus) lines. This sampling strategy allowed us to identify candidates for selectively important genes and investigate the likely timing of selection. Thirty-six genes showed evidence of selection during either domestication or improvement based on multiple criteria, and a sequence-based test of selection on a subset of these loci confirmed this result. In view of what is known about the structure of linkage disequilibrium across the sunflower genome, these genes are themselves likely to have been targeted by selection, rather than being merely linked to the actual targets. While the selection candidates showed a broad range of putative functions, they were enriched for genes involved in amino acid synthesis and protein catabolism. Given that a similar pattern has been detected in maize (Zea mays), this finding suggests that selection on amino acid composition may be a general feature of the evolution of crop plants. In terms of genomic locations, the selection candidates were significantly clustered near quantitative trait loci (QTL) that contribute to phenotypic differences between wild and cultivated sunflower, and specific instances of QTL colocalization provide some clues as to the roles that these genes may have played during sunflower evolution.     

Population genetic evidence for rapid changes in intraspecific diversity and allelic cycling of a specialist defense gene in Zea

Two patterns of plant defense gene evolution are emerging from molecular population genetic surveys. One is that specialist defenses experience stronger selection than generalist defenses. The second is that specialist defenses are more likely to be subject to balancing selection, i.e., evolve in a manner consistent with balanced-polymorphism or trench-warfare models of host-parasite coevolution. Because most of the data of specialist defenses come from Arabidopsis thaliana, we examined the genetic diversity and evolutionary history of three defense genes in two outcrossing species, the autotetraploid Zea perennis and its most closely related extant relative the diploid Z. diploperennis. Intraspecific diversity at two generalist defenses, the protease inhibitors wip1 and mpi, were consistent with a neutral model. Like previously studied genes in these taxa, wip1 and mpi harbored similar levels of diversity in Z. diploperennis and Z. perennis. In contrast, the specialist defense hm2 showed strong although distinctly different departures from a neutral model in the two species. Z. diploperennis appears to have experienced a strong and recent selective sweep. Using a rejection-sampling coalescent method, we estimate the strength of selection on Z. diploperennis hm2 to be approximately 3.0%, which is approximately equal to the strength of selection on tb1 during maize domestication. Z. perennis hm2 harbors three highly diverged alleles, two of which are found at high frequency. The distinctly different patterns of diversity may be due to differences in the phase of host-parasite coevolutionary cycles, although higher hm2 diversity in Z. perennis may also reflect reduced efficacy of selection in the autotetraploid relative to its diploid relative.     

Whole-genome resequencing reveals domestication and signatures of selection in Ujimqin,Sunit, and Wu Ranke Mongolian sheep breeds

ObjectiveThe current study aimed to perform whole-genome resequencing of Chinese indigenous Mongolian sheep breeds including Ujimqin, Sunit, and Wu Ranke sheep breeds (UJMQ, SNT, WRK) and deeply analyze genetic variation, population structure, domestication, and selection for domestication traits among these Mongolian sheep breeds.MethodsBlood samples were collected from a total of 60 individuals comprising 20 WRK, 20 UJMQ, and 20 SNT. For genome sequencing, about 1.5 μg of genomic DNA was used for library construction with an insert size of about 350 bp. Pair-end sequencing were performed on Illumina NovaSeq platform, with the read length of 150 bp at each end. We then investigated the domestication and signatures of selection in these sheep breeds.ResultsAccording to the population and demographic analyses, WRK and SNT populations were very similar, which were different from UJMQ populations. Genome wide association study identified 468 and 779 significant loci from SNT vs UJMQ, and UJMQ vs WRK, respectively. However, only 3 loci were identified from SNT vs WRK. Genomic comparison and selective sweep analysis among these sheep breeds suggested that genes associated with regulation of secretion, metabolic pathways including estrogen metabolism and amino acid metabolism, and neuron development have undergone strong selection during domestication.ConclusionOur findings will facilitate the understanding of Chinese indigenous Mongolian sheep breeds domestication and selection for complex traits and provide a valuable genomic resource for future studies of sheep and other domestic animal breeding.     

Molecular evidence for post-domestication selection in the <Emphasis Type="Italic">Waxy</Emphasis> gene of Chinese waxy maize

Waxy maize was first reported in China in 1909 and is mainly used in food production in Asia. The evidence for strong domestication selection in the Waxy locus of rice and a selective sweep around its genomic region make us to wonder whether there has been similar selection in Waxy in glutinous maize. To address this issue, DNA sequences of Waxy, three flanking genes and an unlinked gene (Adh1) of 30 accessions sampled from Chinese waxy maize accessions, including representative landraces and inbred lines, were determined in this study. Sharp reduction of nucleotide diversity and significant neutrality tests (Tajima’s D and Fu and Li’s F*) were observed in the Waxy locus in Chinese waxy maize but not in nonglutinous maize; comparison with the unlinked gene confirmed that this pattern was different to Waxy. Sequence analysis across a 143 kb genomic segment centered on the Waxy locus revealed patterns consistent with a selective sweep in the upstream region of Waxy. The selective sweep detected based on current limited genomic sequences exceeded over 50 kb, indicating strong selection in this or a bigger region. However, No sweep effect was detected in the repetitive downstream region of Waxy. Phylogenetic analysis indicated that Chinese waxy maize was domesticated from the cultivated flint maize (Zea mays ssp. mays) that was introduced from the new world. At least two independent deletions in exon 7 (30 bp) and 10 (15 bp) were identified in the Chinese accessions respectively. These findings demonstrate a similar pattern of domestication selection in the Waxy genomic region in both glutinous maize and rice, suggesting that this pattern in the rise of glutinous phenotype is likely in other cereal crops.     

An analysis of selection on candidate genes for regulation, mobilization, uptake, and transport of iron in maize

Insufficient iron (Fe) availability, which frequently occurs in soils with high pH levels, can lead to leaf chlorosis, a reduced Fe content in harvest products, and yield reduction in maize. The objectives of this study were (i) to describe patterns of sequence variation of 14 candidate genes for mobilization, uptake, and transport of Fe in maize, as well as regulatory function on these processes; (ii) to examine whether Fe-efficiency is an adaptive trait by determining if these genes were targets of selection during domestication; and (iii) to test if the allele distribution at these candidate genes is different for the different subpopulations of maize. The nucleotide diversity of Mtk was reduced by 78% in maize compared with teosinte. The results of our study revealed for the genes Naat1, Nas1, Nramp3, Mtk, and Ys1 a selective sweep, which suggests that these genes might be important for the fast adaptation of maize to new environments with different Fe availabilities.     

Classic Selective Sweeps Revealed by Massive Sequencing in Cattle

Human driven selection during domestication and subsequent breed formation has likely left detectable signatures within the genome of modern cattle. The elucidation of these signatures of selection is of interest from the perspective of evolutionary biology, and for identifying domestication-related genes that ultimately may help to further genetically improve this economically important animal. To this end, we employed a panel of more than 15 million autosomal SNPs identified from re-sequencing of 43 Fleckvieh animals. We mainly applied two somewhat complementary statistics, the integrated Haplotype Homozygosity Score (iHS) reflecting primarily ongoing selection, and the Composite of Likelihood Ratio (CLR) having the most power to detect completed selection after fixation of the advantageous allele. We find 106 candidate selection regions, many of which are harboring genes related to phenotypes relevant in domestication, such as coat coloring pattern, neurobehavioral functioning and sensory perception including KIT, MITF, MC1R, NRG4, Erbb4, TMEM132D and TAS2R16, among others. To further investigate the relationship between genes with signatures of selection and genes identified in QTL mapping studies, we use a sample of 3062 animals to perform four genome-wide association analyses using appearance traits, body size and somatic cell count. We show that regions associated with coat coloring significantly (P<0.0001) overlap with the candidate selection regions, suggesting that the selection signals we identify are associated with traits known to be affected by selection during domestication. Results also provide further evidence regarding the complexity of the genetics underlying coat coloring in cattle. This study illustrates the potential of population genetic approaches for identifying genomic regions affecting domestication-related phenotypes and further helps to identify specific regions targeted by selection during speciation, domestication and breed formation of cattle. We also show that Linkage Disequilibrium (LD) decays in cattle at a much faster rate than previously thought.     

Molecular evolution and phylogenetic analysis of genes related to cotton fibers development from wild and domesticated cotton species in Gossypium

The domestication of both diploid and tetraploid cotton species was carried out for fiber utilization. To understand the origin and domestication of fibers, 18 genes related to fiber development were individually cloned and sequenced from 22 different cotton species. Their structures, phylogenetic relationship and molecular evolution were further studied. In the orthologous and homeologous loci of the 18 genes, the sequence and structure of 72.22% were conserved and 27.78% were diverse. Tree topologies constructed based on the combined sequences showed that all 13 D-genome species were congruent with Fryxell's subsection taxonomy, the A- and D-subgenomes independently evolved in the allopolyploid after polyploid formation, and Gossypium raimondii had the closest relationship with all allotetraploids of D-subgenomes. The molecular evolutionary rates revealed approximately equivalent rates among different D-genome species, and purifying selection acted on all genes in the wild D-genome species. Among orthologs and homeologs, the D-subgenomes had higher evolutionary rates than the A-subgenomes in tetraploid cotton species, and the cultivars had higher evolutionary rates than either the semi-domesticated or wild species. Our study revealed that human domestication altered the molecular evolutionary pattern of genes related to fiber development, and Gossypium hirsutum endured greater selective pressures than Gossypium barbadense during the domestication process.     

Maize Sh2 gene is constrained by natural selection but escaped domestication

In Zea mays L., we studied the molecular evolution of Shrunken2 (Sh2), a gene that encodes the large subunits of a major enzyme in endosperm starch biosynthesis, ADP-glucose pyrophosphorylase. We compared 4669 bp of the Sh2 coding region on 50 accessions of maize and teosinte. Very few nucleotide polymorphisms were found when compared with other genes in Z. mays, revealing an effect of purifying selection in the whole species that predates domestication. Additionally, the comparison of Sh2 sequences in all Z. mays subspecies and outgroups Z. diploperennis and Tripsacum dactyloides suggests the occurrence of an ancient selective sweep in the Sh2 3' region. The amount and nature of nucleotide diversity are similar in both maize and teosinte, confirming previous results that suggested that Sh2 has not been involved in maize domestication. The very low level of nucleotide diversity as well as the highly conserved protein sequence suggest that natural selection retained effective Sh2 allele(s) long before agriculture started, making human selection inefficient on this gene.     

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.     

Genome-wide footprints of pig domestication and selection revealed through massive parallel sequencing of pooled DNA

Background

Artificial selection has caused rapid evolution in domesticated species. The identification of selection footprints across domesticated genomes can contribute to uncover the genetic basis of phenotypic diversity.

Methodology/Main Findings

Genome wide footprints of pig domestication and selection were identified using massive parallel sequencing of pooled reduced representation libraries (RRL) representing ∼2% of the genome from wild boar and four domestic pig breeds (Large White, Landrace, Duroc and Pietrain) which have been under strong selection for muscle development, growth, behavior and coat color. Using specifically developed statistical methods that account for DNA pooling, low mean sequencing depth, and sequencing errors, we provide genome-wide estimates of nucleotide diversity and genetic differentiation in pig. Widespread signals suggestive of positive and balancing selection were found and the strongest signals were observed in Pietrain, one of the breeds most intensively selected for muscle development. Most signals were population-specific but affected genomic regions which harbored genes for common biological categories including coat color, brain development, muscle development, growth, metabolism, olfaction and immunity. Genetic differentiation in regions harboring genes related to muscle development and growth was higher between breeds than between a given breed and the wild boar.

Conclusions/Significance

These results, suggest that although domesticated breeds have experienced similar selective pressures, selection has acted upon different genes. This might reflect the multiple domestication events of European breeds or could be the result of subsequent introgression of Asian alleles. Overall, it was estimated that approximately 7% of the porcine genome has been affected by selection events. This study illustrates that the massive parallel sequencing of genomic pools is a cost-effective approach to identify footprints of selection.     

The western Mediterranean region provided the founder population of domesticated narrow-leafed lupin

Key message

This study revealed that the western Mediterranean provided the founder population for domesticated narrow-leafed lupin and that genetic diversity decreased significantly during narrow-leafed lupin domestication.

Abstract

The evolutionary history of plants during domestication profoundly shaped the genome structure and genetic diversity of today’s crops. Advances in next-generation sequencing technologies allow unprecedented opportunities to understand genome evolution in minor crops, which constitute the majority of plant domestications. A diverse set of 231 wild and domesticated narrow-leafed lupin (Lupinus angustifolius L.) accessions were subjected to genotyping-by-sequencing using diversity arrays technology. Phylogenetic, genome-wide divergence and linkage disequilibrium analyses were applied to identify the founder population of domesticated narrow-leafed lupin and the genome-wide effect of domestication on its genome. We found wild western Mediterranean population as the founder of domesticated narrow-leafed lupin. Domestication was associated with an almost threefold reduction in genome diversity in domesticated accessions compared to their wild relatives. Selective sweep analysis identified no significant footprints of selection around domestication loci. A genome-wide association study identified single nucleotide polymorphism markers associated with pod dehiscence. This new understanding of the genomic consequences of narrow-leafed lupin domestication along with molecular marker tools developed here will assist plant breeders more effectively access wild genetic diversity for crop improvement.

     

Genetic consequences of selection during the evolution of cultivated sunflower

We mapped quantitative trait loci (QTL) controlling differences in seed oil content and composition between cultivated and wild sunflower and used the results, along with those of a previous study of domestication-related QTL, to guide a genome-wide analysis of genetic variation for evidence of past selection. The effects of the seed oil QTL were almost exclusively in the expected direction with respect to the parental phenotypes. A major, oil-related QTL cluster mapped near a cluster of domestication-related QTL on linkage group six (LG06), the majority of which have previously been shown to have effects that are inconsistent with the parental phenotypes. To test the hypothesis that this region was the target of a past selective sweep, perhaps resulting in the fixation of the antagonistic domestication-related QTL, we analyzed simple sequence repeat (SSR) diversity from 102 markers dispersed throughout the sunflower genome. Our results indicate that LG06 was most likely the target of multiple selective sweeps during the postdomestication era. Strong directional selection in concert with genetic hitchhiking therefore offers a possible explanation for the occurrence of numerous domestication-related QTL with apparently maladaptive phenotypic effects.