Restriction Site Tiling Analysis: accurate discovery and quantitative genotyping of genome-wide polymorphisms using nucleotide arrays

Whole-Genome Sequencing Reveals That Regulatory and Low Pleiotropy Variants Underlie Local Adaptation to Environmental Variability in Purple Sea Urchins

AbstractOrganisms experience environments that vary across both space and time. Such environmental heterogeneity shapes standing genetic variation and may influence species' capacity to adapt to rapid environmental change. However, we know little about the kind of genetic variation that is involved in local adaptation to environmental variability. To address this gap, we sequenced the whole genomes of 140 purple sea urchins (Strongylocentrotus purpuratus) from seven populations that vary in their degree of pH variability. Despite no evidence of global population structure, we found a suite of single-nucleotide polymorphisms (SNPs) tightly correlated with local pH variability (outlier SNPs), which were overrepresented in regions putatively involved in gene regulation (long noncoding RNA and enhancers), supporting the idea that variation in regulatory regions is important for local adaptation to variability. In addition, outliers in genes were found to be (i) enriched for biomineralization and ion homeostasis functions related to low pH response, (ii) less central to the protein-protein interaction network, and (iii) underrepresented among genes highly expressed during early development. Taken together, these results suggest that loci that underlie local adaptation to pH variability in purple sea urchins fall in regions with potentially low pleiotropic effects (based on analyses involving regulatory regions, network centrality, and expression time) involved in low pH response (based on functional enrichment).

Rare genetic variation and balanced polymorphisms are important for survival in global change conditions

Standing genetic variation is important for population persistence in extreme environmental conditions. While some species may have the capacity to adapt to predicted average future global change conditions, the ability to survive extreme events is largely unknown. We used single-generation selection experiments on hundreds of thousands of Strongylocentrotus purpuratus sea urchin larvae generated from wild-caught adults to identify adaptive genetic variation responsive to moderate (pH 8.0) and extreme (pH 7.5) low-pH conditions. Sequencing genomic DNA from pools of larvae, we identified consistent changes in allele frequencies across replicate cultures for each pH condition and observed increased linkage disequilibrium around selected loci, revealing selection on recombined standing genetic variation. We found that loci responding uniquely to either selection regime were at low starting allele frequencies while variants that responded to both pH conditions (11.6% of selected variants) started at high frequencies. Loci under selection performed functions related to energetics, pH tolerance, cell growth and actin/cytoskeleton dynamics. These results highlight that persistence in future conditions will require two classes of genetic variation: common, pH-responsive variants maintained by balancing selection in a heterogeneous environment, and rare variants, particularly for extreme conditions, that must be maintained by large population sizes.

Population structure, genetic connectivity, and adaptation in the Olympia oyster (Ostrea lurida) along the west coast of North America

Effective management of threatened and exploited species requires an understanding of both the genetic connectivity among populations and local adaptation. The Olympia oyster (Ostrea lurida), patchily distributed from Baja California to the central coast of Canada, has a long history of population declines due to anthropogenic stressors. For such coastal marine species, population structure could follow a continuous isolation-by-distance model, contain regional blocks of genetic similarity separated by barriers to gene flow, or be consistent with a null model of no population structure. To distinguish between these hypotheses in O. lurida, 13,424 single nucleotide polymorphisms (SNPs) were used to characterize rangewide population structure, genetic connectivity, and adaptive divergence. Samples were collected across the species range on the west coast of North America, from southern California to Vancouver Island. A conservative approach for detecting putative loci under selection identified 235 SNPs across 129 GBS loci, which were functionally annotated and analyzed separately from the remaining neutral loci. While strong population structure was observed on a regional scale in both neutral and outlier markers, neutral markers had greater power to detect fine-scale structure. Geographic regions of reduced gene flow aligned with known marine biogeographic barriers, such as Cape Mendocino, Monterey Bay, and the currents around Cape Flattery. The outlier loci identified as under putative selection included genes involved in developmental regulation, sensory information processing, energy metabolism, immune response, and muscle contraction. These loci are excellent candidates for future research and may provide targets for genetic monitoring programs. Beyond specific applications for restoration and management of the Olympia oyster, this study lends to the growing body of evidence for both population structure and adaptive differentiation across a range of marine species exhibiting the potential for panmixia. Computational notebooks are available to facilitate reproducibility and future open-sourced research on the population structure of O. lurida.

