Comparative Genome-Wide Alternative Splicing Analysis of Longissimus Dorsi Muscles Between Japanese Black (Wagyu) and Chinese Red Steppes Cattle

Copy number variations within fibroblast growth factor 13 gene influence growth traits and alternative splicing in cattle

Previous researches revealed a copy number variation (CNV) region in the bovine fibroblast growth factor 13 (FGF13) gene. However, its effects remain unknown. This study detected the various copy number types in seven Chinese cattle breeds and analysed their population genetic characteristics and effects on growth traits and transcription levels. Copy number Loss was more frequent in Caoyuan Red cattle and Xianan cattle than in the other breeds. Association analysis between CNV and growth traits of Qinchuan indicated that the CNV was significantly related to chest depth, hip width and hucklebone width (P < 0.05). Additionally, the growth traits of individuals with copy number Loss were significantly inferior to those with copy number Gain or Median (P < 0.05). Besides, we found two splicing isoforms, AS1 and AS2, in FGF13 gene, which resulted from alternative 5' splicing sites of intron 1. These isoforms showed varied expression levels in various tissues. Moreover, CNV was significantly and negatively associated with the mRNA expression of AS1 (r = -0.525, P < 0.05). The CNVs in bovine FGF13 gene negatively regulated growth traits and gene transcription. These observations provide new insights into bovine FGF13 gene, delivering potentially useful information for future Chinese cattle breeding programs.

Comparative Genome-Wide Alternative Splicing Analysis between Preadipocytes and Adipocytes

Alternative splicing is a ubiquitous regulatory mechanism in gene expression that allows a single gene to generate multiple messenger RNAs (mRNAs). Adipocyte development is regulated by many processes, and recent studies have found that splicing factors also play an important role in adipogenic development. In the present study, we further investigated the differences in selective shearing during different periods of adipocyte differentiation. We identified five alternative splicing types including skipped exon, mutually exclusive exon, Alternative 5' splice site, Alternative 3' splice site, and Retained intron, with skipped exons being the most abundant type of selective shearing. In total, 641 differentially expressed selective shearing genes were obtained, enriched in 279 pathways, from which we selected and verified the accuracy of the sequencing results. Overall, RNA-seq revealed changes in the splicing and expression levels of these new candidate genes between precursor adipocytes and adipocytes, suggesting that they may be involved in adipocyte generation and differentiation.

Deciphering the role of alternative splicing as a potential regulator in fat-tail development of sheep: a comprehensive RNA-seq based study

Although research on alternative splicing (AS) has been widely conducted in mammals, no study has investigated the splicing profiles of genes involved in fat-tail formation in sheep. Here, for the first time, a comprehensive study was designed to investigate the profile of AS events and their involvement in fat-tail development of sheep. In total, 45 RNA-Seq samples related to seven different studies, which have compared the fat-tailed vs thin-tailed sheep breeds, were analyzed. Two independent tools, rMATS and Whippet, along with a set of stringent filters were applied to identify differential AS (DAS) events between the breeds per each study. Only DAS events that were detected by both tools as well as in at least three datasets with the same ΔPSI trend (percent spliced in), were considered as the final high-confidence set of DAS genes. Final results revealed 130 DAS skipped exon events (69 negative and 61 positive ΔPSI) belonged to 124 genes. Functional enrichment analysis highlighted the importance of the genes in the underlying molecular mechanisms of fat metabolism. Moreover, protein-protein interaction network analysis revealed that DAS genes are significantly connected. Of DAS genes, five transcription factors were found that were enriched in the biological process associated with lipid metabolism like "Fat Cell Differentiation". Further investigations of the findings along with a comprehensive literature review provided a reliable list of candidate genes that may potentially contribute to fat-tail formation including HSD11B1, SIRT2, STRN3 and TCF7L2. Based on the results, it can be stated that the AS patterns may have evolved, during the evolution of sheep breeds, as another layer of regulation to contribute to biological complexity by reprogramming the gene regulatory networks. This study provided the theoretical basis of the molecular mechanisms behind the sheep fat-tail development in terms of AS.

