The Genome of the Self-Fertilizing Mangrove Rivulus Fish, Kryptolebias marmoratus: A Model for Studying Phenotypic Plasticity and Adaptations to Extreme Environments

Neutral Genetic Diversity in Mixed Mating Systems

BACKGROUND/OBJECTIVES

Systems of reproduction differ with respect to the magnitude of neutral genetic diversity maintained in a population. In particular, the partitioning of reproductive organisms into mating types and regular inbreeding have long been recognized as key factors that influence effective population number. Here, a range of reproductive systems are compared with respect to the maintenance of neutral genetic diversity. This study addresses full gonochorism, full hermaphroditism, androdioecy (male and hermaphroditic reproductives), and gynodioecy (female and hermaphroditic reproductives).

METHODS

Coalescence theory is used to determine the level of diversity maintained under each mating system considered.

RESULTS

For each mating system, the nature of the dependence of the level of neutral diversity on inbreeding depression, sex-specific viability, and other factors is described. In particular, the models account for the effects of sex-specific viability on the evolutionarily stable sex ratio and the collective contribution of each mating type (sex) to the offspring generation.

CONCLUSIONS

Within the context of conservation biology, population genetic and quantitative genetic theory has addressed the determination of the target minimum effective population size. In contrast, this study proposes and explores a summary statistic (a ratio of effective numbers) as a means of characterizing the context in which evolution occurs.

Warming, stochastic diel thermal fluctuations affect physiological performance and gill plasticity in an amphibious mangrove fish

Natural temperature variation in many marine ecosystems is stochastic and unpredictable, and climate change models indicate that this thermal irregularity is likely to increase. Temperature acclimation may be more challenging when conditions are highly variable and stochastic, and there is a need for empirical physiological data in these thermal environments. Using the hermaphroditic, amphibious mangrove rivulus (Kryptolebias marmoratus), we hypothesized that compared with regular, warming diel thermal fluctuations, stochastic warm fluctuations would negatively affect physiological performance. To test this, we acclimated fish to: (1) non-stochastic and (2) stochastic thermal fluctuations with a similar thermal load (27-35°C), and (3) a stable/consistent control temperature at the low end of the cycle (27°C). We determined that fecundity was reduced in both cycles, with reproduction ceasing in stochastic thermal environments. Fish acclimated to non-stochastic thermal cycles had growth rates lower than those of control fish. Exposure to warm, fluctuating cycles did not affect emersion temperature, and only regular diel cycles modestly increased critical thermal tolerance. We predicted that warm diel cycling temperatures would increase gill surface area. Notably, fish acclimated to either thermal cycle had a reduced gill surface area and increased intralamellar cell mass when compared with control fish. This decreased gill surface area with warming contrasts with what is observed for exclusively aquatic fish and suggests a preparatory gill response for emersion in these amphibious fish. Collectively, our data reveal the importance of considering stochastic thermal variability when studying the effects of temperature on fishes.

Genome sequencing, comparative analysis, and gene expression responses of cytochrome P450 genes in Oryzias curvinotus provide insights into environmental adaptation

The mangrove fish (Oryzias curvinotus) serves as a model for researching environmental adaptation and sexual development. To further such research, we sequenced and assembled a high-quality 842 Mb reference genome for O. curvinotus. Comparative genomic analysis revealed 891 expanded gene families, including significantly expanded cytochrome P450 (CYP) detoxification genes known to be involved in xenobiotic defense. We identified 69 O. curvinotus CYPs (OcuCYPs) across 18 families and 10 clans using multiple methods. Extensive RNA-seq and qPCR analysis demonstrated diverse spatiotemporal expression patterns of OcuCYPs by developmental stage, tissue type, sex, and pollutant exposure (17β-estradiol (E2) and testosterone (MT)). Many OcuCYPs exhibited sexual dimorphism in gonads, suggesting reproductive roles in steroidogenesis, while their responsiveness to model toxicants indicates their importance in environmental adaptation through enhanced detoxification. Pathway analysis highlighted expanded CYP genes in arachidonic acid metabolism, drug metabolism, and steroid hormone biosynthesis. This chromosome-level genomic resource provides crucial biological insights to elucidate the functional roles of expanded CYPs in environmental adaptation, sexual development, early life history, and conservation in the anthropogenically impacted mangrove habitats of O. curvinotus. It also enables future ecotoxicology research leveraging O. curvinotus as a pollution sentinel species.