Genome-wide signals of positive selection in strongylocentrotid sea urchins

Background

Comparative genomics studies investigating the signals of positive selection among groups of closely related species are still rare and limited in taxonomic breadth. Such studies show great promise in advancing our knowledge about the proportion and the identity of genes experiencing diversifying selection. However, methodological challenges have led to high levels of false positives in past studies. Here, we use the well-annotated genome of the purple sea urchin, Strongylocentrotus purpuratus, as a reference to investigate the signals of positive selection at 6520 single-copy orthologs from nine sea urchin species belonging to the family Strongylocentrotidae paying careful attention to minimizing false positives.

Results

We identified 1008 (15.5%) candidate positive selection genes (PSGs). Tests for positive selection along the nine terminal branches of the phylogeny identified 824 genes that showed lineage-specific adaptive diversification (1.67% of branch-sites tests performed). Positively selected codons were not enriched at exon borders or near regions containing missing data, suggesting a limited contribution of false positives caused by alignment or annotation errors. Alignments were validated at 10 loci with re-sequencing using Sanger methods. No differences were observed in the rates of synonymous substitution (d_S), GC content, and codon bias between the candidate PSGs and those not showing positive selection. However, the candidate PSGs had 68% higher rates of nonsynonymous substitution (d_N) and 33% lower levels of heterozygosity, consistent with selective sweeps and opposite to that expected by a relaxation of selective constraint. Although positive selection was identified at reproductive proteins and innate immunity genes, the strongest signals of adaptive diversification were observed at extracellular matrix proteins, cell adhesion molecules, membrane receptors, and ion channels. Many candidate PSGs have been widely implicated as targets of pathogen binding, inactivation, mimicry, or exploitation in other groups (notably mammals).

Conclusions

Our study confirmed the widespread action of positive selection across sea urchin genomes and allowed us to reject the possibility that annotation and alignment errors (including paralogs) were responsible for creating false signals of adaptive molecular divergence. The candidate PSGs identified in our study represent promising targets for future research into the selective agents responsible for their adaptive diversification and their contribution to speciation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3944-7) contains supplementary material, which is available to authorized users.

Sensitivity of sea urchin fertilization to pH varies across a natural pH mosaic

In the coastal ocean, temporal fluctuations in pH vary dramatically across biogeographic ranges. How such spatial differences in pH variability regimes might shape ocean acidification resistance in marine species remains unknown. We assessed the pH sensitivity of the sea urchin Strongylocentrotus purpuratus in the context of ocean pH variability. Using unique male-female pairs, originating from three sites with similar mean pH but different variability and frequency of low pH (pHT ≤ 7.8) exposures, fertilization was tested across a range of pH (pHT 7.61-8.03) and sperm concentrations. High fertilization success was maintained at low pH via a slight right shift in the fertilization function across sperm concentration. This pH effect differed by site. Urchins from the site with the narrowest pH variability regime exhibited the greatest pH sensitivity. At this site, mechanistic fertilization dynamics models support a decrease in sperm-egg interaction rate with decreasing pH. The site differences in pH sensitivity build upon recent evidence of local pH adaptation in S. purpuratus and highlight the need to incorporate environmental variability in the study of global change biology.

Breaking RAD: an evaluation of the utility of restriction site-associated DNA sequencing for genome scans of adaptation

