Diversity, distribution, and significance of transposable elements in the genome of the only selfing hermaphroditic vertebrate Kryptolebias marmoratus

Transposable Element Tissue-Specific Response to Temperature Stress in the Stenothermal Fish Puntius tetrazona

Ray-finned fish represent a very interesting group of vertebrates comprising a variety of organisms living in different aquatic environments worldwide. In the case of stenothermal fish, thermal fluctuations are poorly tolerated, thus ambient temperature represents a critical factor. In this paper, we considered the tiger barb Puntius tetrazona, a freshwater fish belonging to the family Cyprinidae, living at 21-28 °C. We analyzed the available RNA-Seq data obtained from specimens exposed at 27 °C and 13 °C to investigate the transcriptional activity of transposable elements (TEs) and genes encoding for proteins involved in their silencing in the brain, gill, and liver. TEs are one of the tools generating genetic variability that underlies biological evolution, useful for organisms to adapt to environmental changes. Our findings highlighted a different response of TEs in the three analyzed tissues. While in the brain and gill, no variation in TE transcriptional activity was observed, a remarkable increase at 13 °C was recorded in the liver. Moreover, the transcriptional analysis of genes encoding proteins involved in TE silencing such as heterochromatin formation, the NuRD complex, and the RISC complex (e.g., AGO and GW182 proteins) highlighted their activity in the hepatic tissue. Overall, our findings suggested that this tissue is a target organ for this kind of stress, since TE activation might regulate the expression of stress-induced genes, leading to a better response of the organism to temperature changes. Therefore, this view corroborates once again the idea of a potential role of TEs in organism rapid adaptation, hence representing a promising molecular tool for species resilience.

Additive and non-additive epigenetic signatures of natural hybridisation between fish species with different mating systems

Hybridization is a major source of evolutionary innovation. In plants, epigenetic mechanisms can help to stabilize hybrid genomes and contribute to reproductive isolation, but the relationship between genetic and epigenetic changes in animal hybrids is unclear. We analysed the relationship between genetic background and methylation patterns in natural hybrids of two genetically divergent fish species with different mating systems, Kryptolebias hermaphroditus (self-fertilizing) and K. ocellatus (outcrossing). Co-existing parental species displayed highly distinct genetic (SNPs) and methylation patterns (37,000 differentially methylated cytosines). Hybrids had predominantly intermediate methylation patterns (88.5% of the sites) suggesting additive effects, as expected from hybridization between genetically distant species. The large number of differentially methylated cytosines between hybrids and parental species (n = 5,800) suggests that hybridization may play a role in increasing genetic and epigenetic variation. Although most of the observed epigenetic variation was additive and had a strong genetic component, we also found a small percentage of non-additive, potentially stochastic, methylation differences that might act as an evolutionary bet-hedging strategy and increase fitness under environmental instability.

Living in Temporary Ponds Loading Giant Genomes: The Neotropical Annual Killifish Genus Austrolebias as New Outstanding Evolutionary Model

The term Annual killifish describes a short-lived and amazing group of vertebrates inhabiting temporary ponds exposed to an extremely variable environment during its short lifespan in South America and Africa, leading to the death of the entire adult population during the dry season. Austrolebias is a specious genus of the family Rivulidae, with ∼58 currently recognized species, extensively distributed in the temperate Neotropical region. Herein, we reviewed different aspects of the evolutionary biology with emphasis on the genome dynamic linked to the burst speciation process in this genus. Austrolebias constitutes an excellent model to study the genomic evolutionary processes underlying speciation events, since all the species of this genus analyzed so far share an unusually large genome size, with an average DNA content of 5.95 ± 0.45 picograms per diploid cell (mean C-value of about 2.98 pg). The drastic nuclear DNA–increasing would be associated with a considerable proportion of transposable elements (TEs) found in the Austrolebias genomes. The genomic proportion of the moderately repetitive DNA in the A. charrua genome represents approximately twice (45%) the amount of the repetitive components of the highly related sympatric and syntopic rivulinae taxon Cynopoecilus melanotaenia (25%), as well as from other rivulids and actinopterygian fish. These events could explain the great genome instability, the high genetic diversity, chromosome variability, as well as the morphological diversity in species of Austrolebias. Thus, species of this genus represent new model systems linking different evolutionary processes: drastic genome increase, massive TEs genomic representation, high chromosome instability, occurrence of natural hybridization between sister species, and burst speciation events.

