Transcriptome-wide polymorphisms of red abalone (Haliotis rufescens) reveal patterns of gene flow and local adaptation

Combination of RNAseq and RADseq to Identify Physiological and Adaptive Responses to Acidification in the Eastern Oyster (Crassostrea virginica)

Ocean acidification (OA) is a major stressor threatening marine calcifiers, including the eastern oyster (Crassostrea virginica). In this paper, we provide insight into the molecular mechanisms associated with resilience to OA, with the dual intentions of probing both acclimation and adaptation potential in this species. C. virginica were spawned, and larvae were reared in control or acidified conditions immediately after fertilization. RNA samples were collected from larvae and juveniles, and DNA samples were collected from juveniles after undergoing OA-induced mortality and used to contrast gene expression (RNAseq) and SNP (ddRADseq) profiles from animals reared under both conditions. Results showed convergence of evidence from both approaches, particularly in genes involved in biomineralization that displayed significant changes in variant frequencies and gene expression levels among juveniles that survived acidification as compared to controls. Downregulated genes were related to immune processes, supporting previous studies demonstrating a reduction in immunity from exposure to OA. Acclimation to OA via regulation of gene expression might confer short-term resilience to immediate threats; however, the costs may not be sustainable, underscoring the importance of selection of resilient genotypes. Here, we identified SNPs associated with survival under OA conditions, suggesting that this commercially and ecologically important species might have the genetic variation needed for adaptation to future acidification. The identification of genetic features associated with OA resilience is a highly-needed step for the development of marker-assisted selection of oyster stocks for aquaculture and restoration activities.

Genomic Signatures of Local Adaptation under High Gene Flow in Lumpfish-Implications for Broodstock Provenance Sourcing and Larval Production

Aquaculture of the lumpfish (Cyclopterus lumpus L.) has become a large, lucrative industry owing to the escalating demand for "cleaner fish" to minimise sea lice infestations in Atlantic salmon mariculture farms. We used over 10K genome-wide single nucleotide polymorphisms (SNPs) to investigate the spatial patterns of genomic variation in the lumpfish along the coast of Norway and across the North Atlantic. Moreover, we applied three genome scans for outliers and two genotype-environment association tests to assess the signatures and patterns of local adaptation under extensive gene flow. With our 'global' sampling regime, we found two major genetic groups of lumpfish, i.e., the western and eastern Atlantic. Regionally in Norway, we found marginal evidence of population structure, where the population genomic analysis revealed a small portion of individuals with a different genetic ancestry. Nevertheless, we found strong support for local adaption under high gene flow in the Norwegian lumpfish and identified over 380 high-confidence environment-associated loci linked to gene sets with a key role in biological processes associated with environmental pressures and embryonic development. Our results bridge population genetic/genomics studies with seascape genomics studies and will facilitate genome-enabled monitoring of the genetic impacts of escapees and allow for genetic-informed broodstock selection and management in Norway.

A Draft Reference Genome Assembly of the Critically Endangered Black Abalone, Haliotis cracherodii

The once abundant black abalone, Haliotis cracherodii, is a large, long-lived grazing marine mollusk that inhabits the rocky intertidal along the coast of California. The species has experienced dramatic declines since the mid-1980s largely due to the fatal bacterial disease called withering syndrome, leading to the collapse of an economically important fishery and to its inclusion into the IUCN listing as a critically endangered species. In some places impacted by the disease, populations of black abalone have declined by more than 90%, prompting population crashes associated with very little recruitment of new individuals and changes to intertidal communities. Habitats that were dominated by crustose coralline algae and bare rock have become dominated instead by fleshy algae and sessile invertebrates. Here, we present the first high-quality black abalone reference genome, assembled with PacBio HiFi long-reads and assembled with Dovetail Omni-C data to generate a scaffold-level assembly. The black abalone reference genome will be an essential resource in understanding the evolutionary history of this species as well as for exploring its current levels of genetic diversity and establishing future management and restoration plans.

Differential gene expression and chemical patterns of an intertidal crab inhabiting a polluted port and an adjacent marine protected area

The acquisition of data to safeguard marine protected areas located close to ports is important in order to develop plans that allow effective protection from pollution as well as sustainable development of the port. The area Secche della Meloria is a Marine Protected Area (MPA-MEL) three miles from Livorno Harbour (LH), which is characterized by a long history of pollution. Here we studied the bioaccumulation and transcriptomic patterns of the marbled crab, Pachygrapsus marmoratus (Fabricius, 1787) (Crustacea; Brachyura, Grapsidae), inhabiting the two selected sites. Results showed that the two crab populations are significantly different in their chemical composition of trace elements and Polyciclic Aromatic Hydrocarbons (PAHs), and gene expression patterns (1280 DEGs). Enrichment analysis indicated that crabs at LH had the highest stress response genes, and they were associated with higher levels of bioaccumulation detected in body tissues. We are confident that the significant differential gene expression profiles observed between crabs, characterized by significant chemical differences, is associated with responses to contaminant exposure.

Long-read RNA sequencing of Pacific abalone Haliotis discus hannai reveals innate immune system responses to environmental stress

Haliotis discus hannai is a commercially important mollusk species, and the abalone aquaculture sector has been jeopardized by deteriorating environmental circumstances such as bacterial infection and thermal stress during the hot summers. However, due to a paucity of genetic information, such as transcriptome resources, our understanding of their stress adaptation is restricted. In this research, using single-molecule long-read (SMRT) sequencing technology, a library composed of ten tissues (i.e., haemocytes, gills, muscle, hepatopancreas, digestive tract, mantle, mucous gland, ovary, testis and head) was constructed and sequenced. In all, 41,855 high-quality unique transcripts, among which 24,778 were successfully annotated. Additionally, 13,463 SSRs, 1,169 transcription factors, and 18,124 lncRNAs were identified in H. discus hannai transcriptome. Furthermore, multiple immune-related transcripts were identified according to KEGG annotation, and a portion of these transcripts were mapped into several classical immune-related pathways, including the PI3K-AKT signaling pathway and Toll-like receptor signaling pathway. Additionally, 24 typical sequences related to the immunity pathway were detected by RT-PCR; the results showed that most of the immune-related genes showed significantly high expression at 72 h after bacterial challenges and thermal stress, especially the expression level of genes in gills was significantly higher than that in haemocytes under V. parahaemolyticus stress at 24 h. At the same time. The analysis of alternative splicing identified several innate immunity-related functions genes, including CD109 and caspase 2. These results suggest that the complex immune system, particularly the powerful innate immunity system, was crucial for H. discus hannai response to numerous environmental challenges.

