Determination of dosage compensation and comparison of gene expression in a triploid hybrid fish

Identification and effective regulation of scarb1 gene involved in pigmentation change in autotetraploid Carassius auratus

The autotetraploid Carassius auratus (4nRR, 4 n=200, RRRR) is derived from whole-genome duplication of Carassius auratus red var. (RCC, 2 n=100, RR). In the current study, we demonstrated that chromatophores and pigment changes directly caused the coloration and variation of 4nRR skin (red in RCC, brownish-yellow in 4nRR). To further explore the molecular mechanisms underlying coloration formation and variation in 4nRR, we performed transcriptome profiling and molecular functional verification in RCC and 4nRR. Results revealed that scarb1, associated with carotenoid metabolism, underwent significant down-regulation in 4nRR. Efficient editing of this candidate pigment gene provided clear evidence of its significant role in RCC coloration. Subsequently, we identified four divergent scarb1 homeologs in 4nRR: two original scarb1 homeologs from RCC and two duplicated ones. Notably, three of these homeologs possessed two highly conserved alleles, exhibiting biased and allele-specific expression in the skin. Remarkably, after precise editing of both the original and duplicated scarb1 homeologs and/or alleles, 4nRR individuals, whether singly or multiply mutated, displayed a transition from brownish-yellow skin to a cyan-gray phenotype. Concurrently, the proportional areas of the cyan-gray regions displayed a gene-dose correlation. These findings illustrate the subfunctionalization of duplicated scarb1, with all scarb1 genes synergistically and equally contributing to the pigmentation of 4nRR. This is the first report concerning the functional differentiation of duplicated homeologs in an autopolyploid fish, substantially enriching our understanding of coloration formation and change within this group of organisms.

Transcriptome profiling explores the immune defence mechanism of triploid Pacific oyster (Crassostrea gigas) blood against Vibrio alginolyticus based on protein interaction networks

Triploid oysters have provided the oyster industry with many benefits, such as fast growth rates, meat quality improvement, and increased oyster production and economic benefits, since the first report on triploid oysters was published. The development of polyploid technology has remarkably increased the output of triploid oysters to meet the increasing demand of consumers for Crassostrea gigas in the past decades. At present, research on triploid oyster has mainly focused on breeding and growth, but studies on the immunity of triploid oysters are limited. According to recent reports, Vibrio alginolyticus is a highly virulent strain that can cause disease and death in shellfish, shrimp, as well as serious economic losses. V. alginolyticus may be a reason why oysters die during summer. Therefore, using V. alginolyticus to explore the resistance and immune defense mechanisms of triploid oysters against pathogens presents practical significance. Transcriptome analysis of gene expression was performed in triploid C. gigas at 12 and 48 h after infection with V. alginolyticus, and the respective 2257 and 191 differentially expressed genes (DEGs) were identified. The results of GO and KEGG enrichment analyses showed that multiple significantly enriched GO terms and KEGG signaling pathways are associated with immunity. A protein-protein interaction network was constructed to investigate the interaction relationship of immune-related genes. Finally, we verified the expression situation of 16 key genes using quantitative RT-PCR. This study is the first to use the PPI network in exploring the immune defense mechanism of triploid C. gigas blood to fill the gap in the immune mechanism of triploid oysters and other mollusks, and provide valuable reference for future triploid farming and pathogen prevention and control.

Proteomics-based molecular and functional characteristic profiling of muscle tissue in Triploid crucian carp

Triploid crucian carp (TCC) is a kind of artificially bred fish with huge economic value to China. It has several excellent characteristics, such as fast growth, strong disease resistance and delicious taste. However, as a regionally specific fish, the underlying molecular mechanisms of these characteristics are largely unknown. In this study, we performed quantitative proteomics on the muscle tissues of TCC and its parents, allotetraploid (♂), red crucian carp (♀) and common carp. Combined with multiple bioinformatic analysis, we found that the taste of TCC can be mainly attributed to umami amino acid-enriched proteins such as PURBA, PVALBI and ATP5F1B, and that its rapid growth can be mainly ascribed to the high expression and regulation of metabolism-related proteins such as NDUFS1, ENO1A and CS. These play significant roles in substrate and energy metabolism, as well as in bias transformation. Subsequently, we identified several proteins, including MDH1AA, GOT1 and DLAT, that may serve as potential regulators of innate immunity by regulating the biosynthesis and transformation of significant antibiotics and antimicrobial peptides. In conclusion, this study can serve as a significant reference for similar investigations and shed light on the molecular and biological functions of individual proteins in TCC muscle tissue.