Understanding how and why populations evolve is of fundamental importance to molecular ecology. Restriction site-associated DNA sequencing (RADseq), a popular reduced representation method, has ushered in a new era of genome-scale research for assessing population structure, hybridization, demographic history, phylogeography and migration. RADseq has also been widely used to conduct genome scans to detect loci involved in adaptive divergence among natural populations. Here, we examine the capacity of those RADseq-based genome scan studies to detect loci involved in local adaptation. To understand what proportion of the genome is missed by RADseq studies, we developed a simple model using different numbers of RAD-tags, genome sizes and extents of linkage disequilibrium (length of haplotype blocks). Under the best-case modelling scenario, we found that RADseq using six- or eight-base pair cutting restriction enzymes would fail to sample many regions of the genome, especially for species with short linkage disequilibrium. We then surveyed recent studies that have used RADseq for genome scans and found that the median density of markers across these studies was 4.08 RAD-tag markers per megabase (one marker per 245 kb). The length of linkage disequilibrium for many species is one to three orders of magnitude less than density of the typical recent RADseq study. Thus, we conclude that genome scans based on RADseq data alone, while useful for studies of neutral genetic variation and genetic population structure, will likely miss many loci under selection in studies of local adaptation.

AFSM sequencing approach: a simple and rapid method for genome-wide SNP and methylation site discovery and genetic mapping

We describe methods for the assessment of amplified-fragment single nucleotide polymorphism and methylation (AFSM) sites using a quick and simple molecular marker-assisted breeding strategy based on the use of two restriction enzyme pairs (EcoRI-MspI and EcoRI-HpaII) and a next-generation sequencing platform. Two sets of 85 adapter pairs were developed to concurrently identify SNPs, indels and methylation sites for 85 lines of cassava population in this study. In addition to SNPs and indels, the simplicity of the AFSM protocol makes it particularly suitable for high-throughput full methylation and hemi-methylation analyses. To further demonstrate the ease of this approach, a cassava genetic linkage map was constructed. This approach should be widely applicable for genetic mapping in a variety of organisms and will improve the application of crop genomics in assisted breeding.

The impact of gene expression variation on the robustness and evolvability of a developmental gene regulatory network

Regulatory interactions buffer development against genetic and environmental perturbations, but adaptation requires phenotypes to change. We investigated the relationship between robustness and evolvability within the gene regulatory network underlying development of the larval skeleton in the sea urchin Strongylocentrotus purpuratus. We find extensive variation in gene expression in this network throughout development in a natural population, some of which has a heritable genetic basis. Switch-like regulatory interactions predominate during early development, buffer expression variation, and may promote the accumulation of cryptic genetic variation affecting early stages. Regulatory interactions during later development are typically more sensitive (linear), allowing variation in expression to affect downstream target genes. Variation in skeletal morphology is associated primarily with expression variation of a few, primarily structural, genes at terminal positions within the network. These results indicate that the position and properties of gene interactions within a network can have important evolutionary consequences independent of their immediate regulatory role.

Signs of adaptation to local pH conditions across an environmental mosaic in the California Current Ecosystem

Little is known about the potential for rapid evolution in natural populations in response to the high rate of contemporary climatic change. Organisms that have evolved in environments that experience high variability across space and time are of particular interest as they may harbor genetic variation that can facilitate evolutionary response to changing conditions. Here we review what is known about genetic capacity for adaptation in the purple sea urchin, Strongylocentrotus purpuratus, a species that has evolved in the upwelling ecosystem of the Northeast Pacific Ocean. We also present new results testing for adaptation to local pH conditions in six populations from Oregon to southern California. We integrate data on 19,493 genetic polymorphisms with data on local pH conditions. We find correlations between allele frequency and rank average time spent at pH <7.8 in 318 single-nucleotide polymorphisms in 275 genes. Two of the genes most correlated with local pH are a protein associated with the cytoskeleton and a proton pump, with functional roles in maintenance of cell volume and with internal regulation of pH, respectively. Across all loci tested, high correlations with local pH were concentrated in genes related to transport of ions, biomineralization, lipid metabolism, and cell-cell adhesion, functional pathways important for maintaining homeostasis at low pH. We identify a set of seven genes as top candidates for rapid evolutionary response to acidification of the ocean. In these genes, the putative low-pH-adapted allele, based on allele frequencies in natural populations, rapidly increases in frequency in purple sea urchin larvae raised at low pH. We also found that populations from localities with high pH show a greater change in allele frequency toward putative low-pH-adapted alleles under experimental acidification, compared with low-pH populations, suggesting that both natural and artificial selection favor the same alleles for response to low pH. These results illustrate that purple sea urchins may be adapted to local pH and suggest that this species may possess the genetic capacity for rapid evolution in response to acidification. This adaptive capacity likely comes from standing genetic variation maintained in nature by balancing selection across the spatial and temporal environmental mosaic that characterizes the California Current Ecosystem.

