The Selection of Reference Genes for Quantitative Real-Time PCR in the Ashidan Yak Mammary Gland During Lactation and Dry Period

Transcriptome-Based Evaluation of Optimal Reference Genes for Quantitative Real-Time PCR in Yak Stomach throughout the Growth Cycle

Efficient nutritional assimilation and energy metabolism in the stomachs of yaks contribute to their adaption to harsh environments. Accurate gene expression profile analysis will help further reveal the molecular mechanism of nutrient and energy metabolism in the yak stomach. RT-qPCR is regarded as an accurate and dependable method for analyzing gene expression. The selection of reference genes is essential to obtain meaningful RT-qPCR results, especially in longitudinal gene expression studies of tissues and organs. Our objective was to select and validate optimal reference genes from across the transcriptome as internal controls for longitudinal gene expression studies in the yak stomach. In this study, 15 candidate reference genes (CRGs) were determined according to transcriptome sequencing (RNA-seq) results and the previous literature. The expression levels of these 15 CRGs were quantified using RT-qPCR in the yak stomach, including the rumen, reticulum, omasum and abomasum at five stages: 0 days, 20 days, 60 days, 15 months and three years old (adult). Subsequently, the expression stabilities of these 15 CRGs were evaluated via four algorithms: geNorm, NormFinder, BestKeeper and the comparative C_T method. Furthermore, RefFinder was employed to obtain a comprehensive ranking of the stability of CRGs. The analysis results indicate that RPS15, MRPL39 and RPS23 are the most stable genes in the yak stomach throughout the growth cycle. In addition, to verify the reliability of the selected CRGs, the relative expression levels of HMGCS2 were quantified via RT-qPCR using the three most stable or the three least stable CRGs. Overall, we recommend combining RPS15, MRPL39 and RPS23 as reference genes for the normalization of RT-qPCR data in the yak stomach throughout the growth cycle.

Copy Number Variations in the MICALL2 and MOGAT2 Genes Are Associated with Ashidan Yak Growth Traits

Copy number variations (CNVs) are a result of genomic rearrangement affecting DNA regions over 1 kb in length, and can include inversions, translocations, deletions, and duplications. The molecule interacting with CasL-like protein 2 (MICALL2) gene is primarily associated with mitochondrial protein targeting and exhibits predicted stress fiber colocalization. The monoacylglycerol O-acyltransferase 2 (MOGAT2) gene encodes an enzyme responsible for catalyzing diacylglycerol synthesis from 2-monoacylglycerol and fatty acyl-CoA. For this study, blood samples were obtained from 315 yaks, and the body weight, body length, withers height, and chest girth of these animals were measured at 6, 12, 18, and 30 months of age. Genomic DNA was harvested from the collected blood samples, and CNVs in these samples were detected by qPCR. The resultant data were compared using ANOVAs, revealing significant associations between MICALL2 gene CNVs and body weight at 6 months of age (p < 0.05), body length and chest girth at 30 months of age (p < 0.05), and withers height at 18 months of age (p < 0.01) in Ashidan yaks. Similarly, MOGAT2 CNVs were significantly associated with body weight at 6 and 30 months of age (p < 0.05), and with withers height at 18 months of age (p < 0.01) in these Ashidan yaks. MICALL2 and MOGAT2 gene expression was further analyzed in yak tissue samples, revealing that MICALL2 was most highly expressed in the adipose tissue, whereas MOGAT2 was most highly expressed in the lung. These results thus confirmed the relationship between CNVs in the MICALL2 and MOGAT2 genes and Ashidan yak growth traits, providing a valuable gene locus that can be leveraged for future marker-assisted yak breeding efforts.

Evaluation of candidate reference genes for quantitative RTqPCR analysis in goldfish (Carassius auratus L.) in healthy and CyHV-2 infected fish

