Differential expression of miRNAs in skeletal muscles of Indian sheep with diverse carcass and muscle traits

Exploring transcriptomic mechanisms underlying pulmonary adaptation to diverse environments in Indian rams

BACKGROUND

The Changthangi sheep thrive at high altitudes in the cold desert regions of Ladakh, India while Muzaffarnagri sheep are well-suited to the low altitude plains of northern India. This study investigates the molecular mechanisms of pulmonary adaptation to diverse environments by analyzing gene expression profiles of lung tissues through RNA sequencing.

METHODS AND RESULTS

Four biological replicates of lung tissue from each breed were utilized to generate the transcriptomic data. Differences in gene expression analysis revealed discrete expression profiles in lungs of each breed. In Changthangi sheep, genes related to immune responses, particularly cytokine signaling, were significantly enriched. Pathway analysis highlighted the activation of NF-kB signaling, a key mediator of inflammation and immune response. Additionally, the gene network analysis indicated a strong association between cytokine signaling, hypoxia-inducible factor (HIF) and NF-kB activation, suggesting a coordinated response to hypoxic stress in lungs of Changthangi sheep. In Muzaffarnagri sheep, the gene expression profiles were enriched for pathways related to energy metabolism, homeostasis and lung physiology. Key pathways identified include collagen formation and carbohydrate metabolism, both of which are crucial for maintaining lung function and structural integrity. Gene network analysis further reinforced this by revealing a strong connection between genes associated with lung structure and function.

CONCLUSIONS

Our findings shed light on the valuable insights into gene expression mechanisms that enable these sheep breeds to adapt to their respective environments and contribute to a better understanding of high altitude adaptation in livestock.

Skeletal muscle transcriptomics of sheep acclimated to cold desert and tropical regions identifies genes and pathways accentuating their diversity

The current study attempts to investigate the differences in gene expression in longissimus thoracis muscles between sheep breeds acclimated to diverse environments. Changthangi sheep inhabits the cold arid plateau of Ladakh, at an altitude above 3000 m with prevalence of rarefied atmosphere. Muzzafarnagri sheep, on the other hand is found in the sub-tropical hot and humid plains at an altitude of about 250 m. Comparative transcriptomics was used to provide a molecular perspective of the differential adaptation of the two breeds. RNA sequencing data was generated from four biological replicates of the longissimus thoracis muscles from both breeds. The common genes expressed in both breeds were involved in muscle contraction and muscle fibre organization. The most significant pathways enriched in Changthangi muscles were glycogen metabolism, reduction of cytosolic Ca++ levels and NFE2L2 regulating anti-oxidant, while those in Muzzafarnagri were extracellular matrix organization and collagen formation. The hub genes identified in Changthangi were involved in hematopoiesis and HIF signaling pathway, suggesting the molecular acclimatization of Changthangi to the high altitude cold desert of Ladakh. The nodal genes discovered in Muzzafarnagri sheep were associated with the extracellular matrix which accentuates its significance in the development, growth and repair of muscles. The observed transcriptomic differences underscore the morphological and adaptive disparity between the two breeds. The candidate genes and pathways identified in this study will form the basis for future research on adaptation to high altitude and body size in small ruminants.

microRNA-1 Regulates Metabolic Flexibility in Skeletal Muscle via Pyruvate Metabolism

MicroRNA-1 (miR-1) is the most abundant miRNA in adult skeletal muscle. To determine the function of miR-1 in adult skeletal muscle, we generated an inducible, skeletal muscle-specific miR-1 knockout (KO) mouse. Integration of RNA-sequencing (RNA-seq) data from miR-1 KO muscle with Argonaute 2 enhanced crosslinking and immunoprecipitation sequencing (AGO2 eCLIP-seq) from human skeletal muscle identified miR-1 target genes involved with glycolysis and pyruvate metabolism. The loss of miR-1 in skeletal muscle induced cancer-like metabolic reprogramming, as shown by higher pyruvate kinase muscle isozyme M2 (PKM2) protein levels, which promoted glycolysis. Comprehensive bioenergetic and metabolic phenotyping combined with skeletal muscle proteomics and metabolomics further demonstrated that miR-1 KO induced metabolic inflexibility as a result of pyruvate oxidation resistance. While the genetic loss of miR-1 reduced endurance exercise performance in mice and in C. elegans, the physiological down-regulation of miR-1 expression in response to a hypertrophic stimulus in both humans and mice causes a similar metabolic reprogramming that supports muscle cell growth. Taken together, these data identify a novel post-translational mechanism of adult skeletal muscle metabolism regulation mediated by miR-1.

Molecular insights into Pashmina fiber production: comparative skin transcriptomic analysis of Changthangi goats and sheep

