Selection of suitable reference genes for normalization of quantitative RT-PCR (RT-qPCR) expression data across twelve tissues of riverine buffaloes (Bubalus bubalis)

Hypoxia related genes modulate in similar fashion in skin fibroblast cells of yak (Bos grunniens) adapted to high altitude and native cows (Bos indicus) adapted to tropical climate during hypoxia stress

The present study was conducted to understand transcriptional response of skin fibroblast of yak (Bos grunniens) and cows of Bos indicus origin to hypoxia stress. Six primary fibroblast cell lines derived from three individuals each of Ladakhi yak (Bos grunniens) and Sahiwal cows (Bos indicus) were exposed to low oxygen concentration for a period of 24 h, 48 h and 72 h. The expression of 10 important genes known to regulate hypoxia response such as HIF1A, VEGFA, EPAS1, ATP1A1, GLUT1, HMOX1, ECE1, TNF-A, GPx and SOD were evaluated in fibroblast cells of Ladakhi yak (LAY-Fb) and Sahiwal cows (SAC-Fb) during pre- and post-hypoxia stress. A panel of 10 reference genes (GAPDH, RPL4, EEF1A1, RPS9, HPRT1, UXT, RPS23, B2M, RPS15, ACTB) were also evaluated for their expression stability to perform accurate normalization. The expression of HIF1A was significantly (p < 0.05) induced in both LAY-Fb (2.29-fold) and SAC-Fb (2.07-fold) after 24 h of hypoxia stress. The angiogenic (VEGFA), metabolic (GLUT1) and antioxidant genes (SOD and GPx) were also induced after 24 h of hypoxia stress. However, EPAS1 and ATP1A1 induced significantly (p < 0.05) after 48 h whereas, ECE1 expression induced significantly (p < 0.05) at 72 h after exposure to hypoxia. The TNF-alpha which is a pro-inflammatory gene induced significantly (p < 0.05) at 24 h in SAC-Fb and at 72 h in LAY-Fb. The induction of hypoxia associated genes indicated the utility of skin derived fibroblast as cellular model to evaluate transcriptome signatures post hypoxia stress in populations adapted to diverse altitudes.

Selection of reference genes for normalizing gene expression data across seasons in spermatozoa of water buffalo (Bubalus bubalis)

Selection of the most stably expressed reference genes is key to monitoring accurate target gene expression across any tissue or cell type. The mRNA in spermatozoa stores valuable information related to changes in spermatogenesis due to variations in environmental conditions, especially during heat stress, which affects various sperm functions. Semen quality in buffalo bulls is significantly influenced by the seasons. In the study, a panel of nine genes was evaluated to identify the most stably expressed internal control gene (ICG) for the normalization of real-time gene expression data generated across various seasons for Murrah buffalo bulls' spermatozoa. Sperm cells were purified from the semen samples collected during different seasons, with temperature-humidity index (THI) ranging from 80.80 ± 1.47 (hot summer) to 55.88 ± 1.98 (winter), using the BoviPure™ gradient purification method. The RNA isolated from the purified spermatozoa fraction was quality checked prior to reverse transcription and subjected to qPCR (quantitative real-time PCR) based expression analysis. An automated 'endoGene' pipeline was employed to apply the geNorm, NormFinder, and BestKeeper algorithms for data analysis. The result indicated that GAPDH and PP1A were the most stably expressed among the gene panel, whereas ATPSF1 and ACTB were the two least stable expressed reference genes. Further, the most suitable ICGs identified were validated by normalization of real time expression data of heat stress and sperm quality genes, HSFY2 and AKAP4, respectively. The genes identified would help in generating the most reliable results for the expression profiling of the genes dictating sperm quality and heat stress cope-up mechanism in buffalo spermatozoa, collected during different seasons.

Sequence characterization and comparative expression profile of buffalo WNT10B gene in adult and fetal tissues