Evolutionary change during experimental ocean acidification

Rising atmospheric carbon dioxide (CO2) conditions are driving unprecedented changes in seawater chemistry, resulting in reduced pH and carbonate ion concentrations in the Earth's oceans. This ocean acidification has negative but variable impacts on individual performance in many marine species. However, little is known about the adaptive capacity of species to respond to an acidified ocean, and, as a result, predictions regarding future ecosystem responses remain incomplete. Here we demonstrate that ocean acidification generates striking patterns of genome-wide selection in purple sea urchins (Strongylocentrotus purpuratus) cultured under different CO2 levels. We examined genetic change at 19,493 loci in larvae from seven adult populations cultured under realistic future CO2 levels. Although larval development and morphology showed little response to elevated CO2, we found substantial allelic change in 40 functional classes of proteins involving hundreds of loci. Pronounced genetic changes, including excess amino acid replacements, were detected in all populations and occurred in genes for biomineralization, lipid metabolism, and ion homeostasis--gene classes that build skeletons and interact in pH regulation. Such genetic change represents a neglected and important impact of ocean acidification that may influence populations that show few outward signs of response to acidification. Our results demonstrate the capacity for rapid evolution in the face of ocean acidification and show that standing genetic variation could be a reservoir of resilience to climate change in this coastal upwelling ecosystem. However, effective response to strong natural selection demands large population sizes and may be limited in species impacted by other environmental stressors.

Going where traditional markers have not gone before: utility of and promise for RAD sequencing in marine invertebrate phylogeography and population genomics

Characterization of large numbers of single-nucleotide polymorphisms (SNPs) throughout a genome has the power to refine the understanding of population demographic history and to identify genomic regions under selection in natural populations. To this end, population genomic approaches that harness the power of next-generation sequencing to understand the ecology and evolution of marine invertebrates represent a boon to test long-standing questions in marine biology and conservation. We employed restriction-site-associated DNA sequencing (RAD-seq) to identify SNPs in natural populations of the sea anemone Nematostella vectensis, an emerging cnidarian model with a broad geographic range in estuarine habitats in North and South America, and portions of England. We identified hundreds of SNP-containing tags in thousands of RAD loci from 30 barcoded individuals inhabiting four locations from Nova Scotia to South Carolina. Population genomic analyses using high-confidence SNPs resulted in a highly-resolved phylogeography, a result not achieved in previous studies using traditional markers. Plots of locus-specific FST against heterozygosity suggest that a majority of polymorphic sites are neutral, with a smaller proportion suggesting evidence for balancing selection. Loci inferred to be under balancing selection were mapped to the genome, where 90% were located in gene bodies, indicating potential targets of selection. The results from analyses with and without a reference genome supported similar conclusions, further highlighting RAD-seq as a method that can be efficiently applied to species lacking existing genomic resources. We discuss the utility of RAD-seq approaches in burgeoning Nematostella research as well as in other cnidarian species, particularly corals and jellyfishes, to determine phylogeographic relationships of populations and identify regions of the genome undergoing selection.

Development under elevated pCO2 conditions does not affect lipid utilization and protein content in early life-history stages of the purple sea urchin, Strongylocentrotus purpuratus

