Hitchhiking mapping reveals a candidate genomic region for natural selection in three-spined stickleback chromosome VIII

Hitchhiking Mapping of Candidate Regions Associated with Fat Deposition in Iranian Thin and Fat Tail Sheep Breeds Suggests New Insights into Molecular Aspects of Fat Tail Selection

Simple Summary

Fatness-related traits are economically very important in sheep production and are associated with serious diseases in humans. Using a denser set of SNP markers and a variety of statistical approaches, our results were able to refine the regions associated with fat deposition and to suggest new insights into molecular aspects of fat tail selection. These results may provide a strong foundation for studying the regulation of fat deposition in sheep and do offer hope that the causal mutations and the mode of inheritance of this trait will soon be discovered by further investigation.

Abstract

The fat tail is a phenotype that divides indigenous Iranian sheep genetic resources into two major groups. The objective of the present study is to refine the map location of candidate regions associated with fat deposition, obtained via two separate whole genome scans contrasting thin and fat tail breeds, and to determine the nature of the selection occurring in these regions using a hitchhiking approach. Zel (thin tail) and Lori-Bakhtiari (fat tail) breed samples that had previously been run on the Illumina Ovine 50 k BeadChip, were genotyped with a denser set of SNPs in the three candidate regions using a Sequenom Mass ARRAY platform. Statistical tests were then performed using different and complementary methods based on either site frequency (F_ST and Median homozygosity) or haplotype (iHS and XP-EHH). The results from candidate regions on chromosome 5 and X revealed clear evidence of selection with the derived haplotypes that was consistent with selection to near fixation for the haplotypes affecting fat tail size in the fat tail breed. An analysis of the candidate region on chromosome 7 indicated that selection differentiated the beneficial alleles between breeds and homozygosity has increased in the thin tail breed which also had the ancestral haplotype. These results enabled us to confirm the signature of selection in these regions and refine the critical intervals from 113 kb, 201 kb, and 2831 kb to 28 kb, 142 kb, and 1006 kb on chromosome 5, 7, and X respectively. These regions contain several genes associated with fat metabolism or developmental processes consisting of TCF7 and PPP2CA (OAR5), PTGDR and NID2 (OAR7), AR, EBP, CACNA1F, HSD17B10,SLC35A2, BMP15, WDR13, and RBM3 (OAR X), and each of which could potentially be the actual target of selection. The study of core haplotypes alleles in our regions of interest also supported the hypothesis that the first domesticated sheep were thin tailed, and that fat tail animals were developed later. Overall, our results provide a comprehensive assessment of how and where selection has affected the patterns of variation in candidate regions associated with fat deposition in thin and fat tail sheep breeds.

Recent selection for self-compatibility in a population of Leavenworthia alabamica

The evolution of self-compatibility (SC) is the first step in the evolutionary transition in plants from outcrossing enforced by self-incompatibility (SI) to self-fertilization. In the Brassicaceae, SI is controlled by alleles of two tightly linked genes at the S-locus. Despite permitting inbreeding, mutations at the S-locus leading to SC may be selected if they provide reproductive assurance and/or gain a transmission advantage in a population when SC plants self- and outcross. Positive selection can leave a genomic signature in the regions physically linked to the focus of selection when selection has occurred recently. From an SC population of Leavenworthia alabamica with a known nonfunctional mutation at the S-locus, we collected sequence data from a ∼690 Kb region surrounding the S-locus, as well as from regions not linked to the S-locus. To test for recent positive selection acting at the S-locus, we examined polymorphism and the site-frequency spectra. Using forward simulations, we demonstrate that recent selection of the strength expected for SC at a locus formerly under balancing selection can generate patterns similar to those seen in our empirical data.

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.

Applications in the search for genomic selection signatures in fish