Detection of genomic variations and selection signatures in Wagyu using whole-genome sequencing data

Wagyu is recognized for producing marbled beef with high nutritional value and flavor. Reportedly, Wagyu has been widely used to improve the meat quality of local breeds around the world. However, studies on the genetic mechanism of meat quality in Wagyu at the whole-genome level are rarely reported. Here, whole-genome sequencing data of 11 Wagyu and 115 other individuals were used to explore the genomic variations and genes under selection pressure in Wagyu. A total of 31 349 non-synonymous variants and 53 102 synonymous variants were identified in Wagyu. The population structure analysis showed that Wagyu had the closest genetic relationship with Mishima-Ushi cattle and was apparently separated from other cattle breeds. Then, composite likelihood ratio (CLR), integrated haplotype score, fixation index and cross-population composite likelihood ratio (XP-CLR) tests were performed to identify the candidate genes under positive selection in Wagyu. In total, 770 regions containing 312 genes were identified by at least three methods. Among them, 97 regions containing 27 genes were detected by all four methods. We specifically illustrate a list of interesting genes, including LRP2BP, GAA, CACNG6, CXADR, GPCPD1, KLF2, KLF13, SOX5, MYBPC1, SLC25A10, ATP8A1 and MYH15, which are associated with lipid metabolism, fat deposition, muscle development, bone development, feed intake and growth traits in Wagyu. This is the first study to explore the genomic variations and selection signatures of Wagyu at the whole-genome level. These results will provide significant help to beef cattle improvement and breeding.

Analysis of genetic diversity and selection signals in Chaling cattle of southern China using whole-genome scan

China has diversified resources of indigenous cattle, which are classified into Northern, Central, and Southern groups according to their geographical distribution. Chaling cattle belong to Southern group. This breed is famous for the production of good quality meat with elite meat grades. To analyze the genetic diversity of Chaling cattle, 20 samples were sequenced using whole-genome resequencing technology, along with 138 published whole-genome sequencing data of Indian indicine cattle, Chinese indicine cattle, East Asian taurine cattle, Eurasian taurine cattle, and European taurine cattle as control. It was found that Chaling cattle originated from Chinese indicine cattle. The genetic diversity of Chaling cattle is higher than that of Indian indicine cattle, East Asian taurine cattle, Eurasian taurine cattle, and European taurine cattle, but lower than that of Chinese indicine cattle and Xiangxi cattle. Annotating the selection signals obtained by composite likelihood ratio, θπ, F_ST , π-ratio, and XP-EHH methods, several genes associated with immunity, heat tolerance, reproduction, growth, and meat quality showed strong selection signals. In general, this study provides a theoretical basis for analyzing the genetic mechanism of Chaling cattle with excellent adaptability, rough feeding tolerance, good immune performance, and good meat quality. This work lays a foundation for genetic breeding of Chaling cattle in future.

Beef Nutritional Characteristics, Fat Profile and Blood Metabolic Markers from Purebred Wagyu, Crossbred Wagyu and Crossbred European Steers Raised on a Fattening Farm in Spain