In this study, Wingless-type MMTV (mouse mammary tumor virus) integration site family member (WNT10B) gene was sequence characterized in the Indian water buffalo. Sequence analysis revealed an open reading frame of 1176 nucleotides in buffalo, encoding 391 amino acids long protein. Nineteen nucleotide variations were observed between cattle and buffalo resulting in six amino acid changes. Phylogenetic analysis showed the clustering of ruminant species together. Real-time expression analysis of WNT10B in tissues collected from different organs of fetal and adult buffalo, revealed, the gene being abundantly expressed in the rumen and liver of the fetus. The fetal ovary, heart, kidney, lung, testis and mammary gland showed moderate expression, while in adult tissues, expression was high in the ovary, testis, brain, kidney, small intestine and liver, whereas lower expression was observed in the adult rumen. Significant differences in WNT10B expression levels were found for the brain, small intestine, testes, kidney, heart, rumen, and ovary when adult and fetal tissues were compared. A moderate level of genetic variation was found between cattle and buffalo WNT10B and expression patterns in a variety of tissues in adult buffalo implies that in addition to possible roles in adipogenesis and hematopoiesis, the WNT10B gene might be playing a significant role in other regulatory pathways as well.

Role of Different Members of the AGPAT Gene Family in Milk Fat Synthesis in Bubalus bubalis

During triacylglycerol synthesis, the acylglycerol-3-phosphate acyltransferase (AGPAT) family catalyzes the conversion of lysophosphatidic acid to phosphatidic acid and the acylation of sn-2 fatty acids. However, the catalytic activity of different AGPAT members is different. Therefore, this study aimed to investigate the mechanism through which different AGPATs affect the efficiency of TAG synthesis and fatty acid composition. The conservation of amino acid sequences and protein domains of the AGPAT family was analyzed, and the functions of AGPAT1, AGPAT3, and AGPAT4 genes in buffalo mammary epithelial cells (BMECs) were studied using RNA interference and gene overexpression. Prediction of the protein tertiary structure of the AGPAT family demonstrated that four conservative motifs (motif1, motif2, motif3, and motif6) formed a hydrophobic pocket in AGPAT proteins, except AGPAT6. According to cytological studies, AGPAT1, AGPAT3, and AGPAT4 were found to promote the synthesis and fatty acid compositions of triacylglycerol, especially UFA compositions of triacylglycerol, by regulating ACSL1, FASN, GPAM, DGAT2, and PPARG gene expression. This study provides new insights into the role of different AGPAT gene family members involved in TAG synthesis, and a reference for improving the fatty acid composition of milk.

Validation of reference gene stability for miRNA quantification by reverse transcription quantitative PCR in the peripheral blood of patients with COVID-19 critical illness

The COVID-19 pandemic has created an urgency to study the host gene response that leads to variable clinical presentations of the disease, particularly the critical illness response. miRNAs have been implicated in the mechanism of host immune dysregulation and thus hold potential as biomarkers and/or therapeutic agents with clinical application. Hence, further analyses of their altered expression in COVID-19 is warranted. An important basis for this is identifying appropriate reference genes for high quality expression analysis studies. In the current report, NanoString technology was used to study the expression of 798 miRNAs in the peripheral blood of 24 critically ill patients, 12 had COVID-19 and 12 were COVID-19 negative. A list of potentially stable candidate reference genes was generated that included ten miRNAs. The top six were analyzed using reverse transcription quantitative polymerase chain reaction (RT-qPCR) in a total of 41 patients so as to apply standard computational algorithms for validating reference genes, namely geNorm, NormFinder, BestKeeper and RefFinder. There was general agreement among all four algorithms in the ranking of four stable miRNAs: miR-186-5p, miR-148b-3p, miR-194-5p and miR-448. A detailed analysis of their output rankings led to the conclusion that miR-186-5p and miR-148b-3p are appropriate reference genes for miRNA expression studies using PaxGene tubes in the peripheral blood of patients critically ill with COVID-19 disease.

Identification of quantitative polymerase chain reaction reference genes suitable for normalising gene expression in the brain of normal and dystrophic mice and dogs

