Comparative transcriptome profiling of longissimus muscle tissues from Qianhua Mutton Merino and Small Tail Han sheep

Integrative Meta-Analysis: Unveiling Genetic Factors in Meat Sheep Growth and Muscular Development through QTL and Transcriptome Studies

OBJECTIVE

The study aimed to investigate the effects of castration on performance, carcass characteristics, and meat quality in sheep, as well as explore the expression of key genes related to metabolic pathways and muscle growth following castration.

METHODS

A meta-analysis approach was utilized to analyze data from multiple studies to compare the performance, carcass characteristics, and meat quality of castrated sheep (wethers) with intact rams. Additionally, protein-protein interaction (PPI) networks, differential gene expression (DEG) interactions, Gene Ontology (GO) terms, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were examined to identify molecular mechanisms associated with fat metabolism and muscle development in sheep tails.

RESULTS

The analysis revealed that castrated sheep (wethers) exhibited improved average daily gain, increased tenderness, lower backfat thickness, and a tendency for greater loin muscle area compared to intact rams. This suggests that castration promotes faster growth and results in leaner carcasses with potentially higher muscle content. Furthermore, the identification of downregulated DEGs like ACLY, SLC27A2, and COL1A1 and upregulated DEGs such as HOXA9, PGM2L1, and ABAT provides insights into the molecular mechanisms underlying fat deposition and muscle development in sheep.

CONCLUSIONS

The findings support the practice of castration in sheep production as it enhances growth performance, leads to leaner carcasses with higher muscle content, and improves meat tenderness. The identified changes in gene expression offer valuable insights for further research into understanding the impact of castration on muscle development and fat metabolism in sheep. This meta-analysis contributes to the knowledge of molecular mechanisms involved in fat deposition in sheep, opening avenues for future investigations in livestock fat metabolism research.

The Integration of Genome-Wide DNA Methylation and Transcriptomics Identifies the Potential Genes That Regulate the Development of Skeletal Muscles in Ducks

DNA methylation is a pivotal epigenetic regulatory mechanism in the development of skeletal muscles. Nonetheless, the regulators responsible for DNA methylation in the development of embryonic duck skeletal muscles remain unknown. In the present study, whole genome bisulfite sequencing (WGBS) and transcriptome sequencing were conducted on the skeletal muscles of embryonic day 21 (E21) and day 28 (E28) ducks. The DNA methylation pattern was found to fall mainly within the cytosine-guanine (CG) context, with high methylation levels in the intron, exon, and promoter regions. Overall, 7902 differentially methylated regions (DMRs) were identified, which corresponded to 3174 differentially methylated genes (DMGs). By using integrative analysis of both WGBS with transcriptomics, we identified 1072 genes that are DMGs that are negatively associated with differentially expressed genes (DEGs). The gene ontology (GO) analysis revealed significant enrichment in phosphorylation, kinase activity, phosphotransferase activity, alcohol-based receptors, and binding to cytoskeletal proteins. The Kyoto Encyclopedia of Genes and Genomes (KEGGs) analysis showed significant enrichment in MAPK signaling, Wnt signaling, apelin signaling, insulin signaling, and FoxO signaling. The screening of enriched genes showed that hyper-methylation inhibited the expression of Idh3a, Got1, Bcl2, Mylk2, Klf2, Erbin, and Klhl38, and hypo-methylation stimulated the expression of Col22a1, Dnmt3b, Fn1, E2f1, Rprm, and Wfikkn1. Further predictions showed that the CpG islands in the promoters of Klhl38, Klf2, Erbin, Mylk2, and Got1 may play a crucial role in regulating the development of skeletal muscles. This study provides new insights into the epigenetic regulation of the development of duck skeletal muscles.

Global DNA Methylation, miRNA, and mRNA Profiles in Sheep Skeletal Muscle Promoted by Hybridization

With the development of high-throughput sequencing technology, several nongenetic variations, including noncoding RNAs such as miRNAs, and DNA methylation, have been found to play an important role in animal muscle development and fat metabolism. In this study, Southdown and Suffolk were selected as male parents for hybridization with Hu sheep (Southdown × Hu (NH), Suffolk × Hu (SH), and Hu × Hu (HH)). RNA sequencing, bisulfite sequencing, and small-RNA sequencing were used to study the methylation patterns and differences in miRNA and mRNA expression in the F1 sheep longissimus dorsi muscle tissue. We identified 765 differentially expressed genes (DEGs), 10,161 differentially methylated regions (DMRs), and 164 differentially expressed miRNAs, which were significantly enriched in AMPK signaling, fatty acid degradation, metabolism, and other related pathways (P < 0.05). In addition, we constructed a DNA methylation-mRNA and miRNA-mRNA coexpression network. A total of 42 common genes were identified from DMRs and DEGs. Importantly, we predicted that 33 differentially expressed miRNAs directly or indirectly targeted the SLC27A6. The data obtained in this study provide useful information and evidence to support further understanding of the miRNA and DNA methylation of key genes regulating muscle growth and fat metabolism in hybrid sheep populations.

Dynamic Changes in the Global Transcriptome of Postnatal Skeletal Muscle in Different Sheep

Sheep growth performance, mainly skeletal muscle growth, provides direct economic benefits to the animal husbandry industry. However, the underlying genetic mechanisms of different breeds remain unclear. We found that the cross-sectional area (CSA) of skeletal muscle in Dorper (D) and binary cross-breeding (HD) was higher than that in Hu sheep (H) from 3 months to 12 months after birth. The transcriptomic analysis of 42 quadriceps femoris samples showed that a total of 5053 differential expression genes (DEGs) were identified. The differences in the global gene expression patterns, the dynamic transcriptome of skeletal muscle development, and the transcriptome of the transformation of fast and slow muscles were explored using weighted correlation network analysis (WGCNA) and allele-specific expression analysis. Moreover, the gene expression patterns of HD were more similar to D rather than H from 3 months to 12 months, which might be the reason for the difference in muscle growth in the three breeds. Additionally, several genes (GNB2L1, RPL15, DVL1, FBXO31, etc.) were identified as candidates related to skeletal muscle growth. These results should serve as an important resource revealing the molecular basis of muscle growth and development in sheep.

Integration of Transcriptomics and Non-Targeted Metabolomics Reveals the Underlying Mechanism of Skeletal Muscle Development in Duck during Embryonic Stage