Ladakh, one of the highest inhabited regions globally, hosts the unique Changthangi goat, renowned for producing Pashmina, the world's most luxurious natural fiber. In comparison, the fiber derived from Changthangi sheep is considered next only to Pashmina. This research endeavors to compare the skin transcriptome profiles of Changthangi goats and Changthangi sheep, aiming to discern the molecular determinants behind the recognition of Changthangi goats as the source of Pashmina. Drawing upon previously conducted studies, a collective of 225 genes correlated with fiber characteristics were extracted from the differentially expressed genes noticed between the two species (p-value of ≤ 0.05 and a log2 fold change of ≥ 1.5). These genes were analyzed using DAVID software to understand their biological functions and to identify enriched KEGG and Reactome pathways. The protein-protein interaction networks were constructed using Cytoscape, cytoHubba, and STRING to focus on key genes and infer their biological significance. Comparative transcriptome analysis revealed significantly higher expression of genes involved in signaling pathways like Wnt, MAPK, PI3K-Akt, Hedgehog, associated with fiber development and quality in Changthangi goats. These pathways play crucial roles in hair follicle (HF) formation, maintenance of epidermal stem cells, and fiber characteristics. Findings also highlight the enrichment of cell adhesion molecules and ECM-receptor interaction, emphasizing their roles in HF structure, growth, and signaling. This investigation offers an in-depth understanding of the molecular intricacies governing Pashmina production in Changthangi goats, providing valuable insights into their unique genetic makeup and underlying mechanisms influencing the exceptional quality of Pashmina fibers.

Identification and validation of stable reference genes for expression profiling of target genes in diverse ovine tissues

Quantitative PCR (qPCR) is a widely-used technique for quantifying the expression of target genes across various tissues, as well as under different pathological and physiological conditions. One of the challenges associated with this method is the need to identify optimal reference genes (RGs) that maintain consistent expression levels under diverse experimental settings, thereby ensuring accurate biological interpretation. In this study, we conducted a thorough analysis of 18 candidate RGs (ACTB, BACH1, B2M, GAPDH, HMBS, HPRT1, PGK1, PPIA, PPIB, RPLP0, RPL19, RPS9, RPS15, RPS28, SDHA, TBP, UXT, and YWHAZ) across 10 ovine tissues (muscle, skin, kidney, liver, intestine, rumen, lung, testis, heart, and spleen) obtained from five individual sheep. We aimed to identify genes with stable expression across these tissues. A literature-based survey helped us shortlist candidate genes representing various functional classes from multiple livestock species. We employed four algorithms: geNorm, NormFinder, BestKeeper, and Delta Ct (ΔCt), to rank these genes based on their stability. A consistent trend in the rankings was observed across these different algorithms. RefFinder was then used for a comprehensive ranking, integrating the outputs from the various methods. ACTB, PPIB, BACH1, and B2M emerged as the most stable RGs, while RPS9, RPS15, and PGK1 displayed variable expression. We validated our findings through qPCR analysis of four target genes (ACTN2, CRYAB, DLK1, and TRIM54) in the skin samples from two different sheep breeds. Based on these results, we recommend ACTB, PPIB, BACH1, and B2M as reliable internal control genes for qPCR experiments involving diverse ovine tissues.

Genetic insights into fiber quality, coat color and adaptation in Changthangi and Muzzafarnagri sheep: A comparative skin transcriptome analysis

Changthangi sheep, which inhabit the high-altitude regions of Ladakh, are known for their fine fiber production and are characterized by grey skin and either black or white coats. In contrast, Muzzafarnagri sheep from the plains of Uttar Pradesh produce coarse wool and have white skin and coats. We conducted comparative global gene expression profiling on four biological replicates of skin from each breed. Notably, our analysis identified 149 up-regulated genes and 2,139 down-regulated genes in Changthangi sheep compared to Muzzafarnagri sheep, with a p-adjusted value (padj) of ≤0.05 and a Log2 fold change of ≥1.5. Gene Ontology analysis of the up-regulated genes revealed an enrichment of terms related to melanin biosynthesis and developmental pigmentation. Additionally, enriched KEGG pathways included tyrosine metabolism and metabolic pathways. Among the melanogenesis-related genes that exhibited higher expression in Changthangi sheep were TYR, TYRP1, DCT, SLC45A2, PMEL, MLANA, and OCA2. These findings confirm melanin's role in both the animals' black coat color and UV protection at high-altitude. Furthermore, we observed more pronounced expression of genes related to fiber quality, namely KRTAP6, KRTAP7, KRTAP13, and KRTAP2, in the fine wool-producing sheep from Ladakh. The results of the RNA sequencing were validated using real-time PCR on 10 genes governing fiber quality and coat color, with ACTB and PPIB serving as reference genes. In conclusion, our comparative skin transcriptome analysis of Changthangi and Muzzafarnagri sheep sheds light on the genetic differences associated with distinct phenotypic traits and environmental adaptability, offering valuable insights into the underlying mechanisms.

