Genomewide investigation of adaptation to harmful algal blooms in common bottlenose dolphins (Tursiops truncatus)

Marine Neurotoxins' Effects on Environmental and Human Health: An OMICS Overview

Harmful algal blooms (HAB), and the consequent release of toxic metabolites, can be responsible for seafood poisoning outbreaks. Marine wildlife can accumulate these toxins throughout the food chain, which presents a threat to consumers’ health. Some of these toxins, such as saxitoxin (STX), domoic acid (DA), ciguatoxin (CTX), brevetoxin (BTX), tetrodotoxin (TTX), and β-N-methylamino-L-alanine (BMAA), cause severe neurological symptoms in humans. Considerable information is missing, however, notably the consequences of toxin exposures on changes in gene expression, protein profile, and metabolic pathways. This information could lead to understanding the consequence of marine neurotoxin exposure in aquatic organisms and humans. Nevertheless, recent contributions to the knowledge of neurotoxins arise from OMICS-based research, such as genomics, transcriptomics, proteomics, and metabolomics. This review presents a comprehensive overview of the most recent research and of the available solutions to explore OMICS datasets in order to identify new features in terms of ecotoxicology, food safety, and human health. In addition, future perspectives in OMICS studies are discussed.

Detecting genetic signals of selection in heavily bottlenecked reindeer populations by comparing parallel founder events

Founder populations are of special interest to both evolutionary and conservation biologists, but the detection of genetic signals of selection in these populations is challenging due to their demographic history. Geographically separated founder populations likely to have been subjected to similar selection pressures provide an ideal but rare opportunity to overcome these challenges. Here we take advantage of such a situation generated when small, isolated founder populations of reindeer were established on the island of South Georgia, and using this system we look for empirical evidence of selection overcoming strong genetic drift. We generated a 70 k ddRADseq single nucleotide polymorphism database for the two parallel reindeer founder populations and screened for signatures of soft sweeps. We find evidence for a genomic region under selection shared among the two populations, and support our findings with Wright-Fisher model simulations to assess the power and specificity of interpopulation selection scans-namely Bayescan, OutFLANK, PCadapt and a newly developed scan called Genome Wide Differentiation Scan (GWDS)-in the context of pairwise source-founder comparisons. Our simulations indicate that loci under selection in small founder populations are most probably detected by GWDS, and strengthen the hypothesis that the outlier region represents a true locus under selection. We explore possible, relevant functional roles for genes in linkage with the detected outlier loci.

How "simple" methodological decisions affect interpretation of population structure based on reduced representation library DNA sequencing: A case study using the lake whitefish

Reduced representation (RRL) sequencing approaches (e.g., RADSeq, genotyping by sequencing) require decisions about how much to invest in genome coverage and sequencing depth, as well as choices of values for adjustable bioinformatics parameters. To empirically explore the importance of these “simple” methodological decisions, we generated two independent sequencing libraries for the same 142 individual lake whitefish (Coregonus clupeaformis) using a nextRAD RRL approach: (1) a larger number of loci at low sequencing depth based on a 9mer (library A); and (2) fewer loci at higher sequencing depth based on a 10mer (library B). The fish were selected from populations with different levels of expected genetic subdivision. Each library was analyzed using the STACKS pipeline followed by three types of population structure assessment (F_ST, DAPC and ADMIXTURE) with iterative increases in the stringency of sequencing depth and missing data requirements, as well as more specific a priori population maps. Library B was always able to resolve strong population differentiation in all three types of assessment regardless of the selected parameters, largely due to retention of more loci in analyses. In contrast, library A produced more variable results; increasing the minimum sequencing depth threshold (-m) resulted in a reduced number of retained loci, and therefore lost resolution at high -m values for F_ST and ADMIXTURE, but not DAPC. When detecting fine population differentiation, the population map influenced the number of loci and missing data, which generated artefacts in all downstream analyses tested. Similarly, when examining fine scale population subdivision, library B was robust to changing parameters but library A lost resolution depending on the parameter set. We used library B to examine actual subdivision in our study populations. All three types of analysis found complete subdivision among populations in Lake Huron, ON and Dore Lake, SK, Canada using 10,640 SNP loci. Weak population subdivision was detected in Lake Huron with fish from sites in the north-west, Search Bay, North Point and Hammond Bay, showing slight differentiation. Overall, we show that apparently simple decisions about library construction and bioinformatics parameters can have important impacts on the interpretation of population subdivision. Although potentially more costly on a per-locus basis, early investment in striking a balance between the number of loci and sequencing effort is well worth the reduced genomic coverage for population genetics studies. More conservative stringency settings on STACKS parameters lead to a final dataset that was more consistent and robust when examining both weak and strong population differentiation. Overall, we recommend that researchers approach “simple” methodological decisions with caution, especially when working on non-model species for the first time.

