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

Mechanisms of Adipose Tissue Metabolism in Naturally Grazing Sheep at Different Growth Stages: Insights from mRNA and miRNA Profiles

Adipose tissue metabolism plays a crucial role in sheep meat quality and the optimization of adipose tissue utilization. To reveal the molecular mechanisms of adipose tissue metabolism during growth in naturally grazing sheep, we investigated the mRNA and miRNA profiles in subcutaneous adipose tissue (SAT) from naturally grazing Sunit sheep at 6, 18, and 30 months of age (Mth-6, Mth-18, and Mth-30). We identified 927 differentially expressed (DE) genes and 134 DE miRNAs in the SAT of sheep at different growth stages. Specifically, the expressions of ACACA, FASN, DGAT2, GPAM, SCD, ELOVL6, HSD17B12, TECR, PKM, TKT, PCK1, CD44, and THBS2S genes were significantly upregulated in Mth-18 and Mth-30 compared to that in Mth-6. These genes promoted fatty acid synthesis, triglyceride synthesis, gluconeogenesis, and extracellular matrix-receptor interaction and decreased glycolysis, leading to increased adipocyte proliferation and fat deposition. Notably, our findings suggested that the reduced activity of the AMPK signaling pathway may be regulated by CAMKK2 and PP2A during sheep growth. Furthermore, our results revealed several DE miRNAs, mml-miR-320b, chi-miR-1388-3p, bta-miR-6715, oar-miR-143, and miR-424, that potentially influence fat metabolism. Overall, this study provides a theoretical basis and new insights into the molecular mechanisms of adipose tissue metabolism during growth in naturally grazing sheep.

Identification of the crucial circ-mi-mRNA interaction networks regulating testicular development and spermatogenesis in ganders

Semen quality has an important impact on the reproductive performance of ganders, and the quantity and quality of spermatozoa in semen are the determinants of semen quality. In our practical work, a small number of azoospermic ganders were observed in adult goose breeding populations, but the underlying regulatory mechanisms remain unknown. In the present study, we firstly compared the morphological and histological differences in the testes of ganders from normozoospermic group (NG) and azoospermic group (AG), and then analyzed the testicular expression patterns of circRNAs, miRNAs, and mRNAs between the two groups by using whole-transcriptome sequencing technology. Results from histomorphological analysis demonstrated that the body weight alone was not accountable for the occurrence of gander azoospermia, and the possible cause might be the observed testicular abnormalities. At the morphological level, the left, right, and bilateral testicular weights, the right and bilateral testicular organ indexes, and the long, short, and dorsoventral diameters of the left, right and bilateral testes were significantly lower in AG than in NG (P < 0.05). At the histological level, most testicular histological parameters, such as the testicular parenchymal area, the diameter of seminiferous tubules, and the number of germ cells, were significantly higher (P < 0.05) in NG than in AG. The RNA-seq results showed that a total of 683 differentially expressed circRNAs (DEcircRNAs), 24 differentially expressed miRNAs (DEmiRNAs), and 1,118 differentially expressed Genes (DEGs) were identified in the gander testes between NG and AG. Subsequent functional enrichment analysis revealed that most of the DEGs and the target genes of DEcircRNAs and DEmiRNAs were significantly enriched in either the biological processes related to male gonad development, spermatid development, and regulation of cell differentiation or the KEGG terms including the MAPK, TGF-beta, Wnt, and cell cycle signaling pathways. By constructing the core ceRNA regulatory networks, several key DEcircRNAs, including 1:98100313|98104995, 1:171413706|171419341, 6:3414226|3418193, and 2:115876735|115880760, were identified to regulate the expression of TGFB2 and BCL2 through interactions with specific miRNAs such as novel-miR-265 and novel-miR-266, and such interactions could play crucial roles in regulating the gander testicular cell apoptosis, proliferation, and spermatogenesis. This study provides novel insights into the function and molecular mechanisms of ceRNAs in regulating the gander testicular development and semen quality.