A high intramuscular fat content characterizes Wagyu (WY) cattle breed. Our objective was to compare beef from WY, WY-by-Angus, or Wangus (WN) steers with European, Angus-by-Charolais-Limousine crossbred steers (ACL), considering metabolic biomarkers pre-slaughtering and nutritional characteristics, including health-related indexes of the lipid fraction. The fattening system with olein-rich diets and no exercise restriction included 82 steers, 24 WY, 29 WN, and 29 ACL. The slaughter ages and weights were (median and interquartile range) 38.4 mo.-old (34.9-40.3 mo.) and 840 kg (785-895 kg) for WY; for WN, 30.6 mo. (26.9-36.5 mo.) and 832 kg (802-875 kg), and for ACL steers, 20.3 mo.-old (19.0-22.7 mo.) and 780 kg (715-852 kg). Blood lipid-related metabolites, except for non-esterified fatty acids (NEFA) and low-density level cholesterol (LDL), were higher in WY and WN than in ACL, while glucose was lower in WY and WN. Leptin was higher in WN than in ACL. Pre-slaughtering values of plasma HDL underscored as a possible metabolic biomarker directly related to beef quality. The amino-acid content in beef did not differ among experimental groups, except for more crude protein in ACL. Compared to ACL, WY steers showed higher intramuscular fat in sirloin (51.5 vs. 21.9%) and entrecote (59.6 vs. 27.6%), more unsaturated fatty acids in entrecote (55.8 vs. 53.0%), and more oleic acid in sirloin (46 vs. 41.3%) and entrecote (47.5 vs. 43.3%). Compared to ACL entrecote, WY and WN showed better atherogenic (0.6 and 0.55 vs. 0.69), thrombogenicity (0.82 and 0.92 vs. 1.1), and hypocholesterolemic/hypercholesterolemic index (1.9 and 2.1 vs. 1.7). Therefore, beef's nutritional characteristics depend on breed/crossbred, slaughtering age and cut, with WY and WN entrecote samples showing a healthier lipid fraction.

Alternative splicing signature of alveolar type II epithelial cells of Tibetan pigs under hypoxia-induced

Alternative splicing (AS) allows the generation of multiple transcript variants from a single gene and affects biological processes by generating protein diversity in organisms. In total, 41,642 AS events corresponding to 9,924 genes were identified, and SE is the most abundant alternatively spliced type. The analysis of functional categories demonstrates that alternatively spliced differentially expressed genes (DEGs) were enriched in the MAPK signaling pathway and hypoxia-inducible factor 1 (HIF-1) signaling pathway. Proteoglycans in cancer between the normoxic (21% O2, TN and LN) and hypoxic (2% O2, TL and LL) groups, such as SLC2A1, HK1, HK2, ENO3, and PFKFB3, have the potential to rapidly proliferate alveolar type II epithelial (ATII) cells by increasing the intracellular levels of glucose and quickly divert to anabolic pathways by glycolysis intermediates under hypoxia. ACADL, EHHADH, and CPT1A undergo one or two AS types with different frequencies in ATII cells between TN and TL groups (excluding alternatively spliced DEGs shared between normoxic and hypoxic groups), and a constant supply of lipids might be obtained either from the circulation or de novo synthesis for better growth of ATII cells under hypoxia condition. MCM7 and MCM3 undergo different AS types between LN and LL groups (excluding alternatively spliced DEGs shared between normoxic and hypoxic groups), which may bind to the amino-terminal PER-SIM-ARNT domain and the carboxyl terminus of HIF-1α to maintain their stability. Overall, AS and expression levels of candidate mRNAs between Tibetan pigs and Landrace pigs revealed by RNA-seq suggest their potential involvement in the ATII cells grown under hypoxia conditions.

Integrative Analysis of Nanopore and Illumina Sequencing Reveals Alternative Splicing Complexity in Pig Longissimus Dorsi Muscle

Alternative splicing (AS) is a key step in the post-transcriptional regulation of gene expression that can affect intramuscular fat (IMF). In this study, longissimus dorsi muscles from 30 pigs in high- and low- IMF groups were used to perform Oxford Nanopore Technologies (ONT) full-length sequencing and Illumina strand-specific RNA-seq. A total of 43,688 full-length transcripts were identified, with 4,322 novel genes and 30,795 novel transcripts. Using AStalavista, a total of 14,728 AS events were detected in the longissimus dorsi muscle. About 17.79% of the genes produced splicing isoforms, in which exon skipping was the most frequent AS event. By analyzing the expression differences of mRNAs and splicing isoforms, we found that differentially expressed mRNAs with splicing isoforms could participate in skeletal muscle development and fatty acid metabolism, which might determine muscle-related traits. SERBP1, MYL1, TNNT3, and TNNT1 were identified with multiple splicing isoforms, with significant differences in expression. AS events occurring in IFI6 and GADD45G may cause significant differences in gene expression. Other AS events, such as ONT.15153.3, may regulate the function of ART1 by regulating the expression of different transcripts. Moreover, co-expression and protein-protein interaction (PPI) analysis indicated that several genes (MRPL27, AAR2, PYGM, PSMD4, SCNM1, and HNRNPDL) may be related to intramuscular fat. The splicing isoforms investigated in our research provide a reference for the study of alternative splicing regulation of intramuscular fat deposition.