Tracking transcriptomic responses to endogenous and exogenous variation in cetaceans in the Southern California Bight

Marine wildlife populations are adapted to survive in highly dynamic environments. However, identifying the effects of endogenous versus exogenous variables on marine mammal physiology remains a substantial challenge in part because of the logistical constraints that limit the collection of physiological data in free-ranging animals. Measuring genome-wide gene expression is one minimally invasive method that can be used to elucidate how free-ranging cetaceans' physiological responses shift with changing environmental conditions or demographic states, i.e. reproductive status and maturity. We identified transcriptomic differences among bottlenose dolphins (Tursiops truncatus) from the Southern California Bight using RNAseq data from the skin of 75 individuals to examine gene expression associated with sex, pregnancy status, sea surface temperature, geographic location and ecotype. We identified transcriptomic variation between two genetically distinct ecotypes as well as variation related to environmental conditions among groups that exhibit little evidence of genetic divergence. Specifically, we found differential expression of genes associated with structural development, cellular starvation and immune response. Sex and pregnancy status explained a small proportion of the observed variation, in contrast to sea surface temperature, which explained a substantial amount of transcriptomic variation. However, these measured variables did not account for all of the differential expression observed between ecotypes and among geographically distinct groups. Additional research is needed to identify other endogenous or exogenous factors that may be contributing to observed transcriptomic differences among ecotypes.

Contaminant Exposure Linked to Cellular and Endocrine Biomarkers in Southern California Bottlenose Dolphins

Cetaceans in the Southern California Bight (SCB) are exposed to high levels of halogenated organic contaminants (HOCs), which have previously been linked to impaired reproductive health and immune responses. We used a combination of molecular tools to examine the potential physiological impacts of HOC exposure in two bottlenose dolphin ( Tursiops truncatus) ecotypes in the SCB. We quantified 25 HOCs in the blubber of 22 biopsies collected from males between 2012 and 2016. We then analyzed genome-wide gene expression in skin using RNA-sequencing and measured blubber testosterone to compare HOC exposure with cellular and endocrine biomarkers. We found high levels of HOCs in both ecotypes with significantly higher total polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs), tris(4-chlorophenyl)methanol (TCPMOH), and chlordane-related compounds in the coastal ecotype versus the offshore ecotype. We found evidence of PBDE bioaccumulation in both ecotypes, however, the pattern of bioaccumulation or endocrine disruption for other HOCs was different between the ecotypes, suggesting potential endocrine disruption in the coastal ecotype. We also observed correlations between HOCs and gene coexpression networks enriched for xenobiotic metabolism, hormone metabolism, and immune response that could indicate cellular effects from HOC exposure. By integrating measurements of HOC load with both transcriptome profiling and endocrine biomarkers, our approach provides insight into HOC exposure and potential impacts on wild cetacean health in southern California.

Genome-wide association study of an unusual dolphin mortality event reveals candidate genes for susceptibility and resistance to cetacean morbillivirus

