RNA sequencing as a powerful tool in searching for genes influencing health and performance traits of horses

RNA sequencing analysis revealed differentially expressed genes and their functional annotation in porcine longissimus dorsi muscle affected by dietary lysine restriction

The objective of this study was to investigate the effects of dietary lysine restriction on the global gene expression profile of skeletal muscle in growing pigs. Twelve crossbred (Yorkshire × Landrace) barrows (initial BW 22.6 ± 2.04 kg) were randomly assigned to two dietary treatments (LDD: a lysine-deficient diet; LAD: a lysine-adequate diet) according to a completely randomized experiment design (n = 6). After feeding for 8 weeks, skeletal muscle was sampled from the longissimus dorsi of individual pigs. The muscle total RNA was isolated and cDNA libraries were prepared for RNA sequencing (RNA-Seq) analysis. The RNA-Seq data obtained was then analyzed using the CLC Genomics Workbench to identify differentially expressed genes (DEGs). A total of 80 genes (padj ≤ 0.05) were differentially expressed in the longissimus dorsi muscle of the pigs fed LDD vs. LAD, of which 46 genes were downregulated and 34 genes were upregulated. Gene Ontology (GO) analysis of the DEGs (padj ≤ 0.05) for functional annotation identified those GO terms that are mostly associated with the molecular functions of structural molecules and metabolic enzymes (e.g., oxidoreductase and endopeptidase), biological process of acute-phase response, and amino acid metabolism including synthesis and degradation in the extracellular matrix region. Collectively, the results of this study have provided some novel insight regarding the molecular mechanisms of muscle growth that are associated with dietary lysine supply.

A one-health approach to identifying and mitigating the impact of endocrine disorders on human and equine athletes

Endocrinopathies affect multiple species in ever-increasing percentages of their populations, creating an opportunity to apply one-health approaches to determining creative preventative measures and therapies in athletes. Obesity and alterations in insulin and glucose dynamics are medical concerns that play a role in whole-body health and homeostasis in both horses and humans. The role and impact of endocrine disorders on the musculoskeletal, cardiovascular, and reproductive systems are of particular interest to the athlete. Elucidation of both physiologic and pathophysiologic mechanisms involved in disease processes, starting in utero, is important for development of prevention and treatment strategies for the health and well-being of all species. This review focuses on the unrecognized effects of endocrine disorders associated with the origins of metabolic disease; inflammation at the intersection of endocrine disease and related diseases in the musculoskeletal, cardiovascular, and reproductive systems; novel interventions; and diagnostics that are informed via multiomic and one-health approaches. Readers interested in further details on specific equine performance conditions associated with endocrine disease are invited to read the companion Currents in One Health by Manfredi et al, JAVMA, February 2023.

Genomic studies of milk-related traits in water buffalo (Bubalus bubalis) based on single-step genomic best linear unbiased prediction and random regression models