Background: In addition to progressive, debilitating muscle degeneration, ~50% of patients with Duchenne muscular dystrophy (DMD) have associated cognitive and behavioural disorders secondary to deficiency of dystrophin protein in the brain. The brain expresses a variety of dystrophin isoforms (Dp427, Dp140 and Dp71) whose functions remain to be fully elucidated. Detailed comparative analysis of gene expression in healthy and dystrophin-deficient brain is fundamental to understanding the functions of each isoform, and the consequences of their deficiency, with animal models representing a key tool in this endeavour. Reverse transcription quantitative real-time PCR (RT-qPCR) is a widely used method to study gene expression. However, accurate quantitative assessment requires normalisation of expression data using validated reference genes. The aim of this study was to identify a panel of suitable reference genes that can be used to normalise gene expression in the brain of healthy and dystrophic dogs and mice. Methods: Using the DE50-MD dog and mdx mouse models of DMD we performed RT-qPCR from fresh frozen brain tissue and employed the geNorm, BestKeeper and Normfinder algorithms to determine the stability of expression of a panel of candidate reference genes across healthy and dystrophic animals, and across different brain regions. Results: We show that SDHA, UBC and YWHAZ are suitable reference genes for normalising gene expression in healthy and dystrophic canine brain, and GAPDH, RPL13A and CYC1 in healthy and dystrophic murine brain. Notably, there was no overlap in the highest performing reference genes between the two species. Conclusions: Our findings suggest that gene expression normalisation is possible across six regions of the canine brain, and three regions of the murine brain. Our results should facilitate future work to study gene expression in the brains of normal and dystrophic dogs and mice and thus decipher the transcriptional consequences of dystrophin deficiency in the brain.

Abundance of Amino Acid Transporters and mTOR Pathway Components in the Gastrointestinal Tract of Lactating Holstein Cows

Data from non-ruminants indicate that amino acid (AA) transport into cells can regulate mTOR pathway activity and protein synthesis. Whether mTOR is expressed in the ruminant gastrointestinal tract (GIT) and how it may be related to AA transporters and the AA concentrations in the tissue is unknown. Ruminal papillae and the epithelia of the duodenum, jejunum, and ileum collected at slaughter from eight clinically healthy Holstein in mid-lactation were used. Metabolites and RNA were extracted from tissue for liquid chromatography–mass spectrometry and RT-qPCR analysis. The glycine and asparagine concentrations in the rumen were greater than those in the intestine (p < 0.05), but the concentrations of other AAs were greater in the small intestine than those in the rumen. Among the 20 AAs identified, the concentrations of glutamate, alanine, and glycine were the greatest. The mRNA abundances of AKT1 and MTOR were greater in the small intestine than those in the rumen (p < 0.05). Similarly, the SLC1A1, SLC6A6, SLC7A8, SLC38A1, SLC38A7, and SLC43A2 mRNA abundances were greater (p < 0.05) in the small intestine than those in the rumen. The mRNA abundances of SLC1A5, SLC3A2, and SLC7A5 were greater in the rumen than those in the small intestine (p < 0.05). Overall, the present study provides fundamental data on the relationship between mTOR pathway components and the transport of AAs in different sections of the gastrointestinal tract.

Identification of stably expressed Internal Control Genes (ICGs) for normalization of expression data in liver of C57BL/6 mice injected with beta casomorphins

In recent years, beta-casomorphin peptides (BCM7/BCM9) derived from the digestion of cow milk have drawn a lot of attention world over because of their proposed impact on human health. In order to evaluate the transcriptional modulation of target genes through RT-qPCR in response to these peptides, availability of appropriate reference or internal control genes (ICGs) will be the key. The present study was planned to identify a panel of stable ICGs in the liver tissue of C57BL/6 mice injected with BCM7/BCM9 cow milk peptides for 3 weeks. A total of ten candidate genes were evaluated as potential ICGs by assessing their expression stability using software suites; geNorm, NormFinder and BestKeeper. The suitability of the identified ICGs was validated by assessing the relative expression levels of target genes, HP and Cu/Zn SOD. Based on geNorm, PPIA and SDHA gene pair was identified to be most stably expressed in liver tissue during the animal trials. Similarly, NormFinder analysis also identified PPIA as the most stable gene. BestKeeper analysis showed crossing point SD value for all the genes in the acceptable range that is closer to 1. Overall, the study identified a panel of stable ICGs for reliable normalization of target genes expression data in mice liver tissues during BCM7/9 peptides trial.

Selection of species specific panel of reference genes in peripheral blood mononuclear cells of native livestock species adapted to trans-Himalayan region of Leh-Ladakh

