Whole transcriptome RNA Sequencing Reveals the Global Molecular Responses and circRNA/lncRNA-miRNA-mRNA ceRNA Regulatory Network in Chicken Fat Deposition

CircRNAs involved in the red light of effect on follicle selection in pigeons

Red light (RL) can enhance egg production in poultry. CircRNAs play a crucial role by serving as transcriptional regulators. However, their role in influencing follicle development in White King pigeons remains unexplored. In this study, 54 paired White King pigeons were chosen and divided into RL and white light (WL) groups, each with 3 subgroups. The egg production of paired pigeons in each replicate was recorded for 45 d, and the characteristics of follicle development were monitored during the laying interval (LI). The granulosa cell layer from follicles of the second-largest follicle (F2) was collected, and high-throughput sequencing was performed to elucidate the molecular mechanism of follicle development in pigeons. The study confirmed that RL enhances egg production in pigeons. Additionally, under RL, the F2 follicle was selected, while under WL, small follicles were kept on the third day (LI3). A total of 5,510 circRNAs were identified across all samples, revealing differentially expressed circRNAs (DECs) in various comparisons: 627 in RF1 vs. WF1, 900 in RF2 vs. WF2, 606 in RF1 vs. RF2, and 937 in WF1 vs. WF2. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that host genes of DECs were enriched in pathways like steroid hormone biosynthesis, oocyte meiosis, GnRH signaling pathway, and apoptosis pathway. Moreover, circRNA_5497, circRNA_2016, and circRNA_3328 were common DECs across 4 groups, sharing miRNA binding sites with follicle selection-associated genes. In conclusion, our findings suggest that RL promotes egg production by stimulating follicle selection during LI, offering insights into the regulatory mechanisms of circRNAs in follicle selection under RL. This knowledge can help enhance the reproductive performance of pigeons.

Molecular genetic foundation of a sex-linked tailless trait in Hongshan chicken by whole genome data analysis

As a Chinese local chicken breed, Hongshan chickens have 2 kinds of tail feather phenotypes, normal and taillessness. Our previous studies showed that taillessness was a sex-linked dominant trait. Abnormal development of the tail vertebrae could be explained this phenomenon in some chicken breeds. However, the number of caudal vertebrae in rumpless Hongshan chickens was normal, so rumplessness in Hongshan chicken was not related to the development of the caudal vertebrae. Afterwards, we found that rumplessness in Hongshan was due to abnormal development of tail feather rather than abnormal development of caudal vertebrae. In order to understand the genetic foundation of the rumplessness of Hongshan chickens, we compared and reanalyzed 2 sets of data in normal and rumpless Hongshan chickens from our previous studies. By joint analysis of genome-wide selection signature analysis and genome-wide association approach, we found that 1 overlapping gene (EDIL3) and 16 peak genes (ENSGALG00000051843, ENSGALG00000053498, ENSGALG00000054800, KIF27, PTPRD, ENSGALG00000047579, ENSGALG00000041052, ARHGEF28, CAMK4, SERINC5, ENSGALG00000050776, ERCC8, MCC, ADAMTS19, ENSGALG00000053322, CHRNA8) located on the Z chromosome was associated with the rumpless trait. The results of this study furtherly revealed the molecular mechanism of the rumpless trait in Hongshan chickens, and identified the candidate genes associated with this trait. Our results will help to improve the shape of chicken tail feathers and to rise individual economic value in some specific market in China.

Comprehensive analysis of changes in expression of lncRNA, microRNA and mRNA in liver tissues of chickens with high or low abdominal fat deposition