Selection signatures are genomic regions harboring DNA sequences functionally involved in the genetic variation of traits subject to selection. Selection signatures have been intensively studied in recent years because of their relevance to evolutionary biology and their potential association with genes that control phenotypes of interest in wild and domestic populations. Selection signature research in fish has been confined to a smaller scale, due in part to the relatively recent domestication of fish species and limited genomic resources such as molecular markers, genetic mapping, DNA sequences, and reference genomes. However, recent genomic technology advances are paving the way for more studies that may contribute to the knowledge of genomic regions underlying phenotypes of biological and productive interest in fish.

Identification and analysis of genome-wide SNPs provide insight into signatures of selection and domestication in channel catfish (Ictalurus punctatus)

Domestication and selection for important performance traits can impact the genome, which is most often reflected by reduced heterozygosity in and surrounding genes related to traits affected by selection. In this study, analysis of the genomic impact caused by domestication and artificial selection was conducted by investigating the signatures of selection using single nucleotide polymorphisms (SNPs) in channel catfish (Ictalurus punctatus). A total of 8.4 million candidate SNPs were identified by using next generation sequencing. On average, the channel catfish genome harbors one SNP per 116 bp. Approximately 6.6 million, 5.3 million, 4.9 million, 7.1 million and 6.7 million SNPs were detected in the Marion, Thompson, USDA103, Hatchery strain, and wild population, respectively. The allele frequencies of 407,861 SNPs differed significantly between the domestic and wild populations. With these SNPs, 23 genomic regions with putative selective sweeps were identified that included 11 genes. Although the function for the majority of the genes remain unknown in catfish, several genes with known function related to aquaculture performance traits were included in the regions with selective sweeps. These included hypoxia-inducible factor 1β· HIFιβ ¨ and the transporter gene ATP-binding cassette sub-family B member 5 (ABCB5). HIF1β· is important for response to hypoxia and tolerance to low oxygen levels is a critical aquaculture trait. The large numbers of SNPs identified from this study are valuable for the development of high-density SNP arrays for genetic and genomic studies of performance traits in catfish.

New susceptibility locus for obesity and dyslipidaemia on chromosome 3q22.3

BACKGROUND

The muscle Ras (MRAS) gene resides on chromosome 3q22.3 and encodes a member of the membrane-associated Ras small GTPase proteins, which function as signal transducers in multiple processes including cell growth and differentiation. Its role in cardiovascular disease is not fully understood yet. In a preliminary study in heterozygous familial hypercholesterolaemia, we identified a locus linking the early onset of coronary artery disease (CAD) to chromosome 3q.22 and elected to sequence the MRAS gene using the MegaBACE DNA analysis system. In the present study, we investigated the association of seven single-nucleotide polymorphisms (SNPs) at this locus with CAD and its dyslipidaemia-related risk traits in 4,650 Saudi angiographed individuals using TaqMan assays by the Applied Biosystems real-time Prism 7900HT Sequence Detection System.

RESULTS

Among the studied SNPs, rs6782181 (p = 0.017) and rs9818870T (p = 0.009) were associated with CAD following adjustment for sex, age and other confounding risk factors. The rs6782181_GG also conferred risk for obesity (1,764 cases vs. 2,586 controls) [1.16(1.03-1.30); p = 0.017], hypercholesterolaemia (1,686 vs. 2,744) [1.23(1.02-1.47); p = 0.019], hypertriglyceridaemia (1,155 vs. 3,496) [1.29(1.01-1.45); p = 0.043] and low high-density lipoprotein-cholesterol (lHDL-chol) levels (1,935 vs. 2,401) [1.15(1.02-1.30); p = 0.023] after adjustment. Additionally, rs253662_(CT+TT) [1.16(1.01-1.32); p = 0.030] was associated with lHDL-chol levels. Interestingly, rs253662 (p = 0.014) and rs6782181 (p = 0.019) were protective against acquiring high low-density lipoprotein-cholesterol (hLDL-chol) levels (p = 0.014), while rs1720819 showed similar effects against CAD (p < 0.0001). More importantly, a 7-mer haplotype, ACCTGAC (χ2 = 7.66; p = 0.0056), constructed from the studied SNPs, its 6-mer derivative CCTGAC (χ2 = 6.90; p = 0.0086) and several other shorter derivatives conferred risk for obesity. hLDL-chol was weakly linked to CTAA (χ2 = 3.79; p = 0.052) and CCT (χ2 = 4.32; p = 0.038), while several other haplotypes were protective against both obesity and hLDL-chol level.