The identification of appropriate references genes is an integral component of any gene expression-based study for getting accuracy and reliability in data interpretation. In this study, we evaluated the expression stability of 10 candidate reference genes (GAPDH, RPL4, EEF1A1, RPS9, HPRT1, UXT, RPS23, B2M, RPS15, ACTB) in peripheral blood mononuclear cells of livestock species that are adapted to high altitude hypoxia conditions of Leh-Ladakh. A total of 37 PBMCs samples from six native livestock species of Leh-Ladakh region such as Ladakhi cattle, Ladakhi yak, Ladakhi donkey, Chanthangi goat, Double hump cattle and Zanskar ponies were included in this study. The commonly used statistical algorithms such as geNorm, Normfinder, BestKeeper and RefFinder were employed to assess the stability of these RGs in all the livestock species. Our study has identified different panel of reference genes in each species; for example, EEF1A1, RPL4 in Ladakhi cattle; GAPDH, RPS9, ACTB in Ladakhi yak; HPRT1, B2M, ACTB in Ladakhi donkey; HPRT1, B2M, ACTB in Double hump camel, RPS9, HPRT1 in Changthangi goat, HPRT1 and ACTB in Zanskar ponies. To the best of our knowledge, this is the first systematic attempt to identify panel of RGs across different livestock species types adapted to high altitude hypoxia conditions. In future, the findings of the present study would be quite helpful in conducting any transcriptional studies to understand the molecular basis of high altitude adaptation of native livestock population of Leh-Ladakh.

Reference gene selection in bovine caruncular epithelial cells under pregnancy-associated hormones exposure

Examination of transcriptional regulation occurring during pregnancy establishment and maintenance requires the identification of endogenous reference genes characterized by high expression stability. Since the expression of some reference genes may be modulated by pregnancy-associated hormones, the goal of our study was to identify suitable reference genes unaffected by hormonal treatment. In our study bovine caruncular epithelial cells were subjected to progesterone, estrogen and prostaglandin F_2α treatment. Ten candidate reference genes (ACTR1A, CNOT11, HDAC1, HPRT1, RPL19, RPS9, SDHA, SUZ12, UXT and ZNF131) were evaluated with the use of four approaches (geNorm, NormFinder, BestKeeper, delta Ct). We found that RPS9 and SUZ12 displayed the highest expression stability in the tested material. Moreover, HPRT1 and SDHA were found inappropriate for RT-qPCR data normalization as they demonstrated the highest expression variability out of all candidates analysed. Hence geNorm calculations shown that the use of just two best-performing genes would be sufficient for obtaining reliable results, we propose that RPS9 and SUZ12 be used as suitable endogenous controls in future studies investigating gene expression in normal and compromised pregnancies.

Identification of Internal Reference Genes in Peripheral Blood Mononuclear Cells of Cattle Populations Adapted to Hot Arid Normoxia and Cold Arid Hypoxia Environments

To estimate gene expression in a reliable manner, quantitative real-time polymerase chain reaction data require normalisation using a panel of stably expressed reference genes (RGs). To date, information on an appropriate panel of RGs in cattle populations reared at cold arid high-altitude hypoxia and hot arid tropical normoxia environments is not available. Therefore, the present study was carried out to identify a panel of stably expressed RGs from 10 candidate genes (GAPDH, RPL4, EEF1A1, RPS9, HPRT1, UXT, HMBS, B2M, RPS15, and ACTB) in peripheral blood mononuclear cells (PBMCs) of cattle populations reared at cold arid high-altitude hypoxia and hot arid normoxia environments. Four different statistical algorithms: geNorm, NormFinder, BestKeeper, and RefFinder were used to assess the stability of these genes. A total of 30 blood samples were collected: six adult heifers each of Ladakhi (LAC) and Holstein Frisian crosses (HFX) and 4 Jersey (JYC) cows from cold arid high-altitude hypoxia environments (group I) and five adult heifers each of Sahiwal (SAC), Karan Fries (KFC), and Holstein Friesian (HFC) cows from hot arid normoxia environments (group II). Combined analysis of group I and group II resulted in identification of a panel of RGs like RPS9, RPS15, and GAPDH that could act as a useful resource to unravel the accurate transcriptional profile of PBMCs from diverse cattle populations adapted to distinct altitudes.

Evaluation of qPCR reference genes for taimen (Hucho taimen) under heat stress