Population transcriptomics reveals the effect of gene flow on the evolution of range limits

One of the most important questions in evolutionary biology is how the spatial distribution of species is limited. Asymmetric gene flow from core populations is suggested to increase the number of poorly adapted immigrants in the populations at the range edge. Genetic load due to migration, i.e., migration load, should prevent adaptation to the local habitat, leading to decreases in distribution range via local extinction or the limiting range expansion. However, few experimental studies have examined the effects of immigration on fitness and natural selection within recipient populations. To investigate the influence of migration load on the evolution of distribution range, we performed field and laboratory observations as well as population transcriptomics for the common river snail, Semisulcospira reiniana. This species meets the conditions that migration from source populations can prevent local adaptation in a sink population because they inhabit the broader range of environments, including middle/upper reaches of a river and estuaries within a single river and they may be more vulnerable to being swept away by water currents due to lowered spontaneous (upward) locomotion activity. We found that river steepness was related to the lower distribution limit of S. reiniana, with a narrower distribution range in the steeper river. Population transcriptomic analysis showed that gene flow was heavily asymmetric from the upstream populations to downstream ones in the steep river, suggesting a greater migration load in the steep river. The number of genes putatively involved in adaptation to the local habitat was lower in the steep river than in the gentle river. Gene expression profiles suggested that individuals achieve better local adaptation in the gentle river. Laboratory experiments suggested that evolutionary differences in salinity tolerance among local populations were only found in the gentle river. Our results consistent with the hypothesis that migration load owing to asymmetric gene flow disturbs local adaptation and restricts the distribution range of river snails.

Genomic Evidence for Speciation with Gene Flow in Broadcast Spawning Marine Invertebrates

How early stages of speciation in free-spawning marine invertebrates proceed is poorly understood. The Western Pacific abalones, Haliotis discus, H. madaka, and H. gigantea, occur in sympatry with shared breeding season and are capable of producing viable F1 hybrids in spite of being ecologically differentiated. Population genomic analyses revealed that although the three species are genetically distinct, there is evidence for historical and ongoing gene flow among these species. Evidence from demographic modeling suggests that reproductive isolation among the three species started to build in allopatry and has proceeded with gene flow, possibly driven by ecological selection. We identified 27 differentiation islands between the closely related H. discus and H. madaka characterized by high FST and dA, but not high dXY values, as well as high genetic diversity in one H. madaka population. These genomic signatures suggest differentiation driven by recent ecological divergent selection in presence of gene flow outside of the genomic islands of differentiation. The differentiation islands showed low polymorphism in H. gigantea, and both high FST, dXY, and dA values between H. discus and H. gigantea, as well as between H. madaka and H. gigantea. Collectively, the Western Pacific abalones appear to occupy the early stages speciation continuum, and the differentiation islands associated with ecological divergence among the abalones do not appear to have acted as barrier loci to gene flow in the younger divergences but appear to do so in older divergences.

Genomic signatures of spatially divergent selection at clownfish range margins

Understanding how evolutionary forces interact to drive patterns of selection and distribute genetic variation across a species' range is of great interest in ecology and evolution, especially in an era of global change. While theory predicts how and when populations at range margins are likely to undergo local adaptation, empirical evidence testing these models remains sparse. Here, we address this knowledge gap by investigating the relationship between selection, gene flow and genetic drift in the yellowtail clownfish, Amphiprion clarkii, from the core to the northern periphery of the species range. Analyses reveal low genetic diversity at the range edge, gene flow from the core to the edge and genomic signatures of local adaptation at 56 single nucleotide polymorphisms in 25 candidate genes, most of which are significantly correlated with minimum annual sea surface temperature. Several of these candidate genes play a role in functions that are upregulated during cold stress, including protein turnover, metabolism and translation. Our results illustrate how spatially divergent selection spanning the range core to the periphery can occur despite the potential for strong genetic drift at the range edge and moderate gene flow from the core populations.

Transcriptome Expression of Biomineralization Genes in Littoraria flava Gastropod in Brazilian Rocky Shore Reveals Evidence of Local Adaptation

Understanding how selection shapes population differentiation and local adaptation in marine species remains one of the greatest challenges in the field of evolutionary biology. The selection of genes in response to environment-specific factors and microenvironmental variation often results in chaotic genetic patchiness, which is commonly observed in rocky shore organisms. To identify these genes, the expression profile of the marine gastropod Littoraria flava collected from four Southeast Brazilian locations in ten rocky shore sites was analyzed. In this first L. flava transcriptome, 250,641 unigenes were generated, and 24% returned hits after functional annotation. Independent paired comparisons between 1) transects, 2) sites within transects, and 3) sites from different transects were performed for differential expression, detecting 8,622 unique differentially expressed genes. Araçá (AR) and São João (SJ) transect comparisons showed the most divergent gene products. For local adaptation, fitness-related differentially expressed genes were chosen for selection tests. Nine and 24 genes under adaptative and purifying selection, respectively, were most related to biomineralization in AR and chaperones in SJ. The biomineralization-genes perlucin and gigasin-6 were positively selected exclusively in the site toward the open ocean in AR, with sequence variants leading to pronounced protein structure changes. Despite an intense gene flow among L. flava populations due to its planktonic larva, gene expression patterns within transects may be the result of selective pressures. Our findings represent the first step in understanding how microenvironmental genetic variation is maintained in rocky shore populations and the mechanisms underlying local adaptation in marine species.

Genomic signals found using RNA sequencing show signatures of selection and subtle population differentiation in walleye (Sander vitreus) in a large freshwater ecosystem

RNA sequencing is an effective approach for studying aquatic species yielding both physiological and genomic data. However, its population genetic applications are not well-characterized. We investigate this possible role for RNA sequencing for population genomics in Lake Winnipeg, Manitoba, Canada, walleye (Sander vitreus). Lake Winnipeg walleye represent the largest component of the second-largest freshwater fishery in Canada. In the present study, large female walleye were sampled via nonlethal gill biopsy over two years at three spawning sites representing a latitudinal gradient in the lake. Genetic variation from sequenced mRNA was analyzed for neutral and adaptive markers to investigate population structure and possible adaptive variation. We find low population divergence (F ST = 0.0095), possible northward gene flow, and outlier loci that vary latitudinally in transcripts associated with cell membrane proteins and cytoskeletal function. These results indicate that Lake Winnipeg walleye may be effectively managed as a single demographically connected metapopulation with contributing subpopulations and suggest genomic differences possibly underlying observed phenotypic differences. Despite its high cost relative to other genotyping methods, RNA sequencing data can yield physiological in addition to genetic information discussed here. We therefore argue that it is useful for addressing diverse molecular questions in the conservation of freshwater species.

