Global analysis of the small RNA transcriptome in different ploidies and genomic combinations of a vertebrate complex--the Squalius alburnoides

Changes in Ploidy Drive Reproduction Transition and Genomic Diversity in a Polyploid Fish Complex

Unisexual animals are commonly found in some polyploid species complexes, and most of these species have had a long evolutionary history. However, their method for avoiding genomic decay remains unclear. The polyploid Carassius complex naturally comprises the sexual amphidiploid C. auratus (crucian carp or goldfish) (AABB) and the gynogenetic amphitriploid C. gibelio (gibel carp) (AAABBB). Recently, we developed a fertile synthetic amphitetraploid (AAAABBBB) male from C. gibelio by incorporating a C. auratus genome. In this study, we generated novel amphitriploids (AAABBB) by backcrossing the amphitetraploid male with the amphidiploid C. auratus. Whole-genome resequencing revealed the genomic changes, including recombination and independent assortment between homologs of C. gibelio and C. auratus. The fertility, sex determination system, oocyte development, and fertilization behaviors of the novel amphitriploids were investigated. Approximately 80% of the novel amphitriploid females recovered the unisexual gynogenesis ability. Intriguingly, two types of primary oocyte (with and without homolog synapsis) were discovered, and their distinct development fates were observed. Type I oocytes entered apoptosis due to improper synaptonemal complex assembly and incomplete double-strand break repair, whereas subsequent type II oocytes bypassed meiosis through an alternative ameiotic pathway to develop into mature eggs. Moreover, gynogenesis was stabilized in their offspring, and a new array of diverse gynogenetic amphitriploid clones was produced. These revealed genomic changes and detailed cytological data provide comprehensive evidence that changes in ploidy drive unisexual and sexual reproduction transition, thereby resulting in genomic diversity and allowing C. gibelio avoid genomic decay.

Cytogenetic mechanisms of unisexuality in rock lizards

Darevskia rock lizards is a unique complex taxa, including more than thirty species, seven of which are parthenogenetic. In mixed populations of Darevskia lizards, tri- and tetraploid forms can be found. The most important issues in the theory of reticulate evolution of Darevskia lizards are the origin of parthenogenetic species and their taxonomic position. However, there is little data on how meiosis proceeds in these species. The present work reports the complex results of cytogenetics in a diploid parthenogenetic species – D. unisexualis. Here we detail the meiotic prophase I progression and the specific features оf mitotic chromosomes organization. The stages of meiosis prophase I were investigated by immunocytochemical analysis of preparations obtained from isolated primary oocytes of D. unisexualis in comparison with maternal species D. raddei nairensis. It has been shown that in D. unisexualis at the leptotene-zygotene stages the axial elements and the synaptonemal complex (SC) form typical “bouquets”. At the pachytene-diplotene stage, 18 autosomal SC-bivalents and thickened asynapted sex Z and w univalents were observed. The presence of SYCP1 protein between the lateral elements of autosomal chromosomes proved the formation of assembled SCs. Comparative genomic hybridization (CGH) on the mitotic metaphase chromosomes of D. unisexualis was carried out using the genomic DNA isolated from the parental species D. raddei nairensis and D. valentini. In the pericentromeric regions of half of the mitotic chromosomes of D. unisexualis, specific regions inherited from maternal species have been found. Following our results, we suggest a model for diploid germ cells formation from diploid oocytes without premeiotic duplication of chromosomes in the oogenesis of diploid parthenogenetic lizards D. unisexualis. Taken as a whole, our findings confirm the hybrid nature of D. unisexualis and shed light on heterozygosity and automixis in diploid parthenogenetic forms.

Integrated Analysis of Transcriptomic, miRNA and Proteomic Changes of a Novel Hybrid Yellow Catfish Uncovers Key Roles for miRNAs in Heterosis