As a powerful and attractive method for detecting gene expression, qRT-PCR has been broadly used in aquaculture research. Understanding the biology of taimen (Hucho taimen) has drawn increasing interest because of its ecological and economic value. Stable reference genes are required for the reliable quantification of gene expression, but such genes have not yet been optimized for taimen. In this study, the stability levels of 10 commonly used candidate reference genes were evaluated using geNorm, NormFinder, BestKeeper, and RefFinder. The expression levels of the 10 genes were detected using 240 samples from 48 experimental groups consisting of 40 individuals treated under four heat-stress conditions (18, 20, 22, and 24 °C) for 24 h and 26 °C for 4, 24, 48, and 72 h. Six tissues (blood, heart, brain, gill, skin, and liver) were collected from each individual. Ribosomal protein S29 (RPS29) and ribosomal protein L19 (RPL19) were the most stable genes among all of the samples, whereas 28S ribosomal RNA (28S rRNA), attachment region binding protein (ARBP), and 18S ribosomal RNA (18S rRNA) were the least stable. These results were verified by an expression analysis of taimen heat-stress genes (heat shock protein 60, hsp60, and heat shock protein 70, hsp70). In conclusion, RPS29 and RPL19 are the optimal reference genes for qRT-PCR analyses of taimen, irrespective of the tissue and experimental conditions. These results allow the reliable study of gene expression in taimen.

Reference gene screening for analyzing gene expression in the heart, liver, spleen, lung and kidney of forest musk deer

The screening of reference genes for real-time quantitative PCR (qPCR) in forest musk deer (FMD) tissue is of great significance to the basic research on FMD. However, there are few reports on the stability analysis of FMD reference genes so far. In this study, We used qPCR to detect the expression levels of 11 reference gene candidates (18S rRNA, beta-actin [ACTB], glyceraldehyde-3-phosphate dehydrogenase [GAPDH], TATA box-binding protein [TBP], hypoxanthine phosphoribosyltransferase 1 [HPRT1], tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta polypeptide [YWHAZ], hydroxymethylbilane synthase [HMBS], eukaryotic translation elongation factor 1 alpha 1 [EEF1A1], succinate dehydrogenase complex flavoprotein subunit A [SDHA], peptidylprolyl isomerase B [PPIB], and ubiquitin C [UBC]) in heart, liver, spleen, lung and kidney of FMD. After removing 18S rRNA on account of its high expression level, geNorm, NormFinder, BestKeeper and ΔCt algorithms were used to evaluate the expression stability of the remaining genes in the five organs, and further comprehensive ranking was calculated by RefFinder. According to the results, the selected reference genes with the most stable expression in the heart of FMD are SDHA and YWHAZ, while in the liver are ACTB and SDHA; in the spleen and lung are YWHAZ and HPRT1; in the kidney are YWHAZ and PPIB. The use of common reference genes in all five organs is not recommended. The analyses showed that tissue is an important variability factor in genes expression stability. Meanwhile, the result can be used as a reference for the selection of reference genes for qPCR in further study.

Molecular characterization and identification of SNPs in ATP1A1 isoform of sodium-phosphate adenosine triphosphatase across diverse breeds of riverine buffaloes (Bubalus bubalis) and Indian native cattle (Bos indicus): A plausible candidate gene for heat tolerance

In the present study, efforts were made to sequence characterize the selected exonic region (18–21) of ATP1A1 gene to identify variations/SNPs in different breeds of Indian riverine buffaloes and native cattle. The sequence characterization of selected intronic/exonic region (17–21) of ATP1A1 gene was carried out in a total of 120 samples which included 6 animals each of 8 buffalo breeds and 72 animals of 12 cattle breeds. Genomic DNA was extracted from the whole blood by enzymatic digestion using proteinase K using phenol:chloroform method.. Three sets of primers were designed using Primer3 software to amplify genomic region from intron 17 to intron 21 of ATP1A1 gene in both cattle and buffaloes. The amplified products were purified by enzymatic method and purified PCR products were sequenced using forward primers in an ABI 3100 Automated DNA Sequencer. The chromatogram of each sequence obtained was checked manually. Base calling was performed with Phred and contig assembly was done via Phrap/Cross_match/Swat tool available in the suite Codon code Aligner v. 3.5.1. The results revealed a total of 26 variations in exons 18–21 of ATP1A1 gene in riverine buffalo. Out of 26 variations, 6 (T27006876C, C27006599T, T27006345C, T27006330C, G27006309T and T27006240C) were distributed across 4 exonic regions and the remaining 20 were located in intronic region whereas in native cattle, only 2 SNPs were identified in exonic regions (18–21). SNP T27007767C was found to be a novel one in 18 intronic region of ATP1A1 gene in Indian cattle breeds. All 7 variations found in exonic region of both buffalo and cattle breeds were synonymous with the predicted changes in amino acids. The variations identified in ATP1A1 gene in the present study could be evaluated in future for their roles in heat tolerance trait in riverine buffaloes and native cattle.