Deep Small RNA Sequencing Reveals Important miRNAs Related to Muscle Development and Intramuscular Fat Deposition in Longissimus dorsi Muscle From Different Goat Breeds

MicroRNAs (miRNAs) are a class of small non-coding RNAs that have been shown to play important post-transcriptional regulatory roles in the growth and development of skeletal muscle tissues. However, limited research into the effect of miRNAs on muscle development in goats has been reported. In this study, Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats with significant phenotype difference in meat production performance were selected and the difference in Longissimus dorsi muscle tissue expression profile of miRNAs between the two goat breeds was then compared using small RNA sequencing. A total of 1,623 miRNAs were identified in Longissimus dorsi muscle tissues of the two goat breeds, including 410 known caprine miRNAs, 928 known species-conserved miRNAs and 285 novel miRNAs. Of these, 1,142 were co-expressed in both breeds, while 230 and 251 miRNAs were only expressed in LC and ZB goats, respectively. Compared with ZB goats, 24 up-regulated miRNAs and 135 miRNAs down-regulated were screened in LC goats. A miRNA-mRNA interaction network showed that the differentially expressed miRNAs would target important functional genes associated with muscle development and intramuscular fat deposition. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the target genes of differentially expressed miRNAs were significantly enriched in Ras, Rap 1, FoxO, and Hippo signaling pathways. This study suggested that these differentially expressed miRNAs may be responsible for the phenotype differences in meat production performance between the two goat breeds, thereby providing an improved understanding of the roles of miRNAs in muscle tissue of goats.

Emerging Roles of Non-Coding RNAs in the Feed Efficiency of Livestock Species

A global population of already more than seven billion people has led to an increased demand for food and water, and especially the demand for meat. Moreover, the cost of feed used in animal production has also increased dramatically, which requires animal breeders to find alternatives to reduce feed consumption. Understanding the biology underlying feed efficiency (FE) allows for a better selection of feed-efficient animals. Non-coding RNAs (ncRNAs), especially micro RNAs (miRNAs) and long non-coding RNAs (lncRNAs), play important roles in the regulation of bio-logical processes and disease development. The functions of ncRNAs in the biology of FE have emerged as they participate in the regulation of many genes and pathways related to the major FE indicators, such as residual feed intake and feed conversion ratio. This review provides the state of the art studies related to the ncRNAs associated with FE in livestock species. The contribution of ncRNAs to FE in the liver, muscle, and adipose tissues were summarized. The research gap of the function of ncRNAs in key processes for improved FE, such as the nutrition, heat stress, and gut–brain axis, was examined. Finally, the potential uses of ncRNAs for the improvement of FE were discussed.

Genome Wide Association Study of Beef Traits in Local Alpine Breed Reveals the Diversity of the Pathways Involved and the Role of Time Stratification

Knowledge of the genetic architecture of key growth and beef traits in livestock species has greatly improved worldwide thanks to genome-wide association studies (GWAS), which allow to link target phenotypes to Single Nucleotide Polymorphisms (SNPs) across the genome. Local dual-purpose breeds have rarely been the focus of such studies; recently, however, their value as a possible alternative to intensively farmed breeds has become clear, especially for their greater adaptability to environmental change and potential for survival in less productive areas. We performed single-step GWAS and post-GWAS analysis for body weight (BW), average daily gain (ADG), carcass fleshiness (CF) and dressing percentage (DP) in 1,690 individuals of local alpine cattle breed, Rendena. This breed is typical of alpine pastures, with a marked dual-purpose attitude and good genetic diversity. Moreover, we considered two of the target phenotypes (BW and ADG) at different times in the individuals' life, a potentially important aspect in the study of the traits' genetic architecture. We identified 8 significant and 47 suggestively associated SNPs, located in 14 autosomal chromosomes (BTA). Among the strongest signals, 3 significant and 16 suggestive SNPs were associated with ADG and were located on BTA10 (50-60 Mb), while the hotspot associated with CF and DP was on BTA18 (55-62 MB). Among the significant SNPs some were mapped within genes, such as SLC12A1, CGNL1, PRTG (ADG), LOC513941 (CF), NLRP2 (CF and DP), CDC155 (DP). Pathway analysis showed great diversity in the biological pathways linked to the different traits; several were associated with neurogenesis and synaptic transmission, but actin-related and transmembrane transport pathways were also represented. Time-stratification highlighted how the genetic architectures of the same traits were markedly different between different ages. The results from our GWAS of beef traits in Rendena led to the detection of a variety of genes both well-known and novel. We argue that our results show that expanding genomic research to local breeds can reveal hitherto undetected genetic architectures in livestock worldwide. This could greatly help efforts to map genomic complexity of the traits of interest and to make appropriate breeding decisions.