Skeletal muscle is an important economic trait in duck breeding; however, little is known about the molecular mechanisms of its embryonic development. Here, the transcriptomes and metabolomes of breast muscle of Pekin duck from 15 (E15_BM), 21 (E21_BM), and 27 (E27_BM) days of incubation were compared and analyzed. The metabolome results showed that the differentially accumulated metabolites (DAMs), including the up-regulated metabolites, l-glutamic acid, n-acetyl-1-aspartylglutamic acid, l-2-aminoadipic acid, 3-hydroxybutyric acid, bilirubin, and the significantly down-regulated metabolites, palmitic acid, 4-guanidinobutanoate, myristic acid, 3-dehydroxycarnitine, and s-adenosylmethioninamine, were mainly enriched in metabolic pathways, biosynthesis of secondary metabolites, biosynthesis of cofactors, protein digestion and absorption, and histidine metabolism, suggesting that these pathways may play important roles in the muscle development of duck during the embryonic stage. Moreover, a total of 2142 (1552 up-regulated and 590 down-regulated), 4873 (3810 up-regulated and 1063 down-regulated), and 2401 (1606 up-regulated and 795 down-regulated) DEGs were identified from E15_BM vs. E21_BM, E15_BM vs. E27_BM and E21_BM vs. E27_BM in the transcriptome, respectively. The significantly enriched GO terms from biological processes were positive regulation of cell proliferation, regulation of cell cycle, actin filament organization, and regulation of actin cytoskeleton organization, which were associated with muscle or cell growth and development. Seven significant pathways, highly enriched by FYN, PTK2, PXN, CRK, CRKL, PAK, RHOA, ROCK, INSR, PDPK1, and ARHGEF, were focal adhesion, regulation of actin cytoskeleton, wnt signaling pathway, insulin signaling pathway, extracellular matrix (ECM)-receptor interaction, cell cycle, and adherens junction, which participated in regulating the development of skeletal muscle in Pekin duck during the embryonic stage. KEGG pathway analysis of the integrated transcriptome and metabolome indicated that the pathways, including arginine and proline metabolism, protein digestion and absorption, and histidine metabolism, were involved in regulating skeletal muscle development in embryonic Pekin duck. These findings suggested that the candidate genes and metabolites involved in crucial biological pathways may regulate muscle development in the Pekin duck at the embryonic stage, and increased our understanding of the molecular mechanisms underlying the avian muscle development.

Research on meat quality of Qianhua Mutton Merino sheep and Small-tail Han sheep

In this study, we analyzed the meat quality of Qianhua Mutton Merino sheep (QHMM). QHMM and Small-tail Han sheep (STH) were selected as the test animals. After slaughtering, the longissimus dorsi muscles were tested to determine the nutrient composition, content of amino acid, fatty acid, and muscle fiber diameter. According to the nutrient test result, the protein content of QHMM was higher than that of STH. However, the fat content was significantly lower (P < 0.05) than STH. Additionally, total amino acid content, total essential amino acid content, total half-essential amino acid content, total non-essential amino acid content, and flavor amino acid content of QHMM were significantly higher (P < 0.05) than those of STH. Moreover, the contents of C22:1 n9, C18:2 n6, C18:3 n6, C22:6 n6, and C10:0 of the muscle in QHMM were significantly higher (P < 0.05) than those of STH. Furthermore, the essential amino acid score, the total unsaturated-to-saturated fatty acid ratio, and the polyunsaturated-to-saturated fatty acid ratio of QHMM were greater than those of STH. Additionally, the muscle fiber diameter of QHMM was considerably lower (P < 0.01) than that of STH. In conclusion, the meat quality of QHMM was better than that of STH.

Co-expression of candidate genes regulating growth performance and carcass traits of Barki lambs in Egypt

Sheep are considered one of the main sources of animal protein in Egypt and the producers of sheep mutton eagers to find biological criteria for selecting fast-growing lambs that reach market weight early. Therefore, the present study aimed to find a link between the expression profile of selected candidate genes with growth performance and carcass traits of Barki lambs. Thirty-eight Barki lambs were kept and fed individually after weaning till 12 months of age and were divided into 3 groups according to growth performance (fast, intermediate, and slow-growing). Three samples were taken from different body tissues (eye muscle, liver, and fat tail) of each group, directly during slaughtering and stored at − 80 °C until RNA isolation. Real-time PCR was used to profile selected candidate genes (RPL7, CTP1, FABP4, ADIPOQ, and CAPN3) and GAPDH was used as a housekeeping gene. The results indicated that the final body weight was significantly (P ≤ 0.05) greater in the fast (49.9 kg) and intermediate (40.7 kg) compared to slow-growing animals (30.8 kg). The hot carcass weight was heavier (P ≤ 0.05) in the fast and intermediate-growing (24.57 and 19.07 kg) than slow-growing lambs (15.10 kg). The blood profiles of T3 and T4 hormones in addition to other parameters such as total protein, total lipids, and calcium level showed no clear variations among different experimental groups. At the molecular level, our data demonstrated upregulation of genes involved in protein biosynthesis (RPL7), fatty acid oxidation (CPT1), and lipolysis (FABP4) in the fast and intermediate-growing lambs in all studied tissues which facilitate protein accretion, energy expenditure, and fatty acid partitioning required for muscle building up. Moreover, the expression profile of the gene involved in muscle development (CAPN3) was increased in fast and intermediate-growing compared to slow-growing lambs in order to support muscle proper development. On the other hand, a candidate gene involved in lipogenesis (ADIPOQ) was expressed similarly in fat and liver tissues; however, its expression was increased in muscles of fast and intermediate-growing lambs compared to slow-growing animals. In conclusion, the current study indicated that the expression profile of genes involved in metabolic activities of liver, muscle, and adipose tissue is linked with the growth performance of lambs although no variations were detected in blood parameters. This provides an evidence for the importance of co-expression of these genes in body tissues to determine the final body weight and carcass characteristics of Barki sheep.

Comparative Transcriptome Analysis Reveals the Mechanism Associated With Dynamic Changes in Meat Quality of the Longissimus Thoracis Muscle in Tibetan Sheep at Different Growth Stages

Tibetan sheep are mainly distributed in the Qinghai–Tibet Plateau. Its meat is not only essential for the local people but also preferred by the non-inhabitant of this plateau also. To investigate the salient development features and molecular mechanism of the meat difference of LT muscle caused by different growth stages in Tibetan sheep, the carcass performance, meat quality, and comparative transcriptome analysis were performed for investigating the potential molecular mechanism of the meat quality difference of the LT muscle caused by four growth stages [4-months old (4 months), 1.5-years old (1.5 years), 3.5-years old (3.5 years), and 6-years old (6 years)] in the Tibetan sheep. The shear force increased with the increase of age (p < 0.05) while the intramuscular fat (IMF) was the highest at 1.5 y. The AMPK signaling pathway was significantly enriched in the four comparative groups. The weighted gene co-expression network analysis (WGCNA) results showed that the hub genes P4HA2, FBXL4, and PPARA were identified to regulate the meat quality. In summary, 1.5 years was found to be the most suitable slaughter age of the Tibetan sheep which ensured better meat tenderness and higher IMF content. Moreover, the genes LIPE, LEP, ADIPOQ, SCD, and FASN may regulate the transformation of the muscle fiber types through the AMPK signaling pathway, further affecting the meat quality.

VMP1 Regulated by chi-miR-124a Effects Goat Myoblast Proliferation, Autophagy, and Apoptosis through the PI3K/ULK1/mTOR Signaling Pathway