Bovine Satellite Cells Isolated after 2 and 5 Days of Tissue Storage Maintain the Proliferative and Myogenic Capacity Needed for Cultured Meat Production

Cultured meat is an emerging alternative food technology which aims to deliver a more ethical, sustainable, and healthy muscle-tissue-derived food item compared to conventional meat. As start-up companies are rapidly forming and accelerating this technology, many aspects of this multi-faceted science have still not been investigated in academia. In this study, we investigated if bovine satellite cells with the ability to proliferate and undergo myogenic differentiation could be isolated after extended tissue storage, for the purpose of increasing the practicality for cultured meat production. Proliferation of bovine satellite cells isolated on the day of arrival or after 2 and 5 days of tissue storage were analyzed by metabolic and DNA-based assays, while their myogenic characteristics were investigated using RT-qPCR and immunofluorescence. Extended tissue storage up to 5 days did not negatively affect proliferation nor the ability to undergo fusion and create myosin heavy chain-positive myotubes. The expression patterns of myogenic and muscle-specific genes were also not affected after tissue storage. In fact, the data indicated a positive trend in terms of myogenic potential after tissue storage, although it was non-significant. These results suggest that the timeframe of which viable myogenic satellite cells can be isolated and used for cultured meat production can be greatly extended by proper tissue storage.

Selection of reference genes for quantitative real-time PCR normalization in European quail tissues

Coturniculture has been standing out as an industrial poultry activity in several countries around the world because of the several adaptive advantages of quails. Research that considers the analysis of gene expression can enhance this activity. This study aimed to analyze the stability of reference genes (RGs) in different tissues of quails (both males and females) for the recommendation of use in gene expression studies by the quantitative reverse transcription-polymerase chain reaction (RT-qPCR). The expression stability of ten RGs (ACTA1, ACTB, B2M, GAPDH, HMBS, SDHA, HPRT1, MRPS27, MRPS30, and RPL5) was analyzed in four tissues (breast muscle, abdominal fat, liver, and intestine), and assessed using the statistical tools geNorm, NormFinder, comparative ΔC_q method, and BestKeeper. The HPRT1 gene was the most stable in all quail tissues tested, followed by MRPS27 and MRPS30 in breast muscle, B2M and RPL5 in abdominal fat, HMBS and B2M in the liver, and RPL5 and HMBS in the intestine. These results may help studies using RT-qPCR assays to assess quail tissues from both sexes because they provide data on the most stable genes, which should be tested as candidate RGs for other experimental conditions.

Fenugreek Stimulates the Expression of Genes Involved in Milk Synthesis and Milk Flow through Modulation of Insulin/GH/IGF-1 Axis and Oxytocin Secretion

We previously demonstrated galactagogue effect of fenugreek in a rat model of lactation challenge, foreshadowing its use in women's breastfeeding management. To assess longitudinal molecular mechanisms involved in milk synthesis/secretion in dams submitted to fenugreek supplementation, inguinal mammary, pituitary glands and plasma were isolated in forty-three rats nursing large 12 pups-litters and assigned to either a control (CTL) or a fenugreek-supplemented (FEN) diet during lactation. RT-PCR were performed at days 12 and 18 of lactation (L12 and L18) and the first day of involution (Inv1) to measure the relative expression of genes related to both milk synthesis and its regulation in the mammary gland and lactogenic hormones in the pituitary gland. Plasma hormone concentrations were measured by ELISA. FEN diet induced 2- to 3-times higher fold change in relative expression of several genes related to macronutrient synthesis (Fasn, Acaca, Fabp3, B4galt1, Lalba and Csn2) and energy metabolism (Cpt1a, Acads) and in IGF-1 receptor in mammary gland, mainly at L12. Pituitary oxytocin expression and plasma insulin concentration (+77.1%) were also significantly increased. Altogether, these findings suggest fenugreek might extend duration of peak milk synthesis through modulation of the insulin/GH/IGF-1 axis and increase milk ejection by activation of oxytocin secretion.