Whole Transcriptome Sequencing Reveals miRNAs and ceRNA Networks in Duck Abdominal Fat Deposition

Excessive deposition of abdominal fat will cause a waste of resources. In order to explore the key miRNAs and circRNA/lncRNA-miRNA-mRNA ceRNA regulatory network involved in regulating abdominal fat deposition, hematoxylin and eosin (H&E) staining was performed on abdominal fat tissues of ducks in the high abdominal fat rate group (HF) and low abdominal fat rate group (LF) at 21 and 42 days of age, and whole transcriptome sequencing was performed on abdominal tissues of ducks in the HF and LF groups at 42 days of age. The results showed that the number of adipocytes in ducks in the HF group was significantly higher than that in the LF group at 21 days of age (p < 0.001), while the number of adipocytes in ducks in the HF group at 42 days of age was significantly lower than that in the LF group (p < 0.001). In addition, transcriptome sequencing screened out a total of 14 differentially expressed miRNAs (10 miRNAs were significantly up-regulated, and 4 miRNAs were significantly down-regulated). By predicting the target genes of these differentially expressed miRNAs, a total of 305 target genes were obtained. Further analysis of miRNA target genes using GO and KEGG functional enrichment analyses revealed that these target genes were significantly enriched in the GnRH signaling pathway, the PPAR signaling pathway, insulin resistance, the mTOR signaling pathway, the AMPK signaling pathway, the FoxO signaling pathway, and other pathways related to adipose development. In addition, miRNA-205-x, miRNA-6529-x, miRNA-194-x, miRNA-215-x, miRNA-3074-x, miRNA-2954-x, novel-m0133-3p, and novel-m0156-5p were found to be important candidate miRNAs for abdominal fat deposition in ducks. These miRNAs were related to the expression of FOXO3, LIFR, Pdk4, PPARA, FBN1, MYH10, Cd44, PRELP, Esrrg, AKT3, and STC2. Based on these eight candidate miRNAs, a ceRNA regulatory network of circRNA/lncRNA-miRNA-mRNA regulating abdominal fat deposition was successfully constructed. The results of this study will provide a useful reference for accelerating the understanding of the molecular mechanism of duck abdominal fat deposition.

Construction of an lncRNA-mediated ceRNA network to investigate the inflammatory regulatory mechanisms of ischemic stroke

Long non-coding RNAs (lncRNAs) are among the most abundant types of non-coding RNAs in the genome and exhibit particularly high expression levels in the brain, where they play crucial roles in various neurophysiological and neuropathological processes. Although ischemic stroke is a complex multifactorial disease, the involvement of brain-derived lncRNAs in its intricate regulatory networks remains inadequately understood. In this study, we established a cerebral ischemia-reperfusion injury model using middle cerebral artery occlusion (MCAO) in male Sprague-Dawley rats. High-throughput sequencing was performed to profile the expression of cortical lncRNAs post-stroke, with subsequent validation using RT-PCR and qRT-PCR. Among the 31,183 lncRNAs detected in the rat cerebral cortex, 551 were differentially expressed between the MCAO and sham-operated groups in the ipsilateral cortex (fold change ≥2.0, P < 0.05). An integrated analysis of the 20 most abundant and significantly differentially expressed lncRNAs (DELs) identified 25 core cytoplasmic DELs, which were used to construct an interaction network based on their targeting relationships. This led to the establishment of a comprehensive lncRNA-miRNA-mRNA regulatory network comprising 12 lncRNAs, 16 sponge miRNAs, and 191 target mRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed that differentially expressed mRNAs (DEmRNAs) were significantly enriched in stroke-related pathways. Our analysis predicted four key lncRNAs, four miRNAs, and eleven crucial mRNAs involved in post-transcriptional regulation through competing endogenous RNA (ceRNA) mechanisms. These molecules were shown to participate extensively in post-stroke processes, including angiogenesis, axonal regeneration, inflammatory responses, microglial activation, blood-brain barrier (BBB) disruption, apoptosis, autophagy, ferroptosis, and thrombocytopenia. These findings highlight the role of lncRNAs as multi-level regulators in the complex network of post-stroke mechanisms, providing novel insights into the pathophysiological processes of stroke.