Heterosis is a complex biological phenomenon in which hybridization produces offspring that exhibit superior phenotypic characteristics compared with the parents. Heterosis is widely utilized in agriculture, for example in fish farming; however, its underlying molecular basis remains elusive. To gain a comprehensive and unbiased molecular understanding of fish heterosis, we analyzed the mRNA, miRNA, and proteomes of the livers of three catfish species, Pelteobagrus fulvidraco, P. vachelli, and their hybrid, the hybrid yellow catfish "Huangyou-1" (P. fulvidraco ♀ × P. vachelli ♂). Using next-generation sequencing and mass spectrometry, we show that the nonadditive, homoeolog expression bias and expression level dominance pattern were readily identified at the transcriptional, post-transcriptional, or protein levels, providing the evidence for the widespread presence of dominant models during hybridization. A number of predicted miRNA-mRNA-protein pairs were found and validated by qRT-PCR and PRM assays. Furthermore, several diverse key pathways were identified, including immune defense, metabolism, digestion and absorption, and cell proliferation and development, suggesting the vital mechanisms involved in the generation of the heterosis phenotype in progenies. We propose that the high parental expression of genes/proteins (growth, nutrition, feeding, and disease resistance) coupled with low parental miRNAs of the offspring, are inherited from the mother or father, thus indicating that the offspring were enriched with the advantages of the father or mother. We provide new and important information about the molecular mechanisms of heterosis, which represents a significant step toward a more complete elucidation of this phenomenon.

MicroRNA Alternations in the Testes Related to the Sterility of Triploid Fish

The knowledge of understanding the molecular traits of the sterile triploid fish is sparse. Herein, we analyzed the microRNA (miRNA) alternations in the testes of the sterile triploid fish produced by crossing the tetraploid fish with the diploid fish, compared with those of tetraploids and diploids used as the controls. A total of 136, 134, and 142 conserved miRNAs and 105, 112, and 119 novel miRNAs were identified in the diploid, triploid, and tetraploid fish, respectively. The genes targeted by the differentially expressed miRNAs were identified and were enriched in the GO term cell surface receptor signaling pathway, cellular process, G-protein coupled receptor signaling pathway, and metabolic process. KEGG pathway enrichment was also assessed to evaluate the target genes with differentially expressed miRNAs and these genes were enriched in four pathways (synthesis and degradation of ketone bodies, pentose and glucuronate interconversions, cyanoamino acid metabolic process, and ascorbate and aldarate metabolism). Nine differentially expressed miRNAs were verified by quantitative real-time PCR analysis (qPCR). The upregulated miRNAs in triploids, including miR-101a, miR-199-5p, miR-214, miR-222, and miR-193a, showed the same results with high-throughput sequencing. Among the selected downregulated miRNAs, miR-7b and miR-153b had significantly lower expression levels in triploids. Dnah3 and Tekt1 genes targeted by miR-199-5p showed lower expression in triploids by qPCR. These verified differentially expressed miRNAs may participate in testicular development and sperm activity by targeting functional genes, which were identified with differential expression in the triploid. This evidence provides insights into the epigenetic regulatory mechanisms of sterility in triploid cyprinids.

Gene copy silencing and DNA methylation in natural and artificially produced allopolyploid fish

Allelic silencing is an important mechanism for coping with gene dosage changes in polyploid organisms that is well known in allopolyploid plants. Only recently, it was shown in the allotriploid fish Squalius alburnoides that this process also occurs in vertebrates. However, it is still unknown whether this silencing mechanism is common to other allopolyploid fish, and which mechanisms might be responsible for allelic silencing. We addressed these questions in a comparative study between Squalius alburnoides and another allopolyploid complex, the Amazon molly (Poecilia formosa). We examined the allelic expression patterns for three target genes in four somatic tissues of natural allo-anorthoploids and laboratory-produced tri-genomic hybrids of S. alburnoides and P. formosa. Also, for both complexes, we evaluated the correlation between total DNA methylation level and the ploidy status and genomic composition of the individuals. We found that allelic silencing also occurs in other allopolyploid organisms besides the single one that was previously known. We found and discuss disparities within and between the two considered complexes concerning the pattern of allele-specific expression and DNA methylation levels. Disparities might be due to intrinsic characteristics of each genome involved in the hybridization process. Our findings also support the idea that long-term evolutionary processes have an effect on the allele expression patterns and possibly also on DNA methylation levels.

Dynamics of Rex3 in the genomes of endangered Iberian Leuciscinae (Teleostei, Cyprinidae) and their natural hybrids

Background

Iberian Leuciscinae are greatly diverse comprising taxa of hybrid origin. With highly conservative karyotypes, Iberian Chondrostoma s.l. have recently demonstrated sub-chromosomal differentiation and rapid genome restructuring in natural hybrids, which was confirmed by ribosomal DNA (rDNA) transposition and/or multiplication. To understand the role of repetitive DNAs in the differentiation of their genomes, a genetic and molecular cytogenetic survey was conducted in Achondrostoma oligolepis, Anaecypris hispanica, Iberochondrostoma lemmingii, I. lusitanicum, Pseudochondrostoma duriense, P. polylepis, Squalius pyrenaicus and hybrids between A. oligolepis x (P. duriense/P. polylepis), representing ‘alburnine’, chondrostomine and Squalius lineages.