Genome of the Rio Pearlfish (Nematolebias whitei), a bi-annual killifish model for Eco-Evo-Devo in extreme environments

The Rio Pearlfish, Nematolebias whitei, is a bi-annual killifish species inhabiting seasonal pools in the Rio de Janeiro region of Brazil that dry twice per year. Embryos enter dormant diapause stages in the soil, waiting for the inundation of the habitat which triggers hatching and commencement of a new life cycle. Rio Pearlfish represents a convergent, independent origin of annualism from other emerging killifish model species. While some transcriptomic datasets are available for Rio Pearlfish, thus far, a sequenced genome has been unavailable. Here, we present a high quality, 1.2 Gb chromosome-level genome assembly, genome annotations, and a comparative genomic investigation of the Rio Pearlfish as representative of a vertebrate clade that evolved environmentally cued hatching. We show conservation of 3D genome structure across teleost fish evolution, developmental stages, tissues, and cell types. Our analysis of mobile DNA shows that Rio Pearlfish, like other annual killifishes, possesses an expanded transposable element profile with implications for rapid aging and adaptation to harsh conditions. We use the Rio Pearlfish genome to identify its hatching enzyme gene repertoire and the location of the hatching gland, a key first step in understanding the developmental genetic control of hatching. The Rio Pearlfish genome expands the comparative genomic toolkit available to study convergent origins of seasonal life histories, diapause, and rapid aging phenotypes. We present the first set of genomic resources for this emerging model organism, critical for future functional genetic, and multiomic explorations of “Eco-Evo-Devo” phenotypes of resilience and adaptation to extreme environments.

Against the Odds: Hybrid Zones between Mangrove Killifish Species with Different Mating Systems

Different mating systems are expected to affect the extent and direction of hybridization. Due to the different levels of sexual conflict, the weak inbreeder/strong outbreeder (WISO) hypothesis predicts that gametes from self-incompatible (SI) species should outcompete gametes from self-compatible (SC) ones. However, other factors such as timing of selfing and unilateral incompatibilities may also play a role on the direction of hybridization. In addition, differential mating opportunities provided by different mating systems are also expected to affect the direction of introgression in hybrid zones involving outcrossers and selfers. Here, we explored these hypotheses with a unique case of recent hybridization between two mangrove killifish species with different mating systems, Kryptolebias ocellatus (obligately outcrossing) and K. hermaphroditus (predominantly self-fertilizing) in two hybrid zones in southeast Brazil. Hybridization rates were relatively high (~20%), representing the first example of natural hybridization between species with different mating systems in vertebrates. All F1 individuals were sired by the selfing species. Backcrossing was small, but mostly asymmetrical with the SI parental species, suggesting pattern commonly observed in plant hybrid zones with different mating systems. Our findings shed light on how contrasting mating systems may affect the direction and extent of gene flow between sympatric species, ultimately affecting the evolution and maintenance of hybrid zones.

Genome-wide DNA methylation of the liver reveals delayed effects of early-life exposure to 17-α-ethinylestradiol in the self-fertilizing mangrove rivulus