Emerging roles of circular RNAs on the regulation of production traits in chicken

Chickens are vital agricultural animals that supply a significant portion of the protein consumed by humans. In society today, enhancing the productive performance of chickens in a safe and efficient manner has become a central focus of research. This performance is determined by various production traits that are primarily influenced by multiple factors, including epigenetics-a critical aspect of gene regulation. Circular RNAs (circRNAs), a unique class of non-coding RNAs, have emerged as key epigenetic regulators. Recent studies have demonstrated that circRNAs are extensively engaged in numerous production traits, which include skeletal muscle formation, fat deposition, ovarian follicle development, liver function, bone development, immunity, and resistance to environmental stress. These processes play crucial roles in determining the overall productivity of chickens. Given the significance of circRNAs in these various traits, this article provides a comprehensive review of the functional circRNAs associated with different traits in chickens, serving as a valuable theoretical reference for future research. Further investigation into the role of circRNAs may reveal novel insights into the molecular mechanisms underlying key economic traits in chickens and pave the way for innovative strategies in molecular breeding aimed at enhancing chicken productive performance.

Characterization of the regulatory network and pathways in duodenum affecting chicken abdominal fat deposition

The excessive accumulation of abdominal fat in chickens has resulted in a reduction in both the feed conversion efficiency and the slaughter yield. To elucidate the regulatory mechanisms and metabolic pathways affecting abdominal fat deposition in the context of broiler breeding, a cohort of 400 Qingyuan partridge chickens with varying abdominal fat deposition was established. Whole transcriptome sequencing analyses were conducted on the duodenum of 20 representative chickens to ascertain the regulatory networks at this vital digestive and absorptive organ. Consequently, 116 differentially expressed genes were identified, exhibiting a trend of increasing or decreasing expression in correlation with the accumulation of abdominal fat. A total of 36 DEmRNAs, 170 DElncRNAs, 92 DEcircRNAs and 88 DEmiRNAs were identified as differentially expressed between chickens with extremely high and low abdominal fat deposition. The functional enrichment analyses demonstrated that the differentially expressed RNA in the duodenum were involved in the regulation of chicken abdominal fat deposition by mediating a series of metabolic pathways, including the Wnt signaling pathway, the PPAR signaling pathway, the Hippo signaling pathway, the FoxO signaling pathway, the MAPK signaling pathway and other signaling pathways that are involved in fatty acid metabolism and degradation. The construction of putative interaction pairs led to the suggestion of two lncRNA-miRNA-mRNA ceRNA networks comprising two mRNAs, two miRNAs, and 29 lncRNAs, as well as two circRNA-lncRNA-miRNA-mRNA ceRNA networks comprising 26 mRNAs, 12 miRNAs, 17 lncRNAs, and nine circRNAs, as core regulatory networks in the duodenum affecting chicken abdominal fat deposition. The aforementioned genes including TMEM150C, REXO1, PIK3C2G, ppp1cb, PARP12, SERPINE2, LRAT, CYP1A1, INSR and APOA4, were proposed as candidate genes, while the miRNAs, including miR-107-y, miR-22-y, miR-25-y, miR-2404-x and miR-16-x, as well as lncRNAs such as ENSGALT00000100291, TCONS_00063508, TCONS_00061201 and TCONS_00079402 were the candidate regulators associated with chicken abdominal fat deposition. The findings of this study provide a theoretical foundation for the molecular mechanisms of mRNAs and non-coding RNAs in duodenal tissues on abdominal fat deposition in chickens.

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.