Genomic selection has been widely implemented in many livestock breeding programs, but it remains incipient in buffalo. Therefore, this study aimed to (1) estimate variance components incorporating genomic information in Murrah buffalo; (2) evaluate the performance of genomic prediction for milk-related traits using single- and multitrait random regression models (RRM) and the single-step genomic best linear unbiased prediction approach; and (3) estimate longitudinal SNP effects and candidate genes potentially associated with time-dependent variation in milk, fat, and protein yields, as well as somatic cell score (SCS) in multiple parities. The data used to estimate the genetic parameters consisted of a total of 323,140 test-day records. The average daily heritability estimates were moderate (0.35 ± 0.02 for milk yield, 0.22 ± 0.03 for fat yield, 0.42 ± 0.03 for protein yield, and 0.16 ± 0.03 for SCS). The highest heritability estimates, considering all traits studied, were observed between 20 and 280 d in milk (DIM). The genetic correlation estimates at different DIM among the evaluated traits ranged from -0.10 (156 to 185 DIM for SCS) to 0.61 (36 to 65 DIM for fat yield). In general, direct selection for any of the traits evaluated is expected to result in indirect genetic gains for milk yield, fat yield, and protein yield but also increase SCS at certain lactation stages, which is undesirable. The predicted RRM coefficients were used to derive the genomic estimated breeding values (GEBV) for each time point (from 5 to 305 DIM). In general, the tuning parameters evaluated when constructing the hybrid genomic relationship matrices had a small effect on the GEBV accuracy and a greater effect on the bias estimates. The SNP solutions were back-solved from the GEBV predicted from the Legendre random regression coefficients, which were then used to estimate the longitudinal SNP effects (from 5 to 305 DIM). The daily SNP effect for 3 different lactation stages were performed considering 3 different lactation stages for each trait and parity: from 5 to 70, from 71 to 150, and from 151 to 305 DIM. Important genomic regions related to the analyzed traits and parities that explain more than 0.50% of the total additive genetic variance were selected for further analyses of candidate genes. In general, similar potential candidate genes were found between traits, but our results suggest evidence of differential sets of candidate genes underlying the phenotypic expression of the traits across parities. These results contribute to a better understanding of the genetic architecture of milk production traits in dairy buffalo and reinforce the relevance of incorporating genomic information to genetically evaluate longitudinal traits in dairy buffalo. Furthermore, the candidate genes identified can be used as target genes in future functional genomics studies.

Development of EST-Molecular Markers from RNA Sequencing for Genetic Management and Identification of Growth Traits in Potato Grouper (Epinephelus tukula)

Simple Summary

The potato grouper is a novel aquaculture species in Taiwan. Due to the lack of genetic information concerning this species, we have developed molecular markers based on transcriptome sequencing and further characterized their association with gene diversity and growth traits of this species. Ultimately, these markers could be utilized as accurate and efficient tools for genetic management and marker-assisted selection of potato grouper with distinct growth traits.

Abstract

The accuracy and efficiency of marker-assisted selection (MAS) has been proven for economically critical aquaculture species. The potato grouper (Epinephelus tukula), a novel cultured grouper species in Taiwan, shows large potential in aquaculture because of its fast growth rate among other groupers. Because of the lack of genetic information for the potato grouper, the first transcriptome and expressed sequence tag (EST)-derived simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers were developed. Initially, the transcriptome was obtained from seven cDNA libraries by using the Illumina platform. De novo transcriptome of the potato grouper yielded 51.34 Gb and 111,490 unigenes. The EST-derived SSR and SNP markers were applied in genetic management, in parentage analysis, and to discover the functional markers of economic traits. The F₁ juveniles were identified as siblings from one pair of parents (80 broodstocks). Fast- and slow-growth individuals were analyzed using functional molecular markers and through their association with growth performance. The results revealed that two SNPs were correlated with growth traits. The transcriptome database obtained in this study and its derived SSR and SNP markers may be applied not only for MAS but also to maintain functional gene diversity in the novel cultured grouper.

Sequencing analysis of matrix metalloproteinase 7-induced genetic changes in Schwann cells

Previous research revealed the positive activity of matrix metalloproteinase 7 (MMP7) on migration and myelin regeneration of Schwann cells (SCs). However, understanding of the molecular changes and biological activities induced by increased amounts of MMP7 in SCs remains limited. To better understand the underlying molecular events, primary SCs were isolated from the sciatic nerve stump of newborn rats and cultured with 10 nM human MMP7 for 24 hours. The results of genetic testing were analyzed at a relatively relaxed threshold value (fold change ≥ 1.5 and P-value < 0.05). Upon MMP7 exposure, 149 genes were found to be upregulated in SCs, whereas 133 genes were downregulated. Gene Ontology analysis suggested that many differentially expressed molecules were related to cellular processes, single-organism processes, and metabolic processes. Kyoto Enrichment of Genes and Genomes pathway analysis further indicated the critical involvement of cell signaling and metabolism in MMP7-induced molecular regulation of SCs. Results of Ingenuity Pathway Analysis (IPA) also revealed that MMP7 regulates biological processes, molecular functions, cellular components, diseases and functions, biosynthesis, material metabolism, cell movement, and axon guidance. The outcomes of further analysis will deepen our comprehension of MMP7-induced biological changes in SCs. This study was approved by the Laboratory Animal Ethics Committee of Nantong University, China (approval No. 20190225-004) on February 27, 2019.