The production of goat meat is determined by the growth speed of muscle fibers, and the autophagy and apoptosis of myoblast cells is a crucial process in the growth of muscle fibers. The rapid growth of muscle fibers occurs from one month old to nine months old in goats; however, the mechanisms of myoblast cells’ autophagy and apoptosis in this process are still unknown. To identify candidate genes and signaling pathway mechanisms involved in myoblast apoptosis and autophagy, we compared the expression characteristics of longissimus dorsi tissues from Wu’an goats—a native goat breed of China—at 1 month old (mon1 group) and 9 months old (mon9 group). Herein, a total of 182 differentially expressed mRNAs (DEGs) in the mon1 vs. mon9 comparison, along with the KEGG enrichments, showed that the PI3K-Akt pathway associated with autophagy and apoptosis was significantly enriched. Among these DEGs, expression of vacuole membrane protein 1 (VMP1)—a key gene for the PI3K-Akt pathway—was significantly upregulated in the older goats relative to the 1-month-old goats. We demonstrated that VMP1 promotes the proliferation and autophagy of myoblasts, and inhibits their apoptosis. The integration analysis of miRNA–mRNA showed that miR-124a was a regulator of VMP1 in muscle tissue, and overexpression and inhibition of miR-124a suppressed the proliferation and autophagy of myoblasts. The PI3K/Akt/mTOR pathway was an important pathway for cell autophagy. Additionally, the activator of the PI3K/Akt/mTOR pathway, the expression of VMP1, and ULK1 were higher than the negative control, and the expression of mTOR was depressed. The expression of VMP1, ULK1, and mTOR was the opposite when the inhibitor was added to the myoblasts. These results show that the PI3K/Akt/mTOR pathway promoted the expression of VMP1 and ULK1. By using adenovirus-mediated apoptosis and proliferation assays, we found that that miR-124a inhibits myoblast proliferation and autophagy, and promotes their apoptosis by targeting VMP1. In conclusion, our results indicated that VMP1 was highly expressed in the LD muscle tissues of nine-month-old goats, and that it was regulated by miR-124a to inhibit myoblast cells’ apoptosis through the PI3K/Akt/mTOR pathway, and to promote proliferation and autophagy. These findings contribute to the understanding of the molecular mechanisms involved in myoblast proliferation, autophagy, and apoptosis.

Profiling and Functional Analysis of mRNAs during Skeletal Muscle Differentiation in Goats

Skeletal myogenesis is a complicated biological event that involves a succession of tightly controlled gene expressions. In order to identify novel regulators of this process, we performed mRNA-Seq studies of goat skeletal muscle satellite cells (MuSCs) cultured under proliferation (GM) and differentiation (DM1/DM5) conditions. A total of 19,871 goat genes were expressed during these stages, 198 of which represented novel transcripts. Notably, in pairwise comparisons at the different stages, 2551 differentially expressed genes (DEGs) were identified (p < 0.05), including 1560 in GM vs. DM1, 1597 in GM vs. DM5, and 959 in DM1 vs. DM5 DEGs. The time-series expression profile analysis clustered the DEGs into eight gene groups, three of which had significantly upregulated and downregulated patterns (p < 0.05). Functional enrichment analysis showed that DEGs were enriched for essential biological processes such as muscle structure development, muscle contraction, muscle cell development, striated muscle cell differentiation, and myofibril assembly, and were involved in pathways such as the MAPK, Wnt and PPAR signaling pathways. Moreover, the expression of eight DEGs (MYL2, DES, MYOG, FAP, PLK2, ADAM, WWC1, and PRDX1) was validated. These findings offer novel insights into the transcriptional regulation of skeletal myogenesis in goats.

Transcriptome Analysis of the Liver and Muscle Tissues of Dorper and Small-Tailed Han Sheep

It is well known that Dorper (DP) is a full-bodied, fast-growing and high dressing percentage breed, while the production performance of Small-tailed Han sheep (STH) is not so excellent, in contrast to DP. Therefore, in this study, a comparative transcriptomic analysis of liver and muscle tissues from DP and STH breeds was carried out to find differentially expressed genes (DEGs) that affect their growth and meat quality traits. The results showed that the total number of DEGs was 2,188 in the two tissues. There were 950, 160 up-regulated and 1,007, 71 down-regulated genes in the liver and muscle, respectively. Several DEGs such as TGFB1, TGFB3, FABP3, LPL may be associated with growth and development in DP. Also, several GO terms were found to be associated with muscle growth and development, such as developmental growth (GO:0048589), and myofibril (GO:0030016). Further validation of eight genes (6 up-regulated, and 2 down-regulated) was performed using quantitative RT-PCR. These findings will provide valuable information for studying growth and development as well as meat quality traits in sheep.

Screening of Genes Related to Growth, Development and Meat Quality of Sahan Crossbred F1 Sheep Based on RNA-Seq Technology

This study aimed to identify genes related to sheep growth, development and meat quality. Small-tailed Han sheep (STH), and small-tailed Han sheep and Suffolk crossbred F1 (STH×SFK), were selected to determine the growth performance, slaughter performance, and meat quality. The longissimus dorsi muscle was selected for transcriptome sequencing, and the target gene was screened based on bioinformatics analysis; real-time fluorescent quantitative PCR (RT-PCR) and western blotting (WB) were conducted to verify the target gene. Locations of genes in tissues were confirmed via immunofluorescence. The results showed that the pre-slaughter live weight, bust circumference, slaughter performance, and marbling score of the STH×SFK population were significantly higher than those of the STH population (P < 0.01). Sequencing results showed that 560 differentially expressed genes (DEGs) were identified in the STH×SFK population, of which 377 exhibited up-regulated and 183 exhibited down-regulated expression levels. GO annotation revealed that DEGs could be classified into 13 cell components, 10 molecular functions, and 22 biological processes. The KEGG enrichment analysis showed that DEGs were mainly enriched in the Rap1 signaling pathway, Ras signaling pathway, and other pathways related to growth and meat quality. Based on the GO and KEGG analyses, four candidate genes related to sheep growth and meat quality, namely myostain (MSTN), interferon-related developmental regulator 1 (IFRD1), peroxisome proliferator activator receptor delta (PPARD), and myosin light chain 2 (MLC2 or MYL2), were screened. The expression levels of genes and proteins were verified via RT-PCR and WB, and the results were consistent with the trend of transcriptome sequencing. Immunofluorescence results showed that IFRD1 was expressed in the cytoplasm and nucleus, and MYL2 was expressed in the cytoplasm. This study revealed the mechanism of gene regulation of sheep growth and development at the molecular level and provided a theoretical basis for studying sheep genetics and breeding.

Circular RNA, microRNA and Protein Profiles of the Longissimus Dorsi of Germany ZIKA and Sichuan White Rabbits

Due to the dietetic properties and remarkable nutritive value of rabbit meat, its industry is increasing rapidly. However, the association between circular RNAs, microRNAs, and proteins and muscle fiber type, and meat quality of rabbit is still unknown. Here, using deep sequencing and iTRAQ proteomics technologies we first identified 3159 circRNAs, 356 miRNAs, and 755 proteins in the longissimus dorsi tissues from Sichuan white (SCWrabs) and Germany great line ZIKA rabbits (ZIKArabs). Next, we identified 267 circRNAs, 3 miRNAs, and 29 proteins differentially expressed in the muscle tissues of SCWrabs and ZIKArabs. Interaction network analysis revealed some key regulation relationships between noncoding RNAs and proteins that might be associated with the muscle fiber type and meat quality of rabbit. Further, miRNA isoforms and gene variants identified in SCWrabs and ZIKArabs revealed some pathways and biological processes related to the muscle development. This is the first study of noncoding RNA and protein profiles for the two rabbit breeds. It provides a valuable resource for future studies in rabbits and will improve our understanding of the molecular regulation mechanisms in the muscle development of livestock. More importantly, the output of our study will benefit the researchers and producers in the rabbit breeding program.

Hepatic transcriptome analysis identifies genes, polymorphisms and pathways involved in the fatty acids metabolism in sheep

