Genetic map for the Atlantic killifish (Fundulus heteroclitus) reference genome assembly (Fundulus_heteroclitus-3.0.2)
Ontology highlight
ABSTRACT:
SUBMITTER: Andrew Whitehead
PROVIDER: S-BSST163 | biostudies-other |
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Killifish (Fundulus heteroclitus) are widely distributed among different aquatic environments where they demonstrate an impressive range of highly-plastic and locally adaptive phenotypes. High-throughput sequencing has begun to unravel the mechanisms and evolutionary history of these interesting features by establishing relationships in the genotype-phenotype map. However, some genotype-phenotype analyses require a higher order of contiguity than what initial scaffolded (fragmented genome assembly where contigs have been assemble into scaffolds) genome assemblies can provide. Here, we used 5,685 high-quality RAD-Seq markers from a single mapping family to order 84% of the scaffolded genome assembly to 24 chromosomes. This serves to: 1) expand the killifish genomic toolkit, 2) estimate genome-wide recombination rates, and 3) compare genome synteny to humans and other fishes. After initially building our map, we found that the selection of thresholds for sequence data filtration highly impacted scaffold placement in the map. We outline each step of the approach that dramatically improved our map to help guide others toward more effective linkage mapping for genome assembly. Our final map supports strong conservation of genomic synteny among closely related fish species and reveals previously described chromosomal rearrangements between more distantly related clades. However, we also commonly found minor scaffold misorientations in F. heteroclitus and in other assemblies, suggesting that further mapping (such as optical mapping) is necessary for finer scale resolution of genome structure. Lastly, we discuss the problems that would be expected from misoriented/unplaced scaffolds and stress the importance of a quality mapped genome as a key feature for further investigating population and comparative genomic questions with F. heteroclitus and other taxa.
| S-EPMC6723127 | biostudies-literature
Project description:The Atlantic killifish (Fundulus heteroclitus) is an ideal model species to study physiological and toxicological adaptations to stressors. Killifish inhabiting the PCB-contaminated Superfund site in New Bedford Harbor, MA (NBH) have evolved resistance to toxicity and activation of the aryl hydrocarbon receptor (AHR) signaling pathway after exposure to PCBs and other AHR agonists. Until recently, a lack of genomic information has limited efforts to understand the molecular mechanisms underlying environmental adaptation to stressors. The advent of high throughput sequencing has facilitated an unbiased assessment of coding as well as non-coding RNAs in any species of interest. Among non-coding RNAs, microRNAs (miRNAs) are important regulators of gene expression and play crucial roles in development and physiology. The objective of this study is to catalog the miRNAs in killifish and determine their expression patterns in the embryos from contaminated (NBH) and pristine (Scorton Creek, MA (SC)) sites. Embryos from NBH and SC were collected daily from 1 to 15 days post-fertilization and RNA from pooled samples from each site was sequenced using SOLiD sequencing. We obtained 7.5 and 11 million raw reads from pooled SC and NBH samples, respectively. Analysis of the sequencing data identified 216 conserved mature miRNA sequences that are expressed during development. Using the draft killifish genome, we retrieved the miRNA precursor sequences. Based on the capacity of these putative precursor sequences to form the characteristic hairpin loop (assessed using RNAfold), we identified 197 conserved miRNA sequences in the genome.
2015-07-20 | GSE70953 | GEO
Project description:The non-migratory killifish Fundulus heteroclitus inhabits clean and polluted environments interspersed throughout its range along the Atlantic coast of North America. Several populations of this species have successfully adapted to environments contaminated with toxic aromatic hydrocarbon pollutants such as polychlorinated biphenyls (PCBs). Previous studies suggest that the mechanism of resistance to these and other "dioxin-like compounds" (DLCs) may involve reduced signaling through the aryl hydrocarbon receptor (AHR) pathway. Here we investigated gene diversity and evidence for positive selection at three AHR-related loci (AHR1, AHR2, AHRR) in F. heteroclitus by comparing alleles from seven locations ranging over 600 km along the northeastern US, including extremely polluted and reference estuaries, with a focus on New Bedford Harbor (MA, USA), a PCB Superfund site, and nearby reference sites.We identified 98 single nucleotide polymorphisms within three AHR-related loci among all populations, including synonymous and nonsynonymous substitutions. Haplotype distributions were spatially segregated and F-statistics suggested strong population genetic structure at these loci, consistent with previous studies showing strong population genetic structure at other F. heteroclitus loci. Genetic diversity at these three loci was not significantly different in contaminated sites as compared to reference sites. However, for AHR2 the New Bedford Harbor population had significant FST values in comparison to the nearest reference populations. Tests for positive selection revealed ten nonsynonymous polymorphisms in AHR1 and four in AHR2. Four nonsynonymous SNPs in AHR1 and three in AHR2 showed large differences in base frequency between New Bedford Harbor and its reference site. Tests for isolation-by-distance revealed evidence for non-neutral change at the AHR2 locus.Together, these data suggest that F. heteroclitus populations in reference and polluted sites have similar genetic diversity, providing no evidence for strong genetic bottlenecks for populations in polluted locations. However, the data provide evidence for genetic differentiation among sites, selection at specific nucleotides in AHR1 and AHR2, and specific AHR2 SNPs and haplotypes that are associated with the PCB-resistant phenotype in the New Bedford Harbor population. The results suggest that AHRs, and especially AHR2, may be important, recurring targets for selection in local adaptation to dioxin-like aromatic hydrocarbon contaminants.
| S-EPMC3899389 | biostudies-other
Project description:BackgroundThe most toxic aromatic hydrocarbon pollutants are categorized as dioxin-like compounds (DLCs) to which extreme tolerance has evolved independently and contemporaneously in (at least) four populations of Atlantic killifish (Fundulus heteroclitus). Surprisingly, the magnitude and phenotype of DLC tolerance is similar among these killifish populations that have adapted to varied, but highly aromatic hydrocarbon-contaminated urban/industrialized estuaries of the US Atlantic coast. Multiple tolerant and neighboring sensitive killifish populations were compared with the expectation that genetic loci associated with DLC tolerance would be revealed.ResultsSince the aryl hydrocarbon receptor (AHR) pathway partly or fully mediates DLC toxicity in vertebrates, single nucleotide polymorphisms (SNPs) from 42 genes associated with the AHR pathway were identified to serve as targeted markers. Wild fish (N = 36/37) from four highly tolerant killifish populations and four nearby sensitive populations were genotyped using 59 SNP markers. Similar to other killifish population genetic analyses, strong genetic differentiation among populations was detected, consistent with isolation by distance models. When DLC-sensitive populations were pooled and compared to pooled DLC-tolerant populations, multi-locus analyses did not distinguish the two groups. However, pairwise comparisons of nearby tolerant and sensitive populations revealed high differentiation among sensitive and tolerant populations at these specific loci: AHR 1 and 2, cathepsin Z, the cytochrome P450s (CYP1A and 3A30), and the NADH dehydrogenase subunits. In addition, significant shifts in minor allele frequency were observed at AHR2 and CYP1A loci across most sensitive/tolerant pairs, but only AHR2 exhibited shifts in the same direction across all pairs.ConclusionsThe observed differences in allelic composition at the AHR2 and CYP1A SNP loci were identified as significant among paired sensitive/tolerant populations of Atlantic killifish with multiple statistical tests. The genetic patterns reported here lend support to the argument that AHR2 and CYP1A play a role in the adaptive response to extreme DLC contamination. Additional functional assays are required to isolate the exact mechanism of DLC tolerance.
| S-EPMC4029433 | biostudies-literature
Project description:Understanding molecular mechanisms of toxicity is facilitated by experimental manipulations, such as disruption of function by gene targeting, that are especially challenging in non-standard model species with limited genomic resources. While loss-of-function approaches have included gene knock-down using morpholino-modified oligonucleotides and random mutagenesis using mutagens or retroviruses, more recent approaches include targeted mutagenesis using zinc finger nuclease (ZFN), transcription activator-like effector nuclease (TALENs) and clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 technology. These latter methods provide more accessible opportunities to explore gene function in non-traditional model species. To facilitate evaluation of toxic mechanisms for important categories of aryl hydrocarbon pollutants, whose actions are known to be receptor mediated, we used ZFN and CRISPR-Cas9 approaches to generate aryl hydrocarbon receptor 2a (AHR2a) and AHR2b gene mutations in Atlantic killifish (Fundulus heteroclitus) embryos. This killifish is a particularly valuable non-traditional model, with multiple paralogs of AHR whose functions are not well characterized. In addition, some populations of this species have evolved resistance to toxicants such as halogenated aromatic hydrocarbons. AHR-null killifish will be valuable for characterizing the role of the individual AHR paralogs in evolved resistance, as well as in normal development. We first used five-finger ZFNs targeting exons 1 and 3 of AHR2a. Subsequently, CRISPR-Cas9 guide RNAs were designed to target regions in exon 2 and 3 of AHR2a and AHR2b. We successfully induced frameshift mutations in AHR2a exon 3 with ZFN and CRISPR-Cas9 guide RNAs, with mutation frequencies of 10% and 16%, respectively. In AHR2b, mutations were induced using CRISPR-Cas9 guide RNAs targeting sites in both exon 2 (17%) and exon 3 (63%). We screened AHR2b exon 2 CRISPR-Cas9-injected embryos for off-target effects in AHR paralogs. No mutations were observed in closely related AHR genes (AHR1a, AHR1b, AHR2a, AHRR) in the CRISPR-Cas9-injected embryos. Overall, our results demonstrate that targeted genome-editing methods are efficient in inducing mutations at specific loci in embryos of a non-traditional model species, without detectable off-target effects in paralogous genes.
| S-EPMC4272816 | biostudies-literature
Project description:Atlantic killifish inhabiting polluted sites along the east coast of the U.S. have evolved resistance to toxic effects of contaminants. One such contaminated site is the Acushnet River estuary, near New Bedford Harbor (NBH), Massachusetts, which is characterized by very high PCB concentrations in the sediments and in the tissues of resident killifish. Though killifish at this site appear to be thriving, the metabolic costs of survival in a highly contaminated environment are not well understood. In this study we compared the hepatic metabolite profiles of resistant (NBH) and sensitive populations (Scorton Creek (SC), Sandwich, MA) using a targeted metabolomics approach in which polar metabolites were extracted from adult fish livers and quantified. Our results revealed differences in the levels of several metabolites between fish from the two sites. The majority of these metabolites are associated with one-carbon metabolism, an important pathway that supports multiple physiological processes including DNA and protein methylation, nucleic acid biosynthesis and amino acid metabolism. We measured the gene expression of DNA methylation (DNA methyltransferase 1, dnmt1) and demethylation genes (Ten-Eleven Translocation (TET) genes) in the two populations, and observed lower levels of dnmt1 and higher levels of TET gene expression in the NBH livers, suggesting possible differences in DNA methylation profiles. Consistent with this, the two populations differed significantly in the levels of 5-methylcytosine and 5-hydroxymethylcytosine nucleotides. Overall, our results suggest that the unique hepatic metabolite signatures observed in NBH and SC reflect the adaptive mechanisms for survival in their respective habitats.
| S-EPMC6246827 | biostudies-literature
Project description:The Atlantic killifish (Fundulus heteroclitus), native to estuarine areas of the Atlantic coast of the United States, has become a valuable ecotoxicological model as a result of its ability to acclimate to rapid environmental changes and adapt to polluted habitats. MicroRNAs (miRNAs) are highly conserved small RNAs that regulate gene expression and play critical roles in stress responses in a variety of organisms. Global miRNA expression in killifish and the potential roles miRNA have in environmental acclimation have yet to be characterized. Accordingly, we profiled miRNA expression in killifish gill for the first time and identified a small group of highly expressed, well-conserved miRNAs as well as 16 novel miRNAs not yet identified in other organisms. Killifish respond to large fluctuations in salinity with rapid changes in gene expression and protein trafficking to maintain osmotic balance, followed by a secondary phase of gene and protein expression changes that enable remodeling of the gills. Arsenic, a major environmental toxicant, was previously shown to inhibit gene expression responses in killifish gill, as well the ability of killifish to acclimate to a rapid increase in salinity. Thus, we examined the individual and combined effects of salinity and arsenic on miRNA expression in killifish gill. Using small RNA sequencing, we identified 270 miRNAs expressed in killifish, and found that miR-135b was differentially expressed in response to arsenic and at 24?h following transfer to salt water. Predicted targets of miR-135b are involved in ion transport, cell motility and migration, GTPase mediated signal transduction and organelle assembly. Consistent with previous studies of these two environmental stressors, we found a significant interaction (i.e., arsenic dependent salinity effect), whereby killifish exposed to arsenic exhibited an opposite response in miR-135b expression at 24?h post hyperosmotic challenge compared to controls. By examining mRNA expression of predicted miRNA targets during salinity acclimation and arsenic exposure, we found that miR-135b targets were significantly more likely to decrease during salinity acclimation than non-targets. Our identification of a significant interaction effect of arsenic and salinity on miR-135b expression supports the hypothesis that arsenic alters upstream regulators of stress response networks, which may adversely affect the killifish response to osmotic stress. The characterization of miRNAs in this ecotoxicological model will be a valuable resource for future studies investigating the role of miRNAs in response to environmental stress.
| S-EPMC6298807 | biostudies-literature
Project description:Exposure of developing fish to polycyclic aromatic hydrocarbons (PAHs) and halogenated aromatic hydrocarbons (HAHs) results in a suite of defects including cardiac malformation, pericardial and yolk sac edema, craniofacial defects, and hemorrhaging. Several populations of Atlantic killifish or mummichog (Fundulus heteroclitus) on the Atlantic coast of the United States are resistant to the developmental and acute toxicity caused by PAHs and HAHs; this has made Fundulus a valuable model for studying aryl hydrocarbon sensitivity and adaptation. In order to further increase the utility of Fundulus, better understanding of the components of the molecular pathways governing aryl hydrocarbon response in Fundulus is required. The aryl hydrocarbon receptor (AHR) is known to mediate many of the toxic responses to PAHs and HAHs. A single AHR has been identified in mammals, but Fundulus has two AHRs and their relative roles are not clear. In the current study, translation-blocking and splice-junction morpholino gene knockdown was used to determine the roles of AHR1 and AHR2 in mediating cardiac teratogenesis induced by beta-naphthoflavone (BNF), benzo[k]fluoranthene (BkF), and 3,3',4,4',5-pentachlorobiphenyl (PCB-126). Here we report that AHR2 and not AHR1 knockdown resulted in rescue of teratogenicity induced by BNF, BkF, and PCB-126. These data demonstrate that AHR2 is the primary mediator of cardiac teratogenesis caused by multiple aryl hydrocarbons in Fundulus and suggest that suppression of the AHR pathway through modulation of AHR2 is a plausible mechanism for PAH resistance in adapted fish. Additionally, this is the first reported use of splice-junction morpholinos in Fundulus.
| S-EPMC2901414 | biostudies-literature
Project description:The Atlantic Wood Industries Superfund site (AWI) on the Elizabeth River in Portsmouth, VA is heavily contaminated with polycyclic aromatic hydrocarbons (PAHs) from a wood treatment facility. Atlantic killifish, or mummichog (Fundulus heteroclitus), at this Superfund site are exposed to very high concentrations of several carcinogens. In this study, we measured PAH concentrations in both fish tissues and sediments. Concurrently, we assessed different aspects of genotoxicity in the killifish exposed in situ. Both sediment and tissue PAH levels were significantly higher in AWI samples, relative to a reference site, but the chemistry profile was different between sediments and tissues. Killifish at AWI exhibited higher levels of DNA damage compared to reference fish, as measured via the flow cytometric method (FCM), and the damage was consistent with sediment PAH concentrations. Covalent binding of benzo[a]pyrene (BaP) metabolites to DNA, as measured via LC-MS/MS adduct detection methods, were also elevated and could be partially responsible for the DNA damage. Using similar LC-MS/MS methods, we found no evidence that oxidative DNA adducts had a role in observed genotoxicity.
| S-EPMC3203518 | biostudies-literature
Project description:Sentinel species such as the Atlantic killifish (Fundulus heteroclitus) living in urban waterways can be used as toxicological models to understand impacts of environmental metabolism disrupting compound (MDC) exposure on both wildlife and humans. Exposure to MDCs is associated with increased risk of metabolic syndrome, including impaired lipid and glucose homeostasis, adipogenesis, appetite control, and basal metabolism. MDCs are ubiquitous in the environment, including in aquatic environments. New Bedford Harbor (NBH), Massachusetts is polluted with polychlorinated biphenyls (PCBs), and, as we show for the first time, tin (Sn). PCBs and organotins are ligands for two receptor systems known to regulate lipid homeostasis, the aryl hydrocarbon receptor (AHR) and the peroxisome proliferator-activated receptors (PPARs), respectively. In the current study, we compared lipid homeostasis in laboratory-reared killifish from NBH (F2) and a reference location (Scorton Creek, Massachusetts; F1 and F2) to evaluate how adaptation to local conditions may influence responses to MDCs. Adult killifish from each population were exposed to 3,3',4,4',5-pentachlorobiphenyl (PCB126, dioxin-like), 2,2',4,4',5,5'-hexachlorobiphenyl (PCB153, non-dioxin-like), or tributyltin (TBT, a PPARγ ligand) by a single intraperitoneal injection and analyzed after 3 days. AHR activation was assessed by measuring cyp1a mRNA expression. Lipid homeostasis was evaluated phenotypically by measuring liver triglycerides and organosomatic indices, and at the molecular level by measuring the mRNA expression of pparg and ppara and a target gene for each receptor. Acute MDC exposure did not affect phenotypic outcomes. However, overall NBH killifish had higher liver triglycerides and adiposomatic indices than SC killifish. Both season and population were significant predictors of the lipid phenotype. Acute MDC exposure altered hepatic gene expression only in male killifish from SC. PCB126 exposure induced cyp1a and pparg, whereas PCB153 exposure induced ppara. TBT exposure did not induce ppar-dependent pathways. Comparison of lipid homeostasis in two killifish populations extends our understanding of how MDCs act on fish and provides a basis to infer adaptive benefits of these differences in the wild.
| S-EPMC6544361 | biostudies-literature