Genetic tools for the study of the mangrove killifish, Kryptolebias marmoratus, an emerging vertebrate model for phenotypic plasticity

Kryptolebias marmoratus (Kmar), a teleost fish of the order Cyprinodontiformes, has a suite of unique phenotypes and behaviors not observed in other fishes. Many of these phenotypes are discrete and highly plastic-varying over time within an individual, and in some cases reversible. Kmar and its interfertile sister species, K. hermaphroditus, are the only known self-fertile vertebrates. This unusual sexual mode has the potential to provide unique insights into the regulation of vertebrate sexual development, and also lends itself to genetics. Kmar is easily adapted to the lab and requires little maintenance. However, its internal fertilization and small clutch size limits its experimental use. To support Kmar as a genetic model, we compared alternative husbandry techniques to maximize recovery of early cleavage-stage embryos. We find that frequent egg collection enhances yield, and that protease treatment promotes the greatest hatching success. We completed a forward mutagenesis screen and recovered several mutant lines that serve as important tools for genetics in this model. Several will serve as useful viable recessive markers for marking crosses. Importantly, the mutant kissylips lays embryos at twice the rate of wild-type. Combining frequent egg collection with the kissylips mutant background allows for a substantial enhancement of early embryo yield. These improvements were sufficient to allow experimental analysis of early development and the successful mono- and bi-allelic targeted knockout of an endogenous tyrosinase gene with CRISPR/Cas9 nucleases. Collectively, these tools will facilitate modern developmental genetics in this fascinating fish, leading to future insights into the regulation of plasticity.

Exploring the cobia (Rachycentron canadum) genome: unveiling putative male heterogametic regions and identification of sex-specific markers

BACKGROUND

Cobia (Rachycentron canadum) is the only member of the Rachycentridae family and exhibits considerable sexual dimorphism in growth rate. Sex determination in teleosts has been a long-standing basic biological question, and the molecular mechanisms of sex determination/differentiation in cobia are completely unknown.

RESULTS

Here, we reported 2 high-quality, chromosome-level annotated male and female cobia genomes with assembly sizes of 586.51 Mb (contig/scaffold N50: 86.0 kb/24.3 Mb) and 583.88 Mb (79.9 kb/22.5 Mb), respectively. Synteny inference among perciform genomes revealed that cobia and the remora Echeneis naucrates were sister groups. Further, whole-genome resequencing of 31 males and 60 females, genome-wide association study, and sequencing depth analysis identified 3 short male-specific regions within a 10.7-kb continuous genomic region on male chromosome 18, which hinted at an undifferentiated sex chromosome system with a putative XX/XY mode of sex determination in cobia. Importantly, the only 2 genes within/between the male-specific regions, epoxide hydrolase 1 (ephx1, renamed cephx1y) and transcription factor 24 (tcf24, renamed ctcf24y), showed testis-specific/biased gene expression, whereas their counterparts cephx1x and ctf24x, located in female chromosome 18, were similarly expressed in both sexes. In addition, male-specific PCR targeting the cephx1y gene revealed that this genomic feature is conserved in cobia populations from Panama, Brazil, Australia, and Japan.

CONCLUSION

The first comprehensive genomic survey presented here is a valuable resource for future studies on cobia population structure and dynamics, conservation, and evolutionary history. Furthermore, it establishes evidence of putative male heterogametic regions with 2 genes playing a potential role in the sex determination of the species, and it provides further support for the rapid evolution of sex-determining mechanisms in teleost fish.

Evidence for selfing in a vertebrate from whole-genome sequencing