CONCLUSION

Our results demonstrate that the genomic locus for the MRAS gene confers risk for CAD, obesity and dyslipidaemia and point to the possible involvement of other genes or regulatory elements at this locus, rather than changes in the M-Ras protein function, in these events.

Outlier Loci and Selection Signatures of Simple Sequence Repeats (SSRs) in Flax (Linum usitatissimum L.)

Genomic microsatellites (gSSRs) and expressed sequence tag-derived SSRs (EST-SSRs) have gained wide application for elucidating genetic diversity and population structure in plants. Both marker systems are assumed to be selectively neutral when making demographic inferences, but this assumption is rarely tested. In this study, three neutrality tests were assessed for identifying outlier loci among 150 SSRs (85 gSSRs and 65 EST-SSRs) that likely influence estimates of population structure in three differentiated flax sub-populations (F_ST = 0.19). Moreover, the utility of gSSRs, EST-SSRs, and the combined sets of SSRs was also evaluated in assessing genetic diversity and population structure in flax. Six outlier loci were identified by at least two neutrality tests showing footprints of balancing selection. After removing the outlier loci, the STRUCTURE analysis and the dendrogram topology of EST-SSRs improved. Conversely, gSSRs and combined SSRs results did not change significantly, possibly as a consequence of the higher number of neutral loci assessed. Taken together, the genetic structure analyses established the superiority of gSSRs to determine the genetic relationships among flax accessions, although the combined SSRs produced the best results. Genetic diversity parameters did not differ statistically (P > 0.05) between gSSRs and EST-SSRs, an observation partially explained by the similar number of repeat motifs. Our study provides new insights into the ability of gSSRs and EST-SSRs to measure genetic diversity and structure in flax and confirms the importance of testing for the occurrence of outlier loci to properly assess natural and breeding populations, particularly in studies considering only few loci.

Divergence hitchhiking and the spread of genomic isolation during ecological speciation-with-gene-flow

In allopatric populations, geographical separation simultaneously isolates the entire genome, allowing genetic divergence to accumulate virtually anywhere in the genome. In sympatric populations, however, the strong divergent selection required to overcome migration produces a genetic mosaic of divergent and non-divergent genomic regions. In some recent genome scans, each divergent genomic region has been interpreted as an independent incidence of migration/selection balance, such that the reduction of gene exchange is restricted to a few kilobases around each divergently selected gene. I propose an alternative mechanism, 'divergence hitchhiking' (DH), in which divergent selection can reduce gene exchange for several megabases around a gene under strong divergent selection. Not all genes/markers within a DH region are divergently selected, yet the entire region is protected to some degree from gene exchange, permitting genetic divergence from mechanisms other than divergent selection to accumulate secondarily. After contrasting DH and multilocus migration/selection balance (MM/SB), I outline a model in which genomic isolation at a given genomic location is jointly determined by DH and genome-wide effects of the progressive reduction in realized migration, then illustrate DH using data from several pairs of incipient species in the wild.

Identifying insecticide resistance genes in mosquito by combining AFLP genome scans and 454 pyrosequencing

