Genome evolution and meiotic maps by massively parallel DNA sequencing: spotted gar, an outgroup for the teleost genome duplication

Multi-genome comparisons reveal gain-and-loss evolution of anti-Mullerian hormone receptor type 2 as a candidate master sex-determining gene in Percidae

BACKGROUND

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii, and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems.

RESULTS

We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex-determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplicates (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been likely lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome 18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variations (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex-determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in the testis than in the ovary.

CONCLUSIONS

Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Genetic adaptations for the oceanic success of fish eggs

Genetic adaptations of organisms living in extreme environments are fundamental to our understanding of where life can evolve. Water is the single limiting parameter in this regard, yet when released in the oceans, the single-celled eggs of marine bony fishes (teleosts) have no means of acquiring it. They are strongly hyposmotic to seawater and lack osmoregulatory systems. Paradoxically, modern teleosts successfully release vast quantities of eggs in the extreme saline environment and recorded the most explosive radiation in vertebrate history. Here, we highlight key genetic adaptations that evolved to solve this paradox by filling the pre-ovulated eggs with water. The degree of water acquisition is uniquely prevalent to marine teleosts, permitting the survival and oceanic dispersal of their eggs.

DNA methylation reprogramming in teleosts

Early embryonic development is crucially important but also remarkably diverse among animal taxa. Axis formation and cell lineage specification occur due to both spatial and temporal control of gene expression. This complex system involves various signaling pathways and developmental genes such as transcription factors as well as other molecular interactants that maintain cellular states, including several types of epigenetic marks. 5mC DNA methylation, the chemical modification of cytosines in eukaryotes, represents one such mark. By influencing the compaction of chromatin (a high-order DNA structure), DNA methylation can either repress or induce transcriptional activity. Mammals exhibit a reprogramming of DNA methylation from the parental genomes in the zygote following fertilization, and later in primordial germ cells (PGCs). Whether these periods of methylation reprogramming are evolutionarily conserved, or an innovation in mammals, is an emerging question. Looking into these processes in other vertebrate lineages is thus important, and teleost fish, with their extensive species richness, phenotypic diversity, and multiple rounds of whole genome duplication, provide the perfect research playground for answering such a question. This review aims to present a concise state of the art of DNA methylation reprogramming in early development in fish by summarizing findings from different research groups investigating methylation reprogramming patterns in teleosts, while keeping in mind the ramifications of the methodology used, then comparing those patterns to reprogramming patterns in mammals.

Population genomic analyses reveal that salinity and geographic isolation drive diversification in a free-living protist

Protists make up the vast diversity of eukaryotic life and play a critical role in biogeochemical cycling and in food webs. Because of their small size, cryptic life cycles, and large population sizes, our understanding of speciation in these organisms is very limited. We performed population genomic analyses on 153 strains isolated from eight populations of the recently radiated dinoflagellate genus Apocalathium, to explore the drivers and mechanisms of speciation processes. Species of this genus inhabit both freshwater and saline habitats, lakes and seas, and are found in cold temperate environments across the world. RAD sequencing analyses revealed that the populations were overall highly differentiated, but morphological similarity was not congruent with genetic similarity. While geographic isolation was to some extent coupled to genetic distance, this pattern was not consistent. Instead, we found evidence that the environment, specifically salinity, is a major factor in driving ecological speciation in Apocalathium. While saline populations were unique in loci coupled to genes involved in osmoregulation, freshwater populations appear to lack these. Our study highlights that adaptation to freshwater through loss of osmoregulatory genes may be an important speciation mechanism in free-living aquatic protists.

Manipulation of Proteostasis Networks in Transgenic ZAAT Zebrafish via CRISPR-Cas9 Gene Editing

The CRISPR-Cas9 genome editing system is used to induce mutations in genes of interest resulting in the loss of functional protein. A transgenic zebrafish α1-antitrypsin deficiency (AATD) model displays an unusual phenotype, in that it lacks the hepatic accumulation of the misfolding Z α1-antitrypsin (ZAAT) evident in human and mouse models. Here we describe the application of the CRISPR-Cas9 system to generate mutant zebrafish with defects in key proteostasis networks likely to be involved in the hepatic processing of ZAAT in this model. We describe the targeting of the atf6a and man1b1 genes as examples.

Multi-genome comparisons reveal gain-and-loss evolution of the anti-Mullerian hormone receptor type 2 gene, an old master sex determining gene, in Percidae

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplications (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome-18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variants (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in testis than ovary. Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Ancient multiplicity in cyclic nucleotide-gated (CNG) cation channel repertoire was reduced in the ancestor of Olfactores before re-expansion by whole genome duplications in vertebrates

Cyclic nucleotide-gated (CNG) cation channels are important heterotetrameric proteins in the retina, with different subunit composition in cone and rod photoreceptor cells: three CNGA3 and one CNGB3 in cones and three CNGA1 and one CNGB1 in rods. CNGA and CNGB subunits form separate subfamilies. We have analyzed the evolution of the CNG gene family in metazoans, with special focus on vertebrates by using sequence-based phylogeny and conservation of chromosomal synteny to deduce paralogons resulting from the early vertebrate whole genome duplications (WGDs). Our analyses show, unexpectedly, that the CNGA subfamily had four sister subfamilies in the ancestor of bilaterians and cnidarians that we named CNGC, CNGD, CNGE and CNGF. Of these, CNGC, CNGE and CNGF were lost in the ancestor of Olfactores while CNGD was lost in the vertebrate ancestor. The remaining CNGA and CNGB genes were expanded by a local duplication of CNGA and the subsequent chromosome duplications in the basal vertebrate WGD events. Upon some losses, this resulted in the gnathostome ancestor having three members in the visual CNGA subfamily (CNGA1-3), a single CNGA4 gene, and two members in the CNGB subfamily (CNGB1 and CNGB3). The nature of chromosomal rearrangements in the vertebrate CNGA paralogon was resolved by including the genomes of a non-teleost actinopterygian and an elasmobranch. After the teleost-specific WGD, additional duplicates were generated and retained for CNGA1, CNGA2, CNGA3 and CNGB1. Furthermore, teleosts retain a local duplicate of CNGB3. The retention of duplicated CNG genes is explained by their subfunctionalisation and photoreceptor-specific expression. In conclusion, this study provides evidence for four previously unknown CNG subfamilies in metazoans and further evidence that the early vertebrate WGD events were instrumental in the evolution of the vertebrate visual and central nervous systems.

Comparative Study of the Molecular Characterization, Evolution, and Structure Modeling of Digestive Lipase Genes Reveals the Different Evolutionary Selection Between Mammals and Fishes

Vertebrates need suitable lipases to digest lipids for the requirement of energy and essential nutrients; however, the main digestive lipase genes of fishes have certain controversies. In this study, two types of digestive lipase genes (pancreatic lipase (pl) and bile salt-activated lipase (bsal)) were identified in mammals and fishes. The neighborhood genes and key active sites of the two lipase genes were conserved in mammals and fishes. Three copies of PL genes were found in mammals, but only one copy of the pl gene was found in most of the fish species, and the pl gene was even completely absent in some fish species (e.g., zebrafish, medaka, and common carp). Additionally, the hydrophobic amino acid residues (Ile and Leu) which are important to pancreatic lipase activity were also absent in most of the fish species. The PL was the main digestive lipase gene in mammals, but the pl gene seemed not to be the main digestive lipase gene in fish due to the absence of the pl gene sequence and the important amino acid residues. In contrast, the bsal gene existed in all fish species, even two to five copies of bsal genes were found in most of the fishes, but only one copy of the BSAL gene was found in mammals. The amino acid residues of bile salt-binding sites and the three-dimensional (3D) structure modeling of Bsal proteins were conserved in most of the fish species, so bsal might be the main digestive lipase gene in fish. The phylogenetic analysis also indicated that pl or bsal showed an independent evolution between mammals and fishes. Therefore, we inferred that the evolutionary selection of the main digestive lipase genes diverged into two types between mammals and fishes. These findings will provide valuable evidence for the study of lipid digestion in fish.

Genome-wide identification and expression patterns of opsin genes during larval development in Chinese perch (Siniperca chuatsi)

Vision is important for fish to forage food and fishes express opsin genes to receive visual signals. Chinese perch (Siniperca chuatsi) larvae prey on other fish species larvae at firstfeeding but donoteat any zooplankton, the expression of opsin genes in S. chuatsilarvae is unknown. In this study, we conducted a whole-genome analysis and demonstrated that S. chuatsihave5cone opsin genes (sws1, sws2Aα, sws2Aβ, rh2and lws)and 2 rod opsin genes (rh1and rh1-exorh). The syntenicanalysisshowedthe flanking genes ofall opsin genes were conserved during fish evolution, but the ancestorof S. chuatsimightlost some opsin gene copies duringtheevolution.The phylogeneticanalysisshowed sws1of S. chuatsiwas closest to those of Lates calcariferwhich had a truncated sws1gene; the sws2Aα, sws2Aβ,lws,rh2,rh1 andrh1-exorh of S. chuatsihad a closer relationship with those of Percomorpha fishes.Importantly, results of in situhybridization showed the sws1 opsingene,which is related to forage zooplankton,had extremely low levelexpression in retinaat early stages.Surprisingly, the rh2 opsin gene had a high level expression at firstfeeding stage. The sws2Aα, sws2Aβand lwshad a little expression at early stages but the lwsshowed a increasing trend with larval development, rh1 opsin gene expression appeared at15 dph. In thisstudy, we found a specialpattern of visual opsin genes expression in S. chuatsi, it might influence the larval first feeding and feeding habit.

Linkage mapping, comparative genome analysis, and QTL detection for growth in a non-model teleost, the meagre Argyrosomus regius, using ddRAD sequencing

Meagre (Argyrosomus regius), is a benthopelagic species rapidly emerging in aquaculture, due to its low food to biomass conversion rate, good fillet yield and ease of production. Tracing a species genomic background along with describing the genetic basis of important traits can greatly influence both conservation strategies and production perspectives. In this study, we employed ddRAD sequencing of 266 fish from six F1 meagre families, to construct a high-density genetic map comprising 4529 polymorphic SNP markers. The QTL mapping analysis provided a genomic appreciation for the weight trait identifying a statistically significant QTL on linkage group 15 (LG15). The comparative genomics analysis with six teleost species revealed an evolutionarily conserved karyotype structure. The synteny observed, verified the already well-known fusion events of the three-spine stickleback genome, reinforced the evidence of reduced evolutionary distance of Sciaenids with the Sparidae family, reflected the evolutionary proximity with Dicentrarchus labrax, traced several putative chromosomal rearrangements and a prominent putative fusion event in meagre’s LG17. This study presents novel elements concerning the genome evolutionary history of a non-model teleost species recently adopted in aquaculture, starts to unravel the genetic basis of the species growth-related traits, and provides a high-density genetic map as a tool that can help to further establish meagre as a valuable resource for research and production.

Population Genomics Analysis with RAD, Reprised: Stacks 2

Restriction enzymes have been one of the primary tools in the population genetics toolkit for 50 years, being coupled with each new generation of technology to provide a more detailed view into the genetics of natural populations. Restriction site-Associated DNA protocols, which joined enzymes with short-read sequencing technology, have democratized the field of population genomics, providing a means to assay the underlying alleles in scores of populations. More than 10 years on, the technique has been widely applied across the tree of life and served as the basis for many different analysis techniques. Here, we provide a detailed protocol to conduct a RAD analysis from experimental design to de novo analysis-including parameter optimization-as well as reference-based analysis, all in Stacks version 2, which is designed to work with paired-end reads to assemble RAD loci up to 1000 nucleotides in length. The protocol focuses on major points of friction in the molecular approaches and downstream analysis, with special attention given to validating experimental analyses. Finally, the protocol provides several points of departure for further analysis.

Presence of two RIG-I-like receptors, MDA5 and LGP2, and their dsRNA binding capacity in a perciform fish, the snakehead Channa argus

Fish retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs) are critical RNA sensors in cytoplasm and are involved in antiviral innate immunity. However, some species of fish lack RIG-I gene, and the function of RLR members in RIG-I-absent fish is poorly understood. In the present study, MDA5, LGP2 and MAVS genes were identified in commercially important snakehead Channa argus. But, RIG-I gene was not found in this fish, and a systematic analysis of RLRs in available genome database of fish indicated the absence of RIG-I in the Acanthomorphata, Clupeiformes and Polypteriformes, suggesting that loss events of RIG-I gene may have occurred independently three times in the evolutionary history of fish. The MDA5, LGP2 and MAVS in snakehead have conserved protein domains and genomic location based on sequence, phylogenetic and syntenic analyses. These genes are constitutively expressed in healthy fish and can be induced by polyinosinic and polycytidylic acid (poly(I:C)) stimulation in vitro. It is further revealed that the snakehead MDA5 and LGP2 have binding capacity with dsRNA, such as poly(I:C), and MDA5 can interact with MAVS, implying the antiviral function of MDA5 in the RIG-I-absent fish.

Molecular identification of Klebsiella pneumoniae and expression of immune genes in infected spotted gar Lepisosteus oculatus

Spotted gar (Lepisosteus oculatus) is a primitive ray-finned fish which has not undergone the third round whole genome duplication and commonly used as a model to study the evolution of immune genes. In this study, a pathogenic strain of Klebsiella pneumoniae (termed KPY01) was isolated from a diseased spotted gar, based on the Gram-stain and phylogenetic analysis of the 16S rDNA and khe genes. Further, the virulence genes and drug resistance genes were determined and drug sensitivity tests were performed to explore the virulence and drug resistance of the KPY01. Putative biosynthetic gene clusters (BGCs) for the biosynthesis of secondary metabolites were predicted using the anti-SMASH5.0 online genome mining platform. Histopathological analysis revealed that the immune cells were significantly decreased in the white pulp of spleen of fish infected with K. pneumonia and tissue inflammation became apparent. Besides, the expression of cytokines including interleukin (il) -8, il-10, il-12a, il-18 and interferon γ (ifn-γ) were shown to be modulated in the spleen, gills and kidney. Our work provides useful information for further investigation on the virulence of K. pneumoniae and host immune responses to K. pneumoniae infection in fish.

Seasonal genotype dynamics of a marine dinoflagellate: Pelagic populations are homogeneous and as diverse as benthic seed banks

Genetic diversity is the basis for evolutionary adaptation and selection under changing environmental conditions. Phytoplankton populations are genotypically diverse, can become genetically differentiated within small spatiotemporal scales and many species form resting stages. Resting stage accumulations in sediments (seed banks) are expected to serve as reservoirs for genetic information, but so far their role in maintaining phytoplankton diversity and in evolution has remained unclear. In this study we used the toxic dinoflagellate Alexandrium ostenfeldii (Dinophyceae) as a model organism to investigate if (i) the benthic seed bank is more diverse than the pelagic population and (ii) the pelagic population is seasonally differentiated. Resting stages (benthic) and plankton (pelagic) samples were collected at a coastal bloom site in the Baltic Sea, followed by cell isolation and genotyping using microsatellite markers (MS) and restriction site associated DNA sequencing (RAD). High clonal diversity (98%–100%) combined with intermediate to low gene diversity (0.58–0.03, depending on the marker) was found. Surprisingly, the benthic and pelagic fractions of the population were equally diverse, and the pelagic fraction was temporally homogeneous, despite seasonal fluctuation of environmental selection pressures. The results of this study suggest that continuous benthic–pelagic coupling, combined with frequent sexual reproduction, as indicated by persistent linkage equilibrium, prevent the dominance of single clonal lineages in a dynamic environment. Both processes harmonize the pelagic with the benthic population and thus prevent seasonal population differentiation. At the same time, frequent sexual reproduction and benthic–pelagic coupling maintain high clonal diversity in both habitats.

The clinical importance of tandem exon duplication-derived substitutions

Most coding genes in the human genome are annotated with multiple alternative transcripts. However, clear evidence for the functional relevance of the protein isoforms produced by these alternative transcripts is often hard to find. Alternative isoforms generated from tandem exon duplication-derived substitutions are an exception. These splice events are rare, but have important functional consequences. Here, we have catalogued the 236 tandem exon duplication-derived substitutions annotated in the GENCODE human reference set. We find that more than 90% of the events have a last common ancestor in teleost fish, so are at least 425 million years old, and twenty-one can be traced back to the Bilateria clade. Alternative isoforms generated from tandem exon duplication-derived substitutions also have significantly more clinical impact than other alternative isoforms. Tandem exon duplication-derived substitutions have >25 times as many pathogenic and likely pathogenic mutations as other alternative events. Tandem exon duplication-derived substitutions appear to have vital functional roles in the cell and may have played a prominent part in metazoan evolution.

Genome-wide identification, evolution of Krüppel-like factors (klfs) and their expressions during GCRV challenge in grass carp (Ctenopharyngodonidella)

The Krüppel-like factors (KLFs) are a family of transcription factors containing three highly conserved tandem zinc finger structures, and each member participates in multiple physiological and pathological processes. The publication of genome sequences and the application of bioinformatics tools have led to the discovery of numerous gene families in fishes. Here, 24 klf genes were re-annotated in grass carp. Subsequently, the number of klf family members were investigated in some representative vertebrate species. Then, a series of bioinformatics analysis showed that grass carp klfs in the same subfamily had similar genome structure patterns and conserved distribution patterns of motifs, which supported their molecular evolutionary relationships. Furthermore, the mRNA expression profiles showed that 24 grass carp klfs were ubiquitously expressed in 11 different tissues, and some of them displayed tissue-enriched expression patterns. Finally, the expressions of the evolutionarily expanded klf members (klf2a, 2b, 2l, 5a, 5b, 5l, 6a, 6b, 7a, 7b, 11a, 11b, 12a, 12b, 15 and 15l) during GCRV infection were also analyzed. The results suggested that grass carp klf genes with common evolutionary sources may share functional diversity and conservation. In conclusion, this study provides preliminary clues for further researches on grass carp klf members and their underlying transcriptional regulatory mechanisms during GCRV infection.

Identification, Expression and Evolution of Short-Chain Dehydrogenases/Reductases in Nile Tilapia (Oreochromis niloticus)

The short-chain dehydrogenases/reductases (SDR) superfamily is involved in multiple physiological processes. In this study, genome-wide identification and comprehensive analysis of SDR superfamily were carried out in 29 animal species based on the latest genome databases. Overall, the number of SDR genes in animals increased with whole genome duplication (WGD), suggesting the expansion of SDRs during evolution, especially in 3R-WGD and polyploidization of teleosts. Phylogenetic analysis indicated that vertebrates SDRs were clustered into five categories: classical, extended, undefined, atypical, and complex. Moreover, tandem duplication of hpgd-a, rdh8b and dhrs13 was observed in teleosts analyzed. Additionally, tandem duplications of dhrs11-a, dhrs7a, hsd11b1b, and cbr1-a were observed in all cichlids analyzed, and tandem duplication of rdh10-b was observed in tilapiines. Transcriptome analysis of adult fish revealed that 93 SDRs were expressed in more than one tissue and 5 in one tissue only. Transcriptome analysis of gonads from different developmental stages showed that expression of 17 SDRs were sexually dimorphic with 11 higher in ovary and 6 higher in testis. The sexually dimorphic expressions of these SDRs were confirmed by in situ hybridization (ISH) and qPCR, indicating their possible roles in steroidogenesis and gonadal differentiation. Taken together, the identification and the expression data obtained in this study contribute to a better understanding of SDR superfamily evolution and functions in teleosts.

RADSex: A computational workflow to study sex determination using restriction site-associated DNA sequencing data

The study of sex determination and sex chromosome organization in nonmodel species has long been technically challenging, but new sequencing methodologies now enable precise and high-throughput identification of sex-specific genomic sequences. In particular, restriction site-associated DNA sequencing (RAD-Seq) is being extensively applied to explore sex determination systems in many plant and animal species. However, software specifically designed to search for and visualize sex-biased markers using RAD-Seq data is lacking. Here, we present RADSex, a computational analysis workflow designed to study the genetic basis of sex determination using RAD-Seq data. RADSex is simple to use, requires few computational resources, makes no prior assumptions about the type of sex-determination system or structure of the sex locus, and offers convenient visualization through a dedicated R package. To demonstrate the functionality of RADSex, we re-analysed a published data set of Japanese medaka, Oryzias latipes, where we uncovered a previously unknown Y chromosome polymorphism. We then used RADSex to analyse new RAD-Seq data sets from 15 fish species spanning multiple taxonomic orders. We identified the sex determination system and sex-specific markers in six of these species, five of which had no known sex-markers prior to this study. We show that RADSex greatly facilitates the study of sex determination systems in nonmodel species thanks to its speed of analyses, low resource usage, ease of application and visualization options. Furthermore, our analysis of new data sets from 15 species provides new insights on sex determination in fish.

A fish with no sex: gonadal and adrenal functions partition between zebrafish NR5A1 co-orthologs

People with NR5A1 mutations experience testicular dysgenesis, ovotestes, or adrenal insufficiency, but we do not completely understand the origin of this phenotypic diversity. NR5A1 is expressed in gonadal soma precursor cells before expression of the sex-determining gene SRY. Many fish have two co-orthologs of NR5A1 that likely partitioned ancestral gene subfunctions between them. To explore ancestral roles of NR5A1, we knocked out nr5a1a and nr5a1b in zebrafish. Single-cell RNA-seq identified nr5a1a-expressing cells that co-expressed genes for steroid biosynthesis and the chemokine receptor Cxcl12a in 1-day postfertilization (dpf) embryos, as does the mammalian adrenal-gonadal (interrenal-gonadal) primordium. In 2dpf embryos, nr5a1a was expressed stronger in the interrenal-gonadal primordium than in the early hypothalamus but nr5a1b showed the reverse. Adult Leydig cells expressed both ohnologs and granulosa cells expressed nr5a1a stronger than nr5a1b. Mutants for nr5a1a lacked the interrenal, formed incompletely differentiated testes, had no Leydig cells, and grew far larger than normal fish. Mutants for nr5a1b formed a disorganized interrenal and their gonads completely disappeared. All homozygous mutant genotypes lacked secondary sex characteristics, including male breeding tubercles and female sex papillae, and had exceedingly low levels of estradiol, 11-ketotestosterone, and cortisol. RNA-seq showed that at 21dpf, some animals were developing as females and others were not, independent of nr5a1 genotype. By 35dpf, all mutant genotypes greatly under-expressed ovary-biased genes. Because adult nr5a1a mutants form gonads but lack an interrenal and conversely, adult nr5a1b mutants lack a gonad but have an interrenal, the adrenal, and gonadal functions of the ancestral nr5a1 gene partitioned between ohnologs after the teleost genome duplication, likely owing to reciprocal loss of ancestral tissue-specific regulatory elements. Identifying such elements could provide hints to otherwise unexplained cases of Differences in Sex Development.

The rise and fall of the ancient northern pike master sex-determining gene

The understanding of the evolution of variable sex determination mechanisms across taxa requires comparative studies among closely related species. Following the fate of a known master sex-determining gene, we traced the evolution of sex determination in an entire teleost order (Esociformes). We discovered that the northern pike (Esox lucius) master sex-determining gene originated from a 65 to 90 million-year-old gene duplication event and that it remained sex linked on undifferentiated sex chromosomes for at least 56 million years in multiple species. We identified several independent species- or population-specific sex determination transitions, including a recent loss of a Y chromosome. These findings highlight the diversity of evolutionary fates of master sex-determining genes and the importance of population demographic history in sex determination studies. We hypothesize that occasional sex reversals and genetic bottlenecks provide a non-adaptive explanation for sex determination transitions.

Genome-wide single nucleotide polymorphism markers reveal population structure and dispersal direction of an expanding nuisance algal bloom species

Species invasion and range expansion are currently under scrutiny due to increasing anthropogenic impact on the natural environment. This is also true for harmful algal blooms, which have been reported to have increased in frequency. However, this research is challenging due to the ephemeral nature, small size and mostly low concentrations of microalgae in the environment. One such species is the nuisance microalga Gonyostomum semen (Raphidophyceae), which has increased in occurrence in northern Europe in recent decades. The question of whether the species has expanded its habitat range or if it was already present in the lakes but was too rare to be detected remains unanswered. The aim of the present study was to determine the genetic structure and dispersal pathways of G. semen using RAD (restriction-site-associated DNA) tag sequencing. For G. semen, which has a huge genome (32 Gbp), we faced particular challenges, but were nevertheless able to recover over 1000 single nucleotide polymorphisms at high coverage. Our data revealed a distinct population genetic structure, demonstrating a divide of western and eastern populations that probably represent different lineages. Despite significant genetic differentiation among lakes, we found only limited isolation-by-distance. While we had expected a pattern of recent expansion northwards, the data demonstrated gene flow from the northeast/east towards the southwest/west. This genetic signature suggests that the observed gene flow may be due to dispersal by autumn migratory birds, which act as dispersal vectors of resistant resting propagules that form at the end of the G. semen blooms.

Zebrafish Kit ligands cooperate with erythropoietin to promote erythroid cell expansion

Kit ligand (Kitlg) is pleiotropic cytokine with a prominent role in vertebrate erythropoiesis. Although the role of Kitlg in this process has not been reported in Danio rerio (zebrafish), in the present study we show that its function is evolutionarily conserved. Zebrafish possess 2 copies of Kitlg genes (Kitlga and Kitlgb) as a result of whole-genome duplication. To determine the role of each ligand in zebrafish, we performed a series of ex vivo and in vivo gain- and loss-of-function experiments. First, we tested the biological activity of recombinant Kitlg proteins in suspension culture from zebrafish whole-kidney marrow, and we demonstrate that Kitlga is necessary for expansion of erythroid progenitors ex vivo. To further address the role of kitlga and kitlgb in hematopoietic development in vivo, we performed gain-of-function experiments in zebrafish embryos, showing that both ligands cooperate with erythropoietin (Epo) to promote erythroid cell expansion. Finally, using the kita mutant (kitab5/b5 or sparse), we show that the Kita receptor is crucial for Kitlga/b cooperation with Epo in erythroid cells. In summary, using optimized suspension culture conditions with recombinant cytokines (Epo, Kitlga), we report, for the first time, ex vivo suspension cultures of zebrafish hematopoietic progenitor cells that can serve as an indispensable tool to study normal and aberrant hematopoiesis in zebrafish. Furthermore, we conclude that, although partial functional diversification of Kit ligands has been described in other processes, in erythroid development, both paralogs play a similar role, and their function is evolutionarily conserved.

The shift from early to late types of ribosomes in zebrafish development involves changes at a subset of rRNA 2'-O-Me sites

During zebrafish development, an early type of rRNA is gradually replaced by a late type that is substantially different in sequence. We applied RiboMeth-seq to rRNA from developmental stages for profiling of 2'-O-Me, to learn if changes in methylation pattern were a component of the shift. We compiled a catalog of 2'-O-Me sites and cognate box C/D guide RNAs comprising 98 high-confidence sites, including 10 sites that were not known from other vertebrates, one of which was specific to late-type rRNA. We identified a subset of sites that changed in methylation status during development and found that some of these could be explained by availability of their cognate SNORDs. Sites that changed during development were enriched in the novel sites revealed in zebrafish. We propose that the early type of rRNA is a specialized form and that its structure and ribose methylation pattern may be an adaptation to features of development, including translation of specific maternal mRNAs.

Gonadotropin-inhibitory hormone (GnIH) distribution in the brain of the ancient fish Atractosteus tropicus (Holostei, Lepisosteiformes)

The Holostei group occupies a critical phylogenetic position as the sister group of the Teleostei. However, little is known about holostean pituitary anatomy or brain distribution of important reproductive neuropeptides, such as the gonadotropin-inhibitory hormone (GnIH). Thus, the present study set out to characterize the structure of the pituitary and to localize GnIH-immunoreactive cells in the brain of Atractosteus tropicus from the viewpoint of comparative neuroanatomy. Juveniles of both sexes were processed for general histology and immunohistochemistry. Based on the differences in cell organization, morphology, and staining properties, the neurohypophysis and three regions in the adenohypophysis were identified: the rostral and proximal pars distalis (PPD) and the pars intermedia. This last region was found to be innervated by the neurohypophysis. This organization, together with the presence of a saccus vasculosus, resembles the general teleost pituitary organization. A vast number of blood vessels were also recognized between the infundibulum floor of the hypothalamus and the PPD, evidencing the characteristic presence of a median eminence and a portal system. However, this well-developed pituitary portal system resembles that of tetrapods. As regards the immunohistochemical localization of GnIH, we found four GnIH-immunoreactive (GnIH-ir) populations in three hypothalamic nuclei (suprachiasmatic, retrotuberal, and tuberal nuclei) and one in the diencephalon (prethalamic nucleus), as well as a few scattered neurons throughout the olfactory bulbs, the telencephalon, and the intersection between them. GnIH-ir fibers showed a widespread distribution over almost all brain regions, suggesting that GnIH function is not restricted to reproduction only. In conclusion, the present study describes, for the first time, the pituitary of A. tropicus and the neuroanatomical localization of GnIH in a holostean fish that exhibits a similar distribution pattern to that of teleosts and other vertebrates, suggesting a high degree of phylogenetic conservation of this system.

Artificial whole genome duplication in paleopolyploid sturgeons yields highest documented chromosome number in vertebrates

Critically endangered sturgeons, having undergone three whole genome duplication events, represent an exceptional example of ploidy plasticity in vertebrates. Three extant ploidy groups, combined with autopolyploidization, interspecific hybridization and the fertility of hybrids are important issues in sturgeon conservation and aquaculture. Here we demonstrate that the sturgeon genome can undergo numerous alterations of ploidy without severe physiological consequences, producing progeny with a range of ploidy levels and extremely high chromosome numbers. Artificial suppression of the first mitotic division alone, or in combination with suppression of the second meiotic division of functionally tetraploid zygotes (4n, C-value = 4.15) of Siberian sturgeon Acipenser baerii and Russian sturgeon A. gueldenstaedtii resulted in progeny of various ploidy levels—diploid/hexaploid (2n/6n) mosaics, hexaploid, octoploid juveniles (8n), and dodecaploid (12n) larvae. Counts between 477 to 520 chromosomes in octoploid juveniles of both sturgeons confirmed the modal chromosome numbers of parental species had been doubled. This exceeds the highest previously documented chromosome count among vertebrates 2n ~ 446 in the cyprinid fish Ptychobarbus dipogon.

Global assessment of organ specific basal gene expression over a diurnal cycle with analyses of gene copies exhibiting cyclic expression patterns

Background

Studying functional divergences between paralogs that originated from genome duplication is a significant topic in investigating molecular evolution. Genes that exhibit basal level cyclic expression patterns including circadian and light responsive genes are important physiological regulators. Temporal shifts in basal gene expression patterns are important factors to be considered when studying genetic functions. However, adequate efforts have not been applied to studying basal gene expression variation on a global scale to establish transcriptional activity baselines for each organ. Furthermore, the investigation of cyclic expression pattern comparisons between genome duplication created paralogs, and potential functional divergence between them has been neglected. To address these questions, we utilized a teleost fish species, Xiphophorus maculatus, and profiled gene expression within 9 organs at 3-h intervals throughout a 24-h diurnal period.

Results

Our results showed 1.3–21.9% of genes in different organs exhibited cyclic expression patterns, with eye showing the highest fraction of cycling genes while gonads yielded the lowest. A majority of the duplicated gene pairs exhibited divergences in their basal level expression patterns wherein only one paralog exhibited an oscillating expression pattern, or both paralogs exhibit oscillating expression patterns, but each gene duplicate showed a different peak expression time, and/or in different organs.

Conclusions

These observations suggest cyclic genes experienced significant sub-, neo-, or non-functionalization following the teleost genome duplication event. In addition, we developed a customized, web-accessible, gene expression browser to facilitate data mining and data visualization for the scientific community.

Directional divergence of Ep300 duplicates in teleosts and its implications

Background

EP300 is a conserved protein in vertebrates, which serves as a key mediator of cellular homeostasis. Mutations and dysregulation of EP300 give rise to severe human developmental disorders and malignancy. Danio rerio is a promising model organism to study EP300 related diseases and drugs; however, the effect of EP300 duplicates derived from teleost-specific whole genome duplication should not just be neglected.

Results

In this study, we obtained EP300 protein sequences of representative teleosts, mammals and sauropsids, with which we inferred a highly supported maximum likelihood tree. We observed that Ep300 duplicates (Ep300a and Ep300b) were widely retained in teleosts and universally expressed in a variety of tissues. Consensus sequences of Ep300a and Ep300b had exactly the same distribution of conserved domains, suggesting that their functions should still be largely overlapped. We analyzed the molecular evolution of Ep300 duplicates in teleosts, using branch-site models, clade models and site models. The results showed that both duplicates were subject to strong positive selection; however, for an extant species, generally at most one copy was under positive selection. At the clade level, there were evident positive correlations between evolutionary rates, the number of positively selected sites and gene expression levels. In Ostariophysi, Ep300a were under stronger positive selection than Ep300b; in Neoteleostei, another species-rich teleost clade, the contrary was the case. We also modeled 3D structures of zf-TAZ domain and its flanking regions of Ep300a and Ep300b of D. rerio and Oryzias latipes and found that in either species the faster evolving copy had more short helixes.

Conclusions

Collectively, the two copies of Ep300 have undoubtedly experienced directional divergence in main teleost clades. The divergence of EP300 between teleosts and mammals should be greater than the divergence between different teleost clades. Further studies are needed to clarify to what extent the EP300 involved regulatory network has diverged between teleosts and mammals, which would also help explain the huge success of teleosts.

Genome of the four-finger threadfin Eleutheronema tetradactylum (Perciforms: Polynemidae)

BACKGROUND

Teleost fish play important roles in aquatic ecosystems and aquaculture. Threadfins (Perciformes: Polynemidae) show a range of interesting biology, and are of considerable importance for both wild fisheries and aquaculture. Additionally, the four-finger threadfin Eleutheronema tetradactylum is of conservation relevance since its populations are considered to be in rapid decline and it is classified as endangered. However, no genomic resources are currently available for the threadfin family Polynemidae.

RESULTS

We sequenced and assembled the first threadfin fish genome, the four-finger threadfin E. tetradactylum. We provide a genome assembly for E. tetradactylum with high contiguity (scaffold N50 = 56.3 kb) and high BUSCO completeness at 96.5%. The assembled genome size of E. tetradactylum is just 610.5 Mb, making it the second smallest perciform genome assembled to date. Just 9.07-10.91% of the genome sequence was found to consist of repetitive elements (standard RepeatMasker analysis vs custom analysis), making this the lowest repeat content identified to date for any perciform fish. A total of 37,683 protein-coding genes were annotated, and we include analyses of developmental transcription factors, including the Hox, ParaHox, and Sox families. MicroRNA genes were also annotated and compared with other chordate lineages, elucidating the gains and losses of chordate microRNAs.

CONCLUSIONS

The four-finger threadfin E. tetradactylum genome presented here represents the first available genome sequence for the ecologically, biologically, and commercially important clade of threadfin fish. Our findings provide a useful genomic resource for future research into the interesting biology and evolution of this valuable group of food fish.

ipyrad: Interactive assembly and analysis of RADseq datasets

SUMMARY

ipyrad is a free and open source tool for assembling and analyzing restriction site-associated DNA sequence datasets using de novo and/or reference-based approaches. It is designed to be massively scalable to hundreds of taxa and thousands of samples, and can be efficiently parallelized on high performance computing clusters. It is available both as a command line interface and as a Python package with an application programming interface, the latter of which can be used interactively to write complex, reproducible scripts and implement a suite of downstream analysis tools.

AVAILABILITY AND IMPLEMENTATION

ipyrad is a free and open source program written in Python. Source code is available from the GitHub repository (https://github.com/dereneaton/ipyrad/), and Linux and MacOS installs are distributed through the conda package manager. Complete documentation, including numerous tutorials, and Jupyter notebooks demonstrating example assemblies and applications of downstream analysis tools are available online: https://ipyrad.readthedocs.io/.

Preferential Mapping of Sex-Biased Differentially-Expressed Genes of Larvae to the Sex-Determining Region of Flathead Grey Mullet (Mugil cephalus)

Flathead gray mullet (Mugil cephalus) is a cosmopolitan mugilid species popular in fishery and aquaculture with an economic preference for all-female population. However, it displays neither sexual dimorphisms nor heteromorphic sex chromosomes. We have previously presented a microsatellite-based linkage map for this species locating a single sex determination region (SDR) on linkage group 9 (LG9) with evidence for XX/XY sex determination (SD) mechanism. In this work, we refine the critical SDR on LG9, and propose positional- and functional- candidate genes for SD. To elucidate the genetic mechanism of SD, we assembled and compared male and female genomic sequences of 19 syntenic genes within the putative SDR on mullet's LG9, based on orthology to tilapia's LG8 (tLG8) physical map. A total of 25 sequence-based markers in 12 genes were developed. For all markers, we observed association with sex in at least one of the two analyzed M. cephalus full-sib families, but not in the wild-type population. Recombination events were inferred within families thus setting the SDR boundaries to a region orthologous to ∼0.9 Mbp with 27 genes on tLG8. As the sexual phenotype is evident only in adults, larvae were assigned into two putative sex-groups according to their paternal haplotypes, following a model of XY/XX SD-system. A total of 107 sex-biased differentially expressed genes in larvae were observed, of which 51 were mapped to tLG8 (48% enrichment), as compared to 5% in random control. Furthermore, 23 of the 107 genes displayed sex-specific expression; and 22 of these genes were positioned to tLG8, indicating 96% enrichment. Of the 27 SDR genes, BCCIP, DHX32A, DOCK1, and FSHR (GTH-RI) are suggested as positional and functional gene candidates for SD.

The Evolution of Duplicated Genes of the Cpi-17/Phi-1 (ppp1r14) Family of Protein Phosphatase 1 Inhibitors in Teleosts

The Cpi-17 (ppp1r14) gene family is an evolutionarily conserved, vertebrate specific group of protein phosphatase 1 (PP1) inhibitors. When phosphorylated, Cpi-17 is a potent inhibitor of myosin phosphatase (MP), a holoenzyme complex of the regulatory subunit Mypt1 and the catalytic subunit PP1. Myosin phosphatase dephosphorylates the regulatory myosin light chain (Mlc2) and promotes actomyosin relaxation, which in turn, regulates numerous cellular processes including smooth muscle contraction, cytokinesis, cell motility, and tumor cell invasion. We analyzed zebrafish homologs of the Cpi-17 family, to better understand the mechanisms of myosin phosphatase regulation. We found single homologs of both Kepi (ppp1r14c) and Gbpi (ppp1r14d) in silico, but we detected no expression of these genes during early embryonic development. Cpi-17 (ppp1r14a) and Phi-1 (ppp1r14b) each had two duplicate paralogs, (ppp1r14aa and ppp1r14ab) and (ppp1r14ba and ppp1r14bb), which were each expressed during early development. The spatial expression pattern of these genes has diverged, with ppp1r14aa and ppp1r14bb expressed primarily in smooth muscle and skeletal muscle, respectively, while ppp1r14ab and ppp1r14ba are primarily expressed in neural tissue. We observed that, in in vitro and heterologous cellular systems, the Cpi-17 paralogs both acted as potent myosin phosphatase inhibitors, and were indistinguishable from one another. In contrast, the two Phi-1 paralogs displayed weak myosin phosphatase inhibitory activity in vitro, and did not alter myosin phosphorylation in cells. Through deletion and chimeric analysis, we identified that the difference in specificity for myosin phosphatase between Cpi-17 and Phi-1 was encoded by the highly conserved PHIN (phosphatase holoenzyme inhibitory) domain, and not the more divergent N- and C- termini. We also showed that either Cpi-17 paralog can rescue the knockdown phenotype, but neither Phi-1 paralog could do so. Thus, we provide new evidence about the biochemical and developmental distinctions of the zebrafish Cpi-17 protein family.

Unravelling the Complex Duplication History of Deuterostome Glycerol Transporters

Transmembrane glycerol transport is an ancient biophysical property that evolved in selected subfamilies of water channel (aquaporin) proteins. Here, we conducted broad level genome (>550) and transcriptome (>300) analyses to unravel the duplication history of the glycerol-transporting channels (glps) in Deuterostomia. We found that tandem duplication (TD) was the major mechanism of gene expansion in echinoderms and hemichordates, which, together with whole genome duplications (WGD) in the chordate lineage, continued to shape the genomic repertoires in craniates. Molecular phylogenies indicated that aqp3-like and aqp13-like channels were the probable stem subfamilies in craniates, with WGD generating aqp9 and aqp10 in gnathostomes but aqp7 arising through TD in Osteichthyes. We uncovered separate examples of gene translocations, gene conversion, and concerted evolution in humans, teleosts, and starfishes, with DNA transposons the likely drivers of gene rearrangements in paleotetraploid salmonids. Currently, gene copy numbers and BLAST are poor predictors of orthologous relationships due to asymmetric glp gene evolution in the different lineages. Such asymmetries can impact estimations of divergence times by millions of years. Experimental investigations of the salmonid channels demonstrated that approximately half of the 20 ancestral paralogs are functional, with neofunctionalization occurring at the transcriptional level rather than the protein transport properties. The combined findings resolve the origins and diversification of glps over >800 million years old and thus form the novel basis for proposing a pandeuterostome glp gene nomenclature.

Reconstruction of the Carbohydrate 6-O Sulfotransferase Gene Family Evolution in Vertebrates Reveals Novel Member, CHST16, Lost in Amniotes

Glycosaminoglycans are sulfated polysaccharide molecules, essential for many biological processes. The 6-O sulfation of glycosaminoglycans is carried out by carbohydrate 6-O sulfotransferases (C6OSTs), previously named Gal/GalNAc/GlcNAc 6-O sulfotransferases. Here, for the first time, we present a detailed phylogenetic reconstruction, analysis of gene synteny conservation and propose an evolutionary scenario for the C6OST family in major vertebrate groups, including mammals, birds, nonavian reptiles, amphibians, lobe-finned fishes, ray-finned fishes, cartilaginous fishes, and jawless vertebrates. The C6OST gene expansion likely started early in the chordate lineage, giving rise to four ancestral genes after the divergence of tunicates and before the emergence of extant vertebrates. The two rounds of whole-genome duplication in early vertebrate evolution (1R/2R) only contributed two additional C6OST subtype genes, increasing the vertebrate repertoire from four genes to six, divided into two branches. The first branch includes CHST1 and CHST3 as well as a previously unrecognized subtype, CHST16 that was lost in amniotes. The second branch includes CHST2, CHST7, and CHST5. Subsequently, local duplications of CHST5 gave rise to CHST4 in the ancestor of tetrapods, and to CHST6 in the ancestor of primates. The teleost-specific gene duplicates were identified for CHST1, CHST2, and CHST3 and are result of whole-genome duplication (3R) in the teleost lineage. We could also detect multiple, more recent lineage-specific duplicates. Thus, the vertebrate repertoire of C6OST genes has been shaped by gene duplications and gene losses at several stages of vertebrate evolution, with implications for the evolution of skeleton, nervous system, and cell-cell interactions.

CRISPR Meets Zebrafish: Accelerating the Discovery of New Therapeutic Targets

Bringing a new drug to the market costs an average of US$2.6 billion and takes more than 10 years from discovery to regulatory approval. Despite the need to reduce cost and time to increase productivity, pharma companies tend to crowd their efforts in the same indications and drug targets. This results in the commercialization of drugs that share the same mechanism of action (MoA) and, in many cases, equivalent efficacies among them-an outcome that helps neither patients nor the balance sheet of the companies trying to bring therapeutics to the same patient population. Indeed, the discovery of new therapeutic targets, based on a deeper understanding of the disease biology, would likely provide more innovative MoAs and potentially greater drug efficacies. It would also bring better chances for identifying appropriate treatments according to the patient's genetic stratification. Nowadays, we count with an enormous amount of unprocessed information on potential disease targets that could be extracted from omics data obtained from patient samples. In addition, hundreds of pharmacological and genetic screenings have been performed to identify innovative drug targets. Traditionally, rodents have been the animal models of choice to perform functional genomic studies. The high experimental cost, combined with the low throughput provided by those models, however, is a bottleneck for discovering and validating novel genetic disease associations. To overcome these limitations, we propose that zebrafish, in conjunction with the use of CRISPR/Cas9 genome-editing tools, could streamline functional genomic processes to bring biologically relevant knowledge on innovative disease targets in a shorter time frame.

Simulation with RADinitio improves RADseq experimental design and sheds light on sources of missing data

Restriction-site associated DNA sequencing (RADseq) has become a powerful and versatile tool in modern population genomics, enabling large-scale evolutionary and genomic analyses in otherwise inaccessible biological systems. With its widespread use, different variants on the protocol have been developed to suit specific experimental needs. Researchers face the challenge of choosing the optimal molecular and sequencing protocols for their reduced representation experimental design, an often-complicated process. Strategic errors can lead to biased data generation that has reduced power to answer biological questions. Here, we present RADinitio, simulation software for the selection and optimization of RADseq experiments via the generation of sequencing data that behave similarly to empirical sources. RADinitio provides an evolutionary simulation of populations, implementation of various RADseq protocols with customizable parameters, and thorough assessment of missing data. We test the efficacy of the software using different RAD protocols across several organisms, highlighting the importance of protocol selection on the magnitude and quality of data acquired. Additionally, we test the effects of RAD library preparation and sequencing on allelic dropout, observing that library preparation and sequencing often contributes more to missing alleles than population-level variation.

Zebrafish models of sarcopenia

Sarcopenia - the accelerated age-related loss of muscle mass and function - is an under-diagnosed condition, and is central to deteriorating mobility, disability and frailty in older age. There is a lack of treatment options for older adults at risk of sarcopenia. Although sarcopenia's pathogenesis is multifactorial, its major phenotypes - muscle mass and muscle strength - are highly heritable. Several genome-wide association studies of muscle-related traits were published recently, providing dozens of candidate genes, many with unknown function. Therefore, animal models are required not only to identify causal mechanisms, but also to clarify the underlying biology and translate this knowledge into new interventions. Over the past several decades, small teleost fishes had emerged as powerful systems for modeling the genetics of human diseases. Owing to their amenability to rapid genetic intervention and the large number of conserved genetic and physiological features, small teleosts - such as zebrafish, medaka and killifish - have become indispensable for skeletal muscle genomic studies. The goal of this Review is to summarize evidence supporting the utility of small fish models for accelerating our understanding of human skeletal muscle in health and disease. We do this by providing a basic foundation of the (zebra)fish skeletal muscle morphology and physiology, and evidence of muscle-related gene homology. We also outline challenges in interpreting zebrafish mutant phenotypes and in translating them to human disease. Finally, we conclude with recommendations on future directions to leverage the large body of tools developed in small fish for the needs of genomic exploration in sarcopenia.

Evolutionary and Molecular Characterization of liver-enriched gene 1

Liver-enriched gene 1 (Leg1) is a newly identified gene with little available functional information. To evolutionarily and molecularly characterize Leg1 genes, a phylogenetic study was first conducted, which indicated that Leg1 is a conserved gene that exists from bacteria to mammals. During the evolution of mammals, Leg1s underwent tandem duplications, which gave rise to Leg1a, Leg1b, and Leg1c clades. Analysis of the pig genome showed the presence of all three paralogs of pig Leg1 genes (pLeg1s), whereas only Leg1a could be found in the human (hLeg1a) or mouse (mLeg1a) genomes. Purifying force acts on the evolution of Leg1 genes, likely subjecting them to functional constraint. Molecularly, pLeg1a and its coded protein, pig LEG1a (pLEG1a), displayed high similarities to its human and mouse homologs in terms of gene organization, expression patterns, and structures. Hence, pLeg1a, hLeg1a, and mLeg1a might preserve similar functions. Additionally, expression analysis of the three Leg1as suggested that eutherian Leg1as might have different functions from those of zebrafish and platypus due to subfunctionalization. Therefore, pLeg1a might provide essential information about eutherian Leg1a. Moreover, a preliminary functional study using RNA-seq suggested that pLeg1a is involved in the lipid homeostasis. In conclusion, our study provides some basic information on the aspects of evolution and molecular function, which could be applied for further validation of Leg1 using pig models.

Genomic Resources for Darters (Percidae: Etheostominae) Provide Insight into Postzygotic Barriers Implicated in Speciation

Comparative genomic approaches are increasingly being used to study the evolution of reproductive barriers in nonmodel species. Although numerous studies have examined prezygotic isolation in darters (Percidae), investigations into postzygotic barriers have remained rare due to long generation times and a lack of genomic resources. Orangethroat and rainbow darters naturally hybridize and provide a remarkable example of male-driven speciation via character displacement. Backcross hybrids suffer from high mortality, which appears to promote behavioral isolation in sympatry. To investigate the genomic architecture of postzygotic isolation, we used Illumina and PacBio sequencing to generate a chromosome-level, annotated assembly of the orangethroat darter genome and high-density linkage maps for orangethroat and rainbow darters. We also analyzed genome-wide RADseq data from wild-caught adults of both species and laboratory-generated backcrosses to identify genomic regions associated with hybrid incompatibles. Several putative chromosomal translocations and inversions were observed between orangethroat and rainbow darters, suggesting structural rearrangements may underlie postzygotic isolation. We also found evidence of selection against recombinant haplotypes and transmission ratio distortion in backcross hybrid genomes, providing further insight into the genomic architecture of genetic incompatibilities. Notably, regions with high levels of genetic divergence between species were enriched for genes associated with developmental and meiotic processes, providing strong candidates for postzygotic isolating barriers. These findings mark significant contributions to our understanding of the genetic basis of reproductive isolation between species undergoing character displacement. Furthermore, the genomic resources presented here will be instrumental for studying speciation in darters, the most diverse vertebrate group in North America.

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

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

New Insights Into the Evolutionary History of Melatonin Receptors in Vertebrates, With Particular Focus on Teleosts

In order to improve our understanding of melatonin signaling, we have reviewed and revised the evolutionary history of melatonin receptor genes (mtnr) in vertebrates. All gnathostome mtnr genes have a conserved gene organization with two exons, except for mtnr1b paralogs of some teleosts that show intron gains. Phylogeny and synteny analyses demonstrate the presence of four mtnr subtypes, MTNR1A, MTNR1B, MTNR1C, MTNR1D that arose from duplication of an ancestral mtnr during the vertebrate tetraploidizations (1R and 2R). In tetrapods, mtnr1d was lost, independently, in mammals, in archosaurs and in caecilian amphibians. All four mtnr subtypes were found in two non-teleost actinopterygian species, the spotted gar and the reedfish. As a result of teleost tetraploidization (3R), up to seven functional mtnr genes could be identified in teleosts. Conservation of the mtnr 3R-duplicated paralogs differs among the teleost lineages. Synteny analysis showed that the mtnr1d was conserved as a singleton in all teleosts resulting from an early loss after tetraploidization of one of the teleost 3R and salmonid 4R paralogs. Several teleosts including the eels and the piranha have conserved both 3R-paralogs of mtnr1a, mtnr1b, and mtnr1c. Loss of one of the 3R-paralogs depends on the lineage: mtnr1ca was lost in euteleosts whereas mtnr1cb was lost in osteoglossomorphs and several ostariophysians including the zebrafish. We investigated the tissue distribution of mtnr expression in a large range of tissues in medaka. The medaka has conserved the four vertebrate paralogs, and these are expressed in brain and retina, and, differentially, in peripheral tissues. Photoperiod affects mtnr expression levels in a gene-specific and tissue-specific manner. This study provides new insights into the repertoire diversification and functional evolution of the mtnr gene family in vertebrates.

Independent Origin of XY and ZW Sex Determination Mechanisms in Mosquitofish Sister Species

Fish are known for the outstanding variety of their sex determination mechanisms and sex chromosome systems. The western (Gambusia affinis) and eastern mosquitofish (G. holbrooki) are sister species for which different sex determination mechanisms have been described: ZZ/ZW for G. affinis and XX/XY for G. holbrooki Here, we carried out restriction-site associated DNA (RAD-) and pool sequencing (Pool-seq) to characterize the sex chromosomes of both species. We found that the ZW chromosomes of G. affinis females and the XY chromosomes of G. holbrooki males correspond to different linkage groups, and thus evolved independently from separate autosomes. In interspecific hybrids, the Y chromosome is dominant over the W chromosome, and X is dominant over Z. In G. holbrooki, we identified a candidate region for the Y-linked melanic pigmentation locus, a rare male phenotype that constitutes a potentially sexually antagonistic trait and is associated with other such characteristics, e.g., large body size and aggressive behavior. We developed a SNP-based marker in the Y-linked allele of GIPC PDZ domain containing family member 1 (gipc1), which was linked to melanism in all tested G. holbrooki populations. This locus represents an example for a color locus that is located in close proximity to a putative sex determiner, and most likely substantially contributed to the evolution of the Y.

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.

Evolution of the Proto Sex-Chromosome in Solea senegalensis

Solea senegalensis is a flatfish belonging to the Soleidae family within the Pleuronectiformes order. It has a karyotype of 2n = 42 (FN = 60; 6M + 4 SM + 8 St + 24 T) and a XX/XY system. The first pair of metacentric chromosomes has been proposed as a proto sex-chromosome originated by a Robertsonian fusion between acrocentric chromosomes. In order to elucidate a possible evolutionary origin of this chromosome 1, studies of genomic synteny were carried out with eight fish species. A total of 88 genes annotated within of 14 BACs located in the chromosome 1 of S. senegalensis were used to elaborate syntenic maps. Six BACs (BAC5K5, BAC52C17, BAC53B20, BAC84K7, BAC56H24, and BAC48P7) were distributed in, at least, 5 chromosomes in the species studied, and a group of four genes from BAC53B20 (grsf1, rufy3, slc4a4 and npffr2) and genes from BAC48K7 (dmrt2, dmrt3, dmrt1, c9orf117, kank1 and fbp1) formed a conserved cluster in all species. The analysis of repetitive sequences showed that the number of retroelements and simple repeat per BAC showed its highest value in the subcentromeric region where 53B20, 16E16 and 48K7 BACs were localized. This region contains all the dmrt genes, which are associated with sex determination in some species. In addition, the presence of a satellite “chromosome Y” (motif length: 860 bp) was detected in this region. These findings allowed to trace an evolutionary trend for the large metacentric chromosome of S. senegalensis, throughout different rearrangements, which could be at an initial phase of differentiation as sex chromosome.

Olfactomedin-like 2 A and B (OLFML2A and OLFML2B) profile expression in the retina of spotted gar (Lepisosteus oculatus) and bioinformatics mining

Olfactomedin-like (OLFML) proteins are members of the olfactomedin domain-containing secreted glycoprotein (OLF) family. OLFML2A and OLFML2B are representative molecules of these glycoproteins. Olfactomedins are critical for the development and functional organization of the nervous system and retina, which is a highly conserved structure in vertebrates, having almost identical anatomical and physiological characteristics in multiple taxa. Spotted gar, a member of the Lepisosteidae family, is a freshwater fish that inhabits rivers, bayous, swamps, and brackish waters. Recently, the complete genome has been sequenced, providing a unique bridge between fish medical models to human biology, making it an excellent animal model. This study was aimed to understanding the evolution OLFML2A and OLFML2B in the retina of spotted gar through looking for the expression of these genes. Spotted gar retina was analyzed with hematoxylin-eosin staining assays to provide an overall view of the retina structure and an immunofluorescence assay to identify OLFML2A and OLFML2B protein expression. A phylogenetic tree was created using the neighbor-joining method. Forces that direct the evolution of the fish genes were tested. Spotted gar retina, as in other vertebrates, is made of several layers. OLFML2A and OLFML2B proteins were detected in the rod and cone photoreceptor layer (PRL), outer nuclear layer (ONL), and inner nuclear layer (INL). Phylogenetic tree analysis confirms the orthology within the OLFML2A gene. Purifying selection is the evolutionary force that directs the OLFML2A genes. OLFML2A genes have a well-conserved function over time and species.

A Hormone That Lost Its Receptor: Anti-Müllerian Hormone (AMH) in Zebrafish Gonad Development and Sex Determination

Fetal mammalian testes secrete Anti-Müllerian hormone (Amh), which inhibits female reproductive tract (Müllerian duct) development. Amh also derives from mature mammalian ovarian follicles, which marks oocyte reserve and characterizes polycystic ovarian syndrome. Zebrafish (Danio rerio) lacks Müllerian ducts and the Amh receptor gene amhr2 but, curiously, retains amh To discover the roles of Amh in the absence of Müllerian ducts and the ancestral receptor gene, we made amh null alleles in zebrafish. Results showed that normal amh prevents female-biased sex ratios. Adult male amh mutants had enormous testes, half of which contained immature oocytes, demonstrating that Amh regulates male germ cell accumulation and inhibits oocyte development or survival. Mutant males formed sperm ducts and some produced a few offspring. Young female mutants laid a few fertile eggs, so they also had functional sex ducts. Older amh mutants accumulated nonvitellogenic follicles in exceedingly large but sterile ovaries, showing that Amh helps control ovarian follicle maturation and proliferation. RNA-sequencing data partitioned juveniles at 21 days postfertilization (dpf) into two groups that each contained mutant and wild-type fish. Group21-1 upregulated ovary genes compared to Group21-2, which were likely developing as males. By 35 dpf, transcriptomes distinguished males from females and, within each sex, mutants from wild types. In adult mutants, ovaries greatly underexpressed granulosa and theca genes, and testes underexpressed Leydig cell genes. These results show that ancestral Amh functions included development of the gonadal soma in ovaries and testes and regulation of gamete proliferation and maturation. A major gap in our understanding is the identity of the gene encoding a zebrafish Amh receptor; we show here that the loss of amhr2 is associated with the breakpoint of a chromosome rearrangement shared among cyprinid fishes.

Analyzing the signaling properties of gar (Lepisosteus oculatus) melanocortin receptors: Evaluating interactions with MRAP1 and MRAP2

RT-PCR analysis of gar pituitary and brain indicated that different combinations of gar melanocortin receptor mRNAs are present in the same tissues with mRNAs for gar mrap1 and gar mrap2. Against this background, an objective of this study was to determine whether the ligand sensitivity for either ACTH or α-MSH was affected when gar (g) melanocortin receptors (Mcrs) were co-expressed with either of the accessory proteins gMrap1 or gMrap2 in Chinese Hamster Ovary cells. The results indicated that gMc2r has an obligatory requirement for co-expression with gMrap1 in order for the receptor to be activated by hACTH(1-24). In addition, activation of gMc2r did not occur when the receptor was expressed alone or co-expressed with gMrap2. Furthermore, co-expression of gMc2r with gMrap1 followed by stimulation with NDP-MSH resulted in a low level of activation (only at 10-7 M and 10-6 M). However, gMc1r, gMc3r, gMc4r, and gMc5r responded to stimulation by NDP-MSH in a more robust manner. Co-expression of gMc1r, gMc3r, gMc4r, and gMc5r with gMRAP1 had no effect on sensitivity to stimulation by NDP-MSH or hACTH(1-24). Co-expression with gMRAP2 had no negative or positive effect on ligand sensitivity for gMc1r, gMc3r, and gMc5r, however this treatment did increase the activation of CHO cells transfected with gMc4r following stimulation with both hACTH(1-24) (p < 0.001), and NDP-MSH (p < 0.001). Co-expression of gMC5R with either gMRAP1 or gMRAP2 increased trafficking of gMC5R to the plasma membrane. These pharmacological observations are compared to the response of melanocortin receptors from other neopterygian fishes, cartilaginous fishes, and tetrapods to stimulation by ACTH(1-24) and forms of α-MSH.

Functional differentiation of three phosphatidylinositol 3-kinase catalytic subunit alpha (PIK3CA) in response to Vibrio anguillarum infection in turbot (Scophthalmus maximus)

PIK3CA has been extensively investigated from its molecular mechanism perspective and association with its mutations in different types of cancers. However, little has been reported regarding the pathological significance of PIK3CA expression in teleost. Here, in our present study, three PIK3CA genes termed SmPIK3CAa, SmPIK3CAb and SmPIK3CA-like were firstly identified in the genome of turbot S. maximus. Although these three genes located in different chromosomes, all of them share the same five domains. Phylogenetic and synteny analysis indicated that SmPIK3CAa, SmPIK3CAb and SmPIK3CA-like were three paralogs that may originate from duplication of the same ancestral PIK3CA gene. Subcellular localization analysis confirmed the cytoplasm distribution of these three paralogs. All three SmPIK3CA were ubiquitously expressed in examined tissues in turbot, with the higher expression levels in immune-related tissues such as blood, spleen, kidney, gills and intestines. Upon Vibrio anguillarum challenge, SmPIK3CAa and SmPIK3CA-like transcripts were significantly induced in spleen, intestine and blood despite of differential expression levels and responsive time points. Additionally, individuals in resistant group showed significantly higher expression level of both two genes than in the susceptible group. Moreover, four SNPs (102, 2530, 3027 and 3060) and one haplotype (Hap2) located in exon region of SmPIK3CA-like were identified and confirmed to be associated with V. anguillarum resistance in turbot by association analysis in different populations. Taken together, these results suggested that functional differentiation occurred in three SmPIK3CA paralogs with Vibrio anguillarum resistance and SmPIK3CAa and SmPIK3CA-like probable play potential roles in innate immune response to pathogenic invasions in turbot.