The accuracy of quantitative real time PCR (RTqPCR) can be attained only when a suitable reference gene is used. The gene expression for a particular gene may vary within different cells at different conditions. Hence, the suitability and stability of various potential reference genes have to be determined for expression studies. In this study, we have examined the potential of four different reference genes including β-Actin (ACTB), 18S ribosomal RNA (18S), glyceraldehyde-3P-dehydrogenase (GAPDH), and elongation factor 1 alpha (EF1A_A) in seven different tissues including gill, liver, kidney, spleen, heart, muscle and intestine of goldfish (Carassius auratus). The housekeeping genes were analyzed from healthy fish and in CyHV-2 challenged fish. Based upon the real time PCR results the gene expression varied among the genes and in tissues. The expression levels of the housekeeping genes were then compared and evaluated with the RefFinder web tool which analyses results using four different algorithms - BestKeeper, delta Ct, geNorm and NormFinder. EF1A_A was ranked to be the best gene in healthy fish by BestKeeper and geNorm analysis. The delta Ct and NormFinder algorithm have found 18S to be a stable gene in healthy fish but 18S was given to be least expressed in challenged fish. ACTB was also given as a stable gene by geNorm analysis in both healthy and challenged fish. Also, in CyHV-2 challenged fish, EF1A_A was identified as the best gene by all the three analysis except by BestKeeper analysis, where it has ranked GADPH as the best housekeeping gene. Expression of the four candidate reference genes differed across all tissue types tested, inferring that a thorough study of the reference genes is necessary for cross tissue comparison. These results can be further used in the immune gene response study of goldfish infected with any viral pathogen to develop better health strategies in the disease management of goldfish aquaculture.

Transcriptome Analysis Reveals the Potential Role of Long Non-coding RNAs in Mammary Gland of Yak During Lactation and Dry Period

The mammary gland is a remarkably dynamic organ of milk synthesis and secretion, and it experiences drastic structural and metabolic changes during the transition from dry periods to lactation, which involves the expression and regulation of numerous genes and regulatory factors. Long non-coding RNA (lncRNA) has considered as a novel type of regulatory factors involved in a variety of biological processes. However, their role in the lactation cycle of yak is still poorly understood. To reveal the involved mechanism, Ribo-zero RNA sequencing was employed to profile the lncRNA transcriptome in mammary tissue samples from yak at two physiological stages, namely lactation (LP) and dry period (DP). Notably, 1,599 lncRNA transcripts were identified through four rigorous steps and filtered through protein-coding ability. A total of 59 lncRNAs showed significantly different expression between two stages. Accordingly, the results of qRT-PCR were consistent with that of the transcriptome data. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that target genes of differentially expressed lncRNAs (DELs) were involved in pathways related to lactation, such as ECM-receptor interaction, PI3K-Akt signaling pathway, biosynthesis of amino acids and focal adhesion etc. Finally, we constructed a lncRNA-gene regulatory network containing some well known candidate genes for milk yield and quality traits. This is the first study to demonstrate a global profile of lncRNA expression in the mammary gland of yak. These results contribute to a valuable resource for future genetic and molecular studies on improving milk yield and quality, and help us to gain a better understanding of the molecular mechanisms underlying lactogenesis and mammary gland development of yak.

Reference gene selection and myosin heavy chain (MyHC) isoform expression in muscle tissues of domestic yak (Bos grunniens)

Domestic yak (Bos grunniens) is the most crucial livestock in the Qinghai-Tibetan Plateau, providing meat and other necessities for local people. The skeletal muscle of adult livestock is composed of muscle fibers, and fiber composition in muscle has influence on meat qualities, such as tenderness, pH, and color. Real-time quantitative polymerase chain reaction (RT-qPCR) is a powerful tool to evaluate the gene expression of muscle fiber, but the normalization of the data depends on the stability of expressed reference genes. Unfortunately, there is no consensus for an ideal reference gene for data normalization in muscle tissues of yak. In this study, we aimed to assess the stability of 14 commonly used candidate reference genes by using five algorithms (GeNorm, NormFinder, BestKeeper, Delat Ct and Refinder). Our results suggested UXT and PRL13A were the most stable reference genes, while the most commonly used reference gene, GAPDH, was most variably expressed across different muscle tissues. We also found that the extensor digitorum lateralis (EDL), trapezius pars thoracica (TPT), and psoas major (PM) muscle had the higher content of type I muscle fibers and the lowest content of type IIB muscle fibers, while gluteobiceps (GB) muscle had the highest content of type IIB muscle fibers. Our study provides the suitable reference genes for accurate analysis of yak muscle fiber composition.

