FishmiRNA: An evolutionarily supported microRNA annotation and expression database for ray-finned fishes

A maternal-to-zygotic-transition gene block on the zebrafish sex chromosome

Wild zebrafish (Danio rerio) have a ZZ/ZW chromosomal sex-determination system with the major sex locus on the right arm of chromosome-4 (Chr4R) near the largest heterochromatic block in the genome, suggesting that Chr4R transcriptomics might differ from the rest of the genome. To test this hypothesis, we conducted an RNA-seq analysis of adult ZW ovaries and ZZ testes in the Nadia strain and identified 4 regions of Chr4 with different gene expression profiles. Unique in the genome, protein-coding genes in a 41.7 Mb section (Region-2) were expressed in testis but silent in ovary. The AB lab strain, which lacks sex chromosomes, verified this result, showing that testis-biased gene expression in Region-2 depends on gonad biology, not on sex-determining mechanism. RNA-seq analyses in female and male brains and livers validated reduced transcripts from Region-2 in somatic cells, but without sex specificity. Region-2 corresponds to the heterochromatic portion of Chr4R and its content of genes and repetitive elements distinguishes it from the rest of the genome. Region-2 lacks protein-coding genes with human orthologs; has zinc finger genes expressed early in zygotic genome activation; has maternal 5S rRNA genes, maternal spliceosome genes, a concentration of tRNA genes, and a distinct set of repetitive elements. The colocalization of (1) genes silenced in ovaries but not in testes that are (2) expressed in embryos briefly at the onset of zygotic genome activation; (3) maternal-specific genes for translation machinery; (4) maternal-specific spliceosome components; and (5) adjacent genes encoding miR-430, which mediates maternal transcript degradation, suggest that this is a maternal-to-zygotic-transition gene regulatory block.

A comprehensive coding and microRNA transcriptome of vertebral bone in postlarvae and juveniles of Senegalese sole (Solea senegalensis)

Understanding vertebral bone development is essential to prevent skeletal malformations in farmed fish related to genetic and environmental factors. This is an important issue in Solea senegalensis, with special impact of spinal anomalies in postlarval and juvenile stages. Vertebral bone transcriptomics in farmed fish mainly comes from coding genes, and barely on miRNA expression. Here, we used RNA-seq of spinal samples to obtain the first comprehensive coding and miRNA transcriptomic repertoire for postlarval and juvenile vertebral bone, covering different vertebral phenotypes and egg-incubation temperatures related to skeleton health in S. senegalensis. Coding genes, miRNA and pathways regulating bone development and growth were identified. Differential transcriptomic profiles and suggestive mRNA-miRNA interactions were found between postlarvae and juveniles. Bone-related genes and functions were associated with the extracellular matrix, development and regulatory processes, calcium binding, retinol and lipid metabolism or response to stimulus, including those revealed by the miRNA targets related to signaling, cellular and metabolic processes, growth, cell proliferation and biological adhesion. Pathway enrichment associated with fish skeleton were identified when comparing postlarvae and juveniles: growth and bone development functions in postlarvae, while actin cytoskeleton, focal adhesion and proteasome related to bone remodeling in juveniles. The transcriptome data disclosed candidate coding and miRNA gene markers related to bone cell processes, references for functional studies of the anosteocytic bone of S. senegalensis. This study establishes a broad transcriptomic foundation to study healthy and anomalous spines under early thermal conditions across life-stages in S. senegalensis, and for comparative analysis of skeleton homeostasis and pathology in fish and vertebrates.

miR-210 promotes immune- and suppresses oocyte meiosis-related genes in the zebrafish ovarian cells

