Implications of Gene Inheritance Patterns on the Heterosis of Abdominal Fat Deposition in Chickens

Synergy of gut microbiota and host genome in driving heterosis expression of chickens

Heterosis has been widely utilized in agricultural production. Despite over a century of extensive research, the underlying mechanisms of heterosis remain elusive. Most hypotheses and research have focused on the genetic basis of heterosis. However, the potential role of gut microbiota in heterosis has been largely ignored. Here, we carefully design a crossbreeding experiment with two distinct broiler breeds and conduct 16S rRNA amplicon and transcriptome sequencing to investigate the synergistic role of gut microbiota and host genes in driving heterosis. We find that the breast muscle weight of hybrids exhibits a high heterosis, 6.28% higher than the mid-parent value. A notable difference is observed in the composition and potential function of cecal microbiota between hybrids and their parents. Over 90% of differentially colonized microbiota and differentially expressed genes exhibit nonadditive patterns. Integrative analyses uncover associations between nonadditive genes and nonadditive microbiota, including a connection between the expression of cellular signaling pathways and metabolism-related genes and the abundance of Odoribacter, Oscillibacter, and Alistipes in hybrids. Moreover, higher abundances of these microbiota are related to better meat yield. In summary, these findings highlight the importance of gut microbiota in heterosis, serving as crucial factors that modulate heterosis expression in chickens.

Identification of potential tissue-specific biomarkers involved in pig fat deposition through integrated bioinformatics analysis and machine learning

Backfat thickness (BT) and intramuscular fat (IMF) content are closely appertained to meat production and quality in pig production. Deposition in subcutaneous adipose (SA) and IMF concerns different genes and regulatory mechanisms. And larger studies with rigorous design should be carried to explore the molecular regulation of fat deposition in different tissues. The purpose of this study is to gain a better understanding of the molecular mechanisms underlying differences in fat deposition among different tissues and identify tissue-specific genes involved in regulating fat deposition. The SA-associated datasets (GSE122349 and GSE145956) and IMF-associated datasets (GSE165613 and GSE207279) were downloaded from the Gene Expression Omnibus (GEO) as the BT and IMF group, respectively. Subsequently, the Robust Rank Aggregation (RRA) algorithm identified 27 down- and 29 up-regulated differentially expressed genes (DEGs) in the BT group. Based on bioinformatics and three machine learning algorithms, four SA deposition-related potential biomarkers, namely ACLY, FASN, ME1, and ARVCF were selected. FASN was evaluated as the most valuable biomarker for the SA mechanism. The 18 down- and 34 up-regulated DEGs in the IMF group were identified, and ACTA2 and HMGCL were screened as the IMF deposition-related candidate core genes, especially the ACTA2 may play the critical role in IMF deposition regulation. Moreover, based on the constructed ceRNA network, we postulated that the role of predicted ceRNA interaction network of XIST, NEAT1/miR-15a-5p, miR-16-5p, miR-424-5p, miR-497-5p/FASN were vital in the SA metabolism, XIST, NEAT1/miR-27a/b-3p, 181a/c-5p/ACTA2 might contribute to the regulation to IMF metabolism, which all gave suggestions in molecular mechanism for regulation of fat deposition. These findings may facilitate advancements in porcine quality at the genetic and molecular levels and assist with human obesity-associated diseases.

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.

Transcriptome analysis of ovarian tissues highlights genes controlling energy homeostasis and oxidative stress as potential drivers of heterosis for egg number and clutch size in crossbred laying hens

