Searching for the bull's eye: agents and targets of selection vary among geographically disparate cyanogenesis clines in white clover (Trifolium repens L.)

The recurrent evolution of adaptive clines within a species can be used to elucidate theselective factors and genetic responses that underlie adaptation. White clover ispolymorphic for cyanogenesis (HCN release with tissue damage), and climate-associatedcyanogenesis clines have evolved throughout the native and introduced species range. Thispolymorphism arises through two independently segregating Mendelian polymorphisms for thepresence/absence of two required components: cyanogenic glucosides and theirhydrolyzing enzyme linamarase. Cyanogenesis is commonly thought to function in herbivoredefense; however, the individual cyanogenic components may also serve other physiologicalfunctions. To test whether cyanogenesis clines have evolved in response to the sameselective pressures acting on the same genetic targets, we examined cyanogenesis clineshape and its environmental correlates in three world regions: southern New Zealand, thecentral United States and the US Pacific Northwest. For some regional comparisons, clineshapes are remarkably similar despite large differences in the spatial scales over whichclines occur (40–1600 km). However, we also find evidence for majordifferences in both the agents and targets of selection among the sampled clines.Variation in cyanogenesis frequency is best predicted using a combination of minimumwinter temperature and aridity variables. Together, our results provide evidence thatrecurrent adaptive clines do not necessarily reflect shared adaptive processes.     

Rapid evolution of an adaptive cyanogenesis cline in introduced North American white clover (Trifolium repens L.)

White clover is polymorphic for cyanogenesis (HCN production after tissue damage), and this herbivore defence polymorphism has served as a classic model for studying adaptive variation. The cyanogenic phenotype requires two interacting biochemical components; the presence/absence of each component is controlled by a simple Mendelian gene (Ac/ac and Li/li). Climate-associated cyanogenesis clines occur in both native (Eurasian) and introduced populations worldwide, with cyanogenic plants predominating in warmer locations. Moreover, previous studies have suggested that epistatic selection may act within populations to maintain cyanogenic (AcLi) plants and acyanogenic plants that lack both components (acli plants) at the expense of plants possessing a single component (Acli and acLi plants). Here, we examine the roles of selection, gene flow and demography in the evolution of a latitudinal cyanogenesis cline in introduced North American populations. Using 1145 plants sampled across a 1650 km transect, we determine the distribution of cyanogenesis variation across the central United States and investigate whether clinal variation is adaptive or an artefact of population introduction history. We also test for the evidence of epistatic selection. We detect a clear latitudinal cline, with cyanogenesis frequencies increasing from 11% to 86% across the transect. Population structure analysis using nine microsatellite loci indicates that the cline is adaptive and not a by-product of demographic history. However, we find no evidence for epistatic selection within populations. Our results provide strong evidence for rapid adaptive evolution in these introduced populations, and they further suggest that the mechanisms maintaining adaptive variation may vary among populations of a species.     

Evolution of the cyclin gene family in plants

Cyclins are key regulators of cell cycle progression. Previous studies have shown that cyclin genes in plants can be divided into 10 groups. However, because those studies only focused on genes from two well-known model plants (i.e., Arabidopsis thaliana (L.) Heynh. and Oryza sativa L.), it remains unclear whether the 10 groups are reasonably defined. In this study, by analyzing the genomes of 10 representative plants (Chlamydomonas reinhardtii P. A. Dang, Physcomitrella patens(Hedw.) Bruch & Schimp., Selaginella moellendorffii Hieron., Picea abies (L.) H. Karst., Amborella trichopoda Baill., A. thaliana, Populus trichocarpa Torr. & A. Gray ex Hook., Vitis vinifera L., O. sativa, and Sorghum bicolor (L.) Moench), we estimated the phylogenetic relationships of plant cyclins and investigated their evolutionary patterns. We confirmed that plant cyclins can be classified into 10 groups, although only eight ancestral genes may have existed in the most recent common ancestor of extant green plants. We also found that, due to the frequent occurrences of gene duplication events, several groups have expanded extensively in seed plants and, particularly, flowering plants, so that multiple genes belonging to different subgroups are present in a species. Reconciliation of the evolutionary histories of these groups and subgroups further led to the identification of evolutionarily highly conserved and rapidly duplicating gene lineages. These results will guide the classification and nomenclature of plant cyclins and help understand the conservativeness and variation in their functions.     

Large chromosome deletions,duplications, and gene conversion events accumulate with age in normal human colon crypts