microRNA-210 (miRNA), a well-documented miRNA, has been implicated in a myriad of biological processes, including responses to hypoxia, angiogenesis, cell proliferation, and male infertility in humans. However, a comprehensive understanding of its functions in fish requires further investigation. This study pursued to elucidate the downstream effect of dre-miR-210-5p on primary ovarian cell culture in zebrafish (Danio rerio), an animal model. A protocol was settled down by incubations with either an miR-210 mimic or a scrambled miRNA in the isolated ovaries. RNA-sequencing analysis identified ∼6000 differentially expressed target genes revealing that downregulated genes were associated with reproduction-related pathways while immune-related pathways displayed an upregulated pattern. To identify molecular markers, predicted target genes were classified into reproduction and immune cell types. These findings underscore the existence of a profound interplay between the reproductive and immune systems, with miR-210 emerging as a pivotal player in orchestrating transcriptomic alterations within fish ovaries.

Taste evolution in an herbivorous drosophilid

Plant secondary metabolites pose a challenge for generalist herbivorous insects because they are not only potentially toxic, they also may trigger aversion. On the contrary, some highly specialized herbivorous insects evolved to use these same compounds as 'token stimuli' for unambiguous determination of their host plants. Two questions that emerge from these observations are how recently derived herbivores evolve to overcome this aversion to plant secondary metabolites and the extent to which they evolve increased attraction to these same compounds. In this study, we addressed these questions by focusing on the evolution of bitter taste preferences in the herbivorous drosophilid Scaptomyza flava, which is phylogenetically nested deep in the paraphyletic Drosophila. We measured behavioral and neural responses of S. flava and a set of non-herbivorous species representing a phylogenetic gradient (S. pallida, S. hsui, and D. melanogaster) towards host- and non-host derived bitter plant compounds. We observed that S. flava evolved a shift in bitter detection, rather than a narrow shift towards glucosinolates, the precursors of mustard-specific defense compounds. In a dye-based consumption assay, S. flava exhibited shifts in aversion toward the non-mustard bitter, plant-produced alkaloids caffeine and lobeline, and reduced aversion towards glucosinolates, whereas the non-herbivorous species each showed strong aversion to all bitter compounds tested. We then examined whether these changes in bitter preferences of S. flava could be explained by changes in sensitivity in the peripheral nervous system and compared electrophysiological responses from the labellar sensilla of S. flava, S. pallida, and D. melanogaster. Using scanning electron microscopy, we also created a map of labellar sensilla in S. flava and S. pallida. We assigned each sensillum to a functional sensilla class based on their morphology and initial response profiles to bitter and sweet compounds. Despite a high degree of conservation in the morphology and spatial placement of sensilla between S. flava and S. pallida, electrophysiological studies revealed that S. flava had reduced sensitivity to glucosinolates to varying degrees. We found this reduction only in I type sensilla. Finally, we speculate on the potential role that evolutionary genetic changes in gustatory receptors between S. pallida and S. flava may play in driving these patterns. Specifically, we hypothesize that the evolution of bitter receptors expressed in I type sensilla may have driven the reduced sensitivity observed in S. flava, and ultimately, its reduced bitter aversion. The S. flava system showcases the importance of reduced aversion to bitter defense compounds in relatively young herbivorous lineages, and how this may be achieved at the molecular and physiological level.

Looking for a needle in a haystack: de novo phenotypic target identification reveals Hippo pathway-mediated miR-202 regulation of egg production

Understanding microRNA (miRNA) functions has been hampered by major difficulties in identifying their biological target(s). Currently, the main limitation is the lack of a suitable strategy to identify biologically relevant targets among a high number of putative targets. Here we provide a proof of concept of successful de novo (i.e. without prior knowledge of its identity) miRNA phenotypic target (i.e. target whose de-repression contributes to the phenotypic outcomes) identification from RNA-seq data. Using the medaka mir-202 knock-out (KO) model in which inactivation leads to a major organism-level reproductive phenotype, including reduced egg production, we introduced novel criteria including limited fold-change in KO and low interindividual variability in gene expression to reduce the list of 2853 putative targets to a short list of 5. We selected tead3b, a member of the evolutionarily-conserved Hippo pathway, known to regulate ovarian functions, due to its remarkably strong and evolutionarily conserved binding affinity for miR-202-5p. Deleting the miR-202-5p binding site in the 3' UTR of tead3b, but not of other Hippo pathway members sav1 and vgll4b, triggered a reduced egg production phenotype. This is one of the few successful examples of de novo functional assignment of a miRNA phenotypic target in vivo in vertebrates.