Fatty acids (FA) in ruminants, especially unsaturated FA (USFA) have important impact in meat quality, nutritional value, and flavour quality of meat, and on consumer’s health. Identification of the genetic factors controlling the FA composition and metabolism is pivotal to select sheep that produce higher USFA and lower saturated (SFA) for the benefit of sheep industry and consumers. Therefore, this study was aimed to investigate the transcriptome profiling in the liver tissues collected from sheep with divergent USFA content in longissimus muscle using RNA deep-sequencing. From sheep (n = 100) population, liver tissues with higher (n = 3) and lower (n = 3) USFA content were analysed using Illumina HiSeq 2500. The total number of reads produced for each liver sample were ranged from 21.28 to 28.51 million with a median of 23.90 million. Approximately, 198 genes were differentially regulated with significance level of p-adjusted value <0.05. Among them, 100 genes were up-regulated, and 98 were down-regulated (p<0.01, FC>1.5) in the higher USFA group. A large proportion of key genes involved in FA biosynthesis, adipogenesis, fat deposition, and lipid metabolism were identified, such as APOA5, SLC25A30, GFPT1, LEPR, TGFBR2, FABP7, GSTCD, and CYP17A. Pathway analysis revealed that glycosaminoglycan biosynthesis- keratan sulfate, adipokine signaling, galactose metabolism, endocrine and other factors-regulating calcium metabolism, mineral metabolism, and PPAR signaling pathway were playing important regulatory roles in FA metabolism. Importantly, polymorphism and association analyses showed that mutation in APOA5, CFHR5, TGFBR2 and LEPR genes could be potential markers for the FA composition in sheep. These polymorphisms and transcriptome networks controlling the FA variation could be used as genetic markers for FA composition-related traits improvement. However, functional validation is required to confirm the effect of these SNPs in other sheep population in order to incorporate them in the sheep breeding program.

Signature selection analysis reveals candidate genes associated with production traits in Iranian sheep breeds

BACKGROUND

Sheep were among the first animals to be domesticated. They are raised all over the world and produce a major scale of animal-based protein for human consumption and play an important role in agricultural economy. Iran is one of the important locations for sheep genetic resources in the world. Here, we compared the Illumina Ovine SNP50 BeadChip data of three Iranian local breeds (Moghani, Afshari and Gezel), as a population that does not undergone artificial breeding programs as yet, and five other sheep breeds namely East Friesian white, East Friesian brown, Lacaune, DorsetHorn and Texel to detect genetic mechanisms underlying economical traits and daptation to harsh environments in sheep.

RESULTS

To identify genomic regions that have been targeted by positive selection, we used fixation index (Fst) and nucleotide diversity (Pi) statistics. Further analysis indicated candidate genes involved in different important traits such as; wool production included crimp of wool (PTPN3, NBEA and KRTAP20-2 genes), fiber diameter (PIK3R4 gene), hair follicle development (LHX2 gene), the growth and development of fiber (COL17A1 gene)), adaptation to hot arid environments (CORIN gene), adaptive in deficit water status (CPQ gene), heat stress (PLCB4, FAM107B, NBEA, PIK3C2B and USP43 genes) in sheep.

CONCLUSIONS

We detected several candidate genes related to wool production traits and adaptation to hot arid environments in sheep that can be applicable for inbreeding goals. Our findings not only include the results of previous researches, but also identify a number of novel candidate genes related to studied traits. However, more works will be essential to acknowledge phenotype- genotype relationships of the identified genes in our study.

Excavation of Genes Related to the Mining of Growth, Development, and Meat Quality of Two Crossbred Sheep Populations Based on Comparative Transcriptomes

Simple Summary

In this study, we measured the performance parameters of two crossbred sheep breeds, using Masson staining of the muscle tissue, and using the Illumina high-throughput sequencing platform to determine the differentially expressed genes (DEGs) in Dorper (DP) × Small-tailed Han (STH) sheep and Mongolia (MG) × Small-tailed Han sheep (STH). New transcripts of the muscle transcriptome were examined for the first time. DP × STH sheep were superior to MG × STH sheep in terms of meat quality and muscle morphology. In addition, 13 DEGs were found to play important roles in growth, development, and meat quality. The findings of this work may provide valuable resources for future research on muscle development in sheep.

Abstract

Crossbreeding can improve production performance and meat quality in sheep. The objective of this study was to look for genes related to sheep growth, development, and muscle. In this study, Dorper (DP) × Small Tailed Han (STH) sheep and Mongolia (MG) × Small-tailed Han (STH) sheep were used to estimate the productive performance and meat quality in a crossbreed. Subsequently, transcriptome analysis and bioinformatic analysis were performed on the Longissimus dorsi muscles of DP × STH and MG × STH sheep to identify differentially expressed genes (DEGs) related to growth, development, and meat quality. The presence of DEGs was confirmed by real-time PCR (qPCR). Productive performance and meat quality of the DP × STH sheep were better than the MG × STH sheep. Compared to DP × STH, a total of 1445 DEGs were identified in MG × STH sheep (1026 DEG were up-regulated and 419 DEG were down-regulated). Of these, 38 DEGs were related to growth, 161 to development, and 43 to muscle. In addition, 13 co-expressed genes (FGFRL1, SIX1, PLCB1, CRYAB, MYL2, ADIPOQ, GPX1, PPARD, GPC1, CDC42, LOC101106246, IGF1, and LARGE) were identified. The expression of DEGs was consistent with the comparative transcriptome analysis. This work provides genetics resources for future research on muscle development in sheep.

Identification of Meat Quality Determining Marker Genes in Fibroblasts of Bovine Muscle Using Transcriptomic Profiling

In the present study, we comparatively analyzed the transcriptomic profiling of fibroblasts derived from two different muscles, biceps femoris and longissimus dorsi with significant difference in the meat quality and tenderness. EBSeq algorithm was applied to analyze the data, and genes were considered to be significantly differentially expressed if the false discovery rate value was <0.05, the P value was <0.01, and the fold change was >0.585. The results revealed that 253 genes were differentially expressed genes (DEGs) (170 genes were upregulated, and 83 were downregulated) and more than 100 DEGs were probably associated with intramuscular fat deposition, tenderness, and toughness, which are driving the meat quality and were involved in biological processes such as collagen synthesis, cell differentiation, and muscle tissue and fiber development; molecular functions such as chemokine activity and collagen activity; cellular components such as cytoplasm and myofibril; and pathways such as collagen signaling and metabolic pathways. A gene-act network and a co-expression network revealed the close relationship between intramuscular fat deposition and meat tenderness. The expressions of 20 DEGs were validated by real-time PCR, and the results suggested that the DEGs are correlated with RNA-seq data and play crucial roles in muscle growth, development processes, toughness, and tenderness of the meat. Together, the genome-wide transcriptome analysis revealed that various genes are responsible for toughness and tenderness variance in the difference muscles of beef.

Comparative analysis of hypothalamus transcriptome between laying hens with different egg-laying rates