Little is known about the types and numbers of mutations that may accumulate in normal human cells with age. Such information would require obtaining enough DNA from a single cell to accurately carry out reliable analysis despite extensive amplification; and complete genomic coverage under these circumstances is difficult. We have compared colon crypts, which are putatively clonal and contain ~2000 cells each, to determine how much somatic genetic variation occurs in vivo (without ex vivo cell culturing). Using high‐density SNP microarrays, we find that chromosome deletions, duplications, and gene conversions were significantly more frequent in colons from the older individuals. These changes affected lengths ranging from 73 kb to 46 Mb. Although detection requires progeny of a single mutant stem cell to reach niche dominance over neighboring stem cells, none of the deletions appear likely to confer a selective advantage. Mutations can become fixed randomly during stem cell evolution through neutral drift in normal human crypts. The fact that chromosomal changes are detected in individual crypts with increasing age suggests that either such changes accumulate with age or single stem cell dominance increases with age, and the former is more likely. This progressive genome‐wide divergence of human somatic cells with age has implications for aging and disease in multicellular organisms.     

Parallel selection on gene copy number variations through evolution of three-spined stickleback genomes

Background

Understanding the genetic basis of adaptive evolution is one of the major goals in evolutionary biology. Recently, it has been revealed that gene copy number variations (GCNVs) constitute significant proportions of genomic diversities within natural populations. However, it has been unclear whether GCNVs are under positive selection and contribute to adaptive evolution. Parallel evolution refers to adaptive evolution of the same trait in related but independent lineages, and three-spined stickleback (Gasterosteus aculeatus) is a well-known model organism. Through identification of genetic variations under parallel selection, i.e., variations shared among related but independent lineages, evidence of positive selection is obtained. In this study, we investigated whole-genome resequencing data from the marine and freshwater groups of three-spined sticklebacks from diverse areas along the Pacific and Atlantic Ocean coastlines, and searched for GCNVs under parallel selection.

Results

We identified 24 GCNVs that showed significant differences in the numbers of mapped reads between the two groups, and this number was significantly larger than that expected by chance. The derived group, i.e., freshwater group, was typically characterized by larger gene-copy numbers, which implied that gene duplications or multiplications helped with adaptation to the freshwater environment. Some of the identified GCNVs were those of multigenic family genes, which is consistent with the theory that fatal effects due to copy-number changes of multigenic family genes tend to be less than those of single-copy genes.

Conclusion

The identification of GCNVs that were likely under parallel selection suggests that contribution of GCNVs should be considered in studies on adaptive evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-735) contains supplementary material, which is available to authorized users.     

Targeted resequencing reveals geographical patterns of differentiation for loci implicated in parallel evolution

Parallel divergence and speciation provide evidence for the role of divergent selection in generating biological diversity. Recent studies indicate that parallel phenotypic divergence may not have the same genetic basis in different geographical locations – ‘outlier loci’ (loci potentially affected by divergent selection) are often not shared among parallel instances of phenotypic divergence. However, limited sharing may be due, in part, to technical issues if false‐positive outliers occur. Here, we test this idea in the marine snail Littorina saxatilis, which has evolved two partly isolated ecotypes (adapted to crab predation vs. wave action) in multiple locations independently. We argue that if the low extent of sharing observed in earlier studies in this system is due to sampling effects, we expect outliers not to show elevated F_ST when sequenced in new samples from the original locations and also not to follow predictable geographical patterns of elevated F_ST. Following a hierarchical sampling design (within vs. between country), we applied capture sequencing, targeting outliers from earlier studies and control loci. We found that outliers again showed elevated levels of F_ST in their original location, suggesting they were not generated by sampling effects. Outliers were also likely to show increased F_ST in geographically close locations, which may be explained by higher levels of gene flow or shared ancestral genetic variation compared with more distant locations. However, in contrast to earlier findings, we also found some outlier types to show elevated F_ST in geographically distant locations. We discuss possible explanations for this unexpected result.     

Similarity-Based Analysis of Allele Frequency Distribution among Multiple Populations Identifies Adaptive Genomic Structural Variants

Structural variants have a considerable impact on human genomic diversity. However, their evolutionary history remains mostly unexplored. Here, we developed a new method to identify potentially adaptive structural variants based on a similarity-based analysis that incorporates genotype frequency data from 26 populations simultaneously. Using this method, we analyzed 57,629 structural variants and identified 576 structural variants that show unusual population differentiation. Of these putatively adaptive structural variants, we further showed that 24 variants are multiallelic and overlap with coding sequences, and 20 variants are significantly associated with GWAS traits. Closer inspection of the haplotypic variation associated with these putatively adaptive and functional structural variants reveals deviations from neutral expectations due to: 1) population differentiation of rapidly evolving multiallelic variants, 2) incomplete sweeps, and 3) recent population-specific negative selection. Overall, our study provides new methodological insights, documents hundreds of putatively adaptive variants, and introduces evolutionary models that may better explain the complex evolution of structural variants.     

