Identification of Candidate Circular RNAs Underlying Intramuscular Fat Content in the Donkey

Muscle transcriptome analysis reveal candidate genes and pathways related to fat and lipid metabolism in Yunling cattle

Yunling cattle (YL) is a recently developed beef breed harboring a quarter of Yunnan ancestral cattle genome, spanning over past 30 years. Compared with Diqing cattle (DQ), a Yunnan native cattle breed, YL presents various advantages, including rapid growth and exquisite meat quality. However, the molecular mechanisms underlying these phenotypic differences are not clearly understood. To further identify the candidate genes responsible for the quality of the meat in the muscle, longissimus dorsi (LD) muscle was used for RNA-Seq analysis. A total of 508 differentially expressed genes (DEGs) were identified in YL (adjusted p-value <0.01 and log2FoldChange >1), of which 243 were up-regulated and 265 were down-regulated. Functional association analysis showed that the identified DEGs mainly enriched the lipid and fat metabolism pathways. Moreover, it was also observed that several fat-related genes were differentially expressed in both cattle breeds, including three up-regulated genes (MOGAT1, ACSM3, PLPP2) and two down-regulated genes (ADIG, GPAT3). In addition, alternative splice analysis was also performed revealing an important 9-11 exon skipping variation of GPAM gene (crucial for beef marbling) in YL, which is three times higher than that in DQ, suggesting that this variation might have played the central role in the 'snow beef' effect in YL. We believe that our results will help in understanding the mechanism of muscle development and promote the further breeding programs in YL cattle.

Functional analysis of differentially expressed circular RNAs in sheep subcutaneous fat

BACKGROUND

Circular RNAs (circRNAs), as important non-coding RNAs (ncRNAs), are involved in many biological activities. However, the exact chemical mechanism behind fat accumulation is unknown. In this paper, we obtained the expression profiles of circRNAs using high-throughput sequencing and investigated their differential expression in subcutaneous fat tissue of Duolang and Small Tail Han sheep.

RESULTS

From the transcriptomic analysis, 141 differentially expressed circRNAs were identified, comprising 61 up-regulated circRNAs and 80 down-regulated circRNAs. These host genes were primarily enriched in the MAPK and AMPK signaling pathways which is closely associated with fat deposition regulation. We identified circRNA812, circRNA91, and circRNA388 as vital genes in fat deposition by miRNA-circRNA target gene prediction. The functional annotation results of target genes of key circRNAs showed that the signaling pathways mainly included PI3K-Akt and AMPK. We constructed the competing endogenous RNA (ceRNA) regulatory network to study the role of circRNAs in sheep lipid deposition, and circRNA812, circRNA91, and circRNA388 can adsorb more miRNAs. NC_040253.1_5757, as the source of miRNA response element (MRE) among the three, may play an important role during the process of sheep fat deposition.

CONCLUSIONS

Our study gives a systematic examination of the circRNA profiles expressed in sheep subcutaneous fat. These results from this study provide some new basis for understanding circRNA function and sheep fat metabolism.

Exploring Genetic Markers: Mitochondrial DNA and Genomic Screening for Biodiversity and Production Traits in Donkeys

Donkeys (Equus asinus) play a pivotal role as essential livestock in arid and semi-arid regions, serving various purposes such as transportation, agriculture, and milk production. Despite their significance, donkey breeding has often been overlooked in comparison to other livestock species, resulting in limited genetic improvement programs. Preserving donkey genetic resources within each country necessitates the establishment of breed conservation programs, focusing on managing genetic diversity among populations. In recent years, significant strides have been made in sequencing and analyzing complete mitochondrial DNA (mtDNA) molecules in donkeys. Notably, numerous studies have honed in on the mitochondrial D-loop region, renowned for its remarkable variability and higher substitution rate within the mtDNA genome, rendering it an effective genetic marker for assessing genetic diversity in donkeys. Furthermore, genetic markers at the RNA/DNA level have emerged as indispensable tools for enhancing production and reproduction traits in donkeys. Traditional animal breeding approaches based solely on phenotypic traits, such as milk yields, weight, and height, are influenced by both genetic and environmental factors. To overcome these challenges, genetic markers, such as polymorphisms, InDel, or entire gene sequences associated with desirable traits in animals, have achieved widespread usage in animal breeding practices. These markers have proven increasingly valuable for facilitating the selection of productive and reproductive traits in donkeys. This comprehensive review examines the cutting-edge research on mitochondrial DNA as a tool for assessing donkey biodiversity. Additionally, it highlights the role of genetic markers at the DNA/RNA level, enabling the informed selection of optimal production and reproductive traits in donkeys, thereby driving advancements in donkey genetic conservation and breeding programs.

