Rapid SNP Discovery and Genetic Mapping Using Sequenced RAD Markers

Single nucleotide polymorphism (SNP) discovery and genotyping are essential to genetic mapping. There remains a need for a simple, inexpensive platform that allows high-density SNP discovery and genotyping in large populations. Here we describe the sequencing of restriction-site associated DNA (RAD) tags, which identified more than 13,000 SNPs, and mapped three traits in two model organisms, using less than half the capacity of one Illumina sequencing run. We demonstrated that different marker densities can be attained by choice of restriction enzyme. Furthermore, we developed a barcoding system for sample multiplexing and fine mapped the genetic basis of lateral plate armor loss in threespine stickleback by identifying recombinant breakpoints in F2 individuals. Barcoding also facilitated mapping of a second trait, a reduction of pelvic structure, by in silico re-sorting of individuals. To further demonstrate the ease of the RAD sequencing approach we identified polymorphic markers and mapped an induced mutation in Neurospora crassa. Sequencing of RAD markers is an integrated platform for SNP discovery and genotyping. This approach should be widely applicable to genetic mapping in a variety of organisms

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

Duplication and Diversification of the Hypoxia-Inducible IGFBP-1 Gene in Zebrafish

Gene duplication is the primary force of new gene evolution. Deciphering whether a pair of duplicated genes has evolved divergent functions is often challenging. The zebrafish is uniquely positioned to provide insight into the process of functional gene evolution due to its amenability to genetic and experimental manipulation and because it possess a large number of duplicated genes.We report the identification and characterization of two hypoxia-inducible genes in zebrafish that are co-ortholgs of human IGF binding protein-1 (IGFBP-1). IGFBP-1 is a secreted protein that binds to IGF and modulates IGF actions in somatic growth, development, and aging. Like their human and mouse counterparts, in adult zebrafish igfbp-1a and igfbp-1b are exclusively expressed in the liver. During embryogenesis, the two genes are expressed in overlapping spatial domains but with distinct temporal patterns. While zebrafish IGFBP-1a mRNA was easily detected throughout embryogenesis, IGFBP-1b mRNA was detectable only in advanced stages. Hypoxia induces igfbp-1a expression in early embryogenesis, but induces the igfbp-1b expression later in embryogenesis. Both IGFBP-1a and -b are capable of IGF binding, but IGFBP-1b has much lower affinities for IGF-I and -II because of greater dissociation rates. Overexpression of IGFBP-1a and -1b in zebrafish embryos caused significant decreases in growth and developmental rates. When tested in cultured zebrafish embryonic cells, IGFBP-1a and -1b both inhibited IGF-1-induced cell proliferation but the activity of IGFBP-1b was significantly weaker.These results indicate subfunction partitioning of the duplicated IGFBP-1 genes at the levels of gene expression, physiological regulation, protein structure, and biological actions. The duplicated IGFBP-1 may provide additional flexibility in fine-tuning IGF signaling activities under hypoxia and other catabolic conditions

Contrasting Patterns of Sequence Evolution at the Functionally Redundant bric à brac Paralogs in Drosophila melanogaster

Genes with overlapping expression and function may gradually diverge despite retaining some common functions. To test whether such genes show distinct patterns of molecular evolution within species, we examined sequence variation at the bric à brac (bab) locus of Drosophila melanogaster. This locus is composed of two anciently duplicated paralogs, bab1 and bab2, which are involved in patterning the adult abdomen, legs, and ovaries. We have sequenced the 148 kb genomic region spanning the bab1 and bab2 genes from 94 inbred lines of D. melanogaster sampled from a single location. Two non-coding regions, one in each paralog, appear to be under selection. The strongest evidence of directional selection is found in a region of bab2 that has no known functional role. The other region is located in the bab1 paralog and is known to contain a cis-regulatory element that controls sex-specific abdominal pigmentation. The coding region of bab1 appears to be under stronger functional constraint than the bab2 coding sequences. Thus, the two paralogs are evolving under different selective regimes in the same natural population, illuminating the different evolutionary trajectories of partially redundant duplicate genes

The genomic basis of adaptive evolution in threespine sticklebacks