Fixation dynamics of beneficial alleles in prokaryotic polyploid chromosomes and plasmids

Theoretical population genetics has been mostly developed for sexually reproducing diploid and for monoploid (haploid) organisms, focusing on eukaryotes. The evolution of bacteria and archaea is often studied by models for the allele dynamics in monoploid populations. However, many prokaryotic organisms harbor multicopy replicons—chromosomes and plasmids—and theory for the allele dynamics in populations of polyploid prokaryotes remains lacking. Here, we present a population genetics model for replicons with multiple copies in the cell. Using this model, we characterize the fixation process of a dominant beneficial mutation at 2 levels: the phenotype and the genotype. Our results show that depending on the mode of replication and segregation, the fixation of the mutant phenotype may precede genotypic fixation by many generations; we term this time interval the heterozygosity window. We furthermore derive concise analytical expressions for the occurrence and length of the heterozygosity window, showing that it emerges if the copy number is high and selection strong. Within the heterozygosity window, the population is phenotypically adapted, while both alleles persist in the population. Replicon ploidy thus allows for the maintenance of genetic variation following phenotypic adaptation and consequently for reversibility in adaptation to fluctuating environmental conditions.     

Ongoing transposition in cell culture reveals the phylogeny of diverse Drosophila S2 sublines

Cultured cells are widely used in molecular biology despite poor understanding of how cell line genomes change in vitro over time. Previous work has shown that Drosophila cultured cells have a higher transposable element content than whole flies, but whether this increase in transposable element content resulted from an initial burst of transposition during cell line establishment or ongoing transposition in cell culture remains unclear. Here, we sequenced the genomes of 25 sublines of Drosophila S2 cells and show that transposable element insertions provide abundant markers for the phylogenetic reconstruction of diverse sublines in a model animal cell culture system. DNA copy number evolution across S2 sublines revealed dramatically different patterns of genome organization that support the overall evolutionary history reconstructed using transposable element insertions. Analysis of transposable element insertion site occupancy and ancestral states support a model of ongoing transposition dominated by episodic activity of a small number of retrotransposon families. Our work demonstrates that substantial genome evolution occurs during long-term Drosophila cell culture, which may impact the reproducibility of experiments that do not control for subline identity.     

Genomic Architecture of Rapid Parallel Adaptation to Fresh Water in a Wild Fish

Rapid adaptation to novel environments may drive changes in genomic regions through natural selection. However, the genetic architecture underlying these adaptive changes is still poorly understood. Using population genomic approaches, we investigated the genomic architecture that underlies rapid parallel adaptation of Coilia nasus to fresh water by comparing four freshwater-resident populations with their ancestral anadromous population. Linkage disequilibrium network analysis and population genetic analyses revealed two putative large chromosome inversions on LG6 and LG22, which were enriched for outlier loci and exhibited parallel association with freshwater adaptation. Drastic frequency shifts and elevated genetic differentiation were observed for the two chromosome inversions among populations, suggesting that both inversions would undergo divergent selection between anadromous and resident ecotypes. Enrichment analysis of genes within chromosome inversions showed significant enrichment of genes involved in metabolic process, immunoregulation, growth, maturation, osmoregulation, and so forth, which probably underlay differences in morphology, physiology and behavior between the anadromous and freshwater-resident forms. The availability of beneficial standing genetic variation, large optimum shift between marine and freshwater habitats, and high efficiency of selection with large population size could lead to the observed rapid parallel adaptive genomic change. We propose that chromosomal inversions might have played an important role during the evolution of rapid parallel ecological divergence in the face of environmental heterogeneity in C. nasus. Our study provides insights into the genomic basis of rapid adaptation of complex traits in novel habitats and highlights the importance of structural genomic variants in analyses of ecological adaptation.     

Adaptive evolution of the insulin gene in caviomorph rodents