Patterns of gene expression in ovaries of sexual vs. asexual lineages of a freshwater snail

Why sexual reproduction is so common when asexual reproduction should be much more efficient and less costly remains an open question in evolutionary biology. Comparisons between otherwise similar sexual and asexual taxa allow us to characterize the genetic architecture underlying asexuality, which can, in turn, illuminate how this reproductive mode transition occurred and the mechanisms by which it is maintained or disrupted. Here, we used transcriptome sequencing to compare patterns of ovarian gene expression between actively reproducing obligately sexual and obligately asexual females from multiple lineages of Potamopyrgus antipodarum, a freshwater New Zealand snail characterized by frequent separate transitions to asexuality and coexistence of otherwise similar sexual and asexual lineages. We also used these sequence data to evaluate whether population history accounts for variation in patterns of gene expression. We found that source population was a major source of gene expression variation, and likely more influential than reproductive mode. This outcome for these common garden-raised snails is strikingly similar to earlier results from field-collected snails. While we did not identify a likely set of candidate genes from expression profiles that could plausibly explain how transitions to asexuality occurred, we identified around 1,000 genes with evidence of differential expression between sexual and asexual reproductive modes, and 21 genes that appear to exhibit consistent expression differences between sexuals and asexuals across genetic backgrounds. This second smaller set of genes provides a good starting point for further exploration regarding a potential role in the transition to asexual reproduction. These results mark the first effort to characterize the causes of asexuality in P. antipodarum, demonstrate the apparently high heritability of gene expression patterns in this species, and hint that for P. antipodarum, transitions to asexuality might not necessarily be strongly associated with broad changes in gene expression.

Heterosis of growth trait regulated by DNA methylation and miRNA in allotriploid fish

Background

Heterosis of growth traits in allotriploid fish has benefited the production of aquaculture for many years, yet its genetic and molecular basis has remained obscure. Now, an allotriploid complex, including two triploids and their diploid inbred parents, has provided an excellent model for investigating the potential regulatory mechanisms of heterosis.

Results

Here, we performed a series of analyses on DNA methylation modification and miRNA expression in combination with gene expression in the allotriploid complex. We first established a model of cis- and trans-regulation related to DNA methylation and miRNA in allotriploids. Then, comparative analyses showed that DNA methylation contributed to the emergence of a dosage compensation effect, which reduced gene expression levels in the triploid to the diploid state. We detected 31 genes regulated by DNA methylation in the subgenomes of the allotriploids. Finally, the patterns of coevolution between small RNAs and their homoeologous targets were classified and used to predict the regulation of miRNA expression in the allotriploids.

Conclusions

Our results uncovered the regulatory network between DNA methylation and miRNAs in allotriploids, which not only helps us understand the regulatory mechanisms of heterosis of growth traits but also benefits the study and application of epigenetics in aquaculture.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-022-00455-6.

Production of donor-derived eggs after ovarian germ cell transplantation into the gonads of adult, germ cell-less, triploid hybrid fish†

In animals, spermatogonial transplantation in sterile adult males is widely developed; however, despite its utility, ovarian germ cell transplantation is not well developed. We previously showed that the interspecific hybrid offspring of sciaenid was a suitable model for germ cell transplantation studies as they have germ cell-less gonads. However, all these gonads have testis-like characteristics. Here, we tested whether triploidization in hybrid embryos could result in germ cell-less ovary development. Gonadal structure dimorphism and sex-specific gene expression patterns were examined in 6-month-old triploid hybrids (3nHybs). Thirty-one percent of 3nHybs had germ cell-less gonads with an ovarian cavity. cyp19a1a and foxl2, ovarian differentiation-related genes, were expressed in these gonads, whereas dmrt1 and vasa were not expressed, suggesting ovary-like germ cell-less gonad development. Some (26%) 3nHybs had testis-like germ cell-less gonads. Ovarian germ cells collected from homozygous green fluorescent protein (GFP) transgenic blue drum (BD) (Nibea mitsukurii) were transplanted into 6-month-old 3nHybs gonads via the urogenital papilla or oviduct. After 9 months, the recipients were crossed with wild type BD. Among the six 3nHyb recipients that survived, one female and one male produced fertile eggs and motile sperm carrying gfp-specific DNA sequences. Progeny tests revealed that all F1 offspring possessed gfp-specific DNA sequences, suggesting that these recipients produced only donor-derived eggs or sperm. Histological observation confirmed donor-derived gametogenesis in the 3nHyb recipients' gonads. Overall, triploidization reduces male-biased sex differentiation in germ cell-less gonads. We report, for the first time, donor-derived egg production in an animal via direct ovarian germ cell transplantation into a germ cell-less ovary.