An empirical comparison of population genetic analyses using microsatellite and SNP data for a species of conservation concern

BACKGROUND

Use of genomic tools to characterize wildlife populations has increased in recent years. In the past, genetic characterization has been accomplished with more traditional genetic tools (e.g., microsatellites). The explosion of genomic methods and the subsequent creation of large SNP datasets has led to the promise of increased precision in population genetic parameter estimates and identification of demographically and evolutionarily independent groups, as well as questions about the future usefulness of the more traditional genetic tools. At present, few empirical comparisons of population genetic parameters and clustering analyses performed with microsatellites and SNPs have been conducted.

RESULTS

Here we used microsatellite and SNP data generated from Gunnison sage-grouse (Centrocercus minimus) samples to evaluate concordance of the results obtained from each dataset for common metrics of genetic diversity (HO, HE, FIS, AR) and differentiation (FST, GST, DJost). Additionally, we evaluated clustering of individuals using putatively neutral (SNPs and microsatellites), putatively adaptive, and a combined dataset of putatively neutral and adaptive loci. We took particular interest in the conservation implications of any differences. Generally, we found high concordance between microsatellites and SNPs for HE, FIS, AR, and all differentiation estimates. Although there was strong correlation between metrics from SNPs and microsatellites, the magnitude of the diversity and differentiation metrics were quite different in some cases. Clustering analyses also showed similar patterns, though SNP data was able to cluster individuals into more distinct groups. Importantly, clustering analyses with SNP data suggest strong demographic independence among the six distinct populations of Gunnison sage-grouse with some indication of evolutionary independence in two or three populations; a finding that was not revealed by microsatellite data.

CONCLUSION

We demonstrate that SNPs have three main advantages over microsatellites: more precise estimates of population-level diversity, higher power to identify groups in clustering methods, and the ability to consider local adaptation. This study adds to a growing body of work comparing the use of SNPs and microsatellites to evaluate genetic diversity and differentiation for a species of conservation concern with relatively high population structure and using the most common method of obtaining SNP genotypes for non-model organisms.

Multi-omics investigations within the Phylum Mollusca, Class Gastropoda: from ecological application to breakthrough phylogenomic studies

Gastropods are the largest and most diverse class of mollusc and include species that are well studied within the areas of taxonomy, aquaculture, biomineralization, ecology, microbiome and health. Gastropod research has been expanding since the mid-2000s, largely due to large-scale data integration from next-generation sequencing and mass spectrometry in which transcripts, proteins and metabolites can be readily explored systematically. Correspondingly, the huge data added a great deal of complexity for data organization, visualization and interpretation. Here, we reviewed the recent advances involving gastropod omics ('gastropodomics') research from hundreds of publications and online genomics databases. By summarizing the current publicly available data, we present an insight for the design of useful data integrating tools and strategies for comparative omics studies in the future. Additionally, we discuss the future of omics applications in aquaculture, natural pharmaceutical biodiscovery and pest management, as well as to monitor the impact of environmental stressors.

Error, noise and bias in de novo transcriptome assemblies

De novo transcriptome assembly is a powerful tool, and has been widely used over the last decade for making evolutionary inferences. However, it relies on two implicit assumptions: that the assembled transcriptome is an unbiased representation of the underlying expressed transcriptome, and that expression estimates from the assembly are good, if noisy approximations of the relative abundance of expressed transcripts. Using publicly available data for model organisms, we demonstrate that, across assembly algorithms and data sets, these assumptions are consistently violated. Bias exists at the nucleotide level, with genotyping error rates ranging from 30% to 83%. As a result, diversity is underestimated in transcriptome assemblies, with consistent underestimation of heterozygosity in all but the most inbred samples. Even at the gene level, expression estimates show wide deviations from map-to-reference estimates, and positive bias at lower expression levels. Standard filtering of transcriptome assemblies improves the robustness of gene expression estimates but leads to the loss of a meaningful number of protein-coding genes, including many that are highly expressed. We demonstrate a computational method, length-rescaled CPM, to partly alleviate noise and bias in expression estimates. Researchers should consider ways to minimize the impact of bias in transcriptome assemblies.

Molecular adaptation of molluscan biomineralisation to high-CO2 oceans - The known and the unknown

High-CO₂ induced ocean acidification (OA) reduces the calcium carbonate (CaCO₃) saturation level (Ω) and the pH of oceans. Consequently, OA is causing a serious threat to several ecologically and economically important biomineralising molluscs. Biomineralisation is a highly controlled biochemical process by which molluscs deposit their calcareous structures. In this process, shell matrix proteins aid the nucleation, growth and assemblage of the CaCO₃ crystals in the shell. These molluscan shell proteins (MSPs) are, ultimately, responsible for determination of the diverse shell microstructures and mechanical strength. Recent studies have attempted to integrate gene and protein expression data of MSPs with shell structure and mechanical properties. These advances made in understanding the molecular mechanism of biomineralisation suggest that molluscs either succumb or adapt to OA stress. In this review, we discuss the fate of biomineralisation process in future high-CO₂ oceans and its ultimate impact on the mineralised shell's structure and mechanical properties from the perspectives of limited substrate availability theory, proton flux limitation model and the omega myth theory. Furthermore, studying the interplay of energy availability and differential gene expression is an essential first step towards understanding adaptation of molluscan biomineralisation to OA, because if there is a need to change gene expression under stressors, any living system would require more energy than usual. To conclude, we have listed, four important future research directions for molecular adaptation of molluscan biomineralisation in high-CO₂ oceans: 1) Including an energy budgeting factor while understanding differential gene expression of MSPs and ion transporters under OA. 2) Unraveling the genetic or epigenetic changes related to biomineralisation under stressors to help solving a bigger picture about future evolution of molluscs, and 3) Understanding Post Translational Modifications of MSPs with and without stressors. 4) Understanding carbon uptake mechanisms across taxa with and without OA to clarify the OA theories on Ω.

Genetic structure of Australian glass shrimp, Paratya australiensis, in relation to altitude