A growing number of recent genomic studies report asexual parthenogenetic reproduction in a wide range of taxa, including vertebrate species from the reptile, bird, and fish lineages. Yet, self-fertilization (selfing) has been recorded only in a single vertebrate, the mangrove killifish Kryptolebias marmoratus In cichlid fishes, sex determination is notably diverse and can be influenced by the environment, and sequential hermaphroditism has been reported for some species. Here, we present evidence for a case of facultative selfing in the cichlid fish Benitochromis nigrodorsalis, which is otherwise known as biparentally reproducing ovophilic mouthbrooder from Western Africa. Our laboratory observations revealed that a wild-caught individual produced repeatedly viable offspring in absence of a mating partner. By analyzing genome-wide single-nucleotide polymorphism (SNP) data, we compare that individual and two of its offspring to shed light on its reproductive mode. First, our results confirm uniparental reproduction. Second, overall heterozygosity is reduced in the offspring compared with outbred individuals. Retained maternal heterozygosity in the offspring is ∼51%, which is close to the theoretically expected value of a heterozygosity reduction of 50% by selfing. Heterozygosity patterns along individual chromosomes do not point to alternative parthenogenetic reproductive mechanisms like automixis by terminal or central fusion. Facultative selfing may represent an adaptive strategy ensuring reproduction when mating partners are absent and, hence, contribute to the cichlids' enormous evolutionary success.

Genome size estimation and its associations with body length, chromosome number and evolution in teleost fishes

Precise estimation of genome size (GS) is vital for various genomic studies, such as deciding genome sequencing depth, genome assembly, biodiversity documentation, evolution, genetic disorders studies, duplication events etc. Animal Genome Size Database provides GS of over 2050 fish species, which ranges from 0.35 pg in pufferfish (Tetraodon nigroviridis) to 132.83 pg in marbled lungfish (Protopterus aethiopicus). The GS of majority of the fishes inhabiting waters of Indian subcontinent are still missing. In present study, we estimated GS of 51 freshwater teleost (31 commercially important, 7 vulnerable and 13 ornamental species) that ranged from 0.58 pg in banded gourami (Trichogaster fasciata) to 1.92 pg in scribbled goby (Awaous grammepomus). Substantial variation in GS was observed within the same fish orders (0.64-1.45 pg in cypriniformes, 0.70-1.41 pg in siluriformes and 0.58-1.92 pg in perciformes). We examined the relationship between the GS, chromosome number and body length across all the fishes. Body length was found to be associated with GS, whereas no relationship was noticed between the GS and the chromosome number. The analysis using ancestral information revealed haploid chromosome number 25, 27 and 24 for the most recent common ancestor of cypriniformes, siluriformes and perciformes, respectively. The study led to generation of new records on GS of 43 fish species and revalidated records for 8 species. The finding is valuable resource for further research in the areas of fish genomics, molecular ecology and evolutionary conservation genetics.

Early-life exposure to methylmercury induces reversible behavioral impairments and gene expression modifications in one isogenic lineage of mangrove rivulus fish Kryptolebias marmoratus

Methylmercury (MeHg) is a ubiquitous bioaccumulative neurotoxicant present in aquatic ecosystems. It is known to alter behaviors, sensory functions and learning abilities in fish and other vertebrates. Developmental and early-life stages exposure to MeHg can lead to brain damage with immediate consequences on larvae behavior, but may also induce long term effects in adults after a detoxification period. However, very little is known about developmental origin of behavioral impairment in adults due to early exposure to MeHg. The aim of this study is to assess whether early-life MeHg exposure induces immediate and/or delayed effects on behaviors, related genes expression and DNA methylation (one of epigenetic mechanisms). To reach this goal, newly hatched larvae of mangrove rivulus fish, Kryptolebias marmoratus, were exposed to two sub-lethal concentrations of MeHg (90 μg/L and 135 µg/L) for 7 days, and immediate and delayed effects were assessed respectively in 7 dph (days post-hatching) and 90 dph fish. This species naturally produces isogenic lineages due to its self-fertilizing reproduction system, which is unique among vertebrates. It allows to study how environment stressors can influence organism's phenotype while minimizing genetic variability. As results, both MeHg exposures are associated with a decreased foraging efficiency and thigmotaxis, and a dose-dependent reduction in larvae locomotor activity. Regarding molecular analysis in larvae whole bodies, both MeHg exposures induced significant decreased expression of DNMT3a, MAOA, MeCP2 and NIPBL, and significant increase of GSS, but none of those genes underwent methylation changes in targeted CpGs. None of significant behavioral and molecular impairments observed in 7-dph larvae were found in 90-dph adults, which highlight a distinction between immediate and delayed effects of developmental MeHg exposure. Our results suggest implications of aminergic system and its neurotransmitters, redox/methylation trade-off and possibly other epigenetic mechanisms in MeHg neurotoxicity underlying behavioral alterations in rivulus.

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.