Comparative transcriptome analysis reveals that PCK1 is a potential gene affecting IMF deposition in buffalo

BACKGROUND

In China, although buffaloes are abundant, beef is mainly obtained from cattle, and this preference is mainly attributed to the low intramuscular fat (IMF) content of buffalo. Genetic factors are an important driver that affects IMF deposition.

RESULTS

To reveal the intrinsic factors responsible for the low IMF content of buffalo, mRNA expression patterns in muscle and adipose tissue between buffalo and cattle were characterized by RNA sequencing analysis. The IMF content in Nanyang cattle was higher than that in Xinyang buffalo. A total of 1566 mRNAs expressed in adipose tissue showed differential expression between the longissimus dorsi muscles of buffalo and cattle. Functional annotation suggested a difference in the glycolysis/gluconeogenesis pathway between the two species. The results of RT-qPCR analysis and gain-of-function experiments confirmed the positive association between the IMF content and phosphoenolpyruvate carboxykinase 1 (PCK1) expression in buffalo. In both mouse C2C12 cells and cultured bovine myocytes, the activity of the PCK1 promoter in buffalo is lower than that in cattle. However, in mouse 3T3-L1 adipocytes and cultured bovine adipocytes, the activity of PCK1 in buffalo promoter is higher than that in cattle.

CONCLUSIONS

These results indicate the important role of PCK1 in buffalo IMF deposition and illustrate the differences between buffalo and cattle promoter activity that drive PCK1 expression. This research helps to establish a foundation for further studies investigating IMF deposition in buffalo.

Catalogue of stage-specific transcripts in Ixodes ricinus and their potential functions during the tick life-cycle

Background

The castor bean tick Ixodes ricinus is an important vector of several clinically important diseases, whose prevalence increases with accelerating global climate changes. Characterization of a tick life-cycle is thus of great importance. However, researchers mainly focus on specific organs of fed life stages, while early development of this tick species is largely neglected.

Methods

In an attempt to better understand the life-cycle of this widespread arthropod parasite, we sequenced the transcriptomes of four life stages (egg, larva, nymph and adult female), including unfed and partially blood-fed individuals. To enable a more reliable identification of transcripts and their comparison in all five transcriptome libraries, we validated an improved-fit set of five I. ricinus-specific reference genes for internal standard normalization of our transcriptomes. Then, we mapped biological functions to transcripts identified in different life stages (clusters) to elucidate life stage-specific processes. Finally, we drew conclusions from the functional enrichment of these clusters specifically assigned to each transcriptome, also in the context of recently published transcriptomic studies in ticks.

Results

We found that reproduction-related transcripts are present in both fed nymphs and fed females, underlining the poorly documented importance of ovaries as moulting regulators in ticks. Additionally, we identified transposase transcripts in tick eggs suggesting elevated transposition during embryogenesis, co-activated with factors driving developmental regulation of gene expression. Our findings also highlight the importance of the regulation of energetic metabolism in tick eggs during embryonic development and glutamate metabolism in nymphs.

Conclusions

Our study presents novel insights into stage-specific transcriptomes of I. ricinus and extends the current knowledge of this medically important pathogen, especially in the early phases of its development.

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.

Inter-genus gene expression analysis in livestock fibroblasts using reference gene validation based upon a multi-species primer set