Evolutionary Dynamics and Consequences of Parthenogenesis in Vertebrates

Parthenogenesis is asexual reproduction without any required participation from males and, as such, is a null model for sexual reproduction. In a comparative context, we can expand our understanding of the evolution and ecology of sex by investigating the consequences of parthenogenesis. In this review, we examine the theoretical predictions of and empirical results on the evolution of asexual reproduction in vertebrates, focusing on recent studies addressing the origins and geographic spread of parthenogenetic lineages and the genomic consequences of an asexual life history. With advances in computational methods and genome technologies, researchers are poised to make rapid and significant progress in studying the origin and evolution of parthenogenesis in vertebrates, thus providing an important perspective on understanding biodiversity patterns of both asexual and sexual populations.

Transcriptome sequencing and histology reveal dosage compensation in the liver of triploid pre-smolt Atlantic salmon

Triploid Atlantic salmon (Salmo salar L.) is seen as one of the best solutions to solve key issues in the salmon farming industry, such as the impact of escapees on wild stocks and pre-harvest sexual maturation. However, the effects of triploidy on salmon smoltification are poorly understood at the molecular level, even though smoltification is a very sensitive period that has a major influence on survival rate and performance of farmed salmon. In this study, we have compared the liver transcriptomes of diploid and triploid Atlantic salmon at three ontogeny stages: fry, parr and smolt. In diploid fish, a total of 2,655 genes were differentially expressed between fry and parr, whereas 506 genes had significantly different transcript levels between parr and smolts. In triploids, 1,507 and 974 genes were differentially expressed between fry and parr, and between parr and smolts, respectively. Most of these genes were down-regulated and 34 genes were differentially expressed between ploidies at the same stage. In both ploidy groups, the top differentially expressed genes with ontogeny stage belonged to common functional categories that can be related to smoltification. Nucleotide and energy metabolism were significantly down-regulated in fry when compared to parr, while immune system processes were significantly down-regulated in parr when compared to smolts. The close resemblance of enriched biological processes and pathways between ploidy groups suggests that triploidy is regulated by genome dosage compensation in Atlantic salmon. Histological analysis revealed that areas of vacuolization (steatosis) were present only in fry and parr stages, in contrast to a compact cellular histology with glycogen granules after smoltification. There was no significant difference in vacuolization between ploidy groups at the fry stage but the liver of diploid parr had a 33.5% higher vacuolization area compared to their triploid counterparts. Taken together, our data provide novel insights into the changes that occur at the molecular and histological level in the liver of both diploid and triploid Atlantic salmon prior to and during smoltification.

The Legacy of Sexual Ancestors in Phenotypic Variability, Gene Expression, and Homoeolog Regulation of Asexual Hybrids and Polyploids

Hybridization and polyploidization are important evolutionary processes whose impacts range from the alteration of gene expression and phenotypic variation to the triggering of asexual reproduction. We investigated fishes of the Cobitis taenia-elongatoides hybrid complex, which allowed us to disentangle the direct effects of both processes, due to the co-occurrence of parental species with their diploid and triploid hybrids. Employing morphological, ecological, and RNAseq approaches, we investigated the molecular determinants of hybrid and polyploid forms. In contrast with other studies, hybridization and polyploidy induced relatively very little transgressivity. Instead, Cobitis hybrids appeared intermediate with a clear effect of genomic dosing when triploids expressed higher similarity to the parent contributing two genome sets. This dosage effect was symmetric in the germline (oocyte gene expression), interestingly though, we observed an overall bias toward C. taenia in somatic tissues and traits. At the level of individual genes, expression-level dominance vastly prevailed over additivity or transgressivity. Also, trans-regulation of gene expression was less efficient in diploid hybrids than in triploids, where the expression modulation of homoeologs derived from the "haploid" parent was stronger than those derived from the "diploid" parent. Our findings suggest that the apparent intermediacy of hybrid phenotypes results from the combination of individual genes with dominant expression rather than from simple additivity. The efficiency of cross-talk between trans-regulatory elements further appears dosage dependent. Important effects of polyploidization may thus stem from changes in relative concentrations of trans-regulatory elements and their binding sites between hybridizing genomes. Links between gene regulation and asexuality are discussed.

