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

Transcriptomic analysis reveals the molecular mechanisms of heterosis in low-temperature tolerance in the hybrids of Argopecten scallops

The F1 hybrid of Argopecten irradians irradians (♀) × A. purpuratus (♂) exhibits significant heterosis in growth performance and mid-parent heterosis in low-temperature tolerance. This study presents a comparative transcriptomic analysis of A. irradians irradians (Ai), A. purpuratus (Ap), and the hybrid A. irradians irradians♀ × A. purpuratus♂ (Aip) to explore the mechanisms underlying low-temperature tolerance heterosis in Aip. A total of 33,376 differentially expressed genes (DEGs) were identified between F1 hybrids and purebreds under cold stress. In Aip, 80.32 % of DEGs exhibited non-additive expression patterns, with over-dominant expression observed in 30.65 % of these genes. Pairwise comparisons among the transcriptomes of Ai, Ap, and Aip revealed 14,959 alternative splicing events, affecting 8169 genes. KEGG pathway analysis indicated substantial enrichment of overlapping genes from common DEGs and non-additively expressed genes (NAGs) in apoptosis, longevity regulation, ABC transporters, and spliceosome pathways. Furthermore, analysis of DEGs, DAGs (Differentially Alternative Splicing genes), and NAGs identified 6 genes undergoing alternative splicing. These pathways and genes may be crucial in Aip's response to low-temperature stress and offer insights for advancing scallop cross-breeding strategies.

Strand-Specific RNA Sequencing Reveals Gene Expression Patterns in F1 Chick Breast Muscle and Liver after Hatching

Heterosis refers to the phenomenon where hybrids exhibit superior performance compared to the parental phenotypes and has been widely utilized in crossbreeding programs for animals and crops, yet the molecular mechanisms underlying this phenomenon remain enigmatic. A better understanding of the gene expression patterns in post-hatch chickens is very important for exploring the genetic basis underlying economically important traits in the crossbreeding of chickens. In this study, breast muscle and liver tissues (n = 36) from full-sib F1 birds and their parental pure lines were selected to identify gene expression patterns and differentially expressed genes (DEGs) at 28 days of age by strand-specific RNA sequencing (ssRNA-seq). This study indicates that additivity is the predominant gene expression pattern in the F1 chicken post-hatch breast muscle (80.6% genes with additivity) and liver (94.2% genes with additivity). In breast muscle, Gene Ontology (GO) enrichment analysis revealed that a total of 11 biological process (BP) terms closely associated with growth and development were annotated in the identified DEG sets and non-additive gene sets, including STAT5A, TGFB2, FGF1, IGF2, DMA, FGF16, FGF12, STAC3, GSK3A, and GRB2. Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation presented that a total of six growth- and development-related pathways were identified, involving key genes such as SLC27A4, GLUL, TGFB2, COX17, and GSK3A, including the PPAR signaling pathway, TGF-beta signaling pathway, and mTOR signaling pathway. Our results may provide a theoretical basis for crossbreeding in domestic animals.

Nonadditive and allele-specific expression of ghrelin in hybrid tilapia

Background

Interspecies hybridization is an important breeding method to generate fishes with heterosis in aquaculture. Using this method, hybrid Nile tilapia (Oreochromis niloticus, ♀) × blue tilapia (Oreochromis aureus, ♂) has been produced and widely farmed due to its growth and appetite superiorities. However, the genetic mechanism of these advanced traits is still not well understood. Ghrelin is a crucial gene that regulates growth and appetite in fishes. In the present study, we focused on the expression characteristics and its regulation of ghrelin in the hybrid.

Results

The tissue distribution analysis showed that ghrelin was predominantly expressed in the stomach in the hybrid. Ghrelin was more highly expressed in the stomach in the hybrid and Nile tilapia, compared to blue tilapia, showing a nonadditive pattern. Two single-nucleotide polymorphism (SNP) sites were identified including T/C and C/G from the second exon in the ghrelin gene from Nile tilapia and blue tilapia. By pyrosequencing based on the SNP sites, the allele-specific expression (ASE) of ghrelin in the hybrid was assayed. The result indicated that ghrelin in the hybrid showed higher maternal allelic transcript ratios. Fasting significantly increased ghrelin overall expression at 4, 8, 12, 24, and 48 h. In addition, higher maternal allelic transcript ratios were not changed in the fasting hybrids at 48 h. The cis and trans effects were determined by evaluating the overall expression and ASE values in the hybrid. The expression of ghrelin was mediated by compensating cis and trans effects in hybrid.