1. The liver of chickens is a dominant lipid biosynthetic tissue and plays a vital role in fat deposition, particularly in the abdomen. To determine the molecular mechanisms involved in its lipid metabolism, the livers of chickens with high (H) or low (L) abdominal fat content were sampled and sequencing on long non-coding RNA (lncRNA), messenger RNA (mRNA) and small RNA (microRNA) was performed.2. In total, 351 expressed protein-coding genes for long non-coding RNA (DEL; 201 upregulated and 150 downregulated), 400 differentially expressed genes (DEG; 223 upregulated and 177 downregulated) and 10 differentially expressed miRNA (DEM; four upregulated and six downregulated) were identified between the two groups. Multiple potential signalling pathways related to lipogenesis and lipid metabolism were identified via pathway enrichment analysis. In addition, 173 lncRNA - miRNA - mRNA interaction regulatory networks were identified, including 30 lncRNA, 27 mRNA and seven miRNA.3. These networks may help regulate lipid metabolism and fat deposition. Five promising candidate genes and two lncRNA may play important roles in the regulation of adipogenesis and lipid metabolism in chickens.

Whole transcriptome sequencing reveals key genes and ceRNA regulatory networks associated with pimpled eggs in hens

Eggshell is one of the most important indicators of egg quality, and due to low shell strength, pimple eggs (PE) are more susceptible to breakage, thus causing huge economic losses to the egg industry. At the current time, the molecular mechanisms that regulate the formation of pimple eggs are poorly understood. In this study, uterine tissues of PE-laying hens (n = 8) and normal egg (NE) -laying hens (n = 8) were analyzed by whole transcriptome sequencing, and a total of 619 differentially expressed mRNAs (DE mRNAs), 122 differentially expressed lncRNAs (DE lncRNAs) and 21 differentially expressed miRNAs (DE miRNAs) were obtained. Based on the targeting relationship among DE mRNAs, DE lncRNAs and DE miRNAs, we constructed a competitive endogenous RNA (ceRNA) network including 12 DE miRNAs, 19 DE lncRNAs, and 128 DE mRNAs. Considering the large amount of information contained in the network, we constructed a smaller ceRNA network to better understand the complex mechanisms of pimple egg formation. The smaller ceRNA network network contains 7 DE lncRNAs (LOC107056551, LOC121109367, LOC121108909, LOC121108862, LOC112530033, LOC121113165, LOC107054145), 5 DE miRNAs (gga-miR-6568-3p, gga-miR-31-5p, gga-miR-18b-3p, gga-miR-1759-3p, gga-miR-12240-3p) and 7 DE mRNAs (CABP1, DNAJC5, HCN3, HPCA, IBSP, KCNT1, OTOP3), and these differentially expressed genes may play key regulatory roles in the formation of pimpled eggs in hens. This study provides the overall expression profiles of mRNAs, lncRNAs and miRNAs in the uterine tissues of hens, which provides a theoretical basis for further research on the molecular mechanisms of pimpled egg formation, and has potential applications in improving eggshell quality.

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.

Analysis of circRNA-miRNA-mRNA regulatory network of embryonic gonadal development in Mulard duck

The aim of the study was to explore the regulatory mechanism of differences in embryonic gonadal development between intergeneric distance hybrid offspring Mulard ducks and parent ducks. The morphological differences gonadal tissues of Muscovy ducks, Pekin ducks and Mulard ducks at 12.5-day embryonic age were observed by sectioning and hematoxylin-eosin (HE) staining. Then followed by transcriptome sequencing to screen for gonadal development-related differentially expressed circRNAs and mRNAs to construct a competitive endogenous RNA (ceRNA) regulatory network. Finally, qRT-PCR and luciferase reporter system were used to verify the sequencing data and targeting relationship of ceRNA pairs. The results showed that the seminiferous tubule lumen of Mulard ducks was not obvious, while there were obvious seminiferous tubules and tubular structures in testis of Pekin ducks and Muscovy ducks, with number and shape indicating maturity. There were 18 upregulated circRNAs and 16 downregulated circRNAs in Mulard ducks and Pekin ducks, respectively, and 39 upregulated circRNAs and 1 downregulated circRNA in Mulard ducks and Muscovy ducks, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis found that genes involves in dorso-ventral axis formation, for example, neurogenic locus notch homolog protein 1 (NOTCH1), were significantly enriched (P < 0.05). The novel_circ_0002265-gga-miR-122-5p-PAFAH1B2 regulatory network was constructed. The qRT-PCR results showed that the sequencing results were reliable. The dual-luciferase reporter assay showed that gga-miR-122-5p exists binding site of circ_0002265 and PAFAH1B2, indicating circ_0002265-gga-miR-122-5p-PAFAH1B2 targeting relationship. In summary, the embryonic gonadal development of intergeneric hybrid Mulard ducks may be regulated by differentially expressed circRNAs and genes, such as novel_circ_0000519, novel_circ_0003537, NOTCH1, FGFR2, PAFAH1B1, and PAFAH1B2, among which circ_0002265-gga-miR-122-5p-PAFAH1B2 may participate in the targeted regulation of gonadal development in Mulard ducks. The findings of this study are helpful for analyzing the mechanism of embryonic gonadal development differences in avians.