A Transient Hermaphroditic Stage in Early Male Gonadal Development in Little Yellow Croaker, Larimichthys polyactis

Animal taxa show remarkable variability in sexual reproduction, where separate sexes, or gonochorism, is thought to have evolved from hermaphroditism for most cases. Hermaphroditism accounts for 5% in animals, and sequential hermaphroditism has been found in teleost. In this study, we characterized a novel form of the transient hermaphroditic stage in little yellow croaker (Larimichthys polyactis) during early gonadal development. The ovary and testis were indistinguishable from 7 to 40 days post-hatching (dph). Morphological and histological examinations revealed an intersex stage of male gonads between 43 and 80 dph, which consist of germ cells, somatic cells, efferent duct, and early primary oocytes (EPOs). These EPOs in testis degenerate completely by 90 dph through apoptosis yet can be rescued by exogenous 17-β-estradiol. Male germ cells enter the mitotic flourishing stage before meiosis is initiated at 180 dph, and they undergo normal spermatogenesis to produce functional sperms. This transient hermaphroditic stage is male-specific, and the ovary development appears to be normal in females. This developmental pattern is not found in the sister species Larimichthys crocea or any other closely related species. Further examinations of serum hormone levels indicate that the absence of 11-ketotestosterone and elevated levels of 17-β-estradiol delineate the male intersex gonad stage, providing mechanistic insights on this unique phenomenon. Our research is the first report on male-specific transient hermaphroditism and will advance the current understanding of fish reproductive biology. This unique gonadal development pattern can serve as a useful model for studying the evolutionary relationship between hermaphroditism and gonochorism, as well as teleost sex determination and differentiation strategies.

Whole-genome sequencing reveals sex determination and liver high-fat storage mechanisms of yellowstripe goby (Mugilogobius chulae)

As a promising novel marine fish model for future research on marine ecotoxicology as well as an animal model of human disease, the genome information of yellowstripe goby (Mugilogobius chulae) remains unknown. Here we report the first annotated chromosome-level reference genome assembly for yellowstripe goby. A 20.67-cM sex determination region was discovered on chromosome 5 and seven potential sex-determining genes were identified. Based on combined genome and transcriptome data, we identified three key lipid metabolic pathways for high-fat accumulation in the liver of yellowstripe goby. The changes in the expression patterns of MGLL and CPT1 at different development stage of the liver, and the expansion of the ABCA1 gene, innate immune gene TLR23, and TRIM family genes may help in balancing high-fat storage in hepatocytes and steatohepatitis. These results may provide insights into understanding the molecular mechanisms of sex determination and high-fat storage in the liver of marine fishes.

Chromonomer: A Tool Set for Repairing and Enhancing Assembled Genomes Through Integration of Genetic Maps and Conserved Synteny

The pace of the sequencing and computational assembly of novel reference genomes is accelerating. Though DNA sequencing technologies and assembly software tools continue to improve, biological features of genomes such as repetitive sequence as well as molecular artifacts that often accompany sequencing library preparation can lead to fragmented or chimeric assemblies. If left uncorrected, defects like these trammel progress on understanding genome structure and function, or worse, positively mislead this research. Fortunately, integration of additional, independent streams of information, such as a marker-dense genetic map and conserved orthologous gene order from related taxa, can be used to scaffold together unlinked, disordered fragments and to restructure a reference genome where it is incorrectly joined. We present a tool set for automating these processes, one that additionally tracks any changes to the assembly and to the genetic map, and which allows the user to scrutinize these changes with the help of web-based, graphical visualizations. Chromonomer takes a user-defined reference genome, a map of genetic markers, and, optionally, conserved synteny information to construct an improved reference genome of chromosome models: a “chromonome”. We demonstrate Chromonomer’s performance on genome assemblies and genetic maps that have disparate characteristics and levels of quality.