Paratya australiensis Kemp (Decapoda: Atyidae) is a widely distributed freshwater shrimp in eastern Australia. The species has been considered as an important stream organism for studying genetics, dispersal, biology, behaviour and evolution in atyids and is a major food source for stream dwelling fishes. Paratya australiensis is a cryptic species complex consisting of nine highly divergent mitochondrial DNA lineages. Previous studies in southeast Queensland showed that “lineage 4” favours upstream sites at higher altitudes, with cooler water temperatures. This study aims to identify putative selection and population structure between high elevation and low elevation populations of this lineage at relatively small spatial scales. Sample localities were selected from three streams: Booloumba Creek, Broken Bridge Creek and Obi Obi Creek in the Conondale Range, southeast Queensland. Six sample localities, consisting of 142 individuals in total were sequenced using double digest Restriction Site Associated DNA-sequencing (ddRAD-seq) technique. Among the 142 individuals, 131 individuals shared 213 loci. Outlier analysis on 213 loci showed that 27 loci were putatively under selection between high elevation and low elevation populations. Outlier analysis on individual streams was also done to test for parallel patterns of adaptation, but there was no evidence of a parallel pattern. Population structure was observed using both the 27 outliers and 186 neutral loci and revealed similar population structure in both cases. Therefore, we cannot differentiate between selection and drift here. The highest genetic differentiation was observed between high elevation and low elevation populations of Booloumba Creek, with small levels of differentiation in the other two streams.

Immunity-related genes and signaling pathways under hypoxic stresses in Haliotis diversicolor: a transcriptome analysis

Due to increased temperatures and aquaculture density, thermal and hypoxia stresses have become serious problems for the aquaculture of abalone Haliotis diversicolor. Stresses lead to immunosuppression, which can cause severe negative impacts on aquaculture farms. To study the mechanism of immunosuppression after hypoxia stress and bacterial challenge, transcriptomes of H. diversicolor hemocytes involved in immunity were profiled. A total of 307,395,572 clean reads were generated and assembled into 99,774 unigenes. KEGG analysis indicated that 225 unigenes with immunologic function were mapped into immune-related pathways. Expression of 41 unigenes measured by quantitative real-time PCR (qRT-PCR) showed consistent results with that of transcriptome analysis. When exposure challenge of Vibrio parahaemolyticus, it is indicated that the PI3K-AKT, MAPK, NF-κB and P53 signal pathways were involved in the hypoxia-induced immunosuppression of H. diversicolor. Furthermore, when the AKT gene (HdAKT) was inhibited by double-stranded RNA (dsRNA), expression levels of HdAKT was lower than the blank and control group in hemocytes at 4 h, 12 h and 24 h (p < 0.05).

Signatures of adaptive divergence among populations of an avian species of conservation concern

Understanding the genetic underpinning of adaptive divergence among populations is a key goal of evolutionary biology and conservation. Gunnison sage-grouse (Centrocercus minimus) is a sagebrush obligate species with a constricted range consisting of seven discrete populations, each with distinctly different habitat and climatic conditions. Though geographically close, populations have low levels of natural gene flow resulting in relatively high levels of differentiation. Here, we use 15,033 SNP loci in genomic outlier analyses, genotype-environment association analyses, and gene ontology enrichment tests to examine patterns of putatively adaptive genetic differentiation in an avian species of conservation concern. We found 411 loci within 5 kbp of 289 putative genes associated with biological functions or pathways that were overrepresented in the assemblage of outlier SNPs. The identified gene set was enriched for cytochrome P450 gene family members (CYP4V2, CYP2R1, CYP2C23B, CYP4B1) and could impact metabolism of plant secondary metabolites, a critical challenge for sagebrush obligates. Additionally, the gene set was also enriched with members potentially involved in antiviral response (DEAD box helicase gene family and SETX). Our results provide a first look at local adaption for isolated populations of a single species and suggest adaptive divergence in multiple metabolic and biochemical pathways may be occurring. This information can be useful in managing this species of conservation concern, for example, to identify unique populations to conserve, avoid translocation or release of individuals that may swamp locally adapted genetic diversity, or guide habitat restoration efforts.

Development of EST-based SSR and SNP markers in Gastrodia elata (herbal medicine) by sequencing, de novo assembly and annotation of the transcriptome

Tianma (Gastrodia elata Blume) has unique biological characteristics and high medicinal value. The wild resource of G. elata is being overutilized and should be conserved as it is already included in the list of endangered species in China. The population size of cultivated G. elata is small because of domestication bottleneck. Therefore, it is of utmost importance to evolve high-quality varieties and conserve wild resources of G. elata. In this study, we sequenced tuber transcriptomes of three major cultivated sub-species of Gastrodia elata, namely G. elata BI. f. elata, G. elata Bl. f. glauca S. Chow, and G. elata Bl. f. Viridis, and obtained about 7.8G clean data. The assembled high-quality reads of three sub-species were clustered into 56,884 unigenes. Of these, 31,224 (54.89%), 25,733 (45.24%), 22,629 (39.78%), and 11,856 (20.84%) unigenes were annotated by Nr, Swiss-Port, Eukaryotic Ortholog Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. Here, a total of 3766 EST-SSRs and 128,921 SNPs were identified from the unigenes. The results not only offer huge number of genes that were responsible for the growth, development, and metabolism of bioactive components, but also a large number of molecular markers were detected for future studies on the conservation genetics and molecular breeding of G. elata.

Transcriptional profiles of early stage red sea urchins (Mesocentrotus franciscanus) reveal differential regulation of gene expression across development

The red sea urchin, Mesocentrotus franciscanus, is an ecologically important kelp forest species that also serves as a valuable fisheries resource. In this study, we have assembled and annotated a developmental transcriptome for M. franciscanus that represents eggs and six stages of early development (8- to 16-cell, morula, hatched blastula, early gastrula, prism and early pluteus). Characterization of the transcriptome revealed distinct patterns of gene expression that corresponded to major developmental and morphological processes. In addition, the period during which maternally-controlled transcription was terminated and the zygotic genome was activated, the maternal-to-zygotic transition (MZT), was found to begin during early cleavage and persist through the hatched blastula stage, an observation that is similar to the timing of the MZT in other sea urchin species. The presented developmental transcriptome will serve as a useful resource for investigating, in both an ecological and fisheries context, how the early developmental stages of this species respond to environmental stressors.

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.

Effects of Predator-Prey Interactions on Predator Traits: Differentiation of Diets and Venoms of a Marine Snail

Species interactions are fundamental ecological forces that can have significant impacts on the evolutionary trajectories of species. Nonetheless, the contribution of predator-prey interactions to genetic and phenotypic divergence remains largely unknown. Predatory marine snails of the family Conidae exhibit specializations for different prey items and intraspecific variation in prey utilization patterns at geographic scales. Because cone snails utilize venom to capture prey and venom peptides are direct gene products, it is feasible to examine the evolution of genes associated with changes in resource utilization. Here, we compared feeding ecologies and venom duct transcriptomes of individuals from three populations of Conus miliaris, a species that exhibits geographic variation in prey utilization and dietary breadth, in order to determine the extent to which dietary differences are correlated with differences in venom composition, and if expanded niche breadth is associated with increased variation in venom composition. While populations showed little to no overlap in resource utilization, taxonomic richness of prey was greatest at Easter Island. Changes in dietary breadth were associated with differences in expression patterns and increased genetic differentiation of toxin-related genes. The Easter Island population also exhibited greater diversity of toxin-related transcripts, but did not show increased variance in expression of these transcripts. These results imply that differences in dietary breadth contribute more to the structural and regulatory differentiation of venoms than differences in diet.