Slaughter performance of the main goose breeds raised commercially in China and nutritional value of the meats of the goose breeds: a systematic review

A number of goose breeds are raised commercially in China. However, the data on the slaughter performance of the goose breeds and the nutritional value of their meats lack a thorough comparative analysis. In this systematic review, the slaughter performance of the goose breeds and nutritional value of their meats were comparatively analyzed to provide an overview of the characteristics of the goose breeds raised commercially in China. Fifteen goose breeds were selected from 27 research articles published up to January 2022 on the slaughter performance of the goose breeds raised commercially in China and their nutrient composition after literature searching, literature screening, variety selection, and data collation. The slaughter indexes of the goose breeds and the basic nutrient composition, amino acid composition, and fatty acid composition of the meats of the goose breeds were standardized using min-max normalization and compared. The results suggest that the slaughter indexes and nutritional indicators of the meats of Yangzhou white goose, Xupu goose, Landaise geese, and Sichuan white goose are more balanced than those of the meats of the other goose breeds. The results of this review can lay the foundation for optimizing the breeding methods of the commercially raised goose breeds and processing methods of the meats of the geese. © 2022 Society of Chemical Industry.

Association of MLL3 and TGF-β signaling gene polymorphisms with the susceptibility and prognostic outcomes of Stanford type B aortic dissection : MLL3 with TGF-β signal pathway association with Stanford type B AD

OBJECTIVE

This study aims to investigate the association of lysine methyltransferase 2 C (MLL3) and transforming growth factor β (TGF-β) signaling-related gene polymorphisms with the susceptibility of Stanford type B aortic dissection (AD) and its clinical prognostic outcomes. The methods involved investigating the MLL3 (rs10244604, rs6963460, rs1137721), TGFβ1 (rs1800469), TGFβ2 (rs900), TGFR1 (rs1626340) and TGFR2 (rs4522809) gene polymorphisms. Logistic regression was performed to investigate the association between 7 single nucleotide gene polymorphisms (SNPs) and Stanford type B aortic dissection. The GMDR software was used to analyze gene-gene and gene-environment interactions. The odds ratio (OR) with a 95% confidence interval (CI) was employed to evaluate the association of genes and Stanford type B AD risk.

RESULTS

Genotypes and allele distributions in the case and control groups showed significant differences (P < 0.05). Logistic regression has shown that the Stanford Type B AD risk was highest in individuals with the rs1137721 CT genotype (OR = 4.33, 95% CI = 1.51-12.40). Additionally, WBC, drinking, hypertension, triglycerides (TG), and low-density lipoprotein (LDL-C) were independent risk factors for Stanford Type B AD. Logistic regression showed that the Stanford Type B AD risk was highest in individuals with the MLL3 (rs1137721)-TT + CT and TGFβ1 (rs4522809)-AA genotype (OR = 6.72, 95% CI = 1.56-29.84), and lowest in those with the MLL3 (rs1137721)-CC and TGFβ1 (rs4522809)-AA + GG genotype (OR = 4.38, 95% CI = 0.92-20.83). However, the 55-month median long-term follow-up did not show statistical significance.

CONCLUSION

Carriers of both TT + CT of MLL3 (rs1137721) and AA of TGFβ1 (rs4522809) polymorphisms may be closely related to the development of Stanford type B AD. MLL3 (rs1137721), WBC, and TG/TC were found to be associated with the morbidity of Stanford type B AD. MLL3 (KMT2C) is associated with the TGF-β signaling pathway protein. The risk of Stanford type B AD is related to the interactions of gene-gene and gene-environment.

Research progress of non-coding RNAs regulation on intramuscular adipocytes in domestic animals