Results

Partial Rex3 sequences evidenced high sequence homology among Leuciscinae (≥98 %) and different fish families (80–95 %) proposing a relatively recent activity of these elements in the species inspected. Low nucleotide substitution rates (<20 %) and intact ORFs suggests that Rex3 may in fact be active in these genomes. The chromosomal distribution of Rex3 retroelement was found highly concentrated at pericentromeric and moderately at subtelomeric blocks, co-localizing with 5S rDNA loci, and correlating with blocks of heterochromatin and C₀t-1 DNA. This accumulation was evident in at least 10 chromosome pairs, a pattern that seemed to be shared among the different species, likely pre-dating their divergence. Nevertheless, species-specific clusters were detected in I. lusitanicum, P. duriense, P. polylepis and S. pyrenaicus demonstrating rapid and independent differentiation. Natural hybrids followed the same patterns of accumulation and association with repetitive sequences. An increased number of Rex3 clusters now associating also with translocated 45S rDNA clusters vouched for other genomic rearrangements in hybrids. Rex3 sequence phylogeny did not agree with its hosts’ phylogeny but the observed distribution pattern is congruent with an evolutionary tendency to protect its activity, a robust regulatory system and/or events of horizontal transfer.

Conclusions

This is the first report directed at retroelement physical mapping in Cyprinidae. It helped outlining conceivable ancestral homologies and recognizing retrotransposon activation in hybrids, being possibly associated with genome diversification within the subfamily. The extensive diversity of Iberian Leuciscinae makes them excellent candidates to explore the processes and mechanisms behind the great plasticity distinguishing vertebrate genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13039-015-0180-1) contains supplementary material, which is available to authorized users.

MicroRNA profiles reveal female allotetraploid hybrid fertility

BACKGROUND

The bisexual fertile tetraploid fish is important in biological evolution. Tetraploid fish fertility is the key factor for stable inheritance. Therefore, elucidating tetraploid fish fertility at the molecular level is essential. MicroRNAs regulate gene expression and are involved in many aspects of gonad development.

METHODS

Total RNA was isolated using TRIzol, followed by constructing small RNA libraries. And then, the qualified libraries were sequenced with the HiSeq 2500 SE50 system. The obtained clean reads were analyzed to identify conserved and novel miRNAs, and evaluate the expression, and also predict the target genes. The differential expressions of miRNAs were confirmed by RT-PCR.

RESULTS

In this study, allotetraploid hybrid fish (4nAT) and diploid red crucian carp (RCC) ovaries were used to compare miRNA profiles. The results indicated that most of the highly expressed miRNAs were closely correlated with ovary maturation, and displayed no significant differences in expression. Moreover, 34 up-regulated and nine down-regulated miRNAs were found in 4nAT. The differentially expressed miRNAs were primarily involved in metabolism, defense mechanisms, and cytoskeleton production.

CONCLUSIONS

This is the first study to provide new epigenetic evidences for tetraploid fish fertility and phenotypic changes as a result of increased ploidy.

Global transcriptional and miRNA insights into bases of heterosis in hybridization of Cyprinidae

Hybrid Megalobrama amblycephala × Culter alburnus represents a population newly formed by interspecific crossing between two different genera. Here we assessed the expression pattern of mRNA and small RNA in newly formed F1, F2 and their progenitors. Large amounts of nonadditively expressed protein-coding genes showed parental expression level dominance (ELD). Interestingly, the ELD pattern could inherit from F1 to F2, which guaranteed a stable appearance in progenies. The ELD-B genes were found to contribute to cell development, while the ELD-T genes were enriched in function of stress and adaptability. microRNAs (miRNA) also had similar expression patterns to genes. A high proportion of miRNAs showed nonadditive expression upon hybridization, and were found to target important genes with diverse roles potentially involved in stress adaption and development. Taken together, the gene and miRNA expression divergence contributes to heterosis in the newly formed hybrid, promising the successful existence of hybrid speciation.

Gene expression dosage regulation in an allopolyploid fish