Infectious diseases are significant demographic and evolutionary drivers of populations, but studies about the genetic basis of disease resistance and susceptibility are scarce in wildlife populations. Cetacean morbillivirus (CeMV) is a highly contagious disease that is increasing in both geographic distribution and incidence, causing unusual mortality events (UME) and killing tens of thousands of individuals across multiple cetacean species worldwide since the late 1980s. The largest CeMV outbreak in the Southern Hemisphere reported to date occurred in Australia in 2013, where it was a major factor in a UME, killing mainly young Indo-Pacific bottlenose dolphins (Tursiops aduncus). Using cases (nonsurvivors) and controls (putative survivors) from the most affected population, we carried out a genome-wide association study to identify candidate genes for resistance and susceptibility to CeMV. The genomic data set consisted of 278,147,988 sequence reads and 35,493 high-quality SNPs genotyped across 38 individuals. Association analyses found highly significant differences in allele and genotype frequencies among cases and controls at 65 SNPs, and Random Forests conservatively identified eight as candidates. Annotation of these SNPs identified five candidate genes (MAPK8, FBXW11, INADL, ANK3 and ACOX3) with functions associated with stress, pain and immune responses. Our findings provide the first insights into the genetic basis of host defence to this highly contagious disease, enabling the development of an applied evolutionary framework to monitor CeMV resistance across cetacean species. Biomarkers could now be established to assess potential risk factors associated with these genes in other CeMV-affected cetacean populations and species. These results could also possibly aid in the advancement of vaccines against morbilliviruses.

Guidelines for planning genomic assessment and monitoring of locally adaptive variation to inform species conservation

Identifying and monitoring locally adaptive genetic variation can have direct utility for conserving species at risk, especially when management may include actions such as translocations for restoration, genetic rescue, or assisted gene flow. However, genomic studies of local adaptation require careful planning to be successful, and in some cases may not be a worthwhile use of resources. Here, we offer an adaptive management framework to help conservation biologists and managers decide when genomics is likely to be effective in detecting local adaptation, and how to plan assessment and monitoring of adaptive variation to address conservation objectives. Studies of adaptive variation using genomic tools will inform conservation actions in many cases, including applications such as assisted gene flow and identifying conservation units. In others, assessing genetic diversity, inbreeding, and demographics using selectively neutral genetic markers may be most useful. And in some cases, local adaptation may be assessed more efficiently using alternative approaches such as common garden experiments. Here, we identify key considerations of genomics studies of locally adaptive variation, provide a road map for successful collaborations with genomics experts including key issues for study design and data analysis, and offer guidelines for interpreting and using results from genomic assessments to inform monitoring programs and conservation actions.

Genome-wide scan reveals signatures of selection related to pollution adaptation in non-model estuarine Atlantic killifish (Fundulus heteroclitus)