Heterosis is the major benefit of crossbreeding and has been exploited in laying hens breeding for a long time. This genetic phenomenon has been linked to various modes of nonadditive gene action. However, the molecular mechanism of heterosis for egg production in laying hens has not been fully elucidated. To fill this research gap, we sequenced mRNAs and lncRNAs of the ovary stroma containing prehierarchical follicles in White Leghorn, Rhode Island Red chickens as well as their reciprocal crossbreds that demonstrated heterosis for egg number and clutch size. We further delineated the modes of mRNAs and lncRNAs expression to identify their potential functions in the observed heterosis. Results showed that dominance was the principal mode of nonadditive expression exhibited by mRNAs and lncRNAs in the prehierarchical follicles of crossbred hens. Specifically, low-parent dominance was the main mode of mRNA expression, while high-parent dominance was the predominant mode of lncRNA expression. Important pathways enriched by genes that showed higher expression in crossbreds compared to either one or both parental lines were cell adhesion molecules, tyrosine and purine metabolism. In contrast, ECM-receptor interaction, focal adhesion, PPAR signaling, and ferroptosis were enriched in genes with lower expression in the crossbred. Protein network interaction identified nonadditively expressed genes including apolipoprotein B (APOB), transferrin, acyl-CoA synthetase medium-chain family member (APOBEC) 3, APOBEC1 complementation factor, and cathepsin S as hub genes. Among these potential hub genes, APOB was the only gene with underdominance expression common to the 2 reciprocal crossbred lines, and has been linked to oxidative stress. LncRNAs with nonadditive expression in the crossbred hens targeted natriuretic peptide receptor 1, epidermal differentiation protein beta, spermatogenesis-associated gene 22, sperm-associated antigen 16, melanocortin 2 receptor, dolichol kinase, glycine amiinotransferase, and prolactin releasing hormone receptor. In conclusion, genes with nonadditive expression in the crossbred may play crucial roles in follicle growth and atresia by improving follicle competence and increasing oxidative stress, respectively. These 2 phenomena could underpin heterosis for egg production in crossbred laying hens.

The mRNA-lncRNA landscape of multiple tissues uncovers key regulators and molecular pathways that underlie heterosis for feed intake and efficiency in laying chickens

BACKGROUND

Heterosis is routinely exploited to improve animal performance. However, heterosis and its underlying molecular mechanism for feed intake and efficiency have been rarely explored in chickens. Feed efficiency continues to be an important breeding goal trait since feed accounts for 60 to 70% of the total production costs in poultry. Here, we profiled the mRNA-lncRNA landscape of 96 samples of the hypothalamus, liver and duodenum mucosa from White Leghorn (WL), Beijing-You chicken (YY), and their reciprocal crosses (WY and YW) to elucidate the regulatory mechanisms of heterosis.

RESULTS

We observed negative heterosis for both feed intake and residual feed intake (RFI) in YW during the laying period from 43 to 46 weeks of age. Analysis of the global expression pattern showed that non-additivity was a major component of the inheritance of gene expression in the three tissues for YW but not for WY. The YW-specific non-additively expressed genes (YWG) and lncRNA (YWL) dominated the total number of non-additively expressed genes and lncRNA in the hypothalamus and duodenum mucosa. Enrichment analysis of YWG showed that mitochondria components and oxidation phosphorylation (OXPHOS) pathways were shared among the three tissues. The OXPHOS pathway was enriched by target genes for YWL with non-additive inheritance of expression in the liver and duodenum mucosa. Weighted gene co-expression network analysis revealed divergent co-expression modules associated with feed intake and RFI in the three tissues from WL, YW, and YY. Among the negatively related modules, the OXPHOS pathway was enriched by hub genes in the three tissues, which supports the critical role of oxidative phosphorylation. Furthermore, protein quantification of ATP5I was highly consistent with ATP5I expression in the liver, which suggests that, in crossbred YW, non-additive gene expression is down-regulated and decreases ATP production through oxidative phosphorylation, resulting in negative heterosis for feed intake and efficiency.

CONCLUSIONS

Our results demonstrate that non-additively expressed genes and lncRNA involved in oxidative phosphorylation in the hypothalamus, liver, and duodenum mucosa are key regulators of the negative heterosis for feed intake and RFI in layer chickens. These findings should facilitate the rational choice of suitable parents for producing crossbred chickens.

MicroRNAs with non-additive expression in the ovary of hybrid hens target genes enriched in key reproductive pathways that may influence heterosis for egg laying traits