Differential responses to ocean acidification between populations of Balanophyllia elegans corals from high and low upwelling environments

Ocean acidification (OA), the global decrease in surface water pH from absorption of anthropogenic CO₂ , may put many marine taxa at risk. However, populations that experience extreme localized conditions, and are adapted to these conditions predicted in the global ocean in 2,100, may be more tolerant to future OA. By identifying locally adapted populations, researchers can examine the mechanisms used to cope with decreasing pH. One oceanographic process that influences pH is wind-driven upwelling. Here we compare two Californian populations of the coral Balanophyllia elegans from distinct upwelling regimes, and test their physiological and transcriptomic responses to experimental seawater acidification. We measured respiration rates, protein and lipid content, and gene expression in corals from both populations exposed to pH levels of 7.8 and 7.4 for 29 days. Corals from the population that experiences lower pH due to high upwelling maintained the same respiration rate throughout the exposure. In contrast, corals from the low upwelling site had reduced respiration rates, protein content and lipid-class content at low pH exposure, suggesting they have depleted their energy reserves. Using RNA-Seq, we found that corals from the high upwelling site upregulated genes involved in calcium ion binding and ion transport, most likely related to pH homeostasis and calcification. In contrast, corals from the low upwelling site downregulated stress response genes at low pH exposure. Divergent population responses to low pH observed in B. elegans highlight the importance of multi-population studies for predicting a species' response to future OA.

Different Transcriptomic Responses to Thermal Stress in Heat-Tolerant and Heat-Sensitive Pacific Abalones Indicated by Cardiac Performance

The Pacific abalone Haliotis discus hannai is one of the most economically important mollusks in China. Even though it has been farmed in southern China for almost 20 years, summer mortality remains the most challengeable problem for Pacific abalone aquaculture recently. Here, we determined the different heat tolerance ability for five selective lines of H. discus hannai by measuring the cardiac performance and Arrhenius breakpoint temperature (ABT). The Red line (RL) and Yangxia line (YL) were determined as the most heat-sensitive and most heat-tolerant line, respectively. Heart rates for RL were significantly lower than those of the YL at the same temperature (p < 0.05). The differentially expressed genes (DEGs), which were enriched in several pathways including cardiac muscle contraction, glutathione metabolism and oxidative phosphorylation, were identified between RL and YL at control temperature (20°C) and heat stress temperature (28.5°C, the ABT of the RL) by RNA-seq method. In the RL, 3370 DEGs were identified between the control and the heat-stress temperature, while only 1351 DEGs were identified in YL between these two temperature tests. Most of these DEGs were enriched in the pathways such as protein processing in endoplasmic reticulum, nucleotide binding and oligomerization domain (NOD) like receptor signaling, and ubiquitin mediated proteolysis. Notably, the most heat-tolerant line YL used an effective heat-protection strategy based on moderate transcriptional changes and regulation on the expression of key genes.

Integrative transcriptome analysis and discovery of genes involving in immune response of hypoxia/thermal challenges in the small abalone Haliotis diversicolor

In recent years, the abalone aquaculture industry has been threatened by the deteriorating environmental conditions, such as hypoxia and thermal stress in the hot summers. It is necessary to investigate the molecular mechanism in response to these environmental challenges, and subsequently understand the immune defense system. In this study, the transcriptome profiles by RNA-seq of hemocytes from the small abalone Haliotis diversicolor after exposure to hypoxia, thermal stress, and hypoxia plus thermal stress were established. A total of 103,703,074 clean reads were obtained and 99,774 unigenes were assembled. Of the 99,774 unigenes, 47,154 and 20,455 had homologous sequences in the Nr and Swiss-Prot protein databases, while 16,944 and 10,840 unigenes could be classified by COG or KEGG databases, respectively. RNAseq analysis revealed that the differentially expressed genes (DEGs) after challenges of hypoxia, thermal stress, or hypoxia plus thermal stress were 24,189, 29,165 and 23,665, among which more than 3000 genes involved in at least 230 pathways, including several classical immune-related pathways. The genes and pathways that were involved in immune response to hypoxia/thermal challenges were identified by transcriptome analysis and further validated by quantitative real-time PCR and RNAi technology. The findings in this study can provide information on H. diversicolor innate immunity to improve the abalone aquaculture industry, and the analysis of the potential immune-related genes in innate immunity signaling pathways and the obtained transcriptome data can provide an invaluable genetic resource for the study of the genome and functional genes.

Pathogens in space: Advancing understanding of pathogen dynamics and disease ecology through landscape genetics

Landscape genetics has provided many insights into how heterogeneous landscape features drive processes influencing spatial genetic variation in free-living organisms. This rapidly developing field has focused heavily on vertebrates, and expansion of this scope to the study of infectious diseases holds great potential for landscape geneticists and disease ecologists alike. The potential application of landscape genetics to infectious agents has garnered attention at formative stages in the development of landscape genetics, but systematic examination is lacking. We comprehensively review how landscape genetics is being used to better understand pathogen dynamics. We characterize the field and evaluate the types of questions addressed, approaches used and systems studied. We also review the now established landscape genetic methods and their realized and potential applications to disease ecology. Lastly, we identify emerging frontiers in the landscape genetic study of infectious agents, including recent phylogeographic approaches and frameworks for studying complex multihost and host-vector systems. Our review emphasizes the expanding utility of landscape genetic methods available for elucidating key pathogen dynamics (particularly transmission and spread) and also how landscape genetic studies of pathogens can provide insight into host population dynamics. Through this review, we convey how increasing awareness of the complementarity of landscape genetics and disease ecology among practitioners of each field promises to drive important cross-disciplinary advances.

Population genomic footprints of environmental pollution pressure in natural populations of the Mediterranean mussel

Bivalve molluscs of the genus Mytilus are considered a model organism in ecotoxicology and are known to be well adapted to marine ecosystems affected by multiple anthropogenic factors, including pollution. In order to assess whether pollution interferes with the reproductive success of Mytilus and affects the diversity within and between populations, we sequenced the transcriptomes of 72 individuals from 9 populations of Mytilus galloprovincialis collected along a ca. 130-km north-south transect on the Western coast of the Iberian Peninsula. We found that polluted areas are acting as a barrier to gene flow, potentially because of the detrimental effect of anthropogenic chemicals on larvae carried from more pristine environments. Furthermore, we observed an increase in genetic diversity in populations from polluted site, which could be indicative of higher mutagenicity driven by the environment. We propose that a microevolutionary perspective is essential for a full characterization of human activities on the dispersal of M. galloprovincialis and that it should be incorporated into management, and conservation plans and policies in the context of the effects of pollution on populations.