Using asexual vertebrates to study genome evolution and animal physiology: Banded (Fundulus diaphanus) x Common Killifish (F. heteroclitus) hybrid lineages as a model system

Wild, asexual, vertebrate hybrids have many characteristics that make them good model systems for studying how genomes evolve and epigenetic modifications influence animal physiology. In particular, the formation of asexual hybrid lineages is a form of reproductive incompatibility, but we know little about the genetic and genomic mechanisms by which this mode of reproductive isolation proceeds in animals. Asexual lineages also provide researchers with the ability to produce genetically identical individuals, enabling the study of autonomous epigenetic modifications without the confounds of genetic variation. Here, we briefly review the cellular and molecular mechanisms leading to asexual reproduction in vertebrates and the known genetic and epigenetic consequences of the loss of sex. We then specifically discuss what is known about asexual lineages of Fundulus diaphanus x F. heteroclitus to highlight gaps in our knowledge of the biology of these clones. Our preliminary studies of F. diaphanus and F. heteroclitus karyotypes from Porter's Lake (Nova Scotia, Canada) agree with data from other populations, suggesting a conserved interspecific chromosomal arrangement. In addition, genetic analyses suggest that: (a) the same major clonal lineage (Clone A) of F. diaphanus x F. heteroclitus has remained dominant over the past decade, (b) some minor clones have also persisted, (c) new clones may have recently formed, and iv) wild clones still mainly descend from F. diaphanus ♀ x F. heteroclitus ♂ crosses (96% in 2017-2018). These data suggest that clone formation may be a relatively rare, but continuous process, and there are persistent environmental or genetic factors causing a bias in cross direction. We end by describing our current research on the genomic causes and consequences of a transition to asexuality and the potential physiological consequences of epigenetic variation.

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 Role of Stochasticity in the Origin of Epigenetic Variation in Animal Populations

Epigenetic mechanisms such as DNA methylation modulate gene expression in a complex fashion are consequently recognized as among the most important contributors to phenotypic variation in natural populations of plants, animals, and microorganisms. Interactions between genetics and epigenetics are multifaceted and epigenetic variation stands at the crossroad between genetic and environmental variance, which make these mechanisms prominent in the processes of adaptive evolution. DNA methylation patterns depend on the genotype and can be reshaped by environmental conditions, while transgenerational epigenetic inheritance has been reported in various species. On the other hand, DNA methylation can influence the genetic mutation rate and directly affect the evolutionary potential of a population. The origin of epigenetic variance can be attributed to genetic, environmental, or stochastic factors. Generally less investigated than the first two components, variation lacking any predictable order is nevertheless present in natural populations and stochastic epigenetic variation, also referred to spontaneous epimutations, can sustain phenotypic diversity. Here, potential sources of such stochastic epigenetic variability in animals are explored, with a focus on DNA methylation. To this day, quantifying the importance of stochasticity in epigenetic variability remains a challenge. However, comparisons between the mutation and the epimutation rates showed a high level of the latter, suggesting a significant role of spontaneous epimutations in adaptation. The implications of stochastic epigenetic variability are multifold: by affecting development and subsequently phenotype, random changes in epigenetic marks may provide additional phenotypic diversity, which can help natural populations when facing fluctuating environments. In isogenic lineages and asexually reproducing organisms, poor or absent genetic diversity can hence be tolerated. Further implication of stochastic epigenetic variability in adaptation is found in bottlenecked invasive species populations and populations using a bet-hedging strategy.

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.

Phylogenomics investigation of sparids (Teleostei: Spariformes) using high-quality proteomes highlights the importance of taxon sampling