Intramuscular fat (IMF) is the main determinant of the economic value of domestic animal meat, and has a vital impact on the sensory quality characteristics, while the content of IMF is mainly determined by the size and number of intramuscular adipocytes. In recent years, due to the development of sequencing technology and omics technology, a large number of non-coding RNAs have been identified in intramuscular adipocytes. Non-coding RNAs are a kind of RNA regulatory factors with biological functions but without translation function, which mainly include microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs). These non-coding RNAs regulate the key genes of intramuscular adipocyte growth and development at post-transcriptional level through a variety of regulatory mechanisms, and affect the number and size of intramuscular adipocytes, thus affecting the content of IMF. Here, the review summarizes the candidate non-coding RNAs (miRNAs, lncRNAs, circRNAs) and genes involved in the regulation of intramuscular adipocytes, the related regulation mechanism and signaling pathways, in order to provide reference for further clarifying the molecular regulation mechanism of non-coding RNAs on intramuscular adipocytes in domestic animals.

In Silico Identification and Characterization of circRNAs as Potential Virulence-Related miRNA/siRNA Sponges from Entamoeba histolytica and Encystment-Related circRNAs from Entamoeba invadens

Ubiquitous eukaryotic non-coding circular RNAs regulate transcription and translation. We have reported full-length intronic circular RNAs (flicRNAs) in Entamoeba histolytica with esterified 3′ss and 5′ss. Their 5′ss GU-rich elements are essential for their biogenesis and their suggested role in transcription regulation. Here, we explored whether exonic, exonic-intronic, and intergenic circular RNAs are also part of the E. histolytica and E. invadens ncRNA RNAome and investigated their possible functions. Available RNA-Seq libraries were analyzed with the CIRI-full software in search of circular exonic RNAs (circRNAs). The robustness of the analyses was validated using synthetic decoy sequences with bona fide back splice junctions. Differentially expressed (DE) circRNAs, between the virulent HM1:IMSS and the nonvirulent Rahman E. histolytica strains, were identified, and their miRNA sponging potential was analyzed using the intaRNA software. Respectively, 188 and 605 reverse overlapped circRNAs from E. invadens and E. histolytica were identified. The sequence composition of the circRNAs was mostly exonic although different to human circRNAs in other attributes. 416 circRNAs from E. histolytica were virulent-specific and 267 were nonvirulent-specific. Out of the common circRNAs, 32 were DE between strains. Finally, we predicted that 8 of the DE circRNAs could function as sponges of the bioinformatically reported miRNAs in E. histolytica, whose functions are still unknown. Our results extend the E. histolytica RNAome and allow us to devise a hypothesis to test circRNAs/miRNAs/siRNAs interactions in determining the virulent/nonvirulent phenotypes and to explore other regulatory mechanisms during amoebic encystment.

Profile analysis and functional modeling identify circular RNAs in nonalcoholic fatty liver disease as regulators of hepatic lipid metabolism

Nonalcoholic fatty liver disease (NAFLD) is the leading cause of chronic liver disease, associated with an outcome of hepatic fibrosis/cirrhosis and hepatocellular carcinoma. However, limited exploration of the underlying mechanisms hinders its prevention and treatment. To investigate the mechanisms of epigenetic regulation in NAFLD, the expression profile of circular RNA (circRNA) of rodents in which NAFLD was induced by a high-fat, high-cholesterol (HFHC) diet was studied. Modeling of the circRNA-microRNA (miRNA) -mRNA regulatory network revealed the functional characteristics of NAFLD-specific circRNAs. The targets and effects in the liver of such NAFLD-specific circRNAs were further assessed. Our results uncovered that the downregulation of 28 annotated circRNAs characterizes HFHC diet-induced NAFLD. Among the downregulated circRNAs, long intergenic non-protein coding RNA, P53 induced transcript (LNCPINT) -derived circRNAs (circ_0001452, circ_0001453, and circ_0001454) targeted both miR-466i-3p and miR-669c-3p. Their deficiency in NAFLD abrogated the circRNA-based inhibitory effect on both miRNAs, which further inactivated the AMPK signaling pathway via AMPK-α1 suppression. Inhibition of the AMPK signaling pathway promotes hepatic steatosis, depending on the transcriptional and translational upregulation of lipogenic genes, such as those encoding sterol regulatory element-binding protein 1 (SREBP1) and fatty acid synthase (FASN) in hepatocytes. The levels of LNCPINT-derived circRNAs displayed a negative association with hepatic triglyceride (TG) concentration. These findings suggest that loss of LNCPINT-derived circRNAs may underlie NAFLD via miR-466i-3p- and miR-669c-3p-dependent inactivation of the AMPK signaling pathway.