Conclusion

In summary, the present lines of evidence showed the nonadditive expression of ghrelin in the hybrid tilapia and its regulation by subgenomes, offering new insight into gene expression characteristics in hybrids.

Transcriptome Analysis of Reciprocal Hybrids Between Crassostrea gigas and C. angulata Reveals the Potential Mechanisms Underlying Thermo-Resistant Heterosis

Heterosis, also known as hybrid vigor, is widely used in aquaculture, but the molecular causes for this phenomenon remain obscure. Here, we conducted a transcriptome analysis to unveil the gene expression patterns and molecular bases underlying thermo-resistant heterosis in Crassostrea gigas ♀ × Crassostrea angulata ♂ (GA) and C. angulata ♀ × C. gigas ♂ (AG). About 505 million clean reads were obtained, and 38,210 genes were identified, of which 3779 genes were differentially expressed between the reciprocal hybrids and purebreds. The global gene expression levels were toward the C. gigas genome in the reciprocal hybrids. In GA and AG, 95.69% and 92.00% of the differentially expressed genes (DEGs) exhibited a non-additive expression pattern, respectively. We observed all gene expression modes, including additive, partial dominance, high and low dominance, and under- and over-dominance. Of these, 77.52% and 50.00% of the DEGs exhibited under- or over-dominance in GA and AG, respectively. The over-dominance DEGs common to reciprocal hybrids were significantly enriched in protein folding, protein refolding, and intrinsic apoptotic signaling pathway, while the under-dominance DEGs were significantly enriched in cell cycle. As possible candidate genes for thermo-resistant heterosis, GRP78, major egg antigen, BAG, Hsp70, and Hsp27 were over-dominantly expressed, while MCM6 and ANAPC4 were under-dominantly expressed. This study extends our understanding of the thermo-resistant heterosis in oysters.

Integration of mRNA and miRNA Profiling Reveals Heterosis in Oreochromis niloticus × O. aureus Hybrid Tilapia

Heterosis is a widespread biological phenomenon in fishes, in which hybrids have superior traits to parents. However, the underlying molecular basis for heterosis remains uncertain. Heterosis in growth and survival rates is apparent in hybrid tilapia (Oreochromis niloticus ♀ × O. aureus ♂). Comparisons of growth and hematological biochemical characteristics and mRNA and miRNA transcriptional analyses were performed in hybrid and parents tilapia stocks to investigate the underlying molecular basis for heterosis. Growth characteristics and hematological glucose and cholesterol parameters were significantly improved in hybrids. Of 3097 differentially expressed genes (DEGs) and 120 differentially expressed miRNAs (DEMs) identified among three stocks (O. niloticus, O. aureus, and hybrids), 1598 DEGs and 62 DEMs were non-additively expressed in hybrids. Both expression level dominance and overdominance patterns occurred for DEGs and DEMs, indicating that dominance and overdominance models are widespread in the transcriptional and post-transcriptional regulation of genes involved in growth, metabolism, immunity, and antioxidant capacity in hybrid tilapia. Moreover, potential negative regulation networks between DEMs and predicted target DEGs revealed that most DEGs from miRNA-mRNA pairs are up-regulated. Dominance and overdominance models in levels of transcriptome and miRNAome facilitate the integration of advantageous parental alleles into hybrids, contributing to heterosis of growth and improved survival. The present study provides new insights into molecular heterosis in hybrid tilapia, advancing our understanding of the complex mechanisms involved in this phenomenon in aquatic animals.

Development of a relevant strategy using de novo transcriptome assembly method for transcriptome comparisons between Muscovy and common duck species and their reciprocal inter-specific mule and hinny hybrids fed ad libitum and overfed

Background