Sparidae (Teleostei: Spariformes) are a family of fish constituted by approximately 150 species with high popularity and commercial value, such as porgies and seabreams. Although the phylogeny of this family has been investigated multiple times, its position among other teleost groups remains ambiguous. Most studies have used a single or few genes to decipher the phylogenetic relationships of sparids. Here, we conducted a thorough phylogenomic analysis using five recently available Sparidae gene-sets and 26 high-quality, genome-predicted teleost proteomes. Our analysis suggested that Tetraodontiformes (puffer fish, sunfish) are the closest relatives to sparids than all other groups used. By analytically comparing this result to our own previous contradicting finding, we show that this discordance is not due to different orthology assignment algorithms; on the contrary, we prove that it is caused by the increased taxon sampling of the present study, outlining the great importance of this aspect in phylogenomic analyses in general.

Divergence, evolution and adaptation in ray-finned fish genomes

With the rapid development of next-generation sequencing technologies and bioinformatics, over 50 ray-finned fish genomes by far have been sequenced with high quality. The genomic work provides abundant genetic resources for deep understanding of divergence, evolution and adaptation in the fish genomes. They are also instructive for identification of candidate genes for functional verification, molecular breeding, and development of novel marine drugs. As an example of other omics data, the Fish-T1K project generated a big database of fish transcriptomes to integrate with these published fish genomes for potential applications. In this review, we highlight the above-mentioned recent investigations and core topics on the ray-finned fish genome research, with a main goal to obtain a deeper understanding of fish biology for theoretical and practical applications.

Genome sequence of walking catfish (Clarias batrachus) provides insights into terrestrial adaptation

Background

Walking catfish (Clarias batrachus) is a freshwater fish capable of air-breathing and locomotion on land. It usually inhabits various low-oxygen habitats, burrows inside the mudflat, and sometimes “walks” to search for suitable environments during summer. It has evolved accessory air-breathing organs for respiring air and corresponding mechanisms to survive in such challenging environments. Thereby, it serves as a great model for understanding adaptations to terrestrial life.

Results

Comparative genomics with channel catfish (Ictalurus punctatus) revealed specific adaptations of C. batrachus in DNA repair, enzyme activator activity, and small GTPase regulator activity. Comparative analysis with 11 non-air-breathing fish species suggested adaptive evolution in gene expression and nitrogenous waste metabolic processes. Further, myoglobin, olfactory receptor related to class A G protein-coupled receptor 1, and sulfotransferase 6b1 genes were found to be expanded in the air-breathing walking catfish genome, with 15, 15, and 12 copies, respectively, compared to non-air-breathing fishes that possess only 1–2 copies of these genes. Additionally, we sequenced and compared the transcriptomes of the gill and the air-breathing organ to characterize the mechanism of aerial respiration involved in elastic fiber formation, oxygen binding and transport, angiogenesis, ion homeostasis and acid-base balance. The hemoglobin genes were expressed dramatically higher in the air-breathing organ than in the gill of walking catfish.

Conclusions

This study provides an important genomic resource for understanding the adaptive mechanisms of walking catfish to terrestrial environments. It is possible that the coupling of enhanced abilities for oxygen storage and oxygen transport through genomic expansion of myoglobin genes and transcriptomic up-regulation of hemoglobin and angiogenesis-related genes are important components of the molecular basis for adaptation of this aquatic species to terrestrial life.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5355-9) contains supplementary material, which is available to authorized users.

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.

Solid-phase enzyme catalysis of DNA end repair and 3' A-tailing reduces GC-bias in next-generation sequencing of human genomic DNA

The use of next-generation sequencing (NGS) has been instrumental in advancing biological research and clinical diagnostics. To fully utilize the power of NGS, complete, uniform coverage of the entire genome is required. In this study, we identified the primary sources of bias observed in sequence coverage across AT-rich regions of the human genome with existing amplification-free DNA library preparation methods. We have found evidence that a major source of bias is the inefficient processing of AT-rich DNA in end repair and 3' A-tailing, causing under-representation of extremely AT-rich regions. We have employed immobilized DNA modifying enzymes to catalyze end repair and 3' A-tailing reactions, to notably reduce the GC bias observed with existing library construction methods.

DNA methylation in adults and during development of the self-fertilizing mangrove rivulus, Kryptolebias marmoratus