The Blood and Muscle Expression Pattern of the Equine TCAP Gene during the Race Track Training of Arabian Horses

Horse musculature has been shaped through evolution by environmental and human factors, which has resulted in several extraordinary adaptations to physical effort. Skeletal muscle plasticity results from the response to mechanical stimulation causing hypertrophy, where sarcomeres increase the muscle's cross-sectional area under the influence of contractile forces. The aim of the present study was the evaluation of transcript abundance of the telethonin (TCAP) gene, which is a part of the sarcomere macromolecular mechanosensory complex in the gluteus medius muscle, and the whole blood of Arabian horses during flat race training. The analysis, performed by quantitative PCR, showed an increase of TCAP transcripts in skeletal muscle. However, in whole blood, the transcript abundance decreased after the first stage of training and further increased after the second phase. The obtained results indicate a lack of similarity of TCAP gene expression in both tissues.

Combined cytotoxicity of aflatoxin B1 and deoxynivalenol to hepatoma HepG2/C3A cells

Co-occurrence of aflatoxin B1 (AFB1) and deoxynivalenol (DON) is one of the most common mycotoxin contaminations in cereal crops and food ingredients. However, the mechanism of their combined toxicity is poorly understood. In the current investigation, the hepatoma HepG2/C3A cell line was used to explore the combined cytotoxicity of AFB1 and DON. The values of IC50, based on a sulforhodamine B (SRB) assay, were 4.5 mM and 18.7 mM for DON and AFB1, respectively. The analysis of cytotoxicity endpoints using the combination index (CI) theory revealed that the changes of mitochondrial membrane permeability and ATP resulted from an additive cytotoxic effect (CI≈1) of AFB1 and DON. However, the endpoints of double strand DNA (ds-DNA), reactive oxygen species (ROS) and cell viability (SRB) were synergistically (CI<1) affected in a dose-dependent manner. RNA-seq analysis demonstrated a number of uniquely expressed genes in their combination (AFB1: 2.5 mM, DON: 0.56 mM), indicating a synergistic interaction between AFB1 and DON at a molecular level. Additional transcriptomics analysis showed that the endoplasmic reticulum stress-associated JNK/p38/MAPK pathway was induced by DON, whereas the p53 signalling pathway was activated by AFB1. The expression profiles of apoptosis-related genes caspase-3, Bax and Bcl-2 suggested a mitochondria-mediated apoptosis pathway that was shared between AFB1 and DON. Thus, different cytotoxicity pathways and their converging at the apoptotic process might be the mechanism of the additive/synergistic cytotoxicity of AFB1 and DON to HepG2/C3A cells.

Transcriptome Analysis in Domesticated Species: Challenges and Strategies

Domesticated species occupy a special place in the human world due to their economic and cultural value. In the era of genomic research, domesticated species provide unique advantages for investigation of diseases and complex phenotypes. RNA sequencing, or RNA-seq, has recently emerged as a new approach for studying transcriptional activity of the whole genome, changing the focus from individual genes to gene networks. RNA-seq analysis in domesticated species may complement genome-wide association studies of complex traits with economic importance or direct relevance to biomedical research. However, RNA-seq studies are more challenging in domesticated species than in model organisms. These challenges are at least in part associated with the lack of quality genome assemblies for some domesticated species and the absence of genome assemblies for others. In this review, we discuss strategies for analyzing RNA-seq data, focusing particularly on questions and examples relevant to domesticated species.