Comparative transcriptome analysis of longissimus dorsi tissues with different intramuscular fat contents from Guangling donkeys

Background

Donkey meat has low fat and high protein contents and is rich in various unsaturated fatty acids and trace elements that are beneficial to human digestion and absorption. IMF (intramuscular fat), also known as marbling, is an important indicator of the lean meat to fat ratio, which directly affects the tenderness and juiciness of the meat. At present, the underlying molecular variations affecting IMF content among donkey breeds are unclear. The Guangling donkey is an indigenous species in China. This study explored candidate regulatory genes that affect IMF content in Guangling donkeys. The IMF content of the longissimus dorsi muscle in 30 Guangling donkeys was measured. Six donkeys of similar age were selected according to age factors and divided into two groups, the high (H) and low (L) fat groups, according to their IMF content.

Results

RNA-seq technology was used to compare the muscle transcriptome between the two groups. More than 75.0% of alternative splicing (AS) events were of the skipped exon (SE) type. A total of 887 novel genes were identified; only 386 novel genes were aligned to the annotation information of various databases. Transcriptomics analysis revealed 167 differentially expressed genes (DEGs), of which 64 were upregulated and 103 were downregulated between the H and L groups. Gene ontology analysis showed that the DEGs were enriched in multiple biological processes and pathways that are related to adipocyte differentiation, lipid synthesis, and neutral lipid metabolism. KEGG pathway analysis suggested that arachidonic acid metabolism, the HIF-1 signalling pathway, fructose and mannose metabolism, glycerophospholipid metabolism, and the AMPK signalling pathway were involved in lipid deposition. In addition, a gene–gene interaction network was constructed that revealed that the DEGs, including SCD, LEPR, CIDEA, DLK1, DGAT2, ITGAL, HMOX1, WNT10B, and DGKA, had significant roles in adipocyte differentiation and adipogenesis. The selected DEGs were further validated by qRT–PCR.

Conclusion

This study improves the in-depth understanding of gene regulation and protein expression regarding IMF deposition and lays a basis for subsequent molecular breeding studies in Guangling donkeys.

Genome-Wide Expression Profiling and Networking Reveals an Imperative Role of IMF-Associated Novel CircRNAs as ceRNA in Pigs

Intramuscular fat (IMF) deposition is a biological process that has a strong impact on the nutritional and sensorial properties of meat, with relevant consequences on human health. Pork loins determine the effects of marbling on the sensory attributes and meat quality properties, which differ among various pig breeds. This study explores the crosstalk of non-coding RNAs with mRNAs and analyzes the potential pathogenic role of IMF-associated competing endogenous RNA (ceRNA) in IMF tissues, which offer a framework for the functional validation of key/potential genes. A high-throughput whole-genome transcriptome analysis of IMF tissues from longissimus dorsi muscles of Large White (D_JN) and Laiwu (L_JN) pigs resulted in the identification of 283 differentially expressed circRNAs (DECs), including two key circRNAs (circRNA-23437, circRNA-08840) with potential binding sites for multiple miRNAs regulating the whole network. The potential ceRNA mechanism identified the DEC target miRNAs-mRNAs involved in lipid metabolism, fat deposition, meat quality, and metabolic syndrome via the circRNA-miRNA-mRNA network, concluding that ssc-mir-370 is the most important target miRNA shared by both key circRNAs. TGM2, SLC5A6, ECI1, FASN, PER1, SLC25A34, SOD1, and COL5A3 were identified as hub genes through an intensive protein-protein interaction (PPI) network analysis of target genes acquired from the ceRNA regulatory network. Functional enrichments, pathway examinations, and qRT-PCR analyses infer their implications in fat/cholesterol metabolism, insulin secretion, and fatty acid biosynthesis. Here, circRNAs and miRNA sequencing accompanied by computational techniques were performed to analyze their expressions in IMF tissues from the longissimus dorsi muscles of two pig breeds. Their target gene evolutionary trajectories, expression profiling, functional enrichments, subcellular localizations, and structural advances with high-throughput protein modeling, following genomic organizations, will provide new insights into the underlying molecular mechanisms of adipocyte differentiation and IMF deposition and a much-needed qualitative framework for future research to improve meat quality and its role as a biomarker to treat lipid metabolic syndromes.