Organisms exposed to endocrine disruptors in early life can show altered phenotype later in adulthood. Although the mechanisms underlying these long-term effects remain poorly understood, an increasing body of evidence points towards the potential role of epigenetic processes. In the present study, we exposed hatchlings of an isogenic lineage of the self-fertilizing fish mangrove rivulus for 28 days to 4 and 120 ng/L of 17-α-ethinylestradiol. After a recovery period of 140 days, reduced representation bisulphite sequencing (RRBS) was performed on the liver in order to assess the hepatic genome-wide methylation landscape. Across all treatment comparisons, a total of 146 differentially methylated fragments (DMFs) were reported, mostly for the group exposed to 4 ng/L, suggesting a non-monotonic effect of EE2 exposure. Gene ontology analysis revealed networks involved in lipid metabolism, cellular processes, connective tissue function, molecular transport and inflammation. The highest effect was reported for nipped-B-like protein B (NIPBL) promoter region after exposure to 4 ng/L EE2 (+ 21.9%), suggesting that NIPBL could be an important regulator for long-term effects of EE2. Our results also suggest a significant role of DNA methylation in intergenic regions and potentially in transposable elements. These results support the ability of early exposure to endocrine disruptors of inducing epigenetic alterations during adulthood, providing plausible mechanistic explanations for long-term phenotypic alteration. Additionally, this work demonstrates the usefulness of isogenic lineages of the self-fertilizing mangrove rivulus to better understand the biological significance of long-term alterations of DNA methylation by diminishing the confounding factor of genetic variability.

Transposable Elements and Stress in Vertebrates: An Overview

Since their identification as genomic regulatory elements, Transposable Elements (TEs) were considered, at first, molecular parasites and later as an important source of genetic diversity and regulatory innovations. In vertebrates in particular, TEs have been recognized as playing an important role in major evolutionary transitions and biodiversity. Moreover, in the last decade, a significant number of papers has been published highlighting a correlation between TE activity and exposition to environmental stresses and dietary factors. In this review we present an overview of the impact of TEs in vertebrate genomes, report the silencing mechanisms adopted by host genomes to regulate TE activity, and finally we explore the effects of environmental and dietary factor exposures on TE activity in mammals, which is the most studied group among vertebrates. The studies here reported evidence that several factors can induce changes in the epigenetic status of TEs and silencing mechanisms leading to their activation with consequent effects on the host genome. The study of TE can represent a future challenge for research for developing effective markers able to detect precocious epigenetic changes and prevent human diseases.

Genomic and physiological mechanisms underlying skin plasticity during water to air transition in an amphibious fish

The terrestrial radiation of vertebrates required changes in skin that resolved the dual demands of maintaining a mechanical and physiological barrier while also facilitating ion and gas transport. Using the amphibious killifish Kryptolebias marmoratus, we found that transcriptional regulation of skin morphogenesis was quickly activated upon air exposure (1 h). Rapid regulation of cell-cell adhesion complexes and pathways that regulate stratum corneum formation was consistent with barrier function and mechanical reinforcement. Unique blood vessel architecture and regulation of angiogenesis likely supported cutaneous respiration. Differences in ionoregulatory transcripts and ionocyte morphology were correlated with differences in salinity acclimation and resilience to air exposure. Evolutionary analyses reinforced the adaptive importance of these mechanisms. We conclude that rapid plasticity of barrier, respiratory and ionoregulatory functions in skin evolved to support the amphibious lifestyle of K. marmoratus; similar processes may have facilitated the terrestrial radiation of other contemporary and ancient fishes.

Environmental enrichment induces intergenerational behavioural and epigenetic effects on fish

Parental effects influence offspring phenotypes through pre- and post-natal routes but little is known about their molecular basis, and therefore their adaptive significance. Epigenetic modifications, which control gene expression without changes in the DNA sequence and are influenced by the environment, may contribute to parental effects. We investigated the effects of environmental enrichment on the behaviour, metabolic rate and brain DNA methylation patterns of parents and offspring of the highly inbreed mangrove killifish (Kryptolebias marmoratus). Parental fish reared in enriched environments had lower cortisol levels, lower metabolic rates and were more active and neophobic than those reared in barren environments. They also differed in 1,854 methylated cytosines (DMCs). Offspring activity and neophobia were determined by the parental environment. Among the DMCs of the parents, 98 followed the same methylation patterns in the offspring, three of which were significantly influenced by parental environments irrespective of their own rearing environment. Our results suggest that parental environment influences the behaviour and, to some extent, the brain DNA methylation patterns of the offspring.