A medaka gonad-specific lncRNA may act as pri-miR-202 to regulate testicular endocrine homeostasis and spermatogenesis

A large number of long non-coding RNAs (lncRNAs) are expressed in animal gonads, but their functions are poorly understood. In this study, a gonad-specific lncRNA, termed lnc4, was identified and characterized in the model fish medaka (Oryzias latipes). The expression pattern and in vitro functional analyses indicated that lnc4 was likely to be a primary transcript of miR-202 (pri-miR-202). Results of single-molecule fluorescence in situ hybridization demonstrated that the precursor miR-202 (pre-miR-202) was highly expressed in the nuclei of testicular somatic cells, including Leydig and Sertoli cells, whereas only a small amount of lnc4 molecules could be detected co-expressed with pre-miR-202 in Sertoli cells due to its low expression level. Deletion of the lnc4 locus led to a significant reduction in testis size and a dramatic decrease in the number of male germ cells, as well as a reduction in sperm viability. Moreover, lnc4 knockout resulted in enhanced synthesis and secretion of testicular somatic cells and accelerated differentiation of immature male germ cells. Taken together, functional studies of lnc4 and its mature transcript miR-202 will contribute to the understanding of the important role of non-coding RNAs in animal or human reproductive disorders.

The Discovery and Characterization of Conserved and Novel miRNAs in the Different Developmental Stages and Organs of Pikeperch (Sander lucioperca)

Micro RNAs (miRNAs) are short non-coding RNAs that act as post-transcriptional gene expression regulators. Genes regulated in vertebrates include those affecting growth and development or stress and immune response. Pikeperch (Sander lucioperca) is a species that is increasingly being considered for farming in recirculation aquaculture systems. We characterized the pikeperch miRNA repertoire to increase the knowledge of the genomic mechanisms affecting performance and health traits by applying small RNA sequencing to different developmental stages and organs. There were 234 conserved and 8 novel miRNA genes belonging to 104 families. A total of 375 unique mature miRNAs were processed from these genes. Many mature miRNAs showed high relative abundances or were significantly more expressed at early developmental stages, like the miR-10 and miR-430 family, let-7, the miRNA clusters 106-25-93, and 17-19-92. Several miRNAs associated with immune responses (e.g., slu-mir-731-5p, slu-mir-2188-5p, and slu-mir-8159-5p) were enriched in the spleen. The mature miRNAs slu-mir-203a-3p and slu-mir-205-5p were enriched in gills. These miRNAs are similarly abundant in many vertebrates, indicating that they have shared regulatory functions. There was also a significantly increased expression of the disease-associated miR-462/miR-731 cluster in response to hypoxia stress. This first pikeperch miRNAome reference resource paves the way for future functional studies to identify miRNA-associated variations that can be utilized in marker-assisted breeding programs.

A maternal-to-zygotic-transition gene block on the zebrafish sex chromosome