Heterosis has been extensively exploited in chicken breeding to improve laying traits in commercial hybrid stock. However, the molecular mechanisms underlying it remains elusive. This study characterizes the miRNAome in the pre-hierarchical follicles of purebred and hybrid laying hens, and investigate the functions of miRNAs with non-additive expression in the pre-hierarchical follicles as they modulate heterosis for egg number and clutch size. To achieve that aim, White Leghorn and Rhode Island Red chicken lines were reciprocally crossed to generate hybrids. The crossbreds demonstrated heterosis for egg number and clutch size, and pre-hierarchical follicles from 4 birds of each genotype were collected at 53 weeks of age. Mode of miRNA expression was characterized after miRNA sequencing. A total of 50 miRNAs including 30 novel ones, were found to exhibit non-additive expression. Dominance was the predominant mode of expression exhibited by majority of the miRNAs. Functional analysis of target genes of the known miRNAs with non-additive expression revealed Gene Ontology terms related to regulation of transcription, metabolic processes and gene expression. KEGG and REACTOME pathways including hedgehog, cellular senescence, wnt, TGF-β, progesterone-mediated oocyte maturation, oocyte meiosis, GnRH signaling, signal transduction and generic transcription, which can be linked to primordial follicle activation, growth and ovulation, were significantly enriched by target genes of miRNAs with non-additive expression. Majority of the genes enriched in these biological pathways were targeted by gga-miR-19a, gga-miR-19b, gga-miR-375, gga-miR-135a, and gga-miR-7 and 7b, thus, revealing their synergistic roles in enhancing processes that could influence heterosis for egg number and clutch size in hybrid hens.

Genetic Basis of Sexual Maturation Heterosis: Insights From Ovary lncRNA and mRNA Repertoire in Chicken

Sexual maturation is fundamental to the reproduction and production performance, heterosis of which has been widely used in animal crossbreeding. However, the underlying mechanism have long remained elusive, despite its profound biological and agricultural significance. In the current study, the reciprocal crossing between White Leghorns and Beijing You chickens were performed to measure the sexual maturation heterosis, and the ovary lncRNAs and mRNAs of purebreds and crossbreeds were profiled to illustrate molecular mechanism of heterosis. Heterosis larger than 20% was found for pubic space and oviduct length, whereas age at first egg showed negative heterosis in both crossbreeds. We identified 1170 known lncRNAs and 1994 putative lncRNAs in chicken ovary using a stringent pipeline. Gene expression pattern showed that nonadditivity was predominant, and the proportion of nonadditive lncRNAs and genes was similar between two crossbreeds, ranging from 44.24% to 49.15%. A total of 200 lncRNAs and 682 genes were shared by two crossbreeds, respectively. GO and KEGG analysis showed that the common genes were significantly enriched in the cell cycle, animal organ development, gonad development, ECM-receptor interaction, calcium signaling pathway and GnRH signaling pathway. Weighted gene co-expression network analysis (WGCNA) identified that 7 out of 20 co-expressed lncRNA-mRNA modules significantly correlated with oviduct length and pubic space. Interestingly, genes harbored in seven modules were also enriched in the similar biological process and pathways, in which nonadditive lncRNAs, such as MSTRG.17017.1 and MSTRG.6475.20, were strongly associated with nonadditive genes, such as CACNA1C and TGFB1 to affect gonad development and GnRH signaling pathway, respectively. Moreover, the results of real-time quantitative PCR (RT-qPCR) correlated well with the transcriptome data. Integrated with positive heterosis of serum GnRH and melatonin content detected in crossbreeds, we speculated that nonadditive genes involved in the GnRH signaling pathway elevated the gonad development, leading to the sexual maturation heterosis. We characterized a systematic landscape of ovary lncRNAs and mRNAs related to sexual maturation heterosis in chicken. The quantitative exploration of hybrid transcriptome changes lays foundation for genetic improvement of sexual maturation traits and provides insights into endocrine control of sexual maturation.

Transcriptome-Wide Analyses Identify Dominant as the Predominantly Non-Conservative Alternative Splicing Inheritance Patterns in F1 Chickens