A novel 27-bp indel in the intron region of the YBX3 gene is associated with growth traits in chickens

1. The DNA/RNA binding protein YBX3 is associated with gene transcription, DNA repair, and the progression of various diseases and is highly conserved from bacteria to humans.2. The following experiment found a 27-bp insertion/deletion polymorphism in the intron region of the YBX3 gene through resequencing. In cross-designed, F₂ resource groups, the indel was significantly associated with broiler weight and body size at 0, 2, 4, 6, 8, 10 and 12 weeks of age and several other traits (semi evisceration weight (SEW), evisceration weight (EW), semi evisceration rate (SER), evisceration rate (ER), head weight (HW), claw weight (CLW), wing weight (DWW), gizzard weight (GW), pancreas weight (PW), chest muscle weight (CMW), leg weight (LW), leg muscle weight (LMW), shedding weight (SW), carcass weight (CW) and pectoral area (PA)) (P<0.05).3. The insertion-insertion (II) genotype was significantly associated with the greatest growth traits and carcass traits, whereas the values associated with the insertion-deletion (ID) genotype were the lowest in the F₂ reciprocal cross chickens.4. The mutation sites were genotyped in 3611 individuals from 13 different chicken breeds and cross-designed F₂ resource groups. The II genotype is the most important in commercial broilers, and the I allele frequency observed in these breeds was relatively high. However, there is still considerable potential in breeding dual-purpose chickens and commercial laying hens.5. The mRNA expression of the YBX3 gene in tissues from different breeds and developmental stages demonstrated that the 27-bp indel may affect the entire development process of poultry by affecting muscle development. These findings are beneficial for elucidating the function of the YBX3 gene and facilitating enhanced reproduction in the chicken industry.

The SLC27A1 Gene and Its Enriched PPAR Pathway Are Involved in the Regulation of Flavor Compound Hexanal Content in Chinese Native Chickens

The role of hexanal in flavor as an indicator of the degree of oxidation of meat products is undeniable. However, the genes and pathways of hexanal formation have not been characterized in detail. In this study, we performed differential gene expression analysis and weighted gene co-expression network analysis (WGCNA) on groups of Tiannong partridge chickens with different relative hexanal content in order to find the genes involved in the formation of hexanal and the specific pathways of hexanal formation. Then we confirmed the relationship of these candidate genes with hexanal using Jingxing Yellow chicken and Wenchang chicken. In this study, WGCNA revealed a module of co-expressed genes that were highly associated with the volatile organic compound hexanal. We also compared transcriptome gene expression data of samples from chicken groups with high and low relative contents of hexanal and identified a total of 651 differentially expressed genes (DEGs). Among them, 356 genes were up regulated, and 295 genes were downregulated. The different biological functions associated with the DEGs, hub genes and hexanal were identified by functional analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations. Among all the hub genes in the significant module identified by WGCNA, more were enriched in the PPAR signaling pathway, the proteasome pathway, etc. Additionally, we found that DEGs and hub genes, including SLC27A1, ACOX3, NR4A1, VEGFA, JUN, EGR1, CACNB1, GADD45A and DUSP1, were co-enriched in the peroxisome proliferator-activated receptor (PPAR) signaling pathway, p53 signaling pathway and mitogen-activated protein kinases (MAPK) signaling pathway, etc. Transcriptome results of the Jingxing Yellow chicken population showed that the SLC27A1 gene was significantly associated with hexanal and enriched in the PPAR pathway. Our study provides a comprehensive insight into the key genes related to hexanal content, and can be further explored by functional and molecular studies.