Concordant Changes in Gene Expression and Nucleotides Underlie Independent Adaptation to Hydrogen-Sulfide-Rich Environments

The colonization of novel environments often involves changes in gene expression, protein coding sequence, or both. Studies of how populations adapt to novel conditions, however, often focus on only one of these two processes, potentially missing out on the relative importance of different parts of the evolutionary process. In this study, our objectives were 1) to better understand the qualitative concordance between conclusions drawn from analyses of differential expression and changes in genic sequence and 2) to quantitatively test whether differentially expressed genes were enriched for sites putatively under positive selection within gene regions. To achieve this, we compared populations of fish (Poecilia mexicana) that have independently adapted to hydrogen-sulfide-rich environments in southern Mexico to adjacent populations residing in nonsulfidic waters. Specifically, we used RNA-sequencing data to compare both gene expression and DNA sequence differences between populations. Analyzing these two different data types led to similar conclusions about which biochemical pathways (sulfide detoxification and cellular respiration) were involved in adaptation to sulfidic environments. Additionally, we found a greater overlap between genes putatively under selection and differentially expressed genes than expected by chance. We conclude that considering both differential expression and changes in DNA sequence led to a more comprehensive understanding of how these populations adapted to extreme environmental conditions. Our results imply that changes in both gene expression and DNA sequence-sometimes at the same loci-may be involved in adaptation.

Genomic signatures of environmental selection despite near-panmixia in summer flounder

Rapid environmental change is altering the selective pressures experienced by marine species. While adaptation to local environmental conditions depends on a balance between dispersal and natural selection across the seascape, the spatial scale of adaptation and the relative importance of mechanisms maintaining adaptation in the ocean are not well understood. Here, using population assignment tests, Approximate Bayesian Computation (ABC), and genome scans with double-digest restriction-site associated DNA sequencing data, we evaluated population structure and locus-environment associations in a commercially important species, summer flounder (Paralichthys dentatus), along the U.S. east coast. Based on 1,137 single nucleotide polymorphisms across 232 individuals spanning nearly 1,900 km, we found no indication of population structure across Cape Hatteras, North Carolina (F ST = 0.0014) or of isolation by distance along the coast using individual relatedness. ABC estimated the probability of dispersal across the biogeographic break at Cape Hatteras to be high (95% credible interval: 7%-50% migration). However, we found 15 loci whose allele frequencies were associated with at least one of four environmental variables. Of those, 11 were correlated with bottom temperature. For summer flounder, our results suggest continued fisheries management as a single population and identify likely response mechanisms to climate change. Broadly speaking, our findings suggest that spatial balancing selection can manifest in adaptive divergence on regional scales in marine fish despite high dispersal, and that these conditions likely result in the widespread distribution of adaptive alleles and a high potential for future genetic adaptation in response to changing environmental conditions. In the context of a rapidly changing world, a landscape genomics perspective offers a useful approach for understanding the causes and consequences of genetic differentiation.

Complex signatures of genomic variation of two non-model marine species in a homogeneous environment

BACKGROUND

Genomic tools are increasingly being used on non-model organisms to provide insights into population structure and variability, including signals of selection. However, most studies are carried out in regions with distinct environmental gradients or across large geographical areas, in which local adaptation is expected to occur. Therefore, the focus of this study is to characterize genomic variation and selective signals over short geographic areas within a largely homogeneous region. To assess adaptive signals between microhabitats within the rocky shore, we compared genomic variation between the Cape urchin (Parechinus angulosus), which is a low to mid-shore species, and the Granular limpet (Scutellastra granularis), a high shore specialist.

RESULTS

Using pooled restriction site associated DNA (RAD) sequencing, we described patterns of genomic variation and identified outlier loci in both species. We found relatively low numbers of outlier SNPs within each species, and identified outlier genes associated with different selective pressures than those previously identified in studies conducted over larger environmental gradients. The number of population-specific outlier loci differed between species, likely owing to differential selective pressures within the intertidal environment. Interestingly, the outlier loci were highly differentiated within the two northernmost populations for both species, suggesting that unique evolutionary forces are acting on marine invertebrates within this region.

CONCLUSIONS

Our study provides a background for comparative genomic studies focused on non-model species, as well as a baseline for the adaptive potential of marine invertebrates along the South African west coast. We also discuss the caveats associated with Pool-seq and potential biases of sequencing coverage on downstream genomic metrics. The findings provide evidence of species-specific selective pressures within a homogeneous environment, and suggest that selective forces acting on small scales are just as crucial to acknowledge as those acting on larger scales. As a whole, our findings imply that future population genomic studies should expand from focusing on model organisms and/or studying heterogeneous regions to better understand the evolutionary processes shaping current and future biodiversity patterns, particularly when used in a comparative phylogeographic context.

Oceanographic variation influences spatial genomic structure in the sea scallop, Placopecten magellanicus

Environmental factors can influence diversity and population structure in marine species and accurate understanding of this influence can both improve fisheries management and help predict responses to environmental change. We used 7163 SNPs derived from restriction site-associated DNA sequencing genotyped in 245 individuals of the economically important sea scallop, Placopecten magellanicus, to evaluate the correlations between oceanographic variation and a previously identified latitudinal genomic cline. Sea scallops span a broad latitudinal area (>10 degrees), and we hypothesized that climatic variation significantly drives clinal trends in allele frequency. Using a large environmental dataset, including temperature, salinity, chlorophyll a, and nutrient concentrations, we identified a suite of SNPs (285-621, depending on analysis and environmental dataset) potentially under selection through correlations with environmental variation. Principal components analysis of different outlier SNPs and environmental datasets revealed similar northern and southern clusters, with significant associations between the first axes of each (R2adj = .66-.79). Multivariate redundancy analysis of outlier SNPs and the environmental principal components indicated that environmental factors explained more than 32% of the variance. Similarly, multiple linear regressions and random-forest analysis identified winter average and minimum ocean temperatures as significant parameters in the link between genetic and environmental variation. This work indicates that oceanographic variation is associated with the observed genomic cline in this species and that seasonal periods of extreme cold may restrict gene flow along a latitudinal gradient in this marine benthic bivalve. Incorporating this finding into management may improve accuracy of management strategies and future predictions.