Validation of Suitable Reference Genes for Gene Expression Studies on Yak Testis Development

Simple Summary

Yak (Bos grunniens) provides life materials for herdsmen in high-plateau areas. Improving their low fertility is necessary to meet the demands of the development of the yak industry. The testis is an important organ for male fertility, its development is controlled by a large number of genes. Using real-time reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) to explore the quantitative expression of genes can provide insights for illuminating the molecular mechanisms of testis development, but the RT-qPCR data are influenced by the stability of reference genes (RGs). Unfortunately, no available RGs can normalize the gene expression in yak testis development. In this study, the expression stability of 13 candidate genes in yak testis at different developmental stages was evaluated using five different pieces of software. The results showed that the TATA box-binding protein (TBP) and ubiquitously expressed transcript protein (UXT) exhibited high stability across various developmental stages, TBP and hydroxymethylbilane synthase (HMBS) were the most stably expressed genes in immature stages, and mitochondrial ribosomal protein L39 (MRPL39) and TBP exhibited the most stable expression across mature stages. This study provided suitable RGs for gene expression studies in yak testis development.

Abstract

Testis has an important function in male reproduction. Its development is regulated by a large number of genes. The real-time reserve transcriptase-quantitative polymerase chain reaction (RT-qPCR) is a useful tool to evaluate the gene expression levels. However, unsuitable reference genes (RGs) can cause the misinterpretation of gene expression levels. Unfortunately, the ideal RGs for yak testis development are yet to be studied. In this study, 13 commonly used RGs were selected to identify the most stable RGs in yak testis at four different developmental stages, including two immature stages (6 months and 18 months) and two mature stages (30 months and 6 years). This study used GeNorm, NormFinder, BestKeeper, ∆Ct, and RefFinder programs to evaluate the stability of 13 candidate genes. The results of RefFinder showed that the stabilities of TATA box-binding protein (TBP) and ubiquitously expressed transcript protein (UXT) were ranked the top two across all developmental stages. TBP and hydroxymethylbilane synthase (HMBS) were stably expressed in immature stages, while mitochondrial ribosomal protein L39 (MRPL39) and TBP had higher stability than other candidate genes in mature stages. This study provided valuable information for gene expression studies to assist further investigation on the molecular mechanisms in underlying yak testis development.

Myostatin (MSTN) Gene Indel Variation and Its Associations with Body Traits in Shaanbei White Cashmere Goat

Myostatin (MSTN) gene, also known as growth differentiation factor 8 (GDF8), is a member of the transforming growth factor-beta super-family and plays a negative role in muscle development. It acts as key points during pre- and post-natal life of amniotes that ultimately determine the overall muscle mass of animals. There are several studies that concentrate on the effect of a 5 bp insertion/deletion (indel) within the 5' untranslated region (5' UTR) of goat MSTN gene in goats. However, almost all sample sizes were below 150 individuals. Only in Boer goats, the sample sizes reached 482. Hence, whether the 5 bp indel was still associated with the growth traits of goats in large sample sizes which were more reliable is not clear. To find an effective and dependable DNA marker for goat rearing, we first enlarged the sample sizes (n = 1074, Shaanbei White Cashmere goat) which would enhance the robustness of the analysis and did the association analyses between the 5 bp indel and growth traits. Results uncovered that the 5 bp indel was significantly related to body height, height at hip cross, and chest width index (p < 0.05). In addition, individuals with DD genotype had a superior growing performance than those with the ID genotype. These findings suggested that the 5 bp indel in MSTN gene are significantly associated with growth traits and the specific genotype might be promising for maker-assisted selection (MAS) of goats.