RNA-Seq Analysis Revealed circRNAs and Genes Associated with Abdominal Fat Deposition in Ducks

Fat deposition is an important factor affecting meat quality and feed conversion efficiency in meat ducks. This study aims to identify key circRNAs and genes affecting abdominal fat deposition. The correlations between abdominal fat and other growth performances were analyzed in 304 F2 generation of Cherry Valley duck Runzhou Crested White ducks, and an RNA-seq analysis of abdominal fat tissues from ducks with high and low rates of abdominal fat was performed. Growth performance results showed that Abdominal fat ratio and Intramuscular fat were significantly higher in the high rates of abdominal fat (HF)group than in the low rates of abdominal fat (LF) group for ducks. RNA-seq analysis of abdominal fat tissue unveiled 85 upregulated and 72 downregulated circRNAs among the differentially expressed ones. Notably, 74 circRNAs displayed more than four-fold differential expression, constituting 47.13% of the differentially expressed genes. Functional enrichment analysis of the differentially expressed circRNA source and target genes indicated that 17 circRNAs might partake in regulating duck abdominal fat production by influencing pathways like PPAR signaling, lipid droplets, and triglyceride metabolism. Lastly, multiple circRNA-microRNA-messenger RNA interaction networks were constructed. The results of this study establish the groundwork for understanding the molecular mechanisms that regulate abdominal fat deposition in ducks, offering a theoretical reference for the selective breeding of high-quality meat-producing ducks.

Identification of key genes affecting sperm motility in chicken based on whole-transcriptome sequencing

Sperm motility is an important index for the evaluation of semen quality. Improving sperm motility is important to improve reproductive performance, promote breeding process, and reduce production cost. However, the molecular mechanisms regulating sperm motility in chickens remain unclear. In this study, histological observation and whole-transcriptome analysis were performed on testicular tissue of chickens with high and low sperm motility. Histological observations showed that roosters with high sperm motility exhibited better semen quality than those with low sperm motility. In addition, the germinal epithelial cells of roosters with low sperm motility were loosely arranged and contained many vacuoles. RNA-seq results revealed the expression of 23,033 mRNAs, 2,893 lncRNAs, and 515 miRNAs in chicken testes. Among them, there were 417 differentially expressed mRNAs (DEmRNAs), 106 differentially expressed lncRNAs (DElncRNAs), and 15 differentially expressed miRNAs (DEmiRNAs) between high and low sperm motility testes. These differentially expressed genes were involved in the G protein-coupled receptor signaling pathway, cilia structure, Wnt signaling, MAPK signaling, GnRH signaling, and mTOR signaling. By integrating the competitive relationships between DEmRNAs, DElncRNAs, and DEmiRNAs, we identified the regulatory pathway of MSTRG.3077.3/MSTRG.9085.1-gga-miR-138-5p-CADM1 and MSTRG.2290.1-gga-miR-142-3p-GNAQ/PPP3CA as crucial in the modulation of chicken sperm motility. This study provides new insights into the function and mechanism of ceRNAs in regulating sperm motility in chicken testes.