In many human-altered ecosystems, organisms are increasingly faced with more diverse and complex environmental stressors and pollutant mixtures, to which the adaptations necessary to survive exposure are likely to be numerous and varied. Improving our understanding of the molecular mechanisms that underlie complex polygenic adaptations in natural settings requires significant toxicological, biochemical, physiological, and genomic data rarely available for non-model organisms. Here, we build upon two decades of study of adaptation to anthropogenic pollutants in a population of Atlantic killifish (Fundulus heteroclitus) that inhabits the creosote-contaminated Atlantic Wood Industries Superfund (AW) site on the Elizabeth River, Virginia in the United States. To better understand the genotypes that underlie previously characterized resistance to PCBs and PAHs, we performed Restriction site-Associated DNA sequencing (RADseq) on killifish from AW and two relatively clean reference sites (King's Creek-KC, and Mains Creek-MC). Across the genome, we analyzed over 83,000 loci and 12,000 single nucleotide polymorphisms (SNPs). Shared across both comparisons of killifish from polluted (AW) and relatively unpolluted (KC and MC) sites, we found eight genomic regions with smoothed FST values significantly (p < 0.001) elevated above background. Using the recently published F. heteroclitus reference genome, we identified candidate genes in these significant regions involved in the AHR pathway (e.g. AIP, ARNT1c), as well as genes relating to cardiac structure and function. These genes represent both previously characterized and potentially novel molecular adaptations involved with various aspects of resistance to these environmental toxins.

Immune effects of the neurotoxins ciguatoxins and brevetoxins

Ciguatoxins (CTXs) and brevetoxins (PbTxs) are phycotoxins that can accumulate along the marine food chain and thus cause seafood poisoning in humans, namely "ciguatera fish poisoning" (CFP) and "neurotoxic shellfish poisoning" (NSP), respectively. CFP is characterized by early gastrointestinal symptoms and typical sensory disorders (paraesthesia, pain, pruritus and cold dysaesthesia), which can persist several weeks and, in some cases, several months or years. NSP is considered a mild form of CFP with similar but less severe symptoms. After inhaled exposure, PbTxs can also cause respiratory tract irritation in healthy subjects and asthma exacerbations in predisposed subjects, whose respiratory functions may be disrupted for several days following PbTx inhalation. Mechanistically, it is well established that CTX- or PbTx-induced disturbances are primarily mainly due to voltage-gated sodium channel activation in sensory and motor peripheral nervous system. However, little is known about the pathophysiology or a potential individual susceptibility to long lasting effects of CFP/NSP. In addition to their action on the nervous system, PbTxs and CTXs were also shown to exert effects on the immune system. However, their role in the pathophysiology of syndromes induced by CTX or PbTx exposure is poorly documented. The aim of this review is to inventory the literature thus far on the inflammatory and immune effects of PbTxs and CTXs.

Gamete compatibility genes in mammals: candidates, applications and a potential path forward

Fertilization represents a critical stage in biology, where successful alleles of a previous generation are shuffled into new arrangements and subjected to the forces of selection in the next generation. Although much research has been conducted on how variation in morphological and behavioural traits lead to variation in fertilization patterns, surprisingly little is known about fertilization at a molecular level, and specifically about how genes expressed on the sperm and egg themselves influence fertilization patterns. In mammals, several genes have been identified whose products are expressed on either the sperm or the egg, and which influence the fertilization process, but the specific mechanisms are not yet known. Additionally, in 2014 an interacting pair of proteins was identified: 'Izumo' on the sperm, and 'Juno' on the egg. With the identification of these genes comes the first opportunity to understand the molecular aspects of fertilization in mammals, and to identify how the genetic characteristics of these genes influence fertilization patterns. Here, we review recent progress in our understanding of fertilization and gamete compatibility in mammals, which should provide a helpful guide to researchers interested in untangling the molecular mechanisms of fertilization and the resulting impacts on population biology and evolutionary processes.

Genomic Methods Take the Plunge: Recent Advances in High-Throughput Sequencing of Marine Mammals

The dramatic increase in the application of genomic techniques to non-model organisms (NMOs) over the past decade has yielded numerous valuable contributions to evolutionary biology and ecology, many of which would not have been possible with traditional genetic markers. We review this recent progression with a particular focus on genomic studies of marine mammals, a group of taxa that represent key macroevolutionary transitions from terrestrial to marine environments and for which available genomic resources have recently undergone notable rapid growth. Genomic studies of NMOs utilize an expanding range of approaches, including whole genome sequencing, restriction site-associated DNA sequencing, array-based sequencing of single nucleotide polymorphisms and target sequence probes (e.g., exomes), and transcriptome sequencing. These approaches generate different types and quantities of data, and many can be applied with limited or no prior genomic resources, thus overcoming one traditional limitation of research on NMOs. Within marine mammals, such studies have thus far yielded significant contributions to the fields of phylogenomics and comparative genomics, as well as enabled investigations of fitness, demography, and population structure. Here we review the primary options for generating genomic data, introduce several emerging techniques, and discuss the suitability of each approach for different applications in the study of NMOs.