Common Pekin and Muscovy ducks and their intergeneric hinny and mule hybrids have different abilities for fatty liver production. RNA-Seq analyses from the liver of these different genetic types fed ad libitum or overfed would help to identify genes with different response to overfeeding between them. However RNA-seq analyses from different species and comparison is challenging. The goal of this study was develop a relevant strategy for transcriptome analysis and comparison between different species.

Results

Transcriptomes were first assembled with a reference-based approach. Important mapping biases were observed when heterologous mapping were conducted on common duck reference genome, suggesting that this reference-based strategy was not suited to compare the four different genetic types. De novo transcriptome assemblies were then performed using Trinity and Oases. Assemblies of transcriptomes were not relevant when more than a single genetic type was considered. Finally, single genetic type transcriptomes were assembled with DRAP in a mega-transcriptome. No bias was observed when reads from the different genetic types were mapped on this mega-transcriptome and differences in gene expression between the four genetic types could be identified.

Conclusions

Analyses using both reference-based and de novo transcriptome assemblies point out a good performance of the de novo approach for the analysis of gene expression in different species. It also allowed the identification of differences in responses to overfeeding between Pekin and Muscovy ducks and hinny and mule hybrids.

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.

The subgenomes show asymmetric expression of alleles in hybrid lineages of Megalobrama amblycephala × Culter alburnus

Hybridization drives rapid speciation by shaping novel genotypic and phenotypic profiles. Genomic incompatibility and transcriptome shock have been observed in hybrids, although this is rarer in animals than in plants. Using the newly sequenced genomes of the blunt snout bream (Megalobrama amblycephala [BSB]) and the topmouth culter (Culter alburnus [TC]), we focused on the sequence variation and gene expression changes in the reciprocal intergeneric hybrid lineages (F₁-F₃) of BSB × TC. A genome-wide transcriptional analysis identified 145-974 expressed recombinant genes in the successive generations of hybrid fish, suggesting the rapid emergence of allelic variation following hybridization. Some gradual changes of gene expression with additive and dominance effects and various cis and trans regulations were observed from F₁ to F₃ in the two hybrid lineages. These asymmetric patterns of gene expression represent the alternative strategies for counteracting deleterious effects of the subgenomes and improving adaptability of novel hybrids. Furthermore, we identified positive selection and additive expression patterns in transforming growth factor, beta 1b (tgfb1b), which may account for the morphological variations of the pharyngeal jaw in the two hybrid lineages. Our current findings provide insights into the evolution of vertebrate genomes immediately following hybridization.

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.

Genetic Effects on the Gut Microbiota Assemblages of Hybrid Fish From Parents With Different Feeding Habits

Gut microbiota play critical roles in host nutrition and metabolism. However, little is known about the genetic effects on the gut microbiota assemblages because a suitable model for investigation is lacking. In the present study, we established the reciprocal hybrid fish lineages derived from the parents with different feeding habits, namely, herbivorous blunt snout bream (Megalobrama amblycephala, BSB, 2n = 48) and carnivorous topmouth culter (Culter alburnus, TC, 2n = 48). We investigated the genetic effects on gut microbiota assemblages by using 16S rRNA gene sequencing. The results showed that the gut characteristics (structure, relative gut length, relative gut mass, and Zihler’s index) differed between the two types of hybrids and the two parents. In particular, a strong correlation between genotype and gut microbial assemblages indicated that host genetic (subgenome) significantly altered the gut microbial communities. In addition, the microbial structures (composition and abundance) in the two types of hybrids were more similar to those in BSB parent (P > 0.05) than to those in TC parent (P < 0.05), and the cellulase contents in the gut (produced by gut microbes) also showed the similar results. The results suggested that the host genomic interaction (mainly subgenome domination) had a sizeable effect on shaping the gut microbiota assemblages in reciprocal hybrid fish. This study enriches our understanding of the relationship between host genetic and gut microbiota assemblages, and provides insight into gut microbiota and metabonomics.

Differentiation in stem/progenitor cells along fetal or adult hepatic stages requires transcriptional regulators independently of oscillations in microRNA expression