The reference genome of the selfing fish Kryptolebias hermaphroditus: Identification of phases I and II detoxification genes

The selfing fish Kryptolebias hermaphroditus has unique reproductive system for self-fertilization, making genetically homozygous offsprings. Here, we report on high density genetic map-based genome assembly for the K. hermaphroditus Panama line (PanRS). The numbers of scaffolds were 5212 and the genome was 683,992,224 bp (N50 = 27.45 Mb). The length of anchored scaffold onto 24 linkage groups was 652,231,070 bp (95.3% of genome) with 0.01% of the gap and 39.33% of GC content and complete Benchmarking Universal Single-Copy Orthologs value was 96.6%. The numbers of annotated genes were 36,756 (average gene length 1368 bp) with the GC content of 54.1%. To examine the difference between the two sister species in the genus Kryptolebias, we compared the genomes of K. hermaphroditus PanRS and Kryptolebias marmoratus PAN line on the composition of transposable elements. To demonstrate applications of genome library, phase I and II detoxification related gene families have been analyzed, and compared the syntenies containing loci of CYP and GST genes on linkage groups. This K. hermaphroditus genome information will be helpful for a better understanding on genome-wide mechanistic view of detoxification and antioxidant-related genes over evolution in the view of fish environmental ecotoxicology.

The genome of the Java medaka (Oryzias javanicus): Potential for its use in marine molecular ecotoxicology

The Java medaka (Oryzias javanicus) is distributed in tropical brackish water and is considered as an ecotoxicological experimental organism for assessing diverse pollutions and global climate change effects in the ocean. In this study, we sequenced and assembled the genome of O. javanicus using the Oxford Nanopore technique and anchored the scaffolds to the 24 genetic linkage map of a sister species Oryzias melastigma. The assembled genome consisted of 773 scaffolds including 24 LG-based scaffolds, and the estimated genome length was 846.3 Mb (N50 = 19.3 Mb), containing 24,498 genes. As detoxification processes are crucial in aquatic organisms, antioxidant-related genes including glutathione S-transferases, superoxide dismutase, catalase, and glutathione peroxidase were identified in this study. In the genome of O. javanicus, a total of 21 GSTs, 4 SODs, 1 CAT, and 7 GPxs were identified and showed high similarities between sister species O. melastigma and Oryzias latipes. In addition, despite having 8 classes of cytosolic GSTs family, medaka showed no presence of GST pi and sigma classes, which are predominantly found in carp and salmon, but not in neoteleostei. This study adds another set to genome-library of Oryzias spp. and is a useful resource for better understanding of the molecular ecotoxicology.

Extensive hybridization and past introgression between divergent lineages in a quasi-clonal hermaphroditic fish: Ramifications for species concepts and taxonomy

Extreme inbreeding is expected to reduce the incidence of hybridization, serving as a prezygotic barrier. Mangrove rivulus is a small killifish that reproduces predominantly by self-fertilization, producing highly homozygous lines throughout its geographic range. The Bahamas and Caribbean are inhabited by two highly diverged phylogeographic lineages of mangrove rivulus, Kryptolebias marmoratus and a 'Central clade' closely related to K. hermaphroditus from Brazil. The two lineages are largely allopatric, but recently were found in syntopy on San Salvador, Bahamas, where a single hybrid was reported. To better characterize the degree of hybridization and the possibility of secondary introgression, here we conducted a detailed genetic analysis of the contact zone on San Salvador. Two mixed populations were identified, one of which contained sexually mature hybrids. The distribution of heterozygosity at diagnostic microsatellite loci in hybrids showed that one of these hybrids was an immediate offspring from the K. marmoratus x Central clade cross, whereas the remaining five hybrids were products of reproduction by self-fertilization for 1-3 generations following the initial cross. Two hybrids had mitochondrial haplotypes of K. marmoratus and the remaining four hybrids had a haplotype of the Central clade, indicating that crosses go in both directions. In hybrids, alleles of parental lineages were represented in equal proportions suggesting lack of recent backcrossing to either of the parental lineages. However, sympatric populations of two lineages were less diverged than allopatric populations, consistent with introgression. Results are discussed in terms of applicability of the biological species concept for isogenic, effectively clonal, organisms.