AFLP-based genome scans are widely used to study the genetics of adaptation and to identify genomic regions potentially under selection. However, this approach usually fails to detect the actual genes or mutations targeted by selection owing to the difficulties of obtaining DNA sequences from AFLP fragments. Here, we combine classical AFLP outlier detection with 454 sequencing of AFLP fragments to obtain sequences from outlier loci. We applied this approach to the study of resistance to Bacillus thuringiensis israelensis (Bti) toxins in the dengue vector Aedes aegypti. A genome scan of Bti-resistant and Bti-susceptible A. aegypti laboratory strains was performed based on 432 AFLP markers. Fourteen outliers were detected using two different population genetic algorithms. Out of these, 11 were successfully sequenced. Three contained transposable elements (TEs) sequences, and the 10 outliers that could be mapped at a unique location in the reference genome were located on different supercontigs. One outlier was in the vicinity of a gene coding for an aminopeptidase potentially involved in Bti toxin-binding. Patterns of sequence variability of this gene showed significant deviation from neutrality in the resistant strain but not in the susceptible strain, even after taking into account the known demographic history of the selected strain. This gene is a promising candidate for future functional analysis.

Monitoring adaptive genetic responses to environmental change

Widespread environmental changes including climate change, selective harvesting and landscape alterations now greatly affect selection regimes for most organisms. How animals and plants can adapt to these altered environments via contemporary evolution is thus of strong interest. We discuss how to use genetic monitoring to study adaptive responses via repeated analysis of the same populations over time, distinguishing between phenotypic and molecular genetics approaches. After describing monitoring designs, we develop explicit criteria for demonstrating adaptive responses, which include testing for selection and establishing clear links between genetic and environmental change. We then review a few exemplary studies that explore adaptive responses to climate change in Drosophila, selective responses to hunting and fishing, and contemporary evolution in Daphnia using resurrected resting eggs. We further review a broader set of 44 studies to assess how well they meet the proposed criteria, and conclude that only 23% fulfill all criteria. Approximately half (43%) of these studies failed to rule out the alternative hypothesis of replacement by a different, better-adapted population. Likewise, 34% of the studies based on phenotypic variation did not test for selection as opposed to drift. These shortcomings can be addressed via improved experimental designs and statistical testing. We foresee monitoring of adaptive responses as a future valuable tool in conservation biology, for identifying populations unable to evolve at sufficiently high rates and for identifying possible donor populations for genetic rescue. Technological advances will further augment the realization of this potential, especially next-generation sequencing technologies that allow for monitoring at the level of whole genomes.

Characterization of a divergent chromosome region in the willow warbler Phylloscopus trochilus using avian genomic resources

Genome scans have made it possible to find outlier markers thought to have been influenced by divergent selection in almost any wild population. However, the lack of genomic information in nonmodel species often makes it difficult to associate these markers with certain genes or chromosome regions. Furthermore, the extent of linkage disequilibrium (LD) in the genome will determine the density of markers required to identify the genes under selection. In this study, we investigated a chromosome region in the willow warbler Phylloscopus trochilus surrounding a single marker previously identified in a genome scan. We first located the marker in the assembled genome of another species, the zebra finch Taeniopygia guttata, and amplified surrounding sequences in Fennoscandian willow warblers. Within an investigated chromosome region of 7.3 Mb as mapped to the zebra finch genome, we observed elevated genetic differentiation between a southern and a northern population across a 2.5-Mb interval comprising numerous coding genes. Within the southern and northern populations, higher values of LD were mostly found between SNPs within the same locus, but extended across distantly situated loci when the analyses were restricted to sampling sites showing intermediate allele frequencies of southern and northern alleles. Our study shows that cross-species genome information is a useful resource to obtain candidate sequences adjacent to outlier markers in nonmodel species.

Genomic scans detect signatures of selection along a salinity gradient in populations of the intertidal seaweed Fucus serratus on a 12 km scale