Combined effects of dosage compensation and incomplete dominance on gene expression in triploid cyprinids

Hybridization and polyploidy are pervasive evolutionary features of flowering plants and frequent among some animal groups, such as fish. These processes always lead to novel genotypes and various phenotypes, including growth heterosis. However, its genetic basis in lower vertebrate is still poorly understood. Here, we conducted transcriptome-level analyses of the allopolyploid complex of Carassius auratus red var. (R) (♀) × Cyprinus carpio L. (C) (♂), including the allodiploid and allotetraploid with symmetric subgenomes, and the two allotriploids with asymmetric subgenomes. The gradual changes of gene silencing and novel gene expression suggested the weakening of the constraint of polymorphic expression in genotypic changes. Then, analyses of the direction and magnitude of homoeolog expression exhibited various asymmetric expression patterns, which supported that R incomplete dominance and dosage compensation were co-regulated in the two triploids. Under these effects, various magnitudes of R-homoeolog expression bias were observed in growth-regulated genes, suggesting that they might contribute to growth heterosis in the two triploids. The determination of R incomplete dominance and dosage compensation, which might be led by asymmetric subgenomes and multiple sets of homologous chromosomes, explained why various expression patterns were shaped and their potential contribution to growth heterosis in the two triploids.

Optimization Techniques to Deeply Mine the Transcriptomic Profile of the Sub-Genomes in Hybrid Fish Lineage

It has been shown that reciprocal cross allodiploid lineage with sub-genomes derived from the cross of Megalobrama amblycephala (BSB) × Culter alburnus (TC) generates the variations in phenotypes and genotypes, but it is still a challenge to deeply mine biological information in the transcriptomic profile of this lineage owing to its genomic complexity and lack of efficient data mining methods. In this paper, we establish an optimization model by non-negative matrix factorization approach for deeply mining the transcriptomic profile of the sub-genomes in hybrid fish lineage. A new so-called spectral conjugate gradient algorithm is developed to solve a sequence of large-scale subproblems such that the original complicated model can be efficiently solved. It is shown that the proposed method can provide a satisfactory result of taxonomy for the hybrid fish lineage such that their genetic characteristics are revealed, even for the samples with larger detection errors. Particularly, highly expressed shared genes are found for each class of the fish. The hybrid progeny of TC and BSB displays significant hybrid characteristics. The third generation of TC-BSB hybrid progeny (BTF3 and TBF3) shows larger trait separation.

Musa balbisiana genome reveals subgenome evolution and functional divergence

Banana cultivars (Musa ssp.) are diploid, triploid and tetraploid hybrids derived from Musa acuminata and Musa balbisiana. We presented a high-quality draft genome assembly of M. balbisiana with 430 Mb (87%) assembled into 11 chromosomes. We identified that the recent divergence of M. acuminata (A-genome) and M. balbisiana (B-genome) occurred after lineage-specific whole-genome duplication, and that the B-genome may be more sensitive to the fractionation process compared to the A-genome. Homoeologous exchanges occurred frequently between A- and B-subgenomes in allopolyploids. Genomic variation within progenitors resulted in functional divergence of subgenomes. Global homoeologue expression dominance occurred between subgenomes of the allotriploid. Gene families related to ethylene biosynthesis and starch metabolism exhibited significant expansion at the pathway level and wide homoeologue expression dominance in the B-subgenome of the allotriploid. The independent origin of 1-aminocyclopropane-1-carboxylic acid oxidase (ACO) homoeologue gene pairs and tandem duplication-driven expansion of ACO genes in the B-subgenome contributed to rapid and major ethylene production post-harvest in allotriploid banana fruits. The findings of this study provide greater context for understanding fruit biology, and aid the development of tools for breeding optimal banana cultivars.