Egg-laying performance is one of the most important economic traits in the poultry industry. Commercial layers can lay one egg almost every day during their peak-laying period. However, many Chinese indigenous chicken breeds show a relatively low egg-laying rate, even during their peak-laying period. To understand what makes the difference in egg production, we compared the hypothalamus transcriptome profiles of Lushi blue-shelled-egg chickens (LBS), a Chinese indigenous breed with low egg-laying rate and Rhode Island Red chickens (RIR), a commercial layer with relatively high egg-laying rate using RNA-seq. A total of 753 differentially expressed genes (DEGs) were obtained. Of these DEGs, 38 genes were enriched in 2 Gene Ontology (GO) terms, namely reproduction term and the reproductive process term, and 6 KEGG pathways, namely Wnt signaling pathway, Oocyte meiosis, GnRH signaling pathway, Thyroid hormone signaling pathway, Thyroid hormone synthesis and MAPK signaling pathway, which have been long known to be involved in egg production regulation. To further determine the core genes from the 38 DEGs, protein-protein interaction (PPI) network, co-expression network and transcriptional regulatory network analyses were carried out. After integrated analysis and experimental validation, 4 core genes including RAC1, MRE11A, MAP7 and SOX5 were identified as the potential core genes that are responsible for the laying-rate difference between the 2 breeds. These findings paved the way for future investigating the mechanism of egg-laying regulation and enriched the chicken reproductive regulation theory.

Skeletal Muscle Transcriptome Analysis of Hanzhong Ma Duck at Different Growth Stages Using RNA-Seq

As one of the most important poultry worldwide, ducks (Anas platyrhynchos) are raised mainly for meat and egg products, and muscle development in ducks is important for meat production. Therefore, an investigation of gene expression in duck skeletal muscle would significantly contribute to our understanding of muscle development. In this study, twenty-four cDNA libraries were constructed from breast and leg muscles of Hanzhong Ma ducks at day 17, 21, 27 of the embryo and postnatal at 6-month-old. High-throughput sequencing and bioinformatics were used to determine the abundances and characteristics of transcripts. A total of 632,172,628 (average 52,681,052) and 637,213,938 (average 53,101,162) reads were obtained from the sequencing data of breast and leg muscles, respectively. Over 71.63% and 77.36% of the reads could be mapped to the Anas platyrhynchos genome. In the skeletal muscle of Hanzhong duck, intron variant (INTRON), synonymous variant (SYNONYMOUS_CODING), and prime 3' UTR variant (UTR_3_PRIME) were the main single nucleotide polymorphisms (SNP) annotation information, and "INTRON", "UTR_3_PRIME", and downstream-gene variant (DOWNSTREAM) were the main insertion-deletion (InDel) annotation information. The predicted number of alternative splicing (AS) in all samples were mainly alternative 5' first exon (transcription start site)-the first exon splicing (TSS) and alternative 3' last exon (transcription terminal site)-the last exon splicing (TTS). Besides, there were 292 to 2801 annotated differentially expressed genes (DEGs) in breast muscle and 304 to 1950 annotated DEGs in leg muscle from different databases. It is worth noting that 75 DEGs in breast muscle and 49 DEGs in leg muscle were co-expressed at all developmental points of comparison, respectively. The RNA-Seq data were confirmed to be reliable by qPCR. The identified DEGs, such as CREBL2, RHEB, GDF6, SHISA2, MYLK2, ACTN3, RYR3, and STMN1, were specially highlighted, indicating their strong associations with muscle development in the Hanzhong Ma duck. KEGG pathway analysis suggested that regulation of actin cytoskeleton, oxidative phosphorylation, and focal adhesion were involved in the development of skeletal muscle. The findings from this study can contribute to future investigations of the growth and development mechanism in duck skeletal muscle.

Comparative Transcriptome Profile Analysis of Longissimus dorsi Muscle Tissues From Two Goat Breeds With Different Meat Production Performance Using RNA-Seq

Carcass weight, meat quality and muscle components are important traits economically and they underpin most of the commercial return to goat producers. In this study, the Longissimus dorsi muscle tissues were collected from five Liaoning cashmere (LC) goats and five Ziwuling black (ZB) goats with phenotypic difference in carcass weight, some meat quality traits and muscle components. The histological quantitative of collagen fibers and the transcriptome profiles in the Longissimus dorsi muscle tissues were investigated using Masson-trichrome staining and RNA-Seq, respectively. The percentage of total collagen fibers in the Longissimus dorsi muscle tissues from ZB goats was less than those from LC goats, suggesting that these ZB goats had more tender meat. An average of 15,919 and 15,582 genes were found to be expressed in Longissimus dorsi muscle tissues from LC and ZB goats, respectively. Compared to LC goats, the expression levels of 78 genes were up-regulated in ZB goats, while 133 genes were down-regulated. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that the differentially expressed genes (DEGs) were significantly enriched in GO terms related to the muscle growth and development and the deposition of intramuscular fat and lipid metabolism, hippo signaling pathway and Jak-STAT signaling pathway. The results provide an improved understanding of the genetic mechanisms regulating meat production performance in goats, and will help us improve the accuracy of selection for meat traits in goats using marker-assisted selection based on these differentially expressed genes obtained.

Genome-wide identification and analysis of long noncoding RNAs in longissimus muscle tissue from Kazakh cattle and Xinjiang brown cattle

Objective

In recent years, lncRNAs have been identified in many species, and some of them have been shown to play important roles in muscle development and myogenesis. However, the differences in lncRNAs between Kazakh cattle and Xinjiang brown cattle remain undefined; therefore, we aimed to confirm whether lncRNAs are differentially expressed in the longissimus dorsi between these two types of cattle and whether differentially expressed lncRNAs regulate muscle differentiation.

Methods

We used RNA-seq technology to identify lncRNAs in longissimus muscles from these cattle. The expression of lncRNAs were analyzed using StringTie (1.3.1) in terms of the FPKM values of the encoding genes. The differential expression of the transcripts in the two samples were analyzed using the DESeq R software package. The resulting FDR was controlled by the Benjamini and Hochberg's approach. KOBAS software was utilized to measure the expression of different genes in KEGG pathways. We randomly selected eight lncRNA genes and validated them by RT-qPCR.

Results

We found that 182 lncRNA transcripts, including 102 upregulated and 80 downregulated transcripts, were differentially expressed between Kazakh cattle and Xinjiang brown cattle. The results of RT-qPCR were consistent with the sequencing results. Enrichment analysis and functional annotation of the target genes revealed that the differentially expressed lncRNAs were associated with the MAPK, Ras and PI3k/Akt signaling pathways. We also constructed a lncRNA/mRNA coexpression network for the PI3k/Akt signaling pathway.

Conclusion

Our study provides insights into cattle muscle-associated lncRNAs and will contribute to a more thorough understanding of the molecular mechanism underlying muscle growth and development in cattle.

Variants in myostatin and MyoD family genes are associated with meat quality traits in Santa Inês sheep

Myostatin and MyoD family genes play vital roles in myogenesis and this study aimed to identify association of variants in MyoD1, MyoG, MyF5, MyF6, and MSTN genes with meat traits in Santa Inês sheep. A dataset with 44 variants and records of seven meat traits in 192 lambs (pH0, pH24, a*, b*, L*, tenderness assessed by shear force, and water-holding capacity) was used. Single-locus and haplotype association analyses were performed, and the significance threshold was established according to Bonferroni's method. Single-locus analysis revealed two associations at a Bonferroni level, where the variant c.935-185C > G in MyoD1 had an additive effect (-4.31 ± 1.08 N) on tenderness, while the variant c.464 + 185G > A in MyoG had an additive effect (-2.86 ± 0.64) on a*. Additionally, the haplotype replacement GT>AC in MSTN was associated with pH0 (1.26 ± 0.31), pH24 (1.07 ± 0.27), a* (-1.40 ± 0.51), and tenderness (3.83 ± 1.22 N), while the replacement GT > AG in MyoD1 was associated with pH0 (1.43 ± 0.26), pH24 (1.25 ± 0.22), b* (-1.06 ± 0.39), and tenderness (-4.13 ± 1.16 N). Our results have demonstrated that some variants in MyoG, MyF6, MyoD1, and MSTN can be associated with physicochemical meat traits in Santa Inês sheep.