Effects of intramuscular fat on meat quality and its regulation mechanism in Tan sheep

Intramuscular fat (IMF) contributes importantly to various aspects of meat quality, and genetic regulation is an effective pathway to improve IMF deposition in sheep. In this study, we systematically explored the effect of IMF content on meat quality in Tan sheep and investigated the regulatory mechanism of flavor precursors metabolism and IMF deposition. The results revealed that IMF significantly affected meat color, total muscle fiber numbers, and muscle fiber types in Tan sheep. Widely-targeted metabolomic analysis showed that amino acids were the main differential flavor precursors between lambs with different IMF content. Importantly, the comparison of fatty acid profiles revealed that saturated fatty acids and monounsaturated fatty acids are beneficial for IMF deposition. Furthermore, integrated analysis between metabolome and transcriptome indicated that MME is a key gene resulting in the reduction of amino acids in lambs with high IMF content; and the joint analysis between fatty acid profiles and transcript profiles showed that ADIPOQ, FABP4, PLIN1, PPARGC1A, SLC2A1 accelerated IMF deposition through positive regulation of saturated fatty acids and monounsaturated fatty acids metabolism. These results revealed key changes in meat quality affected by IMF content and the corresponding genetic mechanism, which may provide a new insight for understanding the IMF differential deposition and for improving meat quality in Tan sheep.

Comparative analysis of differentially abundant proteins between high and low intramuscular fat content groups in donkeys

Intramuscular fat (IMF) is an important regulator that determines meat quality, and its content is closely related to flavor, tenderness, and juiciness. Many studies have used quantitative proteomic analysis to identify proteins associated with meat quality traits in livestock, however, the potential candidate proteins that influence IMF in donkey muscle are not fully understood. In this study, we performed quantitative proteomic analysis, with tandem-mass-tagged (TMT) labeling, with samples from the longissimus dorsi (LD) muscle of the donkey. A total of 585,555 spectra were identified from the six muscle samples used in this study. In total, 20,583 peptides were detected, including 15,279 unique peptides, and 2,540 proteins were identified. We analyzed differentially abundant proteins (DAPs) between LD muscles of donkeys with high (H) and low (L) IMF content. We identified 30 DAPs between the H and L IMF content groups, of which 17 were upregulated and 13 downregulated in the H IMF group. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis of these DAPs revealed many GO terms (e.g., bone morphogenetic protein (BMP) receptor binding) and pathways (e.g., Wnt signaling pathway and Hippo signaling pathway) involved in lipid metabolism and adipogenesis. The construction of protein–protein interaction networks identified 16 DAPs involved in these networks. Our data provide a basis for future investigations into candidate proteins involved in IMF deposition and potential new approaches to improve meat quality in the donkey.

Promoter-Bound Full-Length Intronic Circular RNAs-RNA Polymerase II Complexes Regulate Gene Expression in the Human Parasite Entamoeba histolytica

Ubiquitous eukaryotic non-coding circular RNAs are involved in numerous co- and post-transcriptional regulatory mechanisms. Recently, we reported full-length intronic circular RNAs (flicRNAs) in Entamoeba histolytica, with 3′ss–5′ss ligation points and 5′ss GU-rich elements essential for their biogenesis and their suggested role in transcription regulation. Here, we explored how flicRNAs impact gene expression regulation. Using CLIP assays, followed by qRT-PCR, we identified that the RabX13 control flicRNA and virulence-associated flicRNAs were bound to the HA-tagged RNA Pol II C-terminus domain in E. histolytica transformants. The U2 snRNA was also present in such complexes, indicating that they belonged to transcription initiation/elongation complexes. Correspondingly, inhibition of the second step of splicing using boric acid reduced flicRNA formation and modified the expression of their parental genes and non-related genes. flicRNAs were also recovered from chromatin immunoprecipitation eluates, indicating that the flicRNA-Pol II complex was formed in the promoter of their cognate genes. Finally, two flicRNAs were found to be cytosolic, whose functions remain to be uncovered. Here, we provide novel evidence of the role of flicRNAs in gene expression regulation in cis, apparently in a widespread fashion, as an element bound to the RNA polymerase II transcription initiation complex, in E. histolytica.