Genome and transcriptome evolve separately in recently hybridized Trichosporon fungi

Genome hybridization is an important evolutionary event that gives rise to species with novel capabilities. However, the merging of distinct genomes also brings together incompatible regulatory networks that must be resolved during the course of evolution. Understanding of the early stages of post-hybridization evolution is particularly important because changes in these stages have long-term evolutionary consequences. Here, via comparative transcriptomic analyses of two closely related, recently hybridized Trichosporon fungi, T. coremiiforme and T. ovoides, and three extant relatives, we show that early post-hybridization evolutionary processes occur separately at the gene sequence and gene expression levels but together contribute to the stabilization of hybrid genome and transcriptome. Our findings also highlight lineage-specific consequences of genome hybridization, revealing that the transcriptional regulatory dynamics in these hybrids responded completely differently to gene loss events: one involving both subgenomes and another that is strictly subgenome-specific.

Nonadditive and allele-specific expression of insulin-like growth factor 1 in Nile tilapia (Oreochromis niloticus, ♀) × blue tilapia (O. aureus, ♂) hybrids

Hybrid Nile tilapia (Oreochromis niloticus, ♀) × blue tilapia (O. aureus, ♂) is a widely cultured tilapia variety due to its growth vigor compared to the parent species. As a peptide hormone, insulin-like growth factor 1 (IGF-1) plays a critical role in regulating somatic growth. The present study focuses on the expression characteristics of IGF-1 in hybrid tilapia. The cloned complete open reading frame of IGF-1 in hybrid tilapia is 549 bp in length, encoding a protein of 182 amino acids. The deduced protein is highly similar to that of Nile tilapia and blue tilapia. IGF-1 was found to be primarily expressed in the liver and muscle in the hybrid; lower expression levels were found in other tissues such as the intestine, spleen, and head-kidney. Increased mRNA expression was observed in the liver and muscle of the hybrid compared to Nile tilapia and blue tilapia, indicating a nonadditive expression pattern in the hybrid. An IGF-1 SNP site (397 site: C in Nile tilapia, G in blue tilapia) for differentiating the Nile tilapia or blue tilapia subgenome in hybrids was identified. Pyrosequencing analysis of the liver transcriptome indicated that most of the hybrids (9 of 10 individuals) predominantly expressed the G allele, demonstrating bias of the blue tilapia subgenome. The present study provides novel data indicating, for the first time, overall gene expression of IGF-1 and allele-specific expression in hybrid tilapia.

Comparative transcriptome analysis of diploid and triploid hybrid groupers (Epinephelus coioides♀ × E. lanceolatus♂) reveals the mechanism of abnormal gonadal development in triploid hybrids

In our previous studies, diploid and triploid hybrids have been detected from the hybridization of Epinephelus coioides♀ × E. lanceolatus♂. The triploid groupers have been found to be delayed in gonadal development, but the mechanism remains poorly understood. In this study, we examined the gonadal development, assayed the serum steroid hormone levels, and compared the BP (brain and pituitary) and G (gonad) transcriptomes of 18-month-old diploid and triploid hybrids. The results showed that levels of serum estradiol-17β and testosterone were significantly higher in triploid groupers. The RNA-seq data revealed that 1518 and 14,963 differentially expressed genes were identified in the BP and G transcriptome, respectively. Further analysis revealed that the expression levels of genes involved in the sexual differentiation pathway and sex steroid synthesis pathway are significantly higher in triploid hybrids. Our findings provided a comprehensive insight into a better understanding of the regulatory mechanisms of sterility in triploid hybrid fish.

Comparative ploidy response to experimental hydrogen peroxide exposure in Atlantic salmon (Salmo salar)