Detecting natural selection in wild populations is a central challenge in evolutionary biology and genomic scans are an important means of detecting allele frequencies that deviate from neutral expectations among marker loci. We used nine anonymous and 15 EST-linked microsatellites, 362 AFLP loci, and several neutrality tests, to identify outlier loci when comparing four populations of the seaweed Fucus serratus spaced along a 12km intertidal shore with a steep salinity gradient. Under criteria of at least two significant tests in at least two population pairs, three EST-derived and three anonymous loci revealed putative signatures of selection. Anonymous locus FsB113 was a consistent outlier when comparing least saline to fully marine sites. Locus F37 was an outlier when comparing the least saline to more saline areas, and was annotated as a polyol transporter/putative mannitol transporter - an important sugar-alcohol associated with osmoregulation by brown algae. The remaining loci could not be annotated using six different data bases. Exclusion of microsatellite outlier loci did not change either the degree or direction of differentiation among populations. In one outlier test, the number of AFLP outlier loci increased as the salinity differences between population pairs increased (up to 14); only four outliers were detected with the second test and only one was consistent with both tests. Consistency may be improved with a much more rigorous approach to replication and/or may be dependent upon the class of marker used.

The efficacy of divergence hitchhiking in generating genomic islands during ecological speciation

Genes under divergent selection flow less readily between populations than other loci. This observation has led to verbal "divergence hitchhiking" models of speciation in which decreased interpopulation gene flow surrounding loci under divergent selection can generate large regions of differentiation within the genome (genomic islands). The efficacy of this model in promoting speciation depends on the size of the region affected by divergence hitchhiking. Empirical evidence is mixed, with examples of both large and small genomic islands. To address these empirical discrepancies and to formalize the theory, we present mathematical models of divergence hitchhiking, which examine neutral differentiation around selected sites. For a single locus under selection, regions of differentiation do not extend far along a chromosome away from a selected site unless both effective population sizes and migration rates are low. When multiple loci are considered, regions of differentiation can be larger. However, with many loci under selection, genome-wide divergence occurs and genomic islands are erased. The results show that divergence hitchhiking can generate large regions of differentiation, but that the conditions under which this occurs are limited. Thus, speciation may often require multifarious selection acting on many, isolated and physically unlinked genes. How hitchhiking promotes further adaptive divergence warrants consideration.

A microsatellite-based analysis for the detection of selection on BTA1 and BTA20 in northern Eurasian cattle (Bos taurus) populations

BACKGROUND

Microsatellites surrounding functionally important candidate genes or quantitative trait loci have received attention as proxy measures of polymorphism level at the candidate loci themselves. In cattle, selection for economically important traits is a long-term strategy and it has been reported that microsatellites are linked to these important loci.

METHODS

We have investigated the variation of seven microsatellites on BTA1 (Bos taurus autosome 1) and 16 on BTA20, using bovine populations of typical production types and horn status in northern Eurasia. Genetic variability of these loci and linkage disequilibrium among these loci were compared with those of 28 microsatellites on other bovine chromosomes. Four different tests were applied to detect molecular signatures of selection.

RESULTS

No marked difference in locus variability was found between microsatellites on BTA1, BTA20 and the other chromosomes in terms of different diversity indices. Average D' values of pairwise syntenic markers (0.32 and 0.28 across BTA 1 and BTA20 respectively) were significantly (P < 0.05) higher than for non-syntenic markers (0.15). The Ewens-Watterson test, the Beaumont and Nichol's modified frequentist test and the Bayesian FST-test indicated elevated or decreased genetic differentiation, at SOD1 and AGLA17 markers respectively, deviating significantly (P < 0.05) from neutral expectations. Furthermore, lnRV, lnRH and lnRtheta' statistics were used for the pairwise population comparison tests and were significantly less variable in one population relative to the other, providing additional evidence of selection signatures for two of the 51 loci. Moreover, the three Finnish native populations showed evidence of subpopulation divergence at SOD1 and AGLA17. Our data also indicate significant intergenic linkage disequilibrium around the candidate loci and suggest that hitchhiking selection has played a role in shaping the pattern of observed linkage disequilibrium.

CONCLUSION

Hitchhiking due to tight linkage with alleles at candidate genes, e.g. the POLL gene, is a possible explanation for this pattern. The potential impact of selective breeding by man on cattle populations is discussed in the context of selection effects. Our results also suggest that a practical approach to detect loci under selection is to simultaneously apply multiple neutrality tests based on different assumptions and estimations.