Spatial and temporal expression pattern of sex-related genes in ovo-testis of the self-fertilizing mangrove killifish (Kryptolebias marmoratus)

In vertebrates, sex determination and differentiation comprehend a fine balance between female and male factors, leading the bipotential anlage to develop towards ovary or testis, respectively. Nevertheless, the mangrove killifish, (Kryptolebias marmoratus) a simultaneous hermaphroditic species, could overcome those antagonistic pathways and evolved to develop and maintain reproductively active ovarian and testicular tissues in the same organ. Morphological and mRNA localization analyzes of developing and adult gonads demonstrate that genes related to testis (dmrt1 and amh) and ovary differentiation (foxl2 and sox9a) follow the same expression pattern observed in gonochoristic species, thus functioning as two independent organs. In addition, Amh expression patterns make it a strong candidate for initiation of the formation and maintenance of the testicular tissue in the hermaphroditic gonad. Differently from described so far, foxl3 seems to have an important role in oogenesis as well as spermatogenesis and gonadal structure.

Construction of High-Resolution RAD-Seq Based Linkage Map, Anchoring Reference Genome, and QTL Mapping of the Sex Chromosome in the Marine Medaka Oryzias melastigma

Medaka (Oryzias sp.) is an important fish species in ecotoxicology and considered as a model species due to its biological features including small body size and short generation time. Since Japanese medaka Oryzias latipes is a freshwater species with access to an excellent genome resource, the marine medaka Oryzias melastigma is also applicable for the marine ecotoxicology. In genome era, a high-density genetic linkage map is a very useful resource in genomic research, providing a means for comparative genomic analysis and verification of de novo genome assembly. In this study, we developed a high-density genetic linkage map for O. melastigma using restriction-site associated DNA sequencing (RAD-seq). The genetic map consisted of 24 linkage groups with 2,481 single nucleotide polymorphism (SNP) markers. The total map length was 1,784 cM with an average marker space of 0.72 cM. The genetic map was integrated with the reference-assisted chromosome assembly (RACA) of O. melastigma, which anchored 90.7% of the assembled sequence onto the linkage map. The values of complete Benchmarking Universal Single-Copy Orthologs were similar to RACA assembly but N50 (23.74 Mb; total genome length 779.4 Mb; gap 5.29%) increased to 29.99 Mb (total genome length 778.7 Mb; gap 5.2%). Using MapQTL analysis with SNP markers, we identified a major quantitative trait locus for sex traits on the Om10. The integration of the genetic map with the reference genome of marine medaka will serve as a good resource for studies in molecular toxicology, genomics, CRISPR/Cas9, and epigenetics.

New Insights on the Evolution of Genome Content: Population Dynamics of Transposable Elements in Flies and Humans

Understanding the abundance, diversity, and distribution of TEs in genomes is crucial to understand genome structure, function, and evolution. Advances in whole-genome sequencing techniques, as well as in bioinformatics tools, have increased our ability to detect and analyze the transposable element content in genomes. In addition to reference genomes, we now have access to population datasets in which multiple individuals within a species are sequenced. In this chapter, we highlight the recent advances in the study of TE population dynamics focusing on fruit flies and humans, which represent two extremes in terms of TE abundance, diversity, and activity. We review the most recent methodological approaches applied to the study of TE dynamics as well as the new knowledge on host factors involved in the regulation of TE activity. In addition to transposition rates, we also focus on TE deletion rates and on the selective forces that affect the dynamics of TEs in genomes.

Early-life exposure to the endocrine disruptor 17-α-ethinylestradiol induces delayed effects in adult brain, liver and ovotestis proteomes of a self-fertilizing fish