Marine stickleback fish have colonized and adapted to thousands of streams and lakes formed since the last ice age, providing an exceptional opportunity to characterize genomic mechanisms underlying repeated ecological adaptation in nature. Here we develop a high-quality reference genome assembly for threespine sticklebacks. By sequencing the genomes of twenty additional individuals from a global set of marine and freshwater populations, we identify a genome-wide set of loci that are consistently associated with marine–freshwater divergence. Our results indicate that reuse of globally shared standing genetic variation, including chromosomal inversions, has an important role in repeated evolution of distinct marine and freshwater sticklebacks, and in the maintenance of divergent ecotypes during early stages of reproductive isolation. Both coding and regulatory changes occur in the set of loci underlying marine–freshwater evolution, but regulatory changes appear to predominate in this well known example of repeated adaptive evolution in nature.National Human Genome Research Institute (U.S.)National Human Genome Research Institute (U.S.) (NHGRI CEGS Grant P50-HG002568

Habitat-Specific Morphological Variation among Threespine Sticklebacks (Gasterosteus aculeatus) within a Drainage Basin

Habitat-specific morphological variation, often corresponding to resource specialization, is well documented in freshwater fishes. In this study we used landmark based morphometric analyses to investigate morphological variation among threespine sticklebacks (Gasterosteus aculeatus L.) from four interconnected habitat types within a single lowland drainage basin in eastern England. These included the upper and lower reaches of the river, the estuary, a connected ditch network and a coastal salt marsh. We found significant habitat-specific differences in morphology, with three axes of variation describing differences in orbit diameter, body depth, caudal peduncle shape and pectoral fin positioning as well as variation in relative dorsal and pelvic spine size. Interestingly, the ditch system, an artificial and heavily managed habitat, is populated by sticklebacks with a characteristic morphology, suggesting that human management of habitats can in some circumstances lead to morphological variation among the animals that inhabit them. We discuss the mechanisms that conceivably underlie the observed morphological variation and the further work necessary to identify them. Finally, we consider the implications of habitat-specific body shape variation for the behavioural ecology of this ecologically generalist species

Divergence in Sex Steroid Hormone Signaling between Sympatric Species of Japanese Threespine Stickleback

Sex steroids mediate the expression of sexually dimorphic or sex-specific traits that are important both for mate choice within species and for behavioral isolation between species. We investigated divergence in sex steroid signaling between two sympatric species of threespine stickleback (Gasterosteus aculeatus): the Japan Sea form and the Pacific Ocean form. These sympatric forms diverge in both male display traits and female mate choice behaviors, which together contribute to asymmetric behavioral isolation in sympatry. Here, we found that plasma levels of testosterone and 17β-estradiol differed between spawning females of the two sympatric forms. Transcript levels of follicle-stimulating hormone-β (FSHβ) gene were also higher in the pituitary gland of spawning Japan Sea females than in the pituitary gland of spawning Pacific Ocean females. By contrast, none of the sex steroids examined were significantly different between nesting males of the two forms. However, combining the plasma sex steroid data with testis transcriptome data suggested that the efficiency of the conversion of testosterone into 11-ketotestosterone has likely diverged between forms. Within forms, plasma testosterone levels in males were significantly correlated with male body size, a trait important for female mate choice in the two sympatric species. These results demonstrate that substantial divergence in sex steroid signaling can occur between incipient sympatric species. We suggest that investigation of the genetic and ecological mechanisms underlying divergence in hormonal signaling between incipient sympatric species will provide a better understanding of the mechanisms of speciation in animals

Where the Lake Meets the Sea: Strong Reproductive Isolation Is Associated with Adaptive Divergence between Lake Resident and Anadromous Three-Spined Sticklebacks

Contact zones between divergent forms of the same species are often characterised by high levels of phenotypic diversity over small geographic distances. What processes are involved in generating such high phenotypic diversity? One possibility is that introgression and recombination between divergent forms in contact zones results in greater phenotypic and genetic polymorphism. Alternatively, strong reproductive isolation between forms may maintain distinct phenotypes, preventing homogenisation by gene flow. Contact zones between divergent freshwater-resident and anadromous stickleback (Gasterosteus aculeatus L.) forms are numerous and common throughout the species distribution, offering an opportunity to examine these contrasting hypotheses in greater detail. This study reports on an interesting new contact zone located in a tidally influenced lake catchment in western Ireland, characterised by high polymorphism for lateral plate phenotypes. Using neutral and QTL-linked microsatellite markers, we tested whether the high diversity observed in this contact zone arose as a result of introgression or reproductive isolation between divergent forms: we found strong support for the latter hypothesis. Three phenotypic and genetic clusters were identified, consistent with two divergent resident forms and a distinct anadromous completely plated population that migrates in and out of the system. Given the strong neutral differentiation detected between all three morphotypes (mean F-ST = 0.12), we hypothesised that divergent selection between forms maintains reproductive isolation. We found a correlation between neutral genetic and adaptive genetic differentiation that support this. While strong associations between QTL linked markers and phenotypes were also observed in this wild population, our results support the suggestion that such associations may be more complex in some Atlantic populations compared to those in the Pacific. These findings provide an important foundation for future work investigating the dynamics of gene flow and adaptive divergence in this newly discovered stickleback contact zone