Understanding mechanisms in lineage differentiation is critical for organ development, pathophysiology and oncogenesis. To determine whether microRNAs (miRNA) may serve as drivers or adjuncts in hepatic differentiation, we studied human embryonic stem cell-derived hepatocytes and primary hepatocytes representing fetal or adult stages. Model systems were used for hepatic lineage advancement or regression under culture conditions with molecular assays. Profiles of miRNA in primary fetal and adult hepatocytes shared similarities and distinctions from pluripotent stem cells or stem cell-derived early fetal-like hepatocytes. During phenotypic regression in fetal or adult hepatocytes, miRNA profiles oscillated to regain stemness-associated features that had not been extinguished in stem cell-derived fetal-like hepatocytes. These oscillations in stemness-associated features were not altered in fetal-like hepatocytes by inhibitory mimics for dominantly-expressed miRNA, such as hsa-miR-99b, -100, -214 and -221/222. The stem cell-derived fetal-like hepatocytes were permissive for miRNA characterizing mature hepatocytes, including mimics for hsa-miR-122, -126, -192, -194 and -26b, although transfections of the latter did not advance hepatic differentiation. Examination of genome-wide mRNA expression profiles in stem cell-derived or primary fetal hepatocytes indicated targets of highly abundant miRNA regulated general processes, e.g., cell survival, growth and proliferation, functional maintenance, etc., without directing cell differentiation. Among upstream regulators of gene networks in stem cell-derived hepatocytes included HNF4A, SNAI1, and others, which affect transcriptional circuits directing lineage development or maintenance. Therefore, miRNA expression oscillated in response to microenvironmental conditions, whereas lineage-specific transcriptional regulators, such as HNF4A, were necessary for directing hepatic differentiation. This knowledge will be helpful for understanding the contribution of stem cells in pathophysiological states and oncogenesis, as well as for applications of stem cell-derived hepatocytes.

Analyzing homoeolog expression provides insights into the rediploidization event in gynogenetic hybrids of Carassius auratus red var. × Cyprinus carpio

Rediploidization is considered to be a part of the evolutionary history of allotetraploids, and resulted in the emergence of novel epigenetic regulatory activities. To study the changing patterns of gene expression following the reduction of a genome by 50%, we used RNA-seq and quantitative real-time PCR (qPCR) to investigate total gene expression and homoeolog expression in three hybrids of a C. auratus red var. (2n = 100, ♀) (R) and C. carpio (2n = 100, ♂) (C) (i.e., F1, F18, and G4) and their original parents. A comparison of homoeolog expression between G4 and F18 identified 7 genes (0.22%) that exhibited novel R/C homoeolog expression patterns in G4, while 4 genes (0.12%) were affected by R/C homoeolog silencing. We determined the direction and extent of the homoeolog expression bias (HEB). The C-HEB genes (i.e., nrp1a and igf1) and R-HEB genes (i.e., fgf23 and esm1) provided insights into the effects of the dominance of one parental homoeolog expression on growth regulation. This dominance may contribute to the rapid growth of G4 fish. Our findings may be relevant for clarifying the relationship between growth heterosis and differences in homoeolog expression patterns.

Homoeologue expression insights into the basis of growth heterosis at the intersection of ploidy and hybridity in Cyprinidae

Hybridization and polyploidization are considered important driving forces that form new epigenetic regulations. To study the changing patterns of expression accompanying hybridization and polyploidization, we used RNA-seq and qRT-PCR to investigate global expression and homoeologue expression in diploid and tetraploid hybrids of Carassius auratus red var. (♀) (R) and Cyprinus carpio (♂) (C). By comparing the relative expression levels between the hybrids and their parents, we defined the expression level dominance (ELD) and homoeologue expression bias (HEB) in liver tissue. The results showed that polyploidization contributed to the conversion of homoeologue ELD. In addition, hybridization had more effect on the change in HEB than polyploidization, while polyploidization had more effect on the change of global gene expression than hybridization. Meanwhile, similar expression patterns were found in growth-related genes. The results suggested that hybridization and polyploidization result in differential degrees of maternal HEB in three tissues (liver, muscle and ovary) tested. The results of this study will increase our understanding of the underlying regulation mechanism of rapid growth in diploid hybrids and allotetraploids. The differential degrees of global expression and homoeologue expression contribute to growth heterosis in newly formed hybrids, ensuring the on-going success of allotetraploid speciation.