Migration-Selection Balance Drives Genetic Differentiation in Genes Associated with High-Altitude Function in the Speckled Teal (Anas flavirostris) in the Andes

Local adaptation frequently occurs across populations as a result of migration-selection balance between divergent selective pressures and gene flow associated with life in heterogeneous landscapes. Studying the effects of selection and gene flow on the adaptation process can be achieved in systems that have recently colonized extreme environments. This study utilizes an endemic South American duck species, the speckled teal (Anas flavirostris), which has both high- and low-altitude populations. High-altitude speckled teal (A. f. oxyptera) are locally adapted to the Andean environment and mostly allopatric from low-altitude birds (A. f. flavirostris); however, there is occasional gene flow across altitudinal gradients. In this study, we used next-generation sequencing to explore genetic patterns associated with high-altitude adaptation in speckled teal populations, as well as the extent to which the balance between selection and migration have affected genetic architecture. We identified a set of loci with allele frequencies strongly correlated with altitude, including those involved in the insulin-like signaling pathway, bone morphogenesis, oxidative phosphorylation, responders to hypoxia-induced DNA damage, and feedback loops to the hypoxia-inducible factor pathway. These same outlier loci were found to have depressed gene flow estimates, as well as being highly concentrated on the Z-chromosome. Our results suggest a multifactorial response to life at high altitudes through an array of interconnected pathways that are likely under positive selection and whose genetic components seem to be providing an effective genomic barrier to interbreeding, potentially functioning as an avenue for population divergence and speciation.

Transcriptome sequencing and differential gene expression analysis of the schistosome-transmitting snail Oncomelania hupensis inhabiting hilly and marshland regions

The freshwater snail Oncomelania hupensis is the unique intermediate host of the blood fluke Schistosoma japonicum, which is the major cause of schistosomiasis. The snail inhabits two contrasting environments: the hilly and marshland regions. The hilly snails are smaller in size and have the typical smooth shell, whereas the marshland snails are larger and possess the ribbed shell. To reveal the differences in gene expression between the hilly and marshland snails, a total of six snails, three per environment, were individually examined by RNA sequencing technology. All paired-end reads were assembled into contigs from which 34,760 unigenes were predicted. Based on single nucleotide polymorphisms, principal component analysis and neighbor-joining clustering revealed two distinct clusters of hilly and marshland snails. Analysis of expression changes between environments showed that upregulated genes relating to immunity and development were enriched in hilly snails, while those associated with reproduction were over-represented in marshland snails. Eight differentially expressed genes between the two types of snails were validated by qRT-PCR. Our study identified candidate genes that could be targets for future functional studies, and provided a link between expression profiling and ecological adaptation of the snail that may have implications for schistosomiasis control.

Genomic models predict successful coral adaptation if future ocean warming rates are reduced

Population genomic surveys suggest that climate-associated genetic variation occurs widely across species, but whether it is sufficient to allow population persistence via evolutionary adaptation has seldom been quantified. To ask whether rapid adaptation in reef-building corals can keep pace with future ocean warming, we measured genetic variation at predicted warm-adapted loci and simulated future evolution and persistence in a high-latitude population of corals from Rarotonga, Cook Islands. Alleles associated with thermal tolerance were present but at low frequencies in this cooler, southerly locality. Simulations based on predicted ocean warming in Rarotonga showed rapid evolution of heat tolerance resulting in population persistence under mild warming scenarios consistent with low CO₂ emission plans, RCP2.6 and RCP4.5. Under more severe scenarios, RCP6.0 and RCP8.5, adaptation was not rapid enough to prevent extinction. Population adaptation was faster for models based on smaller numbers of additive loci that determine thermal tolerance and for higher population growth rates. Finally, accelerated migration via transplantation of thermally tolerant individuals (1 to 5%/year) sped adaptation. These results show that cool-water corals can adapt to warmer oceans but only under mild scenarios resulting from international emissions controls. Incorporation of genomic data into models of species response to climate change offers a promising method for estimating future adaptive processes.

Signatures of positive selection and local adaptation to urbanization in white-footed mice (Peromyscus leucopus)

Urbanization significantly alters natural ecosystems and has accelerated globally. Urban wildlife populations are often highly fragmented by human infrastructure, and isolated populations may adapt in response to local urban pressures. However, relatively few studies have identified genomic signatures of adaptation in urban animals. We used a landscape genomic approach to examine signatures of selection in urban populations of white-footed mice (Peromyscus leucopus) in New York City. We analysed 154,770 SNPs identified from transcriptome data from 48 P. leucopus individuals from three urban and three rural populations and used outlier tests to identify evidence of urban adaptation. We accounted for demography by simulating a neutral SNP data set under an inferred demographic history as a null model for outlier analysis. We also tested whether candidate genes were associated with environmental variables related to urbanization. In total, we detected 381 outlier loci and after stringent filtering, identified and annotated 19 candidate loci. Many of the candidate genes were involved in metabolic processes and have well-established roles in metabolizing lipids and carbohydrates. Our results indicate that white-footed mice in New York City are adapting at the biomolecular level to local selective pressures in urban habitats. Annotation of outlier loci suggests selection is acting on metabolic pathways in urban populations, likely related to novel diets in cities that differ from diets in less disturbed areas.

Bias and Correction in RNA-seq Data for Marine Species

RNA-seq is a recently developed approach widely used for transcriptome profiling in biological analyses that use next-generation sequencing technologies. Accurate estimation of gene expression levels is critical for answering biological questions. Here, we show that the commonly used measure of gene expression levels, fragments per kilobase of transcript per million mapped reads (FPKM), is biased in transcript length, GC content, and dinucleotide frequencies in the RNA-seq analysis of marine species. We used a generalized linear model to correct the observed biases of FPKM. We used RNA-seq data sets from eight species obtained by different sequencing methods to evaluate the correction methods. Our work contributes to the understanding of potential technical artifacts in RNA-seq experiments for marine species, and presents a means by which more accurate gene expression measures can be obtained.

Invertebrate population genetics across Earth's largest habitat: The deep-sea floor