Whole-Transcriptome RNA Sequencing Uncovers the Global Expression Changes and RNA Regulatory Networks in Duck Embryonic Myogenesis

Duck meat is pivotal in providing high-quality protein for human nutrition, underscoring the importance of studying duck myogenesis. The regulatory mechanisms governing duck myogenesis involve both coding and non-coding RNAs, yet their specific expression patterns and molecular mechanisms remain elusive. To address this knowledge gap, we performed expression profiling analyses of mRNAs, lncRNAs, circRNAs, and miRNAs involved in duck myogenesis using whole-transcriptome RNA-seq. Our analysis identified 1733 differentially expressed (DE)-mRNAs, 1116 DE-lncRNAs, 54 DE-circRNAs, and 174 DE-miRNAs when comparing myoblasts and myotubes. A GO analysis highlighted the enrichment of DE molecules in the extracellular region, protein binding, and exocyst. A KEGG analysis pinpointed pathways related to ferroptosis, PPAR signaling, nitrogen metabolism, cell cycle, cardiac muscle contraction, glycerolipid metabolism, and actin cytoskeleton. A total of 51 trans-acting lncRNAs, including ENSAPLT00020002101 and ENSAPLT00020012069, were predicted to participate in regulating myoblast proliferation and differentiation. Based on the ceRNAs, we constructed lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA ceRNA networks involving five miRNAs (miR-129-5p, miR-133a-5p, miR-22-3p, miR-27b-3p, and let-7b-5p) that are relevant to myogenesis. Furthermore, the GO and KEGG analyses of the DE-mRNAs within the ceRNA network underscored the significant enrichment of the glycerolipid metabolism pathway. We identified five different DE-mRNAs, specifically ENSAPLG00020001677, ENSAPLG00020002183, ENSAPLG00020005019, ENSAPLG00020010497, and ENSAPLG00020017682, as potential target genes that are crucial for myogenesis in the context of glycerolipid metabolism. These five mRNAs are integral to ceRNA networks, with miR-107_R-2 and miR-1260 emerging as key regulators. In summary, this study provides a valuable resource elucidating the intricate interplay of mRNA-lncRNA-circRNA-miRNA in duck myogenesis, shedding light on the molecular mechanisms that govern this critical biological process.

IMF deposition ceRNA network analysis and functional study of HIF1a in yak

The concentration of intramuscular fat (IMF) is a crucial determinant of yak meat quality. However, the molecular mechanisms that regulate IMF in yak remain largely elusive. In our study, we conducted transcriptome sequencing on the longissimus dorsi muscle tissues of yaks with varying IMF contents. We then filtered differentially expressed genes (DEGs), microRNAs (DEMs), and long non-coding RNAs (DELs) to elucidate potential regulatory pathways of adipogenesis in yaks. Overall, our research sheds light on an array of potential mRNAs and noncoding RNAs implicated in IMF deposition and elaborates on the role of HIF1α in yaks. These findings contribute valuable insights that can serve as a guide for further research into the molecular mechanisms governing IMF deposition.

Integrated Comparative Transcriptome and circRNA-lncRNA-miRNA-mRNA ceRNA Regulatory Network Analyses Identify Molecular Mechanisms Associated with Intramuscular Fat Content in Beef Cattle