In addition to genetic variation, epigenetic mechanisms such as DNA methylation might make important contributions to heritable phenotypic diversity in populations. However, it is often difficult to disentangle the contributions of genetic and epigenetic variation to phenotypic diversity. Here, we investigated global DNA methylation and mRNA expression of the methylation-associated enzymes during embryonic development and in adult tissues of one natural isogenic lineage of mangrove rivulus fish, Kryptolebias marmoratus. Being the best-known self-fertilizing hermaphroditic vertebrate affords the opportunity to work with genetically identical individuals to examine, explicitly, the phenotypic effects of epigenetic variance. Using the LUminometric Methylation Assay (LUMA), we described variable global DNA methylation at CpG sites in adult tissues, which differed significantly between hermaphrodite ovotestes and male testes (79.6% and 87.2%, respectively). After fertilization, an immediate decrease in DNA methylation occurred to 15.8% in gastrula followed by re-establishment to 70.0% by stage 26 (liver formation). Compared to zebrafish, at the same embryonic stages, this reprogramming event seems later, deeper, and longer. Furthermore, genes putatively encoding DNA methyltransferases (DNMTs), Ten-Eleven Translocation (TET), and MeCP2 proteins showed specific regulation in adult gonad and brain, and also during early embryogenesis. Their conserved domains and expression profiles suggest that these proteins play important roles during reproduction and development. This study raises questions about mangrove rivulus' peculiar reprogramming period in terms of epigenetic transmission and physiological adaptation of individuals to highly variable environments. In accordance with the general-purpose genotype model, epigenetic mechanisms might allow for the expression of diverse phenotypes among genetically identical individuals. Such phenotypes might help to overcome environmental challenges, making the mangrove rivulus a valuable vertebrate model for ecological epigenetic studies. The mangrove rivulus, Kryptolebias marmoratus, is the best-known self-fertilizing hermaphroditic vertebrate that allows to work with genetically identical individuals to examine, explicitly, the phenotypic effects of epigenetic variance. The reprogramming event is later, more dramatic and longer than in other described vertebrates. High evolutionary conservation and expression patterns of DNMT, TET, and MeCP2 proteins in K. marmoratus suggest biological roles for each member in gametogenesis and development.

Skeletal stiffening in an amphibious fish out of water is a response to increased body weight

Terrestrial animals must support their bodies against gravity, while aquatic animals are effectively weightless because of buoyant support from water. Given this evolutionary history of minimal gravitational loading of fishes in water, it has been hypothesized that weight-responsive musculoskeletal systems evolved during the tetrapod invasion of land and are thus absent in fishes. Amphibious fishes, however, experience increased effective weight when out of water - are these fishes responsive to gravitational loading? Contrary to the tetrapod-origin hypothesis, we found that terrestrial acclimation reversibly increased gill arch stiffness (∼60% increase) in the amphibious fish Kryptolebias marmoratus when loaded normally by gravity, but not under simulated microgravity. Quantitative proteomics analysis revealed that this change in mechanical properties occurred via increased abundance of proteins responsible for bone mineralization in other fishes as well as in tetrapods. Type X collagen, associated with endochondral bone growth, increased in abundance almost ninefold after terrestrial acclimation. Collagen isoforms known to promote extracellular matrix cross-linking and cause tissue stiffening, such as types IX and XII collagen, also increased in abundance. Finally, more densely packed collagen fibrils in both gill arches and filaments were observed microscopically in terrestrially acclimated fish. Our results demonstrate that the mechanical properties of the fish musculoskeletal system can be fine-tuned in response to changes in effective body weight using biochemical pathways similar to those in mammals, suggesting that weight sensing is an ancestral vertebrate trait rather than a tetrapod innovation.

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 Austrofundulus limnaeus offers insights into extreme vertebrate stress tolerance and embryonic development

Background

The annual killifish Austrofundulus limnaeus inhabits ephemeral ponds in northern Venezuela, South America, and is an emerging extremophile model for vertebrate diapause, stress tolerance, and evolution. Embryos of A. limnaeus regularly experience extended periods of desiccation and anoxia as a part of their natural history and have unique metabolic and developmental adaptations. Currently, there are limited genomic resources available for gene expression and evolutionary studies that can take advantage of A. limnaeus as a unique model system.

Results