A high incidence of selection on physiologically important genes in the three-spined stickleback, Gasterosteus aculeatus

Genome scan approaches to detect footprints of directional selection in the genomes of wild animal and plant populations have become popular tools to study local adaptation and speciation at the molecular level. Most studies thus far have used random molecular markers and found footprints of directional selection at, on average, 5% (range: 1-15%) of the examined loci. We focused on physiologically important genes that exhibit transcriptional responses to specific environmental or developmental conditions and assessed if these genes have been subject to directional selection and are responsible for local adaptation in the three-spined stickleback (Gasterosteus aculeatus). Using microsatellite markers located within or closely linked to (<6 kb) target genes, we investigated footprints of directional selection for 157 genes with known physiological functions in three marine and six freshwater populations. A high incidence (16.6%) of footprints of directional selection for these genes was revealed by four different outlier tests. In a subset of four populations screened with both physiologically important and random genes, footprints of directional selection were more frequent in physiologically important genes (13.4%) as compared with random genes (2.4%). In general, our findings indicate strong selective pressures on physiologically important genes, suggesting that these genes have significant functions in evolutionary adaptation to environmental heterogeneity.

Identification of local- and habitat-dependent selection: scanning functionally important genes in nine-spined sticklebacks (Pungitius pungitius)

Understanding the selective forces promoting adaptive population divergence is a central issue in evolutionary biology. The role of environmental salinity in driving adaptation and evolution in aquatic organisms is still poorly understood. We investigated the relative impacts of habitat type (cf. saltwater vs. freshwater) and geographic area in shaping adaptive population divergence, as well as genes responsible for adaptation to different salinities in nine-spined sticklebacks (Pungitius pungitius). To this end, we employed a hitchhiking mapping approach with 111 microsatellite loci and one insertion/deletion locus including 63 loci situated within or close to genes with important physiological functions such as osmoregulation, growth, and thermal response. Using three pairs of marine and freshwater populations from different geographic areas, we identified several loci showing consistent evidence of being under directional selection in different outlier tests. Analyses of molecular variance at the loci under selection indicated that geographic area rather than habitat type has been acting as a central force in shaping adaptive population divergence. Nevertheless, both outlier tests and a spatial analysis method indicated that two loci (growth hormone receptor 2 and DEAD box polypeptide 56) are involved in adaptation to different habitats, implying that environmental salinity has been affecting them as a selective force. These loci are promising candidates for further investigations focusing on the molecular mechanisms of adaptation to marine and freshwater environments.

Rapid selection against inbreeding in a wild population of a rare frog

Populations that are small and isolated can be threatened through loss of fitness due to inbreeding. Nevertheless, an increased frequency of recessive homozygotes could increase the efficiency of selection against deleterious mutants, thus reducing inbreeding depression. In wild populations, observations of evolutionary changes determined by selection against inbreeding are few. We used microsatellite DNA markers to compare the genetic features of tadpoles immediately after hatch with those of metamorphosing froglets belonging to the same cohort in a small, isolated population of the threatened frog Rana latastei. Within a generation, the inbreeding coefficient (F_IS) decreased: at hatch, F_IS was significantly >0, whereas F_IS was <0 after metamorphosis. Furthermore, heterozygosity increased and allelic frequencies changed over time, resulting in the loss of genotypes at metamorphosis that were present in hatchlings. One microsatellite locus exhibited atypically large F_ST values, suggesting it might be linked to a locus under selection. These results support the hypothesis that strong selection against the most inbred genotypes occurred among early life-history stages in our population. Selective forces can promote changes that can affect population dynamics and should be considered in conservation planning.

Nucleotide diversity and linkage disequilibrium in balsam poplar (Populus balsamifera)