Transcriptome analysis has been used to investigate many economically traits in chickens; however, alternative splicing still lacks a systematic method of study that is able to promote proteome diversity, and fine-tune expression dynamics. Hybridization has been widely utilized in chicken breeding due to the resulting heterosis, but the dynamic changes in alternative splicing during this process are significant yet unclear. In this study, we performed a reciprocal crossing experiment involving the White Leghorn and Cornish Game chicken breeds which exhibit major differences in body size and reproductive traits, and conducted RNA sequencing of the brain, muscle, and liver tissues to identify the inheritance patterns. A total of 40 515 and 42 612 events were respectively detected in the brain and muscle tissues, with 39 843 observed in the liver; 2807, 4242, and 4538 events significantly different between two breeds were identified in the brain, muscle, and liver tissues, respectively. The hierarchical cluster of tissues from different tissues from all crosses, based on the alternative splicing profiles, suggests high tissue and strain specificity. Furthermore, a comparison between parental strains and hybrid crosses indicated that over one third of alternative splicing genes showed conserved patterns in all three tissues, while the second prevalent pattern was non-additive, which included both dominant and transgressive patterns; this meant that the dominant pattern plays a more important role than suppression. Our study provides an overview of the inheritance patterns of alternative splicing in layer and broiler chickens, to better understand post-transcriptional regulation during hybridization.

Genetic basis of negative heterosis for growth traits in chickens revealed by genome-wide gene expression pattern analysis

Background

Heterosis is an important biological phenomenon that has been extensively utilized in agricultural breeding. However, negative heterosis is also pervasively observed in nature, which can cause unfavorable impacts on production performance. Compared with systematic studies of positive heterosis, the phenomenon of negative heterosis has been largely ignored in genetic studies and breeding programs, and the genetic mechanism of this phenomenon has not been thoroughly elucidated to date. Here, we used chickens, the most common agricultural animals worldwide, to determine the genetic and molecular mechanisms of negative heterosis.

Results

We performed reciprocal crossing experiments with two distinct chicken lines and found that the body weight presented widely negative heterosis in the early growth of chickens. Negative heterosis of carcass traits was more common than positive heterosis, especially breast muscle mass, which was over − 40% in reciprocal progenies. Genome-wide gene expression pattern analyses of breast muscle tissues revealed that nonadditivity, including dominance and overdominace, was the major gene inheritance pattern. Nonadditive genes, including a substantial number of genes encoding ATPase and NADH dehydrogenase, accounted for more than 68% of differentially expressed genes in reciprocal crosses (4257 of 5587 and 3617 of 5243, respectively). Moreover, nonadditive genes were significantly associated with the biological process of oxidative phosphorylation, which is the major metabolic pathway for energy release and animal growth and development. The detection of ATP content and ATPase activity for purebred and crossbred progenies further confirmed that chickens with lower muscle yield had lower ATP concentrations but higher hydrolysis activity, which supported the important role of oxidative phosphorylation in negative heterosis for growth traits in chickens.

Conclusions

These findings revealed that nonadditive genes and their related oxidative phosphorylation were the major genetic and molecular factors in the negative heterosis of growth in chickens, which would be beneficial to future breeding strategies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-021-00574-2.

Transcriptomics and proteomics-based analysis of heterosis on main economic traits of silkworm, Bombyx mori

The application of silkworm hybrids have promoted the innovation and development of agricultural technology, but the mechanism of heterosis in silkworm has not been explained clearly. In this study, the heterosis of silkworm in the aspects of body weight, silk gland and cocoon weight was investigated by means of silkworm hybridization and multi-omics approaches, including transcriptome and proteome. The results showed that heterosis of silkworm body weight, silk gland and cocoon weight was overdominant, but only part of genes and proteins were overdominant, and most of genes and proteins were non-additive. Combined analysis obtained six up-regulated genes and four down-regulated genes that were consistent both in transcriptome and proteome. Gene functional enrichment analysis indicated that most up-regulated genes and proteins were mostly related to metabolism, which led to accelerated metabolism and protein synthesis and contributing to improved heterosis. The up-regulation of 6-phosphate glucose dehydrogenase (G6PDH), phosphatidylethanolamine-binding protein (PEBP) and sHSP20.4, which are involved in metabolism, might be related to silk gland heterosis. SIGNIFICANCE: A combination of transcriptomic and proteomic analysis was used to understand the molecular mechanism of silkworm heterosis. We found that the phenotypic traits of silkworm are overdominant, while the analysis of transcriptome and proteome showed that only part of genes and proteins were overdominant, and most of genes and proteins were non-additive. Some of the genes had unique expression in F1, which was speculated that genes under heterozygous condition may result in rearrangement and cause metabolic changes in the hybrids. Those both up-regulated in transcriptomic and proteomic analysis were found to be involved in various metabolic processes, so as to accelerate metabolism and protein synthesis, thus exhibiting heterosis.