How allopolyploids are able not only to cope but profit from their condition is a question that remains elusive, but is of great importance within the context of successful allopolyploid evolution. One outstanding example of successful allopolyploidy is the endemic Iberian cyprinid Squalius alburnoides. Previously, based on the evaluation of a few genes, it was reported that the transcription levels between diploid and triploid S. alburnoides were similar. If this phenomenon occurs on a full genomic scale, a wide functional ‘‘diploidization’’ could be related to the success of these polyploids. We generated RNA-seq data from whole juvenile fish and from adult livers, to perform the first comparative quantitative transcriptomic analysis between diploid and triploid individuals of a vertebrate allopolyploid. Together with an assay to estimate relative expression per cell, it was possible to infer the relative sizes of transcriptomes. This showed that diploid and triploid S. alburnoides hybrids have similar liver transcriptome sizes. This in turn made it valid to directly compare the S. alburnoides RNA-seq transcript data sets and obtain a profile of dosage responses across the S. alburnoides transcriptome. We found that 64% of transcripts in juveniles’ samples and 44% in liver samples differed less than twofold between diploid and triploid hybrids (similar expression). Yet, respectively 29% and 15% of transcripts presented accurate dosage compensation (PAA/PA expression ratio of 1 instead of 1.5). Therefore, an exact functional diploidization of the triploid genome does not occur, but a significant down regulation of gene expression in triploids was observed. However, for those genes with similar expression levels between diploids and triploids, expression is not globally strictly proportional to gene dosage nor is it set to a perfect diploid level. This quantitative expression flexibility may be a strong contributor to overcome the genomic shock, and be an immediate evolutionary advantage of allopolyploids.

Rapid genomic DNA changes in allotetraploid fish hybrids

Rapid genomic change has been demonstrated in several allopolyploid plant systems; however, few studies focused on animals. We addressed this issue using an allotetraploid lineage (4nAT) of freshwater fish originally derived from the interspecific hybridization of red crucian carp (Carassius auratus red var., ♀, 2n=100) × common carp (Cyprinus carpio L., ♂, 2n=100). We constructed a bacterial artificial chromosome (BAC) library from allotetraploid hybrids in the 20th generation (F20) and sequenced 14 BAC clones representing a total of 592.126 kb, identified 11 functional genes and estimated the guanine-cytosine content (37.10%) and the proportion of repetitive elements (17.46%). The analysis of intron evolution using nine orthologous genes across a number of selected fish species detected a gain of 39 introns and a loss of 30 introns in the 4nAT lineage. A comparative study based on seven functional genes among 4nAT, diploid F1 hybrids (2nF1) (first generation of hybrids) and their original parents revealed that both hybrid types (2nF1 and 4nAT) not only inherited genomic DNA from their parents, but also demonstrated rapid genomic DNA changes (homoeologous recombination, parental DNA fragments loss and formation of novel genes). However, 4nAT presented more genomic variations compared with their parents than 2nF1. Interestingly, novel gene fragments were found for the iqca1 gene in both hybrid types. This study provided a preliminary genomic characterization of allotetraploid F20 hybrids and revealed evolutionary and functional genomic significance of allopolyploid animals.

Characterization and expression analysis of microRNAs in the tube foot of sea cucumber Apostichopus japonicus

MicroRNAs (miRNAs) are a class of endogenous non-coding small RNA with average length of 22 nucleotides, participating in the post-transcriptional regulation of gene expression. In this study, we report the identification and characterization of miRNAs in the tube foot of sea cucumber (Apostichopus japonicus) by next generation sequencing with Illumina HiSeq 2000 platform. Through the bioinformatic analysis, we identified 260 conserved miRNAs and six novel miRNAs from the tube foot small RNA transcriptome. Quantitative realtime PCR (qRT-PCR) was performed to characterize the specific expression in the tube foot. The results indicated that four miRNAs, including miR-29a, miR-29b, miR-2005 and miR-278-3p, were significantly up-regulated in the tube foot. The target genes of the four specifically expressed miRNAs were predicted in silico and validated by performing qRT-PCR. Gene ontology (GO) and KEGG pathway analyses with the target genes of these four miRNAs were conducted to further understand the regulatory function in the tube foot. This is the first study to profile the miRNA transcriptome of the tube foot in sea cucumber. This work will provide valuable genomic resources to understand the mechanisms of gene regulation in the tube foot, and will be useful to assist the molecular breeding in sea cucumber.