While research into the growth, survival, nutrition and, more recently, disease susceptibility of triploid Atlantic salmon has expanded, there remains an overall lack of studies assessing the response of triploids to chemical treatments. It is essential that the response of triploids to disease treatments be characterised to validate their suitability for commercial production. This study aimed to investigate and compare the stress and immune responses of triploid and diploid Atlantic salmon following an experimental treatment with hydrogen peroxide (H₂O₂). A dose response test was first undertaken to determine a suitable test dose for both diploid and triploid Atlantic salmon. Following this, diploids and triploids were exposed to H₂O₂ (1800 ppm) for 20 min, as per commercial practices, after which blood glucose and lactate, and plasma cortisol and lysozyme were measured, along with the expression of oxidative stress and immune-related genes. In the first 6 h post-exposure to H₂O₂, comparable mortalities occurred in both diploid and triploid Atlantic salmon. Cortisol, glucose and lactate were not significantly influenced by ploidy suggesting that, physiologically, triploid Atlantic salmon are able to cope with the stress associated with H₂O₂ exposure as well as their diploid counterparts. Exposure to H₂O₂ significantly elevated the expression of cat and sod2 in diploid livers and gr, il1β and crp/sap1b in diploid gills, while it significantly decreased the expression of saa5 and crp/sap1a in diploid gills. In triploids, the expression levels of cat, hsp70, sod1, saa5, crp/sap1a and crp/sap1b in liver was significantly higher in fish exposed to H₂O₂ compared to control fish. The expression of gr, sod1 and il1β in triploid gills was also elevated in response to H₂O₂ exposure. This study represents the first experimental evidence of the effects of H₂O₂ exposure on triploid Atlantic salmon and continues to support their application into commercial production.

Asymmetric expression patterns reveal a strong maternal effect and dosage compensation in polyploid hybrid fish

Background

Hybridization and polyploidization are regarded as the major driving forces in plant speciation, diversification, and ecological adaptation. Our knowledge regarding the mechanisms of duplicated-gene regulation following genomic merging or doubling is primarily derived from plants and is sparse for vertebrates.

Results

We successfully obtained an F1 generation (including allodiploid hybrids and triploid hybrids) from female Megalobrama amblycephala Yih (BSB, 2n = 48) × male Xenocypri davidi Bleeker (YB, 2n = 48). The duplicated-gene expression patterns of the two types of hybrids were explored using RNA-Seq data. In total, 5.44 × 10⁸ (69.32 GB) clean reads and 499,631 assembled unigenes were obtained from the testis transcriptomes. The sequence similarity analysis of 4265 orthologs revealed that the merged genomes were dominantly expressed in different ploidy hybrids. The differentially expressed genes in the two types of hybrids were asymmetric compared with those in both parents. Furthermore, the genome-wide expression level dominance (ELD) was biased toward the maternal BSB genome in both the allodiploid and triploid hybrids. In addition, the dosage-compensation mechanisms that reduced the triploid expression levels to the diploid state were determined in the triploid hybrids.

Conclusions

Our results indicate that divergent genomes undergo strong interactions and domination in allopolyploid offspring. Genomic merger has a greater effect on the gene-expression patterns than genomic doubling. The various expression mechanisms (including maternal effect and dosage compensation) in different ploidy hybrids suggest that the initial genomic merger and doubling play important roles in polyploidy adaptation and evolution.

Electronic supplementary material

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

Dosage compensation in the process of inactivation/reactivation during both germ cell development and early embryogenesis in mouse

Ohno proposed that dosage compensation in mammals evolved as a two-step mechanism involving X-inactivation and X-upregulation. While X-inactivation is well characterized, it remains to further analysis whether upregulation of the single activated X chromosome in mammals occurs. We obtained RNA-seq data, including single-cell RNA-seq data, from cells undergoing inactivation/reactivation in both germ cell development and early embryogenesis stages in mouse and calculated the X: A ratio from the gene expression. Our results showed that the X: A ratio is always 1, regardless of the number of X chromosomes being transcribed for expressed genes. Furthermore, the single-cell RNA-seq data across individual cells of mouse preimplantation embryos of mixed backgrounds indicated that strain-specific SNPs could be used to distinguish transcription from maternal and paternal chromosomes and further showed that when the paternal was inactivated, the average gene dosage of the active maternal X chromosome was increased to restore the balance between the X chromosome and autosomes. In conclusion, our analysis of RNA-seq data (particularly single-cell RNA-seq) from cells undergoing the process of inactivation/reactivation provides direct evidence that the average gene dosage of the single active X chromosome is upregulated to achieve a similar level to that of two active X chromosomes and autosomes present in two copies.