*Current perceptions that poplars have high levels of nucleotide variation, large effective population sizes, and rapid decay of linkage disequilibrium are based primarily on studies from one poplar species, Populus tremula. *We analysed 590 gene fragments (average length 565 bp) from each of 15 individuals from different populations from throughout the range of Populus balsamifera. *Nucleotide diversity (theta(total) = 0.0028, pi = 0.0027) was low compared with other trees and model agricultural systems. Patterns of nucleotide diversity and site frequency spectra were consistent with purifying selection on replacement and intron sites. When averaged across all loci we found no evidence for decay of linkage disequilibrium across 750 bp, consistent with the low estimates of the scaled recombination parameter, rho = 0.0092. *Compared with P. tremula, a well studied congener with a similar distribution, P. balsamifera has low diversity and low effective recombination, both of which indicate a lower effective population size in P. balsamifera. Patterns of diversity and linkage indicate that there is considerable variation in population genomic patterns among poplar species and unlike P. tremula, association mapping techniques in balsam poplar should consider sampling single nucleotide polymorphisms (SNPs) at well-spaced intervals.

Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags

Next-generation sequencing technology provides novel opportunities for gathering genome-scale sequence data in natural populations, laying the empirical foundation for the evolving field of population genomics. Here we conducted a genome scan of nucleotide diversity and differentiation in natural populations of threespine stickleback (Gasterosteus aculeatus). We used Illumina-sequenced RAD tags to identify and type over 45,000 single nucleotide polymorphisms (SNPs) in each of 100 individuals from two oceanic and three freshwater populations. Overall estimates of genetic diversity and differentiation among populations confirm the biogeographic hypothesis that large panmictic oceanic populations have repeatedly given rise to phenotypically divergent freshwater populations. Genomic regions exhibiting signatures of both balancing and divergent selection were remarkably consistent across multiple, independently derived populations, indicating that replicate parallel phenotypic evolution in stickleback may be occurring through extensive, parallel genetic evolution at a genome-wide scale. Some of these genomic regions co-localize with previously identified QTL for stickleback phenotypic variation identified using laboratory mapping crosses. In addition, we have identified several novel regions showing parallel differentiation across independent populations. Annotation of these regions revealed numerous genes that are candidates for stickleback phenotypic evolution and will form the basis of future genetic analyses in this and other organisms. This study represents the first high-density SNP-based genome scan of genetic diversity and differentiation for populations of threespine stickleback in the wild. These data illustrate the complementary nature of laboratory crosses and population genomic scans by confirming the adaptive significance of previously identified genomic regions, elucidating the particular evolutionary and demographic history of such regions in natural populations, and identifying new genomic regions and candidate genes of evolutionary significance.

Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags

Next-generation sequencing technology provides novel opportunities for gathering genome-scale sequence data in natural populations, laying the empirical foundation for the evolving field of population genomics. Here we conducted a genome scan of nucleotide diversity and differentiation in natural populations of threespine stickleback (Gasterosteus aculeatus). We used Illumina-sequenced RAD tags to identify and type over 45,000 single nucleotide polymorphisms (SNPs) in each of 100 individuals from two oceanic and three freshwater populations. Overall estimates of genetic diversity and differentiation among populations confirm the biogeographic hypothesis that large panmictic oceanic populations have repeatedly given rise to phenotypically divergent freshwater populations. Genomic regions exhibiting signatures of both balancing and divergent selection were remarkably consistent across multiple, independently derived populations, indicating that replicate parallel phenotypic evolution in stickleback may be occurring through extensive, parallel genetic evolution at a genome-wide scale. Some of these genomic regions co-localize with previously identified QTL for stickleback phenotypic variation identified using laboratory mapping crosses. In addition, we have identified several novel regions showing parallel differentiation across independent populations. Annotation of these regions revealed numerous genes that are candidates for stickleback phenotypic evolution and will form the basis of future genetic analyses in this and other organisms. This study represents the first high-density SNP-based genome scan of genetic diversity and differentiation for populations of threespine stickleback in the wild. These data illustrate the complementary nature of laboratory crosses and population genomic scans by confirming the adaptive significance of previously identified genomic regions, elucidating the particular evolutionary and demographic history of such regions in natural populations, and identifying new genomic regions and candidate genes of evolutionary significance.