Despite the deep sea being the largest habitat on Earth, there are just 77 population genetic studies of invertebrates (115 species) inhabiting non-chemosynthetic ecosystems on the deep-sea floor (below 200 m depth). We review and synthesize the results of these papers. Studies reveal levels of genetic diversity comparable to shallow-water species. Generally, populations at similar depths were well connected over 100s-1,000s km, but studies that sampled across depth ranges reveal population structure at much smaller scales (100s-1,000s m) consistent with isolation by adaptation across environmental gradients, or the existence of physical barriers to connectivity with depth. Few studies were ocean-wide (under 4%), and 48% were Atlantic-focused. There is strong emphasis on megafauna and commercial species with research into meiofauna, "ecosystem engineers" and other ecologically important species lacking. Only nine papers account for ~50% of the planet's surface (depths below 3,500 m). Just two species were studied below 5,000 m, a quarter of Earth's seafloor. Most studies used single-locus mitochondrial genes revealing a common pattern of non-neutrality, consistent with demographic instability or selective sweeps; similar to deep-sea hydrothermal vent fauna. The absence of a clear difference between vent and non-vent could signify that demographic instability is common in the deep sea, or that selective sweeps render single-locus mitochondrial studies demographically uninformative. The number of population genetics studies to date is miniscule in relation to the size of the deep sea. The paucity of studies constrains meta-analyses where broad inferences about deep-sea ecology could be made.

Calcifying algae maintain settlement cues to larval abalone following algal exposure to extreme ocean acidification

Ocean acidification (OA) increasingly threatens marine systems, and is especially harmful to calcifying organisms. One important question is whether OA will alter species interactions. Crustose coralline algae (CCA) provide space and chemical cues for larval settlement. CCA have shown strongly negative responses to OA in previous studies, including disruption of settlement cues to corals. In California, CCA provide cues for seven species of harvested, threatened, and endangered abalone. We exposed four common CCA genera and a crustose calcifying red algae, Peyssonnelia (collectively CCRA) from California to three pCO₂ levels ranging from 419–2,013 µatm for four months. We then evaluated abalone (Haliotis rufescens) settlement under ambient conditions among the CCRA and non-algal controls that had been previously exposed to the pCO₂ treatments. Abalone settlement and metamorphosis increased from 11% in the absence of CCRA to 45–69% when CCRA were present, with minor variation among CCRA genera. Though all CCRA genera reduced growth during exposure to increased pCO₂, abalone settlement was unaffected by prior CCRA exposure to increased pCO₂. Thus, we find no impacts of OA exposure history on CCRA provision of settlement cues. Additionally, there appears to be functional redundancy in genera of CCRA providing cues to abalone, which may further buffer OA effects.

Epipodial Tentacle Gene Expression and Predetermined Resilience to Summer Mortality in the Commercially Important Greenlip Abalone, Haliotis laevigata

"Summer mortality" is a phenomenon that occurs during warm water temperature spikes that results in the mass mortality of many ecologically and economically important mollusks such as abalone. This study aimed to determine whether the baseline gene expression of abalone before a laboratory-induced summer mortality event was associated with resilience to summer mortality. Tentacle transcriptomes of 35 greenlip abalone (Haliotis laevigata) were sequenced prior to the animals being exposed to an increase in water temperature-simulating conditions which have previously resulted in summer mortality. Abalone derived from three source locations with different environmental conditions were categorized as susceptible or resistant to summer mortality depending on whether they died or survived after the water temperature was increased. We detected two genes showing significantly higher expression in resilient abalone relative to susceptible abalone prior to the laboratory-induced summer mortality event. One of these genes was annotated through the NCBI non-redundant protein database using BLASTX to an anemone (Exaiptasia pallida) Transposon Ty3-G Gag Pol polyprotein. Distinct gene expression signatures were also found between resilient and susceptible abalone depending on the population origin, which may suggest divergence in local adaptation mechanisms for resilience. Many of these genes have been suggested to be involved in antioxidant and immune-related functions. The identification of these genes and their functional roles have enhanced our understanding of processes that may contribute to summer mortality in abalone. Our study supports the hypothesis that prestress gene expression signatures are indicative of the likelihood of summer mortality.

Co-Option and De Novo Gene Evolution Underlie Molluscan Shell Diversity

Molluscs fabricate shells of incredible diversity and complexity by localized secretions from the dorsal epithelium of the mantle. Although distantly related molluscs express remarkably different secreted gene products, it remains unclear if the evolution of shell structure and pattern is underpinned by the differential co-option of conserved genes or the integration of lineage-specific genes into the mantle regulatory program. To address this, we compare the mantle transcriptomes of 11 bivalves and gastropods of varying relatedness. We find that each species, including four Pinctada (pearl oyster) species that diverged within the last 20 Ma, expresses a unique mantle secretome. Lineage- or species-specific genes comprise a large proportion of each species' mantle secretome. A majority of these secreted proteins have unique domain architectures that include repetitive, low complexity domains (RLCDs), which evolve rapidly, and have a proclivity to expand, contract and rearrange in the genome. There are also a large number of secretome genes expressed in the mantle that arose before the origin of gastropods and bivalves. Each species expresses a unique set of these more ancient genes consistent with their independent co-option into these mantle gene regulatory networks. From this analysis, we infer lineage-specific secretomes underlie shell diversity, and include both rapidly evolving RLCD-containing proteins, and the continual recruitment and loss of both ancient and recently evolved genes into the periphery of the regulatory network controlling gene expression in the mantle epithelium.

Highly localized divergence within supergenes in Atlantic cod (Gadus morhua) within the Gulf of Maine

BACKGROUND

Atlantic cod (Gadus morhua), is known to vary genetically across the North Atlantic, Greenland, and Newfoundland. This genetic variation occurs both spatially and temporally through decades of heavy fishing, and is concentrated in three linkage disequilibrium blocks, previously defined by pedigreed linkage mapping analysis. Variation within these genomic regions is correlated with both seawater temperature and behavioral ecotype. The full extent and nature of these linkage groups is important information for interpreting cod genetic structure as a tool for future fisheries management.

RESULTS

We conducted whole genome sequencing for 31 individual cod from three sub-populations in the Gulf of Maine. Across the genome, we found 3,390,654 intermediate to high frequency Single Nucleotide Polymorphisms (SNPs). We show that pairwise linkage analysis among these SNPs is a powerful tool to detect linkage disequilibrium clusters by recovering the three previously detected linkage groups and identifying the 1031 genes contained therein. Across these genes, we found significant population differentiation among spawning groups in the Gulf of Maine and between Georges Bank and Gulf of Maine. Coordinated divergence among these genes and their differentiation at both short and long spatial scales suggests that they are acting as linked supergenes in local adaptation of cod populations.

CONCLUSIONS

Differentiation between SNPs in linkage disequilibrium blocks is the major signal of genetic differentiation between all groups tested within the Gulf of Maine. Our data provide a map of genes contained in these blocks, allowing an enhanced search for neutral genetic structure for demographic inference and fisheries modeling. Patterns of selection and the history of populations may be possible to identify in cod using this description of linkage disequilibrium blocks and future data sets to robustly separate neutral and selected genetic markers.