Insulin is a conservative molecule among mammals, maintaining both its structure and function. Rodents that belong to the Suborder Hystricognathi represent an exception, having a very divergent molecule with unusual physiological properties. In this work, we analyzed the evolutionary pattern of the insulin gene in caviomorph rodents (South American hystricomorph rodents). We found that these rodents have higher rates of nonsynonymous:synonymous substitutions (d(N)/d(S)) than nonhystricomorph rodents and that values are heterogeneous inside the group. We estimated codons under positive selection, specifically the second binding site (A13 and B17) and others related with hexamerization (B18, B20, and B22). In the monomer structure, all selected sites formed a single patch around the second binding site. In the hexamer structure, these amino acids were grouped into three major patches. In this structure, contacts between B chains involved all selected sites (except B18), and between faces in the center of the molecule, all contacts were among selected sites. While there is no clear hypothesis regarding the cause of this drastic change, experimental evidence does show that this group of rodents has some peculiarities in growth function, and, whether coincidental or not, these changes appeared together with important changes in life-history traits.     

A study of the diversity and geographical variation in numbers of leg‐bearing segments in centipedes (Chilopoda: Geophilomorpha) in north‐western Europe

Five species of geophilomorphs, Strigamia maritima (Leach, 1817), Geophilus flavus (De Geer, 1778), Geophilus truncorum Bergsøe and Meinert, 1866, Geophilus proximus C. L. Koch, 1847, and Pachymerium ferrugineum (C. L. Koch, 1835), from various sample sites in north‐western Europe were examined for numbers of leg‐bearing segments. Where data was adequate, mean segment numbers were correlated with latitudinal gradients for each species. Solely in S. maritima and some populations of P. ferrugineum did the results show a highly significant geographic pattern towards more leg‐bearing segments in southern populations of both sexes. As concerns G. proximus, we presented for the first time a remarkable geographic shift towards an increased number of pairs of legs in northern populations. Contrary to the conventional geographic pattern, we found that G. flavus and G. truncorum did not exhibit a north–south increase or decrease in segment number. Although there was no general/universal evidence supporting the occurrence of a significant regression slope between mean segment number and latitude/temperature, more information shows that the overall region‐to‐region segment variation was highly significant in both sexes. In S. maritima and P. ferrugineum mean segment number was correlated with the north–south temperature cline. The same was not observed in G. proximus. Parthenogenesis in G. proximus, and a series of ecological characteristics such as habitat preferences, spatial distribution, and fragmented populations could explain the presence or absence of a geographic patterning of segment number variation along a latitudinal cline. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 899–909.     

Gene copy number variations in adaptive evolution: The genomic distribution of gene copy number variations revealed by genetic mapping and their adaptive role in an undomesticated species,white spruce (Picea glauca)

Gene copy number variation (CNV) has been associated with phenotypic variability in animals and plants, but a genomewide understanding of their impacts on phenotypes is largely restricted to human and agricultural systems. As such, CNVs have rarely been considered in investigations of the genomic architecture of adaptation in wild species. Here, we report on the genetic mapping of gene CNVs in white spruce, which lacks a contiguous assembly of its large genome (~20 Gb), and their relationships with adaptive phenotypic variation. We detected 3,911 gene CNVs including de novo structural variations using comparative genome hybridization on arrays (aCGH) in a large progeny set. We inferred the heterozygosity at CNV loci within parents by comparing haploid and diploid tissues and genetically mapped 82 gene CNVs. Our analysis showed that CNVs were distributed over 10 linkage groups and identified four CNV hotspots that we predict to occur in other species of the Pinaceae. Significant relationships were found between 29 of the gene CNVs and adaptive traits based on regression analyses with timings of bud set and bud flush, and height growth, suggesting a role for CNVs in climate adaptation. The importance of CNVs in adaptive evolution of white spruce was also indicated by functional gene annotations and the clustering of 31% of the mapped adaptive gene CNVs in CNV hotspots. Taken together, these results illustrate the feasibility of studying CNVs in undomesticated species and represent a major step towards a better understanding of the roles of CNVs in adaptive evolution.     

Adaptive molecular evolution of HINTW,a female-specific gene in birds

It is well established that many genes on the male-specific Y chromosome of organisms such as mammals are involved in male reproduction and may evolve rapidly because of positive selection on male reproductive traits. In contrast, very little is known about the function and evolution of W-linked genes restricted to the female genome of organisms with female heterogamety. For birds (males ZZ, females ZW), only one W-linked gene (HINTW) is sufficiently different from its Z-linked homolog to indicate a female-specific function. Here, we report that HINTW shows evidence of adaptive molecular evolution, implying strong positive selection for new functional properties in female birds. Moreover, because HINTW is expressed in the gonads of female birds just before sexual differentiation and is thus a candidate for sex determination, it suggests adaptive evolution related to female development. This provides the first example of Darwinian evolution of a gene restricted to the female genome of any organism. Given that HINTW exists in multiple copies on W, similar to some testis-specific genes amplified on mammalian Y, avian HINTW may thus potentially represent a female parallel to the organization and evolution of Y chromosome genes involved in male reproduction and development.     