Early-life represents a critically sensitive window to endocrine disrupting chemicals, potentially leading to long-term repercussions on the phenotype later in life. The mechanisms underlying this phenomenon, referred to as the Developmental Origins of Health and Disease (DOHaD), are still poorly understood. To gain molecular understanding of these effects, we exposed mangrove rivulus (Kryptolebias marmoratus) for 28 days post hatching (dph) to 4 and 120 ng/L 17-α-ethinylestradiol, a model xenoestrogen. After 28 days, fish were raised for 140 days in clean water and we performed quantitative label-free proteomics on brain, liver and ovotestis of 168 dph adults. A total of 820, 888 and 420 proteins were robustly identified in the brain, liver and ovotestis, respectively. Effects of 17-α-ethinylestradiol were tissue and dose-dependent: a total of 31, 51 and 18 proteins were differentially abundant at 4 ng/L in the brain, liver and ovotestis, respectively, compared to 20, 25 and 39 proteins at 120 ng/L. Our results suggest that estrogen-responsive pathways, such as lipid metabolism, inflammation, and the innate immune system were affected months after the exposure. In addition, the potential perturbation of S-adenosylmethionine metabolism encourages future studies to investigate the role of DNA methylation in mediating the long-term effects of early-life exposures. SIGNIFICANCE: The Developmental Origins of Health and Disease (DOHaD) states that early life stages of humans and animals are sensitive to environmental stressors and can develop health issues later in life, even if the stress has ceased. Molecular mechanisms supporting DOHaD are still unclear. The mangrove rivulus is a new fish model species naturally reproducing by self-fertilization, making it possible to use isogenic lineages in which all individuals are highly homozygous. This species therefore permits to strongly reduce the confounding factor of genetic variability in order to investigate the effects of environmental stress on the phenotype. After characterizing the molecular phenotype of brain, liver and ovotestis, we obtained true proteomic reaction norms of these three organs in adults after early life stages have been exposed to the common endocrine disruptor 17-α-ethinylestradiol (EE2). Our study demonstrates long-term effects of early-life endocrine disruption at the proteomic level in diverse estrogen-responsive pathways 5 months after the exposure. The lowest tested and environmentally relevant concentration of 4 ng/L had the highest impact on the proteome in brain and liver, highlighting the potency of endocrine disruptors at low concentrations.

Whole-genome sequencing reveals the extent of heterozygosity in a preferentially self-fertilizing hermaphroditic vertebrate

The mangrove rivulus, Kryptolebias marmoratus, is one of only two self-fertilizing hermaphroditic fish species and inhabits mangrove forests. While selfing can be advantageous, it reduces heterozygosity and decreases genetic diversity. Studies using microsatellites found that there are variable levels of selfing among populations of K. marmoratus, but overall, there is a low rate of outcrossing and, therefore, low heterozygosity. In this study, we used whole-genome data to assess the levels of heterozygosity in different lineages of the mangrove rivulus and infer the phylogenetic relationships among those lineages. We sequenced whole genomes from 15 lineages that were completely homozygous at microsatellite loci and used single nucleotide polymorphisms (SNPs) to determine heterozygosity levels. More variation was uncovered than in studies using microsatellite data because of the resolution of full genome sequencing data. Moreover, missense polymorphisms were found most often in genes associated with immune function and reproduction. Inferred phylogenetic relationships suggest that lineages largely group by their geographic distribution. The use of whole-genome data provided further insight into genetic diversity in this unique species. Although this study was limited by the number of lineages that were available, these data suggest that there is previously undescribed variation within lineages of K. marmoratus that could have functional consequences and (or) inform us about the limits to selfing (e.g., genetic load, accumulation of deleterious mutations) and selection that might favor the maintenance of heterozygosity. These results highlight the need to sequence additional individuals within and among lineages.

The genome of the marine medaka Oryzias melastigma