Genomic hotspots for adaptation: the population genetics of Müllerian mimicry in Heliconius erato

Wing pattern evolution in Heliconius butterflies provides some of the most striking examples of adaptation by natural selection. The genes controlling pattern variation are classic examples of Mendelian loci of large effect, where allelic variation causes large and discrete phenotypic changes and is responsible for both convergent and highly divergent wing pattern evolution across the genus. We characterize nucleotide variation, genotype-by-phenotype associations, linkage disequilibrium (LD), and candidate gene expression patterns across two unlinked genomic intervals that control yellow and red wing pattern variation among mimetic forms of Heliconius erato. Despite very strong natural selection on color pattern, we see neither a strong reduction in genetic diversity nor evidence for extended LD across either patterning interval. This observation highlights the extent that recombination can erase the signature of selection in natural populations and is consistent with the hypothesis that either the adaptive radiation or the alleles controlling it are quite old. However, across both patterning intervals we identified SNPs clustered in several coding regions that were strongly associated with color pattern phenotype. Interestingly, coding regions with associated SNPs were widely separated, suggesting that color pattern alleles may be composed of multiple functional sites, conforming to previous descriptions of these loci as "supergenes." Examination of gene expression levels of genes flanking these regions in both H. erato and its co-mimic, H. melpomene, implicate a gene with high sequence similarity to a kinesin as playing a key role in modulating pattern and provides convincing evidence for parallel changes in gene regulation across co-mimetic lineages. The complex genetic architecture at these color pattern loci stands in marked contrast to the single casual mutations often identified in genetic studies of adaptation, but may be more indicative of the type of genetic changes responsible for much of the adaptive variation found in natural populations.

Identifying footprints of directional and balancing selection in marine and freshwater three-spined stickleback (Gasterosteus aculeatus) populations

Natural selection is expected to leave an imprint on the neutral polymorphisms at the adjacent genomic regions of a selected gene. While directional selection tends to reduce within-population genetic diversity and increase among-population differentiation, the reverse is expected under balancing selection. To identify targets of natural selection in the three-spined stickleback (Gasterosteus aculeatus) genome, 103 microsatellite and two indel markers including expressed sequence tags (EST) and quantitative trait loci (QTL)-associated loci, were genotyped in four freshwater and three marine populations. The results indicated that a high proportion of loci (14.7%) might be affected by balancing selection and a lower proportion (2.8%) by directional selection. The strongest signatures of directional selection were detected in a microsatellite locus and two indel markers located in the intronic regions of the Eda-gene coding for the number of lateral plates. Yet, other microsatellite loci previously found to be informative in QTL-mapping studies revealed no signatures of selection. Two novel microsatellite loci (Stn12 and Stn90) located in chromosomes I and VIII, respectively, showed signals of directional selection and might be linked to genomic regions containing gene(s) important for adaptive divergence. Although the coverage of the total genomic content was relatively low, the predominance of balancing selection signals is in agreement with the contention that balancing, rather than directional selection is the predominant mode of selection in the wild.

Selection and sticklebacks

Over the last decade, there has been increasing circumstantial evidence for the action of natural selection in the genome, arising largely from molecular genetic surveys of large numbers of markers. In nonmodel organisms without densely mapped markers, a frequently used method is to identify loci that have unusually high or low levels of genetic differentiation, or low genetic diversity relative to other populations. The paper by Mäkinen et al. (2008a) in this issue of Molecular Ecology reports the results of a survey of microsatellite allele frequencies at more than 100 loci in seven populations of the three-spined stickleback (Gasterosteus aculeatus). They show that a microsatellite locus and two indel markers located within the intron of the Eda gene, known to control the number of lateral plates in the stickleback (Fig. 1), tend to be much more highly genetically differentiated than other loci, a finding that is consistent with the action of local selection. They identify a further two independent candidates for local selection, and, most intriguingly, they further suggest that up to 15% of their loci may provide evidence of balancing selection.