Quantitative reverse transcription PCR (RT-qPCR) remains as an accurate approach for gene expression analysis but requires labor-intensive validation of reference genes using species-specific primers. To ease such demand, the aim was to design and test a multi-species primer set to validate reference genes for inter-genus RT-qPCR gene expression analysis. Primers were designed for ten housekeeping genes using transcript sequences of various livestock species. All ten gene transcripts were detected by RT-PCR in Bos taurus (cattle), Bubalus bubalis (buffaloes), Capra hircus (goats), and Ovis aries (sheep) cDNA. Primer efficiency was attained for eight reference genes using B. taurus—O. aries fibroblast cDNA (95.54–98.39%). The RT-qPCR data normalization was carried out for B. taurus vs. O. aries relative gene expression using Bestkeeper, GeNorm, Norm-finder, Delta CT method, and RefFinder algorithms. Validation of inter-genus RT-qPCR showed up-regulation of TLR4 and ZFX gene transcripts in B. taurus fibroblasts, irrespectively of normalization conditions (two, three, or four reference genes). In silico search in mammalian transcriptomes showed that the multi-species primer set is expected to amplify transcripts of at least two distinct loci in 114 species, and 79 species would be covered by six or more primers. Hence, a multi-species primer set allows for inter-genus gene expression analysis between O. aries and B. taurus fibroblasts and further reveals species-specific gene transcript abundance of key transcription factors.

Identification of qPCR reference genes suitable for normalizing gene expression in the mdx mouse model of Duchenne muscular dystrophy

The mdx mouse is the most widely-used animal model of the human disease Duchenne muscular dystrophy, and quantitative PCR analysis of gene expression in the muscles of this animal plays a key role in the study of pathogenesis and disease progression and in evaluation of potential therapeutic interventions. Normalization to appropriate stably-expressed reference genes is essential for accurate quantitative measurement, but determination of such genes is challenging: healthy and dystrophic muscles present very different transcriptional environments, further altering with disease progression and muscle use, raising the possibility that no single gene or combination of genes may be stable under all experimental comparative scenarios. Despite the pedigree of this animal model, this problem remains unaddressed. The aim of this work was therefore to comprehensively assess reference gene suitability in the muscles of healthy and dystrophic mice, identifying reference genes appropriate for specific experimental comparisons, and determining whether an essentially universally-applicable set of genes exists. Using a large sample collection comprising multiple muscles (including the tibialis anterior, diaphragm and heart muscles) taken from healthy and mdx mice at three disease-relevant ages, and a panel of sixteen candidate reference genes (FBXO38, FBXW2, MON2, ZFP91, HTATSF1, GAPDH, ACTB, 18S, CDC40, SDHA, RPL13a, CSNK2A2, AP3D1, PAK1IP1, B2M and HPRT1), we used the geNorm, BestKeeper and Normfinder algorithms to identify genes that were stable under multiple possible comparative scenarios. We reveal that no single gene is stable under all conditions, but a normalization factor derived from multiple genes (RPL13a, CSNK2A2, AP3D1 and the widely-used ACTB) appears suitable for normalizing gene expression in both healthy and dystrophic mouse muscle regardless of muscle type or animal age. We further show that other popular reference genes, including GAPDH, are markedly disease- or muscle-type correlated. This study demonstrates the importance of empirical reference gene identification, and should serve as a valuable resource for investigators wishing to study gene expression in mdx mice.

Identification of suitable reference genes for real-time qPCR in homocysteine-treated human umbilical vein endothelial cells

The imbalance in homocysteine (Hcy) metabolism has been implicated in the pathogenesis of human diseases, including cardiovascular and neurodegenerative disorders. When attempting to identify gene expression profiles using quantitative real-time reverse transcription polymerase chain reaction (RT-qPCR), the selection of suitable reference genes is important. Here, the expression levels of 10 commonly used reference genes were assessed for normalization of RT-qPCR in Hcy-treated human umbilical vein endothelial cells (HUVECs) and control cells. The suitability of eight selected candidate genes was comparatively analyzed across the tested samples and separately ranked by four programs, geNorm, NormFinder, BestKeeper, and the ΔCt method. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was the most stable gene in the final ranking using the RankAggreg package. Surprisingly, the β-actin (ACTB) levels decreased significantly in Hcy-treated HUVECs compared with control HUVECs (P<0.05), and further study indicated that Hcy suppressed the expression of ACTB by upregulating the miR-145-5p level in Hcy-treated HUVECs. Our data suggest that GAPDH can be used as a reliable reference gene, while ACTB cannot; normalization of gene expression in RT-qPCR experiments in Hcy-treated HUVECs. The data, which identifies a suitable reference gene in Hcy-treated HUVECs, will contribute to the design of an effective and accurate method for quantitation of gene expression.