Identification and Characterization of Circular RNAs in Association With the Deposition of Intramuscular Fat in Aohan Fine-Wool Sheep

Aohan fine-wool sheep (AFWS) is a high-quality fine-wool sheep breed that supplies wool and meat. Research is needed on the molecular mechanism behind intramuscular fat (IMF) deposition that greatly improves mutton quality. The widely expressed non-coding RNA is physiologically used in roles such as competitive endogenous RNA (ceRNA) that includes circular RNAs (circRNAs). Although circRNAs were studied in many fields, little research was devoted to IMF in sheep. We used the longissimus dorsi muscle of 2 and 12-month-old AWFS as research material to identify circRNAs related to IMF deposition in these sheep by RNA-seq screening for differentially expressed circRNAs in the two age groups. A total of 11,565 candidate circRNAs were identified, of which the 104 differentially expressed circRNAs in the two age groups were analyzed. Enrichment analysis was performed using Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes. The enriched pathways included lipid transport (GO:0006869), negative regulation of canonical Wnt signaling pathway (GO:0090090), fat digestion and absorption (ko04975), and sphingolipid metabolism (ko00600). The differentially expressed circRNAs included ciRNA455, circRNA9086, circRNA7445, circRNA4557, and others. The source genes involved in these pathways might regulate IMF deposition. We used the TargetScan and miRanda software for interaction analysis, and a network diagram of circRNA-miRNA interactions was created. CircRNA455-miR-127, circRNA455-miR-29a, circRNA455-miR-103, circRNA4557-mir149-5p, and circRNA2440-mir-23a might be involved in the IMF deposition process. The targeting relationship of circRNA4557-miR-149-5p was verified by a dual-luciferase reporter assay. The RT-qPCR results of seven randomly selected circRNAs were consistent with the sequencing results. This study provides additional information on circRNA regulation of IMF deposition in AFWS and is a useful resource for future research on this sheep breed.

Dynamic Expression and Regulatory Network of Circular RNA for Abdominal Preadipocytes Differentiation in Chicken (Gallus gallus)

Circular RNA (circRNA), as a novel endogenous biomolecule, has been emergingly demonstrated to play crucial roles in mammalian lipid metabolism and obesity. However, little is known about their genome-wide identification, expression profile, and function in chicken adipogenesis. In present study, the adipogenic differentiation of chicken abdominal preadipocyte was successfully induced, and the regulatory functional circRNAs in chicken adipogenesis were identified from abdominal adipocytes at different differentiation stages using Ribo-Zero RNA-seq. A total of 1,068 circRNA candidates were identified and mostly derived from exons. Of these, 111 differentially expressed circRNAs (DE-circRNAs) were detected, characterized by stage-specific expression, and enriched in several lipid-related pathways, such as Hippo signaling pathway, mTOR signaling pathway. Through weighted gene co-expression network analyses (WGCNA) and K-means clustering analyses, two DE-circRNAs, Z:35565770|35568133 and Z:54674624|54755962, were identified as candidate regulatory circRNAs in chicken adipogenic differentiation. Z:35565770|35568133 might compete splicing with its parental gene, ABHD17B, owing to its strictly negative co-expression. We also constructed competing endogenous RNA (ceRNA) network based on DE-circRNA, DE-miRNA, DE-mRNAs, revealing that Z:54674624|54755962 might function as a ceRNA to regulate chicken adipogenic differentiation through the gga-miR-1635-AHR2/IRF1/MGAT3/ABCA1/AADAC and/or the novel_miR_232-STAT5A axis. Translation activity analysis showed that Z:35565770|35568133 and Z:54674624|54755962 have no protein-coding potential. These findings provide valuable evidence for a better understanding of the specific functions and molecular mechanisms of circRNAs underlying avian adipogenesis.