We describe the first draft genome sequence of A. limnaeus. The genome was assembled de novo using a merged assembly strategy and was annotated using the NCBI Eukaryotic Annotation Pipeline. We show that the assembled genome has a high degree of completeness in genic regions that is on par with several other teleost genomes. Using RNA-seq and phylogenetic-based approaches, we identify several candidate genes that may be important for embryonic stress tolerance and post-diapause development in A. limnaeus. Several of these genes include heat shock proteins that have unique expression patterns in A. limnaeus embryos and at least one of these may be under positive selection.

Conclusion

The A. limnaeus genome is the first South American annual killifish genome made publicly available. This genome will be a valuable resource for comparative genomics to determine the genetic and evolutionary mechanisms that support the unique biology of annual killifishes. In a broader context, this genome will be a valuable tool for exploring genome-environment interactions and their impacts on vertebrate physiology and evolution.

Electronic supplementary material

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

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.

Evolution of the sex ratio and effective number under gynodioecy and androdioecy

We address the evolution of effective number of individuals under androdioecy and gynodioecy. We analyze dynamic models of autosomal modifiers of weak effect on sex expression. In our zygote control models, the sex expressed by a zygote depends on its own genotype, while in our maternal control models, it depends on the genotype of its maternal parent. Our analysis unifies full multi-dimensional local stability analysis with the Li-Price equation, which for all its heuristic appeal, describes evolutionary change over a single generation. We define a point in the neighborhood of a fixation state from which a single-generation step indicates the asymptotic behavior of the frequency of a modifier allele initiated at an arbitrary point near the fixation state. A concept of heritability appropriate for the evolutionary modification of sex emerges from the Li-Priceframework. We incorporate our theoretical analysis into our previously-developed Bayesian inference framework to develop a new method for inferring the viability of gonochores (males or females) relative to hermaphrodites. Applying this approach to microsatellite data derived from natural populations of the gynodioecious plant Schiedea salicaria and the androdioecious killifish Kryptolebias marmoratus, we find that while female and hermaphrodite S. salicaria appear to have similar viabilities, male K. marmoratus appear to survive to reproductive age at less than half the rate of hermaphrodites.

A global community effort to decipher the unique biology of annual killifish

Over the past 50 years, annual killifishes arose as alternative model organisms for studies of vertebrate biology. The annual fish offers exceptional advantages for studies of genetics, genomics, developmental biology, population dynamics, ecology, biogeography, and evolution. They inhabit extremely variable freshwater environments in Africa and South America, have a short lifespan and a set of unique and fascinating developmental characteristics. Embryos survive within the dry substrate during the dry season, whereas the adult population dies. Thus, the survival of the populations is entirely dependent on the buried embryos that hatch the next rainy season. Although Old and New World species share similarities in their life cycle, they also have different adaptive responses associated with climate-related selective pressures. Therefore, contrasting different species from these areas is essential to understand unique adaptations to heterogeneous environment. A network of laboratories (United States, Czech Republic, Italy, Brazil, Chile, and Uruguay) is working and collaborating on many aspects of the biology of annual fishes. Participating researchers share projects and cross-training undergraduate and graduate students. These efforts resulted in two International Symposia (2010 and 2015) that took place in Montevideo and an international book. Herein, we summarize the progress made by this global community of scientists. Developmental Dynamics 246:807-811, 2017. © 2017 Wiley Periodicals, Inc.

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.

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

The Kryptolebias marmoratus is unique because it is the only self-fertilizing hermaphroditic vertebrate, known to date. It primarily reproduces by internal self-fertilization in a mixed ovary/testis gonad. Here, we report on a high-quality genome assembly for the K. marmoratus South Korea (SK) strain highlighting the diversity and distribution of transposable elements (TEs). We find that K. marmoratus genome maintains number and composition of TEs. This can be an important genomic attribute promoting genome recombination in this selfing fish, while, in addition to a mixed mating strategy, it may also represent a mechanism contributing to the evolutionary adaptation to ecological pressure of the species. Future work should help clarify this point further once genomic information is gathered for other taxa of the family Rivulidae that do not self-fertilize. We provide a valuable genome resource that highlights the potential impact of TEs on the genome evolution of a fish species with an uncommon life cycle.