Inheritance patterns of leukocyte gene expression under heat stress in F1 hybrid cattle and their parents

Crossbreeding capitalizes on heterosis effects and results in increased performance of crossbred animals. Dominance hypothesis and overdominance hypothesis are 2 common models proposed to explain heterosis. Differential gene expression between parents and hybrids is hypothesized to be responsible for heterosis. This study aimed to investigate the heat tolerance and inheritance patterns of leukocyte transcriptomics in F₁ hybrid cattle (Angus males × Droughtmaster females) and their parents Red Angus (AN) and Droughtmaster (DR) under heat stress. According to the respiratory rate and heat tolerance coefficient index, DR was better adapted to heat stress than AN. The physiological responses to heat stress of F₁ hybrids were similar to AN. We identified 802 differentially expressed genes in leukocytes between AN and DR under heat stress using mRNA sequencing. Compared with AN, upregulated genes in DR were enriched in biological processes of response to stress, external and chemical stimulus, and cytokine, cell surface receptor signaling pathway, and cardiovascular system development. In contrast, upregulated genes in AN were enriched in B cell activation and regulation of B cell activation. Gene expression levels can be inherited additively or nonadditively and are classified into additive (35%), dominance (44%), and overdominance and underdominance (18%) modes in F₁ hybrids and their parents. Inheritance patterns of gene expression showed that 97% (249/255) of the dominant genes were classified as paternal AN dominant in hybrids. The paternal imprinted PEG10 gene and its regulatory transcription factor MYC showed an AN dominant expression pattern. The MYC interacted with most AN dominant genes. These transcriptomic analyses revealed that DR and AN had specific cellular and humoral immunity and cardiovascular systems development function under heat stress. Inheritance pattern analyses from gene expression partly explained phenotypic differences between parents and F₁ hybrids. The paternal imprinted PEG10 gene interaction with transcription factor MYC may contribute to explaining paternal dominant gene expression in hybrids.

Hybrids generated by crossing elite laying chickens exhibited heterosis for clutch and egg quality traits

Crossbreeding advantage in hybrids compared with their parents, termed heterosis, has been exhaustively exploited in chicken breeding over the last century. Reports for crossbreeding of elite laying chickens covering rearing and laying period remain infrequent. In this study, resource populations of Rhode Island Red (RIR) and White Leghorn (WL) pure-bred chickens were reciprocally crossed to generate 4 distinct groups that were evaluated for prelaying growth, egg production, and egg quality. Birds monitored for prelaying growth consists of 105 (RIR), 131 (WL), 207 (RIR × WL) and 229 (WL × RIR), and 30 pullets from each group were evaluated. Egg laying records were collected from 102, 89, 147, and 191 hens in the 4 populations, respectively. In addition, expression of 5 candidate genes for egg production in the ovarian follicles was measured by RT-qPCR. Results showed that BW of hatched chicks in the WL line was higher than the other populations. However, the 2 crossbreds grew faster than WL purebred throughout the prelaying period. Low to medium heterosis was observed for BW and body length before the onset of lay. White Leghorn and the hybrids commenced laying earlier than RIR pullets and egg production traits were favorable in the crossbreds compared with purebreds. Heterosis for egg number and clutch size was moderate in WL × RIR but low in RIR × WL hens. Expression of antimullerian hormone gene was high in WL and RIR × WL hybrids, suggesting WL parent-specific enhancing dominant expression. Shell weight was higher in the crossbreds than purebreds at 52 wk of age, but RIR hens laid eggs with higher shell ratio than the other populations (P < 0.05). Conversely, WL and the hybrids had higher eggshell strength than RIR birds (P < 0.05). Eggshell strength was the only egg quality trait that showed heterosis above 10% in WL × RIR hybrids at 32 and 52 wk of age. White Leghorn × RIR hens demonstrated higher percent heterosis for economic traits than birds of the reciprocal hybrid. This means that RIR breed is a better dam than a sire line for growth, egg laying, and egg quality traits.