Novel method for analysis of allele specific expression in triploid Oryzias latipes reveals consistent pattern of allele exclusion

Assessing allele-specific gene expression (ASE) on a large scale continues to be a technically challenging problem. Certain biological phenomena, such as X chromosome inactivation and parental imprinting, affect ASE most drastically by completely shutting down the expression of a whole set of alleles. Other more subtle effects on ASE are likely to be much more complex and dependent on the genetic environment and are perhaps more important to understand since they may be responsible for a significant amount of biological diversity. Tools to assess ASE in a diploid biological system are becoming more reliable. Non-diploid systems are, however, not uncommon. In humans full or partial polyploid states are regularly found in both healthy (meiotic cells, polynucleated cell types) and diseased tissues (trisomies, non-disjunction events, cancerous tissues). In this work we have studied ASE in the medaka fish model system. We have developed a method for determining ASE in polyploid organisms from RNAseq data and we have implemented this method in a software tool set. As a biological model system we have used nuclear transplantation to experimentally produce artificial triploid medaka composed of three different haplomes. We measured ASE in RNA isolated from the livers of two adult, triploid medaka fish that showed a high degree of similarity. The majority of genes examined (82%) shared expression more or less evenly among the three alleles in both triploids. The rest of the genes (18%) displayed a wide range of ASE levels. Interestingly the majority of genes (78%) displayed generally consistent ASE levels in both triploid individuals. A large contingent of these genes had the same allele entirely suppressed in both triploids. When viewed in a chromosomal context, it is revealed that these genes are from large sections of 4 chromosomes and may be indicative of some broad scale suppression of gene expression.

Epigenetic control of mobile DNA as an interface between experience and genome change

Mobile DNA in the genome is subject to RNA-targeted epigenetic control. This control regulates the activity of transposons, retrotransposons and genomic proviruses. Many different life history experiences alter the activities of mobile DNA and the expression of genetic loci regulated by nearby insertions. The same experiences induce alterations in epigenetic formatting and lead to trans-generational modifications of genome expression and stability. These observations lead to the hypothesis that epigenetic formatting directed by non-coding RNA provides a molecular interface between life history events and genome alteration.

Chromosomal polymorphism in two species of Hypancistrus (Siluriformes: Loricariidae): an integrative approach for understanding their biodiversity

Structural chromosome changes are widely described in different vertebrate groups and generate genetic, phenotypic and behavioral diversity. During the evolution of loricariids, several rearrangements (fissions, fusions, inversions) seem to have occurred. Hypancistrus, tribe Ancistrini, are highly demanded for fishkeeping around the world. In this tribe, the diploid chromosome number 2n = 52 is considered a synapomorphy, and paracentric-type inversions appear to be involved in the chromosomal evolution of the tribe. The present study investigated the karyotypes of H. zebra and H. cf. debilittera using cytogenetic, classical and molecular tools, as well as DNA barcoding. Data reveal that, although diploid number in both species corroborates the proposed synapomorphy for the tribe, there is a complex karyotype dynamics, reflected in the intense chromosomal polymorphism, resulting from rearrangements involving ribosomal regions (5S and 18S rDNA), which are suggested to be paracentric inversions. Besides, DNA barcode confirms reciprocal monophyletism between the species, validating the existence of two species, only. This scenario, coupled with genomic instability caused by exogenous sequences such as Rex-3 retrotransposons and the species' sedentary lifestyle, which helps the fast polymorphism fixation, may reflect different phenotypic color patterns in natural populations, as observed in H. cf. debilittera.

Genetic and epigenetic changes involving (retro)transposons in animal hybrids and polyploids

Transposable elements (TEs) are discrete genetic units that have the ability to change their location within chromosomal DNA, and constitute a major and rapidly evolving component of eukaryotic genomes. They can be subdivided into 2 distinct types: retrotransposons, which use an RNA intermediate for transposition, and DNA transposons, which move only as DNA. Rapid advances in genome sequencing significantly improved our understanding of TE roles in genome shaping and restructuring, and studies of transcriptomes and epigenomes shed light on the previously unknown molecular mechanisms underlying genetic and epigenetic TE controls. Knowledge of these control systems may be important for better understanding of reticulate evolution and speciation in the context of bringing different genomes together by hybridization and perturbing the established regulatory balance by ploidy changes. See also sister article focusing on plants by Bento et al. in this themed issue.