Wild zebrafish (Danio rerio) have a ZZ/ZW chromosomal sex determination system with the major sex locus on the right arm of chromosome-4 (Chr4R) near the largest heterochromatic block in the genome, suggesting the hypothesis that the Chr4R transcriptome might be different from the rest of the genome. We conducted an RNA-seq analysis of adult ZW ovaries and ZZ testes and identified four regions of Chr4 with different gene expression profiles. Unique in the genome, protein-coding genes in a 41.7 Mb section (Region-2) were expressed in testis but silent in ovary. The AB lab strain, which lacks sex chromosomes, verified this result, showing that testis-biased gene expression in Region-2 depends on gonad biology, not on sex-determining mechanism. RNA-seq analyses in female and male brain and liver validated few transcripts from Region-2 in somatic cells, but without sex-specificity. Region-2 corresponds to the heterochromatic portion of Chr4R and its content of genes and repetitive elements distinguishes it from the rest of the genome. In Region-2, protein-coding genes lack human orthologs; it has zinc finger genes expressed early in zygotic genome activation; it has maternal 5S rRNA genes, maternal spliceosome genes, a concentration of tRNA genes, and an distinct set of repetitive elements. The colocalization of 1) genes silenced in ovaries but not in testes that are 2) expressed in embryos briefly at the onset of zygotic genome activation; 3) maternal-specific genes for translation machinery; 4) maternal-specific spliceosome components; and 4) adjacent genes encoding miR-430, which mediates maternal transcript degradation, suggest that this is a Maternal-to-Zygotic-Transition Gene Regulatory Block.

Deciphering sex-specific miRNAs as heat-recorders in zebrafish

In the last decade, a plethora of microRNAs (miRNAs) has been reported in a wide variety of physiological processes, including reproduction, in many aquatic organisms. However, miRNAome alterations occurred by environmental cues due to water temperature increment have not yet been elucidated. With the aim to identify epigenetic regulations mediated by miRNAs in the gonads in a climate change scenario, the animal model zebrafish (Danio rerio) were subjected to high temperatures during sex differentiation, a treatment that results in male-skewed sex ratios in the adulthood. Once the fish reached adulthood, gonads were sequenced by high-throughput technologies and a total of 23 and 1 differentially expressed miRNAs in ovaries and testes, respectively, were identified two months after the heat treatment. Most of these heat-recorder miRNAs were involved in human sex-related cancer and about 400 predicted-target genes were obtained, some with reproduction-related functions. Their synteny in the zebrafish genome was, for more than half of the predicted target genes, in the chromosomes 7, 2, 4, 3 and 11 in the ovaries, chromosome 4 being the place where the sex-associated-region (sar) is localized in wild zebrafish. Further, spatial localization in the gonads of two selected heat-recorder miRNAs (miR-122-5p and miR-146-5p) showed exclusive expression in the ovarian germ cells. The present study expands the catalog of sex-specific miRNAs and deciphers, for the first time, thermosensitive miRNAs in the zebrafish gonads that might be used as potential epimarkers to predict environmental past events.

An ancient truncated duplication of the anti‐Mullerian hormone receptor type 2 gene is a potential conserved master sex determinant in the Pangasiidae catfish family

The evolution of sex determination (SD) in teleosts is amazingly dynamic, as reflected by the variety of different master sex-determining genes identified. Pangasiids are economically important catfishes in South Asian countries, but little is known about their SD system. Here, we generated novel genomic resources for 12 Pangasiids and characterized their SD system. Based on a Pangasianodon hypophthalmus chromosome-scale genome assembly, we identified an anti-Müllerian hormone receptor type Ⅱ gene (amhr2) duplication, which was further characterized as being sex-linked in males and expressed only in testes. These results point to a Y chromosome male-specific duplication (amhr2by) of the autosomal amhr2a. Sequence annotation revealed that the P. hypophthalmus Amhr2by is truncated in its N-terminal domain, lacking the cysteine-rich extracellular part of the receptor that is crucial for ligand binding, suggesting a potential route for its neofunctionalization. Reference-guided assembly of 11 additional Pangasiids, along with sex-linkage studies, revealed that this truncated amhr2by duplication is a male-specific conserved gene in Pangasiids. Reconstructions of the amhr2 phylogeny suggested that amhr2by arose from an ancient duplication/insertion event at the root of the Siluroidei radiation that is dated to ~100 million years ago. Together these results bring multiple lines of evidence supporting that amhr2by is an ancient and conserved master sex-determining gene in Pangasiids, a finding that highlights the recurrent use of the transforming growth factor β pathway, which is often used for the recruitment of teleost master SD genes, and provides another empirical case towards firther understanding of dynamics of SD systems.