Family‐assisted inference of the genetic architecture of major histocompatibility complex variation

With their direct link to individual fitness, genes of the major histocompatibility complex (MHC) are a popular system to study the evolution of adaptive genetic diversity. However, owing to the highly dynamic evolution of the MHC region, the isolation, characterization and genotyping of MHC genes remain a major challenge. While high‐throughput sequencing technologies now provide unprecedented resolution of the high allelic diversity observed at the MHC, in many species, it remains unclear (i) how alleles are distributed among MHC loci, (ii) whether MHC loci are linked or segregate independently and (iii) how much copy number variation (CNV) can be observed for MHC genes in natural populations. Here, we show that the study of allele segregation patterns within families can provide significant insights in this context. We sequenced two MHC class I (MHC‐I) loci in 1267 European barn owls (Tyto alba), including 590 offspring from 130 families using Illumina MiSeq technology. Coupled with a high per‐individual sequencing coverage (~3000×), the study of allele segregation patterns within families provided information on three aspects of the architecture of MHC‐I variation in barn owls: (i) extensive sharing of alleles among loci, (ii) strong linkage of MHC‐I loci indicating tandem architecture and (iii) the presence of CNV in the barn owl MHC‐I. We conclude that the additional information that can be gained from high‐coverage amplicon sequencing by investigating allele segregation patterns in families not only helps improving the accuracy of MHC genotyping, but also contributes towards enhanced analyses in the context of MHC evolutionary ecology.     

Different Sources of Allelic Variation Drove Repeated Color Pattern Divergence in Cichlid Fishes

The adaptive radiations of East African cichlid fish in the Great Lakes Victoria, Malawi, and Tanganyika are well known for their diversity and repeatedly evolved phenotypes. Convergent evolution of melanic horizontal stripes has been linked to a single locus harboring the gene agouti-related peptide 2 (agrp2). However, where and when the causal variants underlying this trait evolved and how they drove phenotypic divergence remained unknown. To test the alternative hypotheses of standing genetic variation versus de novo mutations (independently originating in each radiation), we searched for shared signals of genomic divergence at the agrp2 locus. Although we discovered similar signatures of differentiation at the locus level, the haplotypes associated with stripe patterns are surprisingly different. In Lake Malawi, the highest associated alleles are located within and close to the 5′ untranslated region of agrp2 and likely evolved through recent de novo mutations. In the younger Lake Victoria radiation, stripes are associated with two intronic regions overlapping with a previously reported cis-regulatory interval. The origin of these segregating haplotypes predates the Lake Victoria radiation because they are also found in more basal riverine and Lake Kivu species. This suggests that both segregating haplotypes were present as standing genetic variation at the onset of the Lake Victoria adaptive radiation with its more than 500 species and drove phenotypic divergence within the species flock. Therefore, both new (Lake Malawi) and ancient (Lake Victoria) allelic variation at the same locus fueled rapid and convergent phenotypic evolution.     

Genomics of local adaptation with gene flow

Gene flow is a fundamental evolutionary force in adaptation that is especially important to understand as humans are rapidly changing both the natural environment and natural levels of gene flow. Theory proposes a multifaceted role for gene flow in adaptation, but it focuses mainly on the disruptive effect that gene flow has on adaptation when selection is not strong enough to prevent the loss of locally adapted alleles. The role of gene flow in adaptation is now better understood due to the recent development of both genomic models of adaptive evolution and genomic techniques, which both point to the importance of genetic architecture in the origin and maintenance of adaptation with gene flow. In this review, we discuss three main topics on the genomics of adaptation with gene flow. First, we investigate selection on migration and gene flow. Second, we discuss the three potential sources of adaptive variation in relation to the role of gene flow in the origin of adaptation. Third, we explain how local adaptation is maintained despite gene flow: we provide a synthesis of recent genomic models of adaptation, discuss the genomic mechanisms and review empirical studies on the genomics of adaptation with gene flow. Despite predictions on the disruptive effect of gene flow in adaptation, an increasing number of studies show that gene flow can promote adaptation, that local adaptations can be maintained despite high gene flow, and that genetic architecture plays a fundamental role in the origin and maintenance of local adaptation with gene flow.