Marine medaka (Oryzias melastigma) is considered to be a useful fish model for marine and estuarine ecotoxicology studies and has good potential for field-based population genomics because of its geographical distribution in Asian estuarine and coastal areas. In this study, we present the first whole-genome draft of O. melastigma. The genome assembly consists of 8,602 scaffolds (N50 = 23.737 Mb) and a total genome length of 779.4 Mb. A total of 23,528 genes were predicted, and 12,670 gene families shared with three teleost species (Japanese medaka, mangrove killifish and zebrafish) were identified. Genome analyses revealed that the O. melastigma genome is highly heterozygous and contains a large number of repeat sequences. This assembly represents a useful genomic resource for fish scientists.

Identification of 74 cytochrome P450 genes and co-localized cytochrome P450 genes of the CYP2K, CYP5A, and CYP46A subfamilies in the mangrove killifish Kryptolebias marmoratus

BACKGROUND

The mangrove killifish Kryptolebias marmoratus is the only vertebrate that reproduces by self-fertilizing and is an important model species in genetics and marine ecotoxicology. Using whole-genome and transcriptome sequences, we identified all members of the cytochrome P450 (CYP) family in this model teleost and compared them with those of other teleosts.

RESULTS

A total of 74 cytochrome P450 genes and one pseudogene were identified in K. marmoratus. Phylogenetic analysis indicated that the CYP genes in clan 2 were most expanded, while synteny analysis with other species showed orthologous relationships of CYP subfamilies among teleosts. In addition to the CYP2K expansions, five tandem duplicated gene copies of CYP5A were observed. These features were unique to K. marmoratus.

CONCLUSIONS

These results shed a light on CYP gene evolution, particularly the co-localized CYP2K, CYP5A, and CYP46A subfamilies in fish. Future studies of CYP expression could identify specific endogenous and exogenous environmental factors that triggered the evolution of tandem CYP duplication in K. marmoratus.

Genome-wide identification of 31 cytochrome P450 (CYP) genes in the freshwater rotifer Brachionus calyciflorus and analysis of their benzo[α]pyrene-induced expression patterns

While marine invertebrate cytochrome P450 (CYP) genes and their roles in detoxification mechanisms have been studied, little information is available regarding freshwater rotifer CYPs and their functions. Here, we used genomic sequences and RNA-seq databases to identify 31 CYP genes in the freshwater rotifer Brachionus calyciflorus. The 31 Bc-CYP genes with a few tandem duplications were clustered into CYP 2, 3, 4, mitochondrial, and 46 clans with two marine rotifers Brachionus plicatilis and Brachionus koreanus. To understand the molecular responses of these 31 Bc-CYP genes, we also examined their expression patterns in response to benzo[α]pyrene (B[α]P). Three Bc-CYP genes (Bc-CYP3044B3, Bc-CYP3049B4, Bc-CYP3049B6) were significantly upregulated (P<0.05) in response to B[α]P, suggesting that these CYP genes can be involved in detoxification in response to B[α]P exposure. These genes might be useful as biomarkers of B[α]P exposure in B. calyciflorus. Overall, our findings expand the repertoire of known CYPs and shed light on their potential roles in xenobiotic detoxification in rotifers.

Genome-wide identification of 52 cytochrome P450 (CYP) genes in the copepod Tigriopus japonicus and their B[α]P-induced expression patterns

Cytochrome P450s (CYPs) are enzymes with a heme-binding domain that are found in all living organisms. CYP enzymes have important roles associated with detoxification of xenobiotics and endogenous compounds (e.g. steroids, fatty acids, and hormones). Although CYP enzymes have been reported in several invertebrates, including insects, little is known about copepod CYPs. Here, we identified the entire repertoire of CYP genes (n=52) from whole genome and transcriptome sequences of the benthic copepod Tigriopus japonicus, including a tandem duplication (CYP3026A3, CYP3026A4, CYP3026A5), and examined patterns of gene expression over various developmental stages and in response to benzo[α]pyrene (B[α]P) exposure. Through phylogenetic analysis, the 52 T. japonicus CYP genes were assigned to five distinct clans: CYP2 (22 genes), CYP3 (19 genes), CYP4 (two genes), CYP20 (one gene), and mitochondrial (eight genes). Developmental stage and gender-specific expression patterns of the 52 T. japonicus CYPs were analyzed. CYP3022A1 was constitutively expressed during all developmental stages. CYP genes in clans 2 and 3 were induced in response to B[α]P, suggesting that these differentially modulated CYP transcripts are likely involved in defense against exposure to B[α]P and other pollutants. This study enhances our understanding of the repertoire of CYP genes in copepods and of their potential role in development and detoxification in copepods.