Ocean acidification (OA) is expected to have a major impact on marine species, particularly during early life-history stages. These effects appear to be species-specific and may include reduced survival, altered morphology, and depressed metabolism. However, less information is available regarding the bioenergetics of development under elevated CO(2) conditions. We examined the biochemical and morphological responses of Strongylocentrotus purpuratus during early development under ecologically relevant levels of pCO(2) (365, 1030, and 1450 μatm) that may occur during intense upwelling events. The principal findings of this study were (1) lipid utilization rates and protein content in S. purpuratus did not vary with pCO(2); (2) larval growth was reduced at elevated pCO(2) despite similar rates of energy utilization; and (3) relationships between egg phospholipid content and larval length were found under control but not high pCO(2) conditions. These results suggest that this species may either prioritize endogenous energy toward development and physiological function at the expense of growth, or that reduced larval length may be strictly due to higher costs of growth under OA conditions. This study highlights the need to further expand our knowledge of the physiological mechanisms involved in OA response in order to better understand how present populations may respond to global environmental change.

Genome-wide polymorphisms show unexpected targets of natural selection

Natural selection can act on all the expressed genes of an individual, leaving signatures of genetic differentiation or diversity at many loci across the genome. New power to assay these genome-wide effects of selection comes from associating multi-locus patterns of polymorphism with gene expression and function. Here, we performed one of the first genome-wide surveys in a marine species, comparing purple sea urchins, Strongylocentrotus purpuratus, from two distant locations along the species' wide latitudinal range. We examined 9112 polymorphic loci from upstream non-coding and coding regions of genes for signatures of selection with respect to gene function and tissue- and ontogenetic gene expression. We found that genetic differentiation (F(ST)) varied significantly across functional gene classes. The strongest enrichment occurred in the upstream regions of E3 ligase genes, enzymes known to regulate protein abundance during development and environmental stress. We found enrichment for high heterozygosity in genes directly involved in immune response, particularly NALP genes, which mediate pro-inflammatory signals during bacterial infection. We also found higher heterozygosity in immune genes in the southern population, where disease incidence and pathogen diversity are greater. Similar to the major histocompatibility complex in mammals, balancing selection may enhance genetic diversity in the innate immune system genes of this invertebrate. Overall, our results show that how genome-wide polymorphism data coupled with growing databases on gene function and expression can combine to detect otherwise hidden signals of selection in natural populations.

Populations of Symbiodinium muscatinei show strong biogeographic structuring in the intertidal anemone Anthopleura elegantissima

Among temperate cnidarian symbioses, the partnership between the intertidal anemone Anthopleura elegantissima and its dinoflagellate and chlorophyte symbionts is one of the most well characterized. Biogeographic, reciprocal transplant, and physiological studies have convincingly demonstrated a relationship between environmental factors such as temperature and irradiance and the distribution of symbionts from both algal phyla. However, little is known about the fine-scale diversity or biogeographic distribution within symbiont lineages of this anemone. We used sequence information from the mitochondrial cytochrome b and chloroplast 23S ribosomal genes and restriction fragment length polymorphism data from the 18S nuclear ribosomal gene to characterize the Symbiodinium populations in tentacles clipped from 105 anemones at 14 sites along the entire California coast, spanning about 1200 km. Our results show the presence of at least three primary biogeographic regions with breaks around Cape Mendocino and Monterey Bay, each dominated by a different Symbiodinium muscatinei genotype. Sharp clines suggest limited gene flow between adjacent regions. Few sampling locations or individual anemones showed symbiont diversity at either organellar locus within the limits of our detection method, while sequence analysis of cloned nr18S polymerase chain reaction product suggests that nuclear pseudogenes may underlie intra-host diversity observed at that locus.

Local adaptation in marine invertebrates

Local adaptation in the sea was regarded historically as a rare phenomenon that was limited to a handful of species with exceptionally low dispersal potential. However, a growing body of experimental studies indicates that adaptive differentiation occurs in numerous marine invertebrates in response to selection imposed by strong gradients (and more complex mosaics) of abiotic and biotic conditions. Moreover, a surprisingly high proportion of the marine invertebrates known or suspected of exhibiting local adaptation are species with planktonic dispersal. Adaptive divergence among populations can occur over a range of spatial scales, including those that are fine-grained (i.e., meters to kilometers), reflecting a balance between scales of gene flow and selection. Addressing the causes and consequences of adaptive genetic differentiation among invertebrate populations promises to advance community ecology, climate change research, and the effective management of marine ecosystems.