Comparative Transcriptome Analysis Identifying the Different Molecular Genetic Markers Related to Production Performance and Meat Quality in Longissimus Dorsi Tissues of MG × STH and STH Sheep

Crossbred sheep have many prominent traits, such as excellent production performance and high-quality meat, when compared to local sheep breeds. However, the genetic molecular markers related to these characteristics remain unclear. The crossbred MG × STH (small-tailed Han sheep (STH) × Mongolian sheep (MG)) breed and the STH breed were selected to measure production performance and meat quality. We used 14 indexes of production performance and meat quality, which in the MG × STH population showed significant differences compared to the STH breed. Subsequently, the longissimusdorsi from the two sheep were subjected to comparative transcriptomic analyses to identify differentially expressed genes (DEGs) related to production performance and meat quality. A total of 874 DEGs were identified between the two sheep groups. A total of 110 unique DEGs related to sheep production performance and meat quality were selected as the candidate DEGs. We found 6 production-performance-related and 30 meat-quality-related DEGs through a correlation analysis, including SPARC, ACVRL1, FNDC5 and FREM1. The expression levels of 11 DEGs were validated by real-time PCR, and the results were in accordance with the results of the comparative transcriptomic and correlation analyses. These results will assist in understanding sheep heterosis and molecular marker-assisted selection.

Transcriptome profile of skeletal muscle at different developmental stages in Large White and Mashen pigs

From the perspectives of promoting individual growth and development, increasing pork yield, and improving feed utilization, it is desirable to screen candidate genes underlying pig muscle growth and regulation. In this study, we investigated transcriptome differences at 1, 90, and 180 d of age in Large White and Mashen pigs, characterized differentially expressed genes (DEGs), and screened candidate genes affecting skeletal muscle growth and development. RNA-seq was applied to analyze the transcriptome of the longissimus dorsi (LD) in the two breeds. In LD samples from the two breeds at three growth stages, 7215, 6332, 237, 3935, 3404, and 846 DEGs were obtained for L01 vs. L90, L01 vs. L180, L90 vs. L180, MS01 vs. MS90, MS01 vs. MS180, and MS90 vs. MS180, respectively. Significant tendencies in DEG expression could be grouped into eight profiles. Based on the functional analysis of DEGs, 16 candidate genes related to skeletal muscle growth and development were identified, including PCK2, GNAS, ADCY2, PRKAB1, PRKAB2, PRKAG1, PRKAG2, PHKA1, PHKA2, PHKG1, PHKG2, ITPR3, IGF1R, FGFR4, FGF1, and FGF18. The results of this study thus provide a theoretical basis for the mechanisms and candidate genes underlying skeletal muscle development in pigs.

A Transcriptomic Study of the Tail Fat Deposition in Two Types of Hulun Buir Sheep According to Tail Size and Sex

Simple Summary

Based on tail types, Hulun Buir sheep were divided into two lines including small and big fat-tailed, but these two lines have similar genetic background. In this study, we investigated the morphology and transcription level differences of tail fat between these two lines. The RNA-seq analyses indicated several differentially expressed genes when compared between sexes or two tail sizes. Interestingly, we also found an obvious sex difference in the fat metabolism in Hulun Buir sheep. Two different co-expression networks were only shown either in male or in female sheep. Our findings will provide theoretical background in understanding the genetic mechanism of fat deposition in sheep.

Abstract

Hulun Buir sheep of similar genetic background were divided into two lines based on tail types: Small- and big fat-tailed. To explore the molecular mechanism of fat deposition in sheep tails, we firstly evaluated the morphology and transcription level differences of tail fat between these two lines. RNA-Seq technology was used to identify differentially expressed genes (DEGs) in phenotypic extremes of tail sizes. Five comparisons were performed taking into account two factors, sex and tail type. We screened out 373 DEGs between big-tailed and small-tailed Hulun Buir sheep, and 775 and 578 DEGs between two types of tails in male and female sheep, respectively. The results showed an obvious sex difference in the fat metabolism in sheep based on gene ontology (GO), pathway, and network analyses. Intriguingly, there were two different co-expression networks only respectively shown in male and female sheep, which were insulin-related network acting on upstream pathways and PPARG-related network effect in downstream pathways. Furthermore, these two networks were linked by a classic pathway of regulating adipogenesis. This is the first study to investigate the sex differences of fat metabolism in domestic animals, and it demonstrates a new experimental way to study fat metabolism. Our findings will provide theoretical background in understanding the tail-size phenotype in sheep and can be exploited in breeding small-tailed sheep.

NK, S. S, Fairoze MN, Kaur M, Sharma A, Girdhar Y, M. SR, Devatkal SK, Ahlawat S, Vijh RK, S. MS. Transcriptome profiling of longissimus thoracis muscles identifies highly connected differentially expressed genes in meat type sheep of India

This study describes the muscle transcriptome profile of Bandur breed, a consumer favoured, meat type sheep of India. The transcriptome was compared to the less desirable, unregistered local sheep population, in order to understand the molecular factors related to muscle traits in Indian sheep breeds. Bandur sheep have tender muscles and higher backfat thickness than local sheep. The longissimus thoracis transcriptome profiles of Bandur and local sheep were obtained using RNA sequencing (RNA Seq). The animals were male, non-castrated, with uniform age and reared under similar environment, as well as management conditions. We could identify 568 significantly up-regulated and 538 significantly down-regulated genes in Bandur sheep (p≤0.05). Among these, 181 up-regulated and 142 down-regulated genes in Bandur sheep, with a fold change ≥1.5, were considered for further analysis. Significant Gene Ontology terms for the up-regulated dataset in Bandur sheep included transporter activity, substrate specific transmembrane, lipid and fatty acid binding. The down-regulated activities in Bandur sheep were mainly related to RNA degradation, regulation of ERK1 and ERK2 cascades and innate immune response. The MAPK signaling pathway, Adipocytokine signaling pathway and PPAR signaling pathway were enriched for Bandur sheep. The highly connected genes identified by network analysis were CNOT2, CNOT6, HSPB1, HSPA6, MAP3K14 and PPARD, which may be important regulators of energy metabolism, cellular stress and fatty acid metabolism in the skeletal muscles. These key genes affect the CCR4-NOT complex, PPAR and MAPK signaling pathways. The highly connected genes identified in this study, form interesting candidates for further research on muscle traits in Bandur sheep.

Integrative microRNA-mRNA Analysis of Muscle Tissues in Qianhua Mutton Merino and Small Tail Han Sheep Reveals Key Roles for oar-miR-655-3p and oar-miR-381-5p

The Qianhua Mutton Merino (QHMM) is a new variety of sheep (Ovis aries) with improved meat performance compared with the traditional Small Tail Han (STH) sheep variety. We recently reported the transcriptome profiling of longissimus muscle tissues between QHMM and STH sheep. In the present study, we aimed to evaluate key micro (mi)RNA-mRNA networks associated with sheep muscle growth and development. We used miRNA sequencing to obtain longissimus muscle miRNA profiles from QHMM and STH sheep. We identified a total of 153 known sheep miRNAs, of which 4 were differentially expressed (DE) between the 2 sheep varieties. We combined these results with mRNA library data to build an miRNA-mRNA network, including 26 target genes of the 4 DE miRNAs. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that 26 target genes were significantly enriched in 86 biological processes, including muscle organogenesis, myoblast migration, cell proliferation, and adipose tissue development, and in 9 metabolic pathways, including carbohydrate, nucleotide, and amino acid metabolic pathways. oar-miR-655-3p and its target gene ACSM3 and oar-miR-381-5p and its target gene ABAT were selected for subsequent analysis based on GO and KEGG analyses. The binding sites of oar-miR-655-3p with ACSM3 and oar-miR-381-5p with ABAT were validated by a dual-luciferase reporter gene detection system. This represents the first integrative analysis of miRNA-mRNA networks in QHMM and STH muscles and suggests that DE miRNAs, especially oar-miR-655-3p and oar-miR-381-5p, play crucial roles in muscle growth and development.

Selection of candidate genes affecting meat quality and preliminary exploration of related molecular mechanisms in the Mashen pig

Objective

The aim of this study was to select the candidate genes affecting meat quality and preliminarily explore the related molecular mechanisms in the Mashen pig.

Methods

The present study explored genetic factors affecting meat quality in the Mashen pig using RNA sequencing (RNA-Seq). We sequenced the transcriptomes of 180-day-old Mashen and Large White pigs using longissimus dorsi to select differentially expressed genes (DEGs).

Results

The results indicated that a total of 425 genes were differentially expressed between Mashen and Large White pigs. A gene ontology enrichment analysis revealed that DEGs were mainly enriched for biological processes associated with metabolism and muscle development, while a Kyoto encyclopedia of genes and genomes analysis showed that DEGs mainly participated in signaling pathways associated with amino acid metabolism, fatty acid metabolism, and skeletal muscle differentiation. A MCODE analysis of the protein-protein interaction network indicated that the four identified subsets of genes were mainly associated with translational initiation, skeletal muscle differentiation, amino acid metabolism, and oxidative phosphorylation pathways.

Conclusion

Based on the analysis results, we selected glutamic-oxaloacetic transaminase 1, malate dehydrogenase 1, pyruvate dehydrogenase 1, pyruvate dehydrogenase kinase 4, and activator protein-1 as candidate genes affecting meat quality in pigs. A discussion of the related molecular mechanisms is provided to offer a theoretical basis for future studies on the improvement of meat quality in pigs.

Comparative transcriptome analysis reveals regulators mediating breast muscle growth and development in three chicken breeds

Objective: The goal of this study was to investigate the mechanisms of muscle growth and development of three chicken breeds. Participants: Eighteen chickens, including three different breeds with different growth speeds (White Broiler, Daheng, and Commercial Layers of Roman), were used. Methods: Total RNA from breast muscle of these chickens was subjected to a gene expression microarray. Differentially expressed genes (DEGs) were screened and functional enrichment analysis was performed using DAVID. Seven DEGs were confirmed by quantitative reverse transcription PCR. Results: Overall, 8,398 DEGs were found among the different lines. The DEGs between each two lines that were unique for a developmental stage were greater than those that were common during all stages. Functional analysis revealed that DEGs across the entire developmental process were primarily involved in positive cell proliferation, growth, cell differentiation, and developmental processes. Genes involved in muscle regulation, muscle construction, and muscle cell differentiation were upregulated in the faster-growing breed compared to the slower-growing breed. DEGs including myosin heavy chain 15 (MYH15), myozenin 2 (MYOZ2), myosin-binding protein C (MYBPC3), insulin-like growth factor 2 (IGF2), apoptosis regulator (BCL-2), AP-1 transcription factor subunit (JUN), and AP-1 transcription factor subunit (FOS) directly regulated muscle growth or were in the center of the protein-protein interaction network. Pathways, including the extracellular matrix (ECM)-receptor interaction, mitogen-activated protein kinase (MAPK) signaling pathway, and focal adhesion, were the most enriched DEGs between lines or within lines under different developmental stages. Conclusions: Genes involved in muscle construction and cell differentiation were differentially expressed among the three breeds.

De Novo Transcriptome Sequencing Analysis of Goose (Anser anser) Embryonic Skin and the Identification of Genes Related to Feather Follicle Morphogenesis at Three Stages of Development

The objective of this study was to evaluate the changes in the goose embryo transcriptome during feather development. RNA-Sequencing (RNA-Seq) was used to find the transcriptome profiles of feather follicles from three stages of embryonic dorsal skin at embryonic day 13, 18, and 28 (E13, E18, E28). The results showed that 3001, 6634, and 13,780 genes were differently expressed in three stages. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that differentially expressed genes (DEGs) in E13 vs. E18 were significantly mapped into the GO term of extracellular structure organization and the pathway of extracellular matrix (ECM)-receptor interaction. In E18 vs. E28, the top significantly mapped into GO term was the single-organism developmental process; the pathway was also the ECM-receptor interaction. DEGs in E13 vs. E28 were significantly mapped into the GO term of the multicellular organismal process and the pathway of cell adhesion molecules. Subsequently, the union of DEGs was categorized by succession cluster into eight profiles, which were then grouped into four ideal profiles. Lastly, the seven genes spatio-temporal expression pattern was confirmed by real-time PCR. Our findings advocate that interleukin 20 receptor subunit alpha (IL20RA), interleukin 6 receptor (IL6R), interleukin 1 receptor type 1 (IL-1R1), Wnt family member 3A (WNT3A), insulin-like growth factor binding protein 3 (IGFBP3), bone morphogenetic protein 7 (BMP7), and secreted-frizzled related protein 2 (SFRP2) might possibly play vital roles in skin and feather follicle development and growth processes.

Comparison of transcriptomic landscapes of different lamb muscles using RNA-Seq

Transcriptome deep sequencing is a powerful tool for exploring the genetic architecture of complex traits. Gene expression patterns may explain a high degree of the observed phenotypic differences in histochemical and metabolic parameters related to meat quality among different muscles. In this study, we sequenced by RNA-Seq the whole transcriptome of nine lamb muscles: Semimembranosus (SM), Semitendinosus (ST), Cranial gluteobiceps, Gluteus medius (GM), Rectus femoris, Supraspinatus (SS), Longissimus lumborum (LL), Adductor and Psoas major. Significant gene expression differences were detected between almost all pairwise comparisons, being more pronounced between SS and ST, SM and LL, and ST and GM. These differences can be explained in terms of ATPase and glycolytic activities, muscle fiber typing and oxidative score, clustering muscles as fast glycolytic, intermediate or slow oxidative. ST showed up-regulation of gene pathways related to carbohydrate metabolism, energy generation and protein turnover as expected from a fast white muscle. SS showed myosin isoforms typical of slow muscles and high expression of genes related to calcium homeostasis and vascularization. SM, LL and GM showed in general intermediate gene expression patterns. Several novel transcripts were detected, mostly related to muscle contraction and structure, oxidative metabolism, lipid metabolism and protein phosphorylation. Expression profiles were consistent with previous histochemical and metabolic characterization of these muscles. Up-regulation of ion transport genes may account for significant differences in water holding capacity. High expression of genes related to cell adhesion, cytoskeleton organization, extracellular matrix components and protein phosphorylation may be related to meat yellowness and lower tenderness scores. Differential expression of genes related to glycolytic activity and lactic acid generation among fast, intermediate and slow muscles may explain the detected final meat pH differences. These results reveal new candidate genes associated with lamb meat quality, and give a deeper insight into the genetic architecture of these complex traits.

Sheep skeletal muscle transcriptome analysis reveals muscle growth regulatory lncRNAs

As widely distributed domestic animals, sheep are an important species and the source of mutton. In this study, we aimed to evaluate the regulatory lncRNAs associated with muscle growth and development between high production mutton sheep (Dorper sheep and Qianhua Mutton Merino sheep) and low production mutton sheep (Small-tailed Han sheep). In total, 39 lncRNAs were found to be differentially expressed. Using co-expression analysis and functional annotation, 1,206 co-expression interactions were found between 32 lncRNAs and 369 genes, and 29 of these lncRNAs were found to be associated with muscle development, metabolism, cell proliferation and apoptosis. lncRNA–mRNA interactions revealed 6 lncRNAs as hub lncRNAs. Moreover, three lncRNAs and their associated co-expressed genes were demonstrated by cis-regulatory gene analyses, and we also found a potential regulatory relationship between the pseudogene lncRNA LOC101121401 and its parent gene FTH1. This study provides a genome-wide resolution of lncRNA and mRNA regulation in muscles from mutton sheep.

Dynamic transcriptomic analysis in hircine longissimus dorsi muscle from fetal to neonatal development stages

Muscle growth and development from fetal to neonatal stages consist of a series of delicately regulated and orchestrated changes in expression of genes. In this study, we performed whole transcriptome profiling based on RNA-Seq of caprine longissimus dorsi muscle tissue obtained from prenatal stages (days 45, 60, and 105 of gestation) and neonatal stage (the 3-day-old newborn) to identify genes that are differentially expressed and investigate their temporal expression profiles. A total of 3276 differentially expressed genes (DEGs) were identified (Q value < 0.01). Time-series expression profile clustering analysis indicated that DEGs were significantly clustered into eight clusters which can be divided into two classes (Q value < 0.05), class I profiles with downregulated patterns and class II profiles with upregulated patterns. Based on cluster analysis, GO enrichment analysis found that 75, 25, and 8 terms to be significantly enriched in biological process (BP), cellular component (CC), and molecular function (MF) categories in class I profiles, while 35, 21, and 8 terms to be significantly enriched in BP, CC, and MF in class II profiles. KEGG pathway analysis revealed that DEGs from class I profiles were significantly enriched in 22 pathways and the most enriched pathway was Rap1 signaling pathway. DEGs from class II profiles were significantly enriched in 17 pathways and the mainly enriched pathway was AMPK signaling pathway. Finally, six selected DEGs from our sequencing results were confirmed by qPCR. Our study provides a comprehensive understanding of the molecular mechanisms during goat skeletal muscle development from fetal to neonatal stages and valuable information for future studies of muscle development in goats.

A high resolution atlas of gene expression in the domestic sheep (Ovis aries)

Sheep are a key source of meat, milk and fibre for the global livestock sector, and an important biomedical model. Global analysis of gene expression across multiple tissues has aided genome annotation and supported functional annotation of mammalian genes. We present a large-scale RNA-Seq dataset representing all the major organ systems from adult sheep and from several juvenile, neonatal and prenatal developmental time points. The Ovis aries reference genome (Oar v3.1) includes 27,504 genes (20,921 protein coding), of which 25,350 (19,921 protein coding) had detectable expression in at least one tissue in the sheep gene expression atlas dataset. Network-based cluster analysis of this dataset grouped genes according to their expression pattern. The principle of 'guilt by association' was used to infer the function of uncharacterised genes from their co-expression with genes of known function. We describe the overall transcriptional signatures present in the sheep gene expression atlas and assign those signatures, where possible, to specific cell populations or pathways. The findings are related to innate immunity by focusing on clusters with an immune signature, and to the advantages of cross-breeding by examining the patterns of genes exhibiting the greatest expression differences between purebred and crossbred animals. This high-resolution gene expression atlas for sheep is, to our knowledge, the largest transcriptomic dataset from any livestock species to date. It provides a resource to improve the annotation of the current reference genome for sheep, presenting a model transcriptome for ruminants and insight into gene, cell and tissue function at multiple developmental stages.

Genome-Wide Analysis Reveals Extensive Changes in LncRNAs during Skeletal Muscle Development in Hu Sheep

As an important type of noncoding RNA molecules, long non-coding RNAs (lncRNAs) act as versatile players in various biological processes. However, little is known about lncRNA regulators during sheep muscle growth. To explore functional lncRNAs during sheep muscle growth, we systematically investigated lncRNAs using strand-specific Ribo-Zero RNA sequencing at three key developmental stages in Hu sheep. A total of 6924 lncRNAs were obtained, and the differentially expressed lncRNAs and genes were screened from (control vs. experiment) fetus vs. lamb, lamb vs. adult, and fetus vs. adult comparisons, respectively. The quantitative real-time polymerase chain reaction (qRT-PCR) analysis results correlated well with the sequencing data. Moreover, functional annotation analysis based on the Gene Ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) databases showed that the target genes of the differentially expressed lncRNAs were significantly enriched in organ morphogenesis, skeletal system development as well as response to stimulus and some other terms related to muscle. Furthermore, a co-expression network of the differentially expressed target genes and lncRNAs was constructed and well-known muscle growth regulators such as retrotransposon-like 1 and Junctophilin-2 were included. Finally, we investigated the expression profiles of seven lncRNAs and their target genes, and found that they played vital roles in muscle growth. This study extends the sheep muscle lncRNA database and provides novel candidate regulators for future genetic and molecular studies on sheep muscle growth, which is helpful for optimizing the production of mutton.

Transcriptome Analysis of Three Sheep Intestinal Regions reveals Key Pathways and Hub Regulatory Genes of Large Intestinal Lipid Metabolism

The large intestine, also known as the hindgut, is an important part of the animal digestive system. Recent studies on digestive system development in ruminants have focused on the rumen and the small intestine, but the molecular mechanisms underlying sheep large intestine metabolism remain poorly understood. To identify genes related to intestinal metabolism and to reveal molecular regulation mechanisms, we sequenced and compared the transcriptomes of mucosal epithelial tissues among the cecum, proximal colon and duodenum. A total of 4,221 transcripts from 3,254 genes were identified as differentially expressed transcripts. Between the large intestine and duodenum, differentially expressed transcripts were found to be significantly enriched in 6 metabolism-related pathways, among which PPAR signaling was identified as a key pathway. Three genes, CPT1A, LPL and PCK1, were identified as higher expression hub genes in the large intestine. Between the cecum and colon, differentially expressed transcripts were significantly enriched in 5 lipid metabolism related pathways, and CEPT1 and MBOAT1 were identified as hub genes. This study provides important information regarding the molecular mechanisms of intestinal metabolism in sheep and may provide a basis for further study.