Nine co-localized cytochrome P450 genes of the CYP2N, CYP2AD, and CYP2P gene families in the mangrove killifish Kryptolebias marmoratus genome: Identification and expression in response to B[α]P, BPA, OP, and NP

The CYP2 genes are the largest and most diverse cytochrome P450 (CYP) subfamily in vertebrates. We have identified nine co-localized CYP2 genes (∼55kb) in a new cluster in the genome of the highly resilient ecotoxicological fish model Kryptolebias marmoratus. Molecular characterization, temporal and tissue-specific expression pattern, and response to xenobiotics of these genes were examined. The CYP2 gene clusters were characterized and designated CYP2N22-23, CYP2AD12, and CYP2P16-20. Gene synteny analysis confirmed that the cluster in K. marmoratus is similar to that found in other teleost fishes, including zebrafish. A gene duplication event with diverged catalytic function was observed in CYP2AD12. Moreover, a high level of divergence in expression was observed among the co-localized genes. Phylogeny of the cluster suggested an orthologous relationship with similar genes in zebrafish and Japanese medaka. Gene expression analysis showed that CYP2P19 and CYP2N20 were consecutively expressed throughout embryonic development, whereas CYP2P18 was expressed in all adult tissues, suggesting that members of each CYP2 gene family have different physiological roles even though they are located in the same cluster. Among endocrine-disrupting chemicals (EDCs), benzo[α]pyrene (B[α]P) induced expression of CYP2N23, bisphenol A (BPA) induced CYP2P18 and CYP2P19, and 4-octylphenol (OP) induced CYP2AD12, but there was no significant response to 4-nonylphenol (NP), implying differential catalytic roles of the enzyme. In this paper, we identify and characterize a CYP2 gene cluster in the mangrove killifish K. marmoratus with differing catalytic roles toward EDCs. Our findings provide insights on the roles of nine co-localized CYP2 genes and their catalytic functions for better understanding of chemical-biological interactions in fish.

Adverse effects of BDE-47 on in vivo developmental parameters, thyroid hormones, and expression of hypothalamus-pituitary-thyroid (HPT) axis genes in larvae of the self-fertilizing fish Kryptolebias marmoratus

2,2',4,4'-tetrabromodiphenylether (BDE-47) is known to have the potential to disrupt the thyroid endocrine system in fishes due to its structural similarity to the thyroid hormones triiodothyronine (T3) and thyroxine (T4). However, the effects of BDE-47 on thyroid function in fishes remain unclear. In this study, abnormal development (e.g. deformity, hemorrhaging) and an imbalance in thyroid hormone (TH) homeostasis was shown in the early developmental stages of the mangrove killifish Kryptolebias marmoratus in response to BDE-47 exposure. To examine the thyroid endocrinal effect of BDE-47 exposure in mangrove killifish K. marmoratus larvae, transcript levels of genes involved in TH homeostasis and hypothalamus-pituitary-thyroid (HPT) axis-related genes were measured. The expression of thyroid hormone metabolism-related genes (e.g. deiodinases, UGT1ab) and HPT axis-related genes was up-regulated and there were significant changes in TH levels (P < 0.05) in response to BDE-47 exposure. This study provides insights into the regulation of TH homeostasis at the transcriptional level and provides a better understanding on the potential impacts of BDE-47 on the thyroid endocrine system of fishes.