miRNA transcriptome and myofiber characteristics of lamb skeletal muscle during hypertrophic growth1

Postnatal muscle growth is achieved through hypertrophy of the muscle fibers and is impacted by the activity of satellite cells, the quiescent muscle stem cell. Several miRNAs are preferentially expressed in skeletal muscle and could provide a mechanism for increasing muscle hypertrophy through satellite cell proliferation and/or differentiation. The objectives of this study were to: 1) Characterize the miRNA transcriptome of the longissimus thoracis et lumborum muscle at several developmental timepoints [gestational d 85 (PN1), 110 (PN2), 133 (PN3), postnatal d 42 (PW1), 65 (PW2), 243 (MAT)] during muscle hypertrophy in lambs, and 2) examine miR-29a, identified in sequencing to be differentially regulated across development, loss of function on satellite cell proliferation and differentiation. Muscle fiber characteristics showed drastic increases (p < 0.0001) in fiber size and alterations in muscle fiber type occur during pre and postnatal development. miRNA sequencing comparisons were performed in developmental order (PN1 vs. PN2, PN2 vs. PN3, PN3 vs. PW1, PW1 vs. PW2, PW2 vs. MAT). There were 184 differentially expressed (P _adj < 0.05) miRNA, 142 unique miRNA, from all 5 comparisons made. The transitional stage (PN3 vs. PW1) had the largest number (115) of differentially expressed miRNA. Inhibition of miR-29a in satellite cell culture increased (p < 0.05) cell proliferation and differentiation capacity. Characterization of the miRNA transcriptome provides valuable insights into the miRNA involved in muscle fiber hypertrophy and the potential importance of the transitional period.

Muscle transcriptome provides the first insight into the dynamics of gene expression with progression of age in sheep

The dynamic synergy of genes and pathways in muscles in relation to age affects the muscle characteristics. Investigating the temporal changes in gene expression will help illustrate the molecular mechanisms underlying muscle development. Here we report the gene expression changes in skeletal muscles through successive age groups in Bandur, a meat type sheep of India. RNA sequencing data was generated from the longissimus thoracis muscles from four age groups, ranging from lamb to adult. Analysis of 20 highest expressed genes common across the groups revealed muscle protein, phosphorylation, acetylation, metal binding and transport as significant functions. Maximum differentiation was observed after 2.5–3 years on transition from lambs to adult. Transcriptional regulation by the TFAP2 transcription factors, IL-6 signaling and PI3K/AKT signaling pathways were enriched in younger animals. The gene-protein network demarcated key interactive genes involved in muscle development and proliferation that can be used as candidates for future research on improvement of muscle characteristics.

Potential Novel Biomarkers for Mastitis Diagnosis in Sheep

Simple Summary

Inflammation of the mammary gland (mastitis) is an important disease of dairy sheep. Mastitis management depends mainly on the diagnosis. Conventional diagnostic methods including somatic cell count, California Mastitis Test, and microbial culture have limitations. Therefore researchers are looking for new diagnostic biomarkers of mastitis including specific proteins produced by the liver in case of disease (acute phase proteins), unique genetic sequences (miRNAs), or antimicrobial peptides produced by immune cells during inflammation (cathelicidines).

Abstract

This review aims to characterize promising novel markers of ovine mastitis. Mastitis is considered as one of the primary factors for premature culling in dairy sheep and has noticeable financial, productional, and animal welfare-related implications. Furthermore, clinical, and subclinical mammary infections negatively affect milk yield and alter the milk composition, thereby leading to lowered quality of dairy products. It is, therefore, crucial to control and prevent mastitis through proper diagnosis, treatment or culling, and appropriate udder health management particularly at the end of the lactation period. The clinical form of mastitis is characterized by abnormalities in milk and mammary gland tissue alteration or systemic symptoms consequently causing minor diagnostic difficulties. However, to identify ewes with subclinical mastitis, laboratory diagnostics is crucial. Mastitis control is primarily dependent on determining somatic cell count (SCC) and the California Mastitis Test (CMT), which aim to detect the quantity of cells in the milk sample. The other useful diagnostic tool is microbial culture, which complements SCC and CMT. However, all mentioned diagnostic methods have their limitations and therefore novel biomarkers of ovine subclinical mastitis are highly desired. These sensitive indicators include acute-phase proteins, miRNA, and cathelicidins measurements, which could be determined in ovine serum and/or milk and in the future may become useful in early mastitis diagnostics as well as a preventive tool. This may contribute to increased detection of ovine mammary gland inflammation in sheep, especially in subclinical form, and consequently improves milk quality and quantity.