Intramuscular fat content (IMF), one of the most important carcass traits in beef cattle, is controlled by complex regulatory factors. At present, molecular mechanisms involved in regulating IMF and fat metabolism in beef cattle are not well understood. Our objective was to integrate comparative transcriptomic and competing endogenous RNA (ceRNA) network analyses to identify candidate messenger RNAs (mRNAs) and regulatory RNAs involved in molecular regulation of longissimus dorsi muscle (LDM) tissue for IMF and fat metabolism of 5 beef cattle breeds (Angus, Chinese Simmental, Luxi, Nanyang, and Shandong Black). In total, 34 circRNAs, 57 lncRNAs, 15 miRNAs, and 374 mRNAs were identified by integrating gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Furthermore, 7 key subnets with 16 circRNAs, 43 lncRNAs, 7 miRNAs, and 237 mRNAs were detected through clustering analyses, whereas GO enrichment analysis of identified RNAs revealed 48, 13, and 28 significantly enriched GO terms related to IMF in biological process, molecular function, and cellular component categories, respectively. The main metabolic-signaling pathways associated with IMF and fat metabolism that were enriched included metabolic, calcium, cGMP-PKG, thyroid hormone, and oxytocin signaling pathways. Moreover, MCU, CYB5R1, and BAG3 genes were common among the 10 comparative groups defined as important candidate marker genes for fat metabolism in beef cattle. Contributions of transcriptome profiles from various beef breeds and a competing endogenous RNA (ceRNA) regulatory network underlying phenotypic differences in IMF provided novel insights into molecular mechanisms associated with meat quality.

CircDOCK7 facilitates the proliferation and adipogenic differentiation of chicken abdominal preadipocytes through the gga-miR-301b-3p/ACSL1 axis

BACKGROUND

Abdominal fat deposition depends on both the proliferation of preadipocytes and their maturation into adipocytes, which is a well-orchestrated multistep process involving many regulatory molecules. Circular RNAs (circRNAs) have emergingly been implicated in mammalian adipogenesis. However, circRNA-mediated regulation in chicken adipogenesis remains unclear. Our previous circRNA sequencing data identified a differentially expressed novel circRNA, 8:27,886,180|27,889,657, during the adipogenic differentiation of chicken abdominal preadipocytes. This study aimed to investigate the regulatory role of circDOCK7 in the proliferation and adipogenic differentiation of chicken abdominal preadipocytes, and explore its molecular mechanisms of competing endogenous RNA underlying chicken adipogenesis.

RESULTS

Our results showed that 8:27,886,180|27,889,657 is an exonic circRNA derived from the head-to-tail splicing of exons 19-22 of the dedicator of cytokinesis 7 (DOCK7) gene, abbreviated as circDOCK7. CircDOCK7 is mainly distributed in the cytoplasm of chicken abdominal preadipocytes and is stable because of its RNase R resistance and longer half-life. CircDOCK7 is significantly upregulated in the abdominal fat tissues of fat chickens compared to lean chickens, and its expression gradually increases during the proliferation and adipogenic differentiation of chicken abdominal preadipocytes. Functionally, the gain- and loss-of-function experiments showed that circDOCK7 promoted proliferation, G0/G1- to S-phase progression, and glucose uptake capacity of chicken abdominal preadipocytes, in parallel with adipogenic differentiation characterized by remarkably increased intracellular lipid droplet accumulation and triglyceride and acetyl coenzyme A content in differentiated chicken abdominal preadipocytes. Mechanistically, a pull-down assay and a dual-luciferase reporter assay confirmed that circDOCK7 interacted with gga-miR-301b-3p, which was identified as an inhibitor of chicken abdominal adipogenesis. Moreover, the ACSL1 gene was demonstrated to be a direct target of gga-miR-301b-3p. Chicken ACSL1 protein is localized in the endoplasmic reticulum and mitochondria of chicken abdominal preadipocytes and acts as an adipogenesis accelerator. Rescue experiments showed that circDOCK7 could counteract the inhibitory effects of gga-miR-301b-3p on ACSL1 mRNA abundance as well as the proliferation and adipogenic differentiation of chicken abdominal preadipocytes.

CONCLUSIONS

CircDOCK7 serves as a miRNA sponge that directly sequesters gga-miR-301b-3p away from the ACSL1 gene, thus augmenting adipogenesis in chickens. These findings may elucidate a new regulatory mechanism underlying abdominal fat deposition in chickens.