Reference gene selection for gene expression analysis of oocytes collected from dairy cattle and buffaloes during winter and summer

Heat-Stress Impacts on Developing Bovine Oocytes: Unraveling Epigenetic Changes, Oxidative Stress, and Developmental Resilience

Extreme temperature during summer may lead to heat stress in cattle and compromise their productivity. It also poses detrimental impacts on the developmental capacity of bovine budding oocytes, which halt their fertility. To mitigate the adverse effects of heat stress, it is necessary to investigate the mechanisms through which it affects the developmental capacity of oocytes. The primary goal of this study was to investigate the impact of heat stress on the epigenetic modifications in bovine oocytes and embryos, as well as on oocyte developmental capacity, reactive oxygen species, mitochondrial membrane potential, apoptosis, transzonal projections, and gene expression levels. Our results showed that heat stress significantly reduced the expression levels of the epigenetic modifications from histone H1, histone H2A, histone H2B, histone H4, DNA methylation, and DNA hydroxymethylation at all stages of the oocyte and embryo. Similarly, heat stress significantly reduced cleavage rate, blastocyst rate, oocyte mitochondrial-membrane potential level, adenosine-triphosphate (ATP) level, mitochondrial DNA copy number, and transzonal projection level. It was also found that heat stress affected mitochondrial distribution in oocytes and significantly increased reactive oxygen species, apoptosis levels and mitochondrial autophagy levels. Our findings suggest that heat stress significantly impacts the expression levels of genes related to oocyte developmental ability, the cytoskeleton, mitochondrial function, and epigenetic modification, lowering their competence during the summer season.

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.

Oocyte competence and gene expression in parthenogenetic produced embryos from repeat breeder and normally fertile buffaloes (Bubalus bubalis) raised in sub-humid tropical climate

The reproductive management of the buffalo species still faces several unresolved problems, which directly affect the productivity of the herd, one of them being the presence of repeat breeder females. Given this scenario, this study aimed to verify the developmental competence of oocytes obtained from repeat breeder females and submitted to parthenogenetic activation. In addition, embryo gene expression was compared to normally fertile females. Murrah buffaloes were divided into two groups: repeat breeder (RB, n = 8) and normally fertile or control (CR, n = 7). Cumulus-oocyte complexes (COCs) were aspirated by transvaginal ovum pick-up from estrus synchronized females. The COCs were submitted to IVM for 24 h, and subsequently, the oocytes were activated using ionomycin, followed by 6-DMAP. Afterwards, the presumptive parthenotes were cultured for six or seven days in a microenvironment of 5 % CO2, 5 % O2, and 90 % N2 at 38.5 °C. The expression of OCT4, GLUT1, BCL2 and TFAM genes from blastocysts was evaluated. The overall COCs recovery rate was 70.9 % (190/268). The maturation (57.8 vs 71.1), cleavage (45.2 vs 62.2) and blastocyst (30.1 vs 45.9) rates did not differ (P > 0.05) between RB and CR females, respectively. Similarly, no significant difference (P > 0.05) was observed for the expression of studied genes in both RB and CR females. In conclusion, oocytes obtained from RB were as developmentally competent as those collected from CR females, with similar energy metabolism and in vitro development capacity. Thus, the low fertility rate of repeat breeder buffaloes, when compared to normal cyclic females, must be due to subsequent events to the blastocyst stage.

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.

Micro-vibration results in vitro-derived bovine blastocysts with greater cryotolerance, epigenetic abnormalities, and a massive transcriptional change

Much effort has been employed to improve the quality of embryos obtained by in vitro production (IVP) given the relevance of this technology to current livestock systems. In this context, dynamic IVP systems have proved beneficial to the embryo once they mimic fluid flows and mechanical forces resulting from the movement of ciliated cells and muscle contraction in the reproductive tract. In the present study, we sought to confirm these initial findings as well as assess potential molecular consequences to the embryo by applying micro-vibration (45 Hz for 5 s once per 60 min) during both oocyte maturation and embryo culture in cattle. As a result, micro-vibration led to lower incidence of apoptosis in blastocysts following vitrification-thawing. Further analyses revealed epigenetic and transcriptional changes in blastocysts derived from the micro-vibration treatment, with a total of 502 differentially expressed genes. Enrichment analyses linked differentially expressed genes to 'Oxidative phosphorylation', 'Cytokine-cytokine receptor interaction', and 'Signaling pathways regulating pluripotency of stem cells'. Yet, a meta-analysis indicated that the transcriptional changes induced by micro-vibration were not toward that of in vivo-derived embryos. In conclusion, micro-vibration increases the cryoresistance of bovine embryos, but caution should be taken given the unclear consequences of epigenetic and transcriptional abnormalities induced by the treatment.

Autophagy deficiency abolishes liver mitochondrial DNA segregation

Mutations in the mitochondrial genome (mtDNA) are ubiquitous in humans and can lead to a broad spectrum of disorders. However, due to the presence of multiple mtDNA molecules in the cell, co-existence of mutant and wild-type mtDNAs (termed heteroplasmy) can mask disease phenotype unless a threshold of mutant molecules is reached. Importantly, the mutant mtDNA level can change across lifespan as mtDNA segregates in an allele- and cell-specific fashion, potentially leading to disease. Segregation of mtDNA is mainly evident in hepatic cells, resulting in an age-dependent increase of mtDNA variants, including non-synonymous potentially deleterious mutations. Here we modeled mtDNA segregation using a well-established heteroplasmic mouse line with mtDNA of NZB/BINJ and C57BL/6N origin on a C57BL/6N nuclear background. This mouse line showed a pronounced age-dependent NZB mtDNA accumulation in the liver, thus leading to enhanced respiration capacity per mtDNA molecule. Remarkably, liver-specific atg7 (autophagy related 7) knockout abolished NZB mtDNA accumulat ion, resulting in close-to-neutral mtDNA segregation through development into adulthood. prkn (parkin RBR E3 ubiquitin protein ligase) knockout also partially prevented NZB mtDNA accumulation in the liver, but to a lesser extent. Hence, we propose that age-related liver mtDNA segregation is a consequence of macroautophagic clearance of the less-fit mtDNA. Considering that NZB/BINJ and C57BL/6N mtDNAs have a level of divergence comparable to that between human Eurasian and African mtDNAs, these findings have potential implications for humans, including the safe use of mitochondrial replacement therapy.Abbreviations: Apob: apolipoprotein B; Atg1: autophagy-related 1; Atg7: autophagy related 7; Atp5a1: ATP synthase, H+ transporting, mitochondrial F1 complex, alpha subunit 1; BL6: C57BL/6N mouse strain; BNIP3: BCL2/adenovirus E1B interacting protein 3; FCCP: carbonyl cyanide 4-(trifluoromethoxy)phenylhydrazone; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; MAP1LC3A: microtubule-associated protein 1 light chain 3 alpha; MAP1LC3B: microtubule-associated protein 1 light chain 3 beta; mt-Atp8: mitochondrially encoded ATP synthase 8; MT-CO1: mitochondrially encoded cytochrome c oxidase I; MT-CO2: mitochondrially encoded cytochrome c oxidase II; mt-Co3: mitochondrially encoded cytochrome c oxidase III; mt-Cytb: mitochondrially encoded cytochrome b; mtDNA: mitochondrial DNA; MUL1: mitochondrial ubiquitin ligase activator of NFKB 1; nDNA: nuclear DNA; Ndufa9: NADH:ubiquinone oxireductase subunit A9; NDUFB8: NADH:ubiquinone oxireductase subunit B8; Nnt: nicotinamide nucleotide transhydrogenase; NZB: NZB/BINJ mouse strain; OXPHOS: oxidative phosphorylation; PINK1: PTEN induced putative kinase 1; Polg2: polymerase (DNA directed), gamma 2, accessory subunit; Ppara: peroxisome proliferator activated receptor alpha; Ppia: peptidylprolyl isomerase A; Prkn: parkin RBR E3 ubiquitin protein ligase; P10: post-natal day 10; P21: post-natal day 21; P100: post-natal day 100; qPCR: quantitative polymerase chain reaction; Rpl19: ribosomal protein L19; Rps18: ribosomal protein S18; SD: standard deviation; SEM: standard error of the mean; SDHB: succinate dehydrogenase complex, subunit B, iron sulfur (Ip); SQSTM1: sequestosome 1; Ssbp1: single-stranded DNA binding protein 1; TFAM: transcription factor A, mitochondrial; Tfb1m: transcription factor B1, mitochondrial; Tfb2m: transcription factor B2, mitochondrial; TOMM20: translocase of outer mitochondrial membrane 20; UQCRC2: ubiquinol cytochrome c reductase core protein 2; WT: wild-type.

Analysis of stability of reference genes for qPCR in bovine preadipocytes during proliferation and differentiation in vitro

The stability of internal reference genes is crucial to the reliability of gene expression results using real-time fluorescence quantitative PCR (qRT-PCR). Inappropriate reference genes may lead to inaccurate results or even wrong conclusions. This study aims to identify stable reference genes for analyzing the expression of proliferation-related and differentiation-inducing genes in bovine primary preadipocytes (BPPs) in vitro. In this study, the stability of 16 candidate internal reference genes (GAPDH, ACTB, PPIA, LRP10, HPRT1, YWHAZ, B2M, TBP, EIF3K, RPS9, UXT, 18S rRNA, RPLP0, MARVELD, EMD and RPS15A) for qRT-PCR at proliferation and differentiation stages of BPPs was investigated by three different algorithms (geNorm, NormFinder and BestKeeper). The expression of two marker genes, PCNA and LPL, was used to determine the validity of the candidate reference genes (RGs) at the proliferation and differentiation stages, respectively. The results showed that GAPDH and RPS15A were the most stable RGs in the proliferation of bovine primary preadipocyte, while PPIA was the least stable internal reference gene. RPLP0 and EIF3K were the most stable RGs in the differentiation induction of bovine primary preadipocyte, while GAPDH was the least stable internal reference gene. This study of RGs laid the foundation for subsequent research into the mechanism of proliferation and differentiation of BPPs in vitro using qRT-PCR.

Screening and validation of reference genes for qRT-PCR of bovine skeletal muscle-derived satellite cells

The accuracy of sixteen commonly used internal reference genes was assessed in skeletal muscle-derived satellite cells of Qinchuan cattle at different stages of proliferation and induction of differentiation to determine the most suitable ones. Quantitative real-time PCR and three commonly used algorithmic programs, GeNorm, NormFinder and BestKeeper, were used to evaluate the stability of expression of the candidate internal reference genes (GAPDH, ACTB, PPIA, LRP10, HPRT1, YWHAZ, B2M, TBP, EIF3K , RPS9, UXT, 18S rRNA, RPLP0, MARVELD, EMD and RPS15A) in skeletal muscle-derived satellite cells at 0, 12, 24, 36 and 48 h of growth and after differentiation for 0, 2, 4, 6 and 8 days. The expression of two satellite cell marker genes, CCNA2 and MYF5, was used for validation analysis. The results of the software analyses showed that GAPDH and RPS15A were the most stable reference gene combinations during in vitro proliferation of bovine skeletal muscle-derived satellite cells, RPS15A and RPS9 were the most stable reference gene combinations during in vitro induction of differentiation of the cells, and PPIA was the least stable reference gene during proliferation and differentiation and was not recommended. This study lays the foundation for the selection of reference genes for qRT-PCR during the proliferation and induction of differentiation of bovine skeletal muscle-derived satellite cells.

The microenvironment of ovarian follicles in fertile dairy cows is associated with high oocyte quality

We hypothesised that heifers and cows with positive genetic merit for fertility would have a follicular microenvironment that resulted in better quality oocytes. To test this, we compared cumulus cell-oocyte complexes (COC) and follicular fluid from preovulatory follicles of 36 Holstein-Friesian nulliparous heifers and 50 primiparous lactating cows with either positive (POS, +5%) or negative (NEG, -5%) fertility breeding values (FertBV). Established gene markers of oocyte quality were measured in individual cumulus cell masses and oocytes, and concentrations of amino acids, steroids, and metabolites were quantified in corresponding follicular fluid and plasma. The timing of visually detectable oestrus in NEG FertBV heifers was inconsistent with their stage of COC maturation. Retrospective analyses of oestrous activity data indicated that NEG FertBV heifers were sampled earlier. Their recovered COC were morphologically less mature and exhibited differential expression of genes that are associated with follicular maturation (lower levels of BMPR2) and protein processing (higher levels of HSP90B1). Despite consistent sampling times being achieved in the lactating cows, lower concentrations of serine, proline, methionine, isoleucine, and non-esterified fatty acids were present in follicular fluid from POS FertBV cows. This was associated with higher expression of gene biomarkers of good oocyte quality (VCAN, PDE8A) in COC recovered from POS FertBV cows. This study supports our hypothesis that the follicular microenvironment in lactating dairy cows with high genetic merit leads to COC with higher metabolic rates and oocytes of superior quality. Moreover, an additional stressor such as lactation is required for this difference to be pronounced.

Genome-Wide Identification of Reference Genes for Reverse-Transcription Quantitative PCR in Goat Rumen

Simple Summary

The rumen plays an essential role as a digestive organ and serves as the primary site of energy substrate absorption for the productive ruminants. Understanding gene expression profiles is necessary to explore the intrinsic regulatory mechanisms of rumen development in goats. The selection of suitable reference genes (RGs) was the primary assay before the real-time quantitative PCR (RT-qPCR). We identified sixteen genome-wide candidate RGs for normalization of gene expression assessments in goat rumen tissues. We demonstrate that the RGs selected (RPS4X and RPS6) were more stably expressed than the commonly used HKGs (ACTB and GAPDH) in goat rumen tissues, suggesting that the ribosomal protein gene family may be another source for the RG pool.

Abstract

As the largest chamber of the ruminant stomach, the rumen not only serves as the principal absorptive surface and nutrient transport pathway from the lumen into the animal, but also plays an important short-chain fatty acid (SCFA) metabolic role in addition to protective functions. Accurate characterization of the gene expression profiles of genes of interest is essential to the exploration of the intrinsic regulatory mechanisms of rumen development in goats. Thus, the selection of suitable reference genes (RGs) is an important prerequisite for real-time quantitative PCR (RT-qPCR). In the present study, 16 candidate RGs were identified from our previous transcriptome sequencing of caprine rumen tissues. The quantitative expressions of the candidate RGs were measured using the RT-qPCR method, and the expression stability of the RGs was assessed using the geNorm, NormFinder, and BestKeeper programs. GeNorm analysis showed that the M values were less than 0.5 for all the RGs except GAPT4, indicating that they were stably expressed in the rumen tissues throughout development. RPS4X and RPS6 were the two most stable RGs. Furthermore, the expressions of two randomly selected target genes (IGF1 and TOP2A), normalized by the selected most stable RGs (RPS4X and RPS6), were consistent with the results of RNA sequencing, while the use of GAPDH and ACTB as RGs resulted in altered profiles. Overall, RPS4X and RPS6 showed the highest expression stability and the lowest coefficients of variation, and could be used as the optimal reference combination for quantifying gene expression in rumen tissues via RT-qPCR analysis.

Cattle In Vitro Induced Pluripotent Stem Cells Generated and Maintained in 5 or 20% Oxygen and Different Supplementation

The event of cellular reprogramming into pluripotency is influenced by several factors, such as in vitro culture conditions (e.g., culture medium and oxygen concentration). Herein, bovine iPSCs (biPSCs) were generated in different levels of oxygen tension (5% or 20% of oxygen) and supplementation (bFGF or bFGF + LIF + 2i-bFL2i) to evaluate the efficiency of pluripotency induction and maintenance in vitro. Initial reprogramming was observed in all groups and bFL2i supplementation initially resulted in a superior number of colonies. However, bFL2i supplementation in low oxygen led to a loss of self-renewal and pluripotency maintenance. All clonal lines were positive for alkaline phosphatase; they expressed endogenous pluripotency-related genes SOX2, OCT4 and STELLA. However, expression was decreased throughout the passages without the influence of oxygen tension. GLUT1 and GLUT3 were upregulated by low oxygen. The biPSCs were immunofluorescence-positive stained for OCT4 and SOX2 and they formed embryoid bodies which differentiated in ectoderm and mesoderm (all groups), as well as endoderm (one line from bFL2i in high oxygen). Our study is the first to compare high and low oxygen environments during and after induced reprogramming in cattle. In our conditions, a low oxygen environment did not favor the pluripotency maintenance of biPSCs.

Use of four genes in exosomes as biomarkers for the identification of lung adenocarcinoma and lung squamous cell carcinoma

The determination of biomarkers in the blood specific for lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) is crucial for the selection of effective treatment strategies and the prediction of prognosis. The purpose of the present study was to analyze the differentially expressed genes (DEGs) in LUSC and LUAD from The Cancer Genome Atlas (TCGA) database. In order to identify the potential biomarkers for non-small cell lung cancer (NSCLC) for clinical diagnosis, bioinformatics was used to analyze the DEGs of two subtypes of NSCLC, LUAD and LUSC. Exosomes were isolated from the serum of patients with LUAD or LUSC and identified using transmission electron microscopy, nanoparticle tracking analysis and western blot analysis. A total of four differential exosomal mRNAs were selected for validation with serum samples from 70 patients with NSCLC via reverse transcription-quantitative polymerase chain reaction. Receiver operating characteristic curves were established to evaluate the clinical diagnostic value of four DEGs for patients with LUAD and LUSC. The analysis based on TCGA data revealed the DEGs in LUSC and LUAD: A total of 1,619 genes were differentially expressed in patients with LUSC and LUAD. DEGs analyzed by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses revealed that inflammation-related signaling pathways, such as complement pathways, and multiple autoimmune diseases, such as systemic lupus erythematosus and asthma were mainly enriched in LUAD. The cell cycle, Hippo signaling pathway, Rap1 signaling pathway and Wnt signaling pathway were the main signaling pathways enriched in LUSC. The combination of tumor protein P63 (TP63), keratin 5 (KRT5), CEA cell adhesion molecule 6 (CEACAM6) and surfactant protein B (SFTPB) improved the specificity and sensitivity in the diagnosis of different lung cancer subtypes. Exosomal TP63, KRT5, CEACAM6 and SFTPB mRNAs can thus be used as biomarkers to differentiate between LUSC and LUAD, and may provide a novel strategy for their differential diagnosis and treatment.

Impact of an acute heat shock during in vitro maturation on interleukin 6 and its associated receptor component transcripts in bovine cumulus-oocyte complexes

An acute heat stress event after the LH surge increased interleukin 6 (IL6) levels in the follicular fluid of the ovulatory follicle in hyperthermic cows. To examine direct consequences of a physiologically-relevant elevated temperature (41.0°C) on the cumulus-oocyte complex (COC), IL6 transcript abundance and related receptor components were evaluated throughout in vitro maturation. Heat-induced increases in IL6 were first noted at 4 hours of in vitro maturation (hIVM); peak levels occurred at 4.67 versus 6.44 hIVM for 41.0 and 38.5°C COCs, respectively (SEM = 0.23; P < 0.001). Peak IL6ST levels occurred at 6.95 versus 8.29 hIVM for 41.0 and 38.5°C, respectively (SEM = 0.23; P < 0.01). Transcript for LIF differed over time (P < 0.0001) but was not affected by 41.0°C exposure. Blastocyst development after performing IVF was not affected by 41.0°C exposure for 4 or 6 h. When limiting analysis to when IL6 was temporally produced, progesterone levels were only impacted by time and temperature (no interaction). Heat-induced shift in the temporal production of IL6 and IL6ST along with its impact on progesterone likely cooperate in heat-induced hastening of meiotic progression described by others.

Mice born to females with oocyte-specific deletion of mitofusin 2 have increased weight gain and impaired glucose homeostasis

Offspring born to obese and diabetic mothers are prone to metabolic diseases, a phenotype that has been linked to mitochondrial dysfunction and endoplasmic reticulum (ER) stress in oocytes. In addition, metabolic diseases impact the architecture and function of mitochondria-ER contact sites (MERCs), changes which associate with mitofusin 2 (MFN2) repression in muscle, liver and hypothalamic neurons. MFN2 is a potent modulator of mitochondrial metabolism and insulin signaling, with a key role in mitochondrial dynamics and tethering with the ER. Here, we investigated whether offspring born to mice with MFN2-deficient oocytes are prone to obesity and diabetes. Deletion of Mfn2 in oocytes resulted in a profound transcriptomic change, with evidence of impaired mitochondrial and ER function. Moreover, offspring born to females with oocyte-specific deletion of Mfn2 presented increased weight gain and glucose intolerance. This abnormal phenotype was linked to decreased insulinemia and defective insulin signaling, but not mitochondrial and ER defects in offspring liver and skeletal muscle. In conclusion, this study suggests a link between disrupted mitochondrial/ER function in oocytes and increased risk of metabolic diseases in the progeny. Future studies should determine whether MERC architecture and function are altered in oocytes from obese females, which might contribute toward transgenerational transmission of metabolic diseases.

TATA box binding protein and ribosomal protein 4 are suitable reference genes for normalization during quantitative polymerase chain reaction study in bovine mesenchymal stem cells

OBJECTIVE

Quantitative polymerase chain reaction (qPCR) has been extensively used in the field of mesenchymal stem cell (MSC) research to elucidate their characteristics and clinical potential by normalization of target genes against reference genes (RGs), which are believed to be stably expressed irrespective of various experimental conditions. However, the expression of RGs is also variable depending on the experimental conditions, which may lead to false or contradictory conclusions upon normalization. Due to the current lack of information for a clear list of stable RGs in bovine MSCs, we conducted this study to identify suitable RGs in bovine MSCs.

METHODS

The cycle threshold values of ten traditionally used RGs (18S ribosomal RNA [18S], beta-2-microglobulin [B2M], H2A histone family, member Z [H2A], peptidylprolyl isomerase A [PPIA], ribosomal protein 4 [RPL4], succinate dehydrogenase complex, subunit A [SDHA], beta actin [ACTB], glyceraldehyde-3-phosphate dehydrogenase [GAPDH], TATA box binding protein [TBP], and hypoxanthine phosphoribosyltrasnfrase1 [HPRT1]) in bovine bone marrow-derived MSCs (bBMMSCs) were validated for their stabilities using three types of RG evaluation algorithms (geNorm, Normfinder, and Bestkeeper). The effect of validated RGs was then verified by normalization of lineage-specific genes (fatty acid binding protein 4 [FABP4] and osteonectin [ON]) expressions during differentiations of bBMMSCs or POU class 5 homeobox 1 (OCT4) expression between bBMMSCs and dermal skins.

RESULTS

Based on the results obtained for the three most stable RGs from geNorm (TBP, RPL4, and H2A), Normfinder (TBP, RPL4, and SDHA), and Bestkeeper (TBP, RPL4, and SDHA), it was comprehensively determined that TBP and RPL4 were the most stable RGs in bBMMSCs. However, traditional RGs were suggested to be the least stable (18S) or moderately stable (GAPDH and ACTB) in bBMMSCs. Normalization of FABP4 or ON against TBP, RPL4, and 18S presented significant differences during differentiation of bBMMSCs. However, although significantly low expression of OCT4 was detected in dermal skins compared to that in bBMMSCs when TBP and RPL4 were used in normalization, normalization against 18S exhibited no significance.

CONCLUSION

This study proposes that TBP and RPL4 were suitable as stable RGs for qPCR study in bovine MSCs.

Frameshift Variant in Novel Adenosine-A1-Receptor Homolog Associated With Bovine Spastic Syndrome/Late-Onset Bovine Spastic Paresis in Holstein Sires

Since their first description almost 100 years ago, bovine spastic paresis (BSP) and bovine spastic syndrome (BSS) are assumed to be inherited neuronal-progressive diseases in cattle. Affected animals are characterized by (frequent) spasms primarily located in the hind limbs, accompanied by severe pain symptoms and reduced vigor, thus initiating premature slaughter or euthanasia. Due to the late onset of BSP and BSS and the massively decreased lifespan of modern cattle, the importance of these diseases is underestimated. In the present study, BSP/BSS-affected German Holstein breeding sires from artificial insemination centers were collected and pedigree analysis, genome-wide association studies, whole genome resequencing, protein-protein interaction network analysis, and protein-homology modeling were performed to elucidate the genetic background. The analysis of 46 affected and 213 control cattle revealed four significantly associated positions on chromosome 15 (BTA15), i.e., AC_000172.1:g.83465449A>G (-log₁₀P = 19.17), AC_000172.1:g.81871849C>T (-log₁₀P = 8.31), AC_000172.1:g.81872621A>T (-log₁₀P = 6.81), and AC_000172.1:g.81872661G>C (-log₁₀P = 6.42). Two additional loci were significantly associated located on BTA8 and BTA19, i.e., AC_000165.1:g.71177788T>C and AC_000176.1:g.30140977T>G, respectively. Whole genome resequencing of five affected individuals and six unaffected relatives (two fathers, two mothers, a half sibling, and a full sibling) belonging to three different not directly related families was performed. After filtering, a homozygous loss of function variant was identified in the affected cattle, causing a frameshift in the so far unknown gene locus LOC100848076 encoding an adenosine-A1-receptor homolog. An allele frequency of the variant of 0.74 was determined in 3,093 samples of the 1000 Bull Genomes Project.

Identification of reliable reference genes for expression studies in maternal reproductive tissues and foetal tissues of pregnant cows

The relationship between the conceptus and the maternal uterine environment is crucial for the successful establishment and maintenance of pregnancy in cattle. Gene expression analysis of the conceptus and maternal reproductive tissues is a favourable method to assess the embryonic maternal interaction. The reliability of the commonly used method reverse transcription-quantitative polymerase chain reaction (RT-qPCR) depends on proper normalization to stable reference genes (RGs). The objective of this study was to determine the expression stability of 10 potential RGs in maternal reproductive tissues and foetal tissues, and to analyse the effect of RG selection on the calculation of the relative expression of target genes. The expression stability of 10 potential RGs was analysed in eight different tissues from three pregnant dairy cows. Three programs-GeNorm, NormFinder and Bestkeeper-were used to identify the best RGs. According to all three programs, the most stable RG was CNOT11, whereas the least stable RGs were GAPDH and HPRT1. GeNorm analysis showed that a combination of five RGs (SDHA, PPIA, CNOT11, RPS9 and RPL19) was necessary for appropriate data normalization. However, NormFinder analysis indicated that the combination of CNOT11 and PPIA was the most suitable. When target genes were normalized to these RGs, the relative expression of the Radical S-adenosyl methionine domain containing 2 gene was not affected by the choice of RGs, whereas a large difference was observed in the expression profile of the Nuclear erythroid2-related factor 2 gene between the most stable and least stable RGs. The results indicate that careful selection of RGs is crucial under different conditions, especially for target genes with relatively small fold changes. Furthermore, the results provide useful information for the selection of RGs for evaluating genes affecting bovine reproduction.

Catalytic inhibition of H3K9me2 writers disturbs epigenetic marks during bovine nuclear reprogramming

Orchestrated events, including extensive changes in epigenetic marks, allow a somatic nucleus to become totipotent after transfer into an oocyte, a process termed nuclear reprogramming. Recently, several strategies have been applied in order to improve reprogramming efficiency, mainly focused on removing repressive epigenetic marks such as histone methylation from the somatic nucleus. Herein we used the specific and non-toxic chemical probe UNC0638 to inhibit the catalytic activity of the histone methyltransferases EHMT1 and EHMT2. Either the donor cell (before reconstruction) or the early embryo was exposed to the probe to assess its effect on developmental rates and epigenetic marks. First, we showed that the treatment of bovine fibroblasts with UNC0638 did mitigate the levels of H3K9me2. Moreover, H3K9me2 levels were decreased in cloned embryos regardless of treating either donor cells or early embryos with UNC0638. Additional epigenetic marks such as H3K9me3, 5mC, and 5hmC were also affected by the UNC0638 treatment. Therefore, the use of UNC0638 did diminish the levels of H3K9me2 and H3K9me3 in SCNT-derived blastocysts, but this was unable to improve their preimplantation development. These results indicate that the specific reduction of H3K9me2 by inhibiting EHMT1/2 during nuclear reprogramming impacts the levels of H3K9me3, 5mC, and 5hmC in preimplantation bovine embryos.

Validation of Real-time PCR Reference Genes of Muscle Metabolism in Harvested Spiny-Cheek Crayfish (Faxonius limosus) Exposed to Seasonal Variation

Real-time quantitative reverse transcription PCR (RT-qPCR) is a sensitive and broadly used technique of assessing gene activity. To obtain a reliable result, stably expressed reference genes are essential for normalization of transcripts in various samples. To our knowledge, this is the first systematic analysis of reference genes for normalization of RT-qPCR data in spiny-cheek crayfish (Faxonius limosus). In this study, expression of five candidate reference genes (actb, β-actin; gapdh, glyceraldehyde-3-phosphate dehydrogenase; eif, eukaryotic translation initiation factor 5a; ef-1α, elongation factor-1α; and tub, α-tubulin) in muscle samples from male and female F. limosus in spring and autumn was analyzed. Additionally, the most stable reference genes were used for accurate normalization of five target genes, i.e., tnnc, troponin c; ak, arginine kinase; fr, ferritin; ccbp-23, crustacean calcium-binding protein 23; and actinsk8, skeletal muscle actin 8. Results obtained using the geNorm and NormFinder algorithms showed high consistency, and differences in the activity of the selected actb with eif genes were successfully identified. The spring and autumn activities of the target genes (except ak) in the muscle tissue of males and females differed significantly, showing that both sexes are immensely involved in an array of breeding behaviors in spring, and females intensively recover in the autumn season. Characterization of first reference genes in spiny-cheek crayfish will facilitate more accurate and reliable expression studies in this key species.

Use of embryo transfer to alleviate infertility caused by heat stress

Heat stress (HS) has a pronounced deleterious effect on fertility in dairy herds throughout the world, especially in hot and humid summer months in tropical and subtropical areas. Summer HS reduces feed intake and increases negative energy balance, induces changes in ovarian follicular dynamics, reduces estrus detection rates and alters oviductal function leading to fertilization failure and early embryonic death. Furthermore, oocytes harvested from lactating cows during summer HS have a decreased ability to develop to the blastocyst stage after in vitro fertilization when compared with oocytes harvested during winter. The present manuscript describes the detrimental effect of HS on reproduction, with emphasis on preovulatory oocytes and carry over effects of HS on embryo development and P/AI. Embryo transfer (ET) has been an effective tool to reestablish fertility during HS because it bypasses the damage to the oocyte and early embryo caused by hyperthermia. Therefore, a management strategy to maintain increased fertility throughout the year would be to produce embryos during the cooler months, when oocyte quality is greater, and use them to produce pregnancies during the periods of HS, when oocyte quality is compromised. However, this strategy only can be implemented using cryopreserved embryos, what is still limiting. During the warmer months, the use of heifers or non-lactating cows as oocyte or embryo donors may facilitate embryo production, mainly because of the lesser deleterious effects of HS comparing to lactating cows. Also, genetic selection of donors for thermoregulation ability is one potential strategy to mitigate effects of HS and increase embryo production during the warmer months. These alternatives enable the transference of fresh embryos with more efficiency during HS periods. Additionally, the application of timed ET protocols, which avoid the need for estrus detection in recipients, has facilitated management and improved the efficiency of ET programs during HS.

Mitofusin 1 is required for oocyte growth and communication with follicular somatic cells

Mitochondrial function, largely regulated by the dynamics of this organelle, is inextricably linked to the oocyte health. In comparison with most somatic cells, mitochondria in oocytes are smaller and rounder in appearance, suggesting limited fusion. The functional implications of this distinct morphology, and how changes in the mitochondrial shape translate to mitochondrial function in oogenesis is little understood. We, therefore, asked whether the pro-fusion proteins mitofusins 1 (MFN1) and 2 (MFN2) are required for the oocyte development. Here we show that oocyte-specific deletion of Mfn1, but not Mfn2, prevents the oocyte growth and ovulation due to a block in folliculogenesis. We pinpoint the loss of oocyte growth and ovulation to impaired PI3K-Akt signaling and disrupted oocyte-somatic cell communication. In support, the double loss of Mfn1 and Mfn2 partially rescues the impaired PI3K-Akt signaling and defects in oocyte development secondary to the single loss of Mfn1. Together, this work demonstrates that the mitochondrial function influences the cellular signaling during the oocyte development, and highlights the importance of distinct, nonredundant roles of MFN1 and MFN2 in oogenesis.

Hepatic mRNA expression of enzymes associated with progesterone metabolism and its impact on ovarian and endocrine responses in Nelore (Bos indicus) and Holstein (Bos taurus) heifers with differing feed intakes

The aim of this study was to evaluate circulating progesterone concentration (P4), LH pulsatility and ovarian follicular dynamics in Nelore (B. indicus) and Holstein (B. taurus) heifers under high (HDMI) and low (LDMI) dry matter/energy intakes. In addition, the effects of dry matter/energy intake and breed on hepatic expression of six genes associated with P4 metabolism (AKR1C4, AKR1D1, CYP3A4, CYP2C19, SRD5A1, and SRD5A3) was evaluated. Heifers received an intravaginal P4 device (1 g), 2 mg of estradiol benzoate (EB) i.m. and 500 μg of PGF_2α at the begging of the synchronization protocol (D0). Eight days later, the P4 device was removed and all heifers received 1 mg of EB 24h later. Regardless of dry matter/energy intake, the number of recruited follicles was greater in Nelore than in Holstein heifers. In contrast, the maximum diameter of the dominant follicle was greater in Holstein than in Nelore heifers. Circulating P4 concentrations were greater in Nelore than in Holstein from D2 to D9, and in heifers receiving LDMI than those receiving HDMI from D1 to D8 of hormonal protocol. In addition, Holstein heifers had greater LH pulsatility and area under the curve of LH peaks compared to Nelore heifers. However, no effects were observed for LH values between feed intake levels. Interestingly, Holstein heifers had higher expression of SRD5A1, AKR1C4, AKR1D1 than Nelore heifers; whereas, for Nelore heifers, only the expression of CYP3A4 was higher compared to Holstein heifers. In conclusion, there are important differences in the follicular dynamics, circulating P4 and LH pulsatility concentrations that need to be considered during synchronization protocols for Nelore and Holstein breeds. More importantly, these differences appear to be at least partially modulated by the level of feed intake and the contrasting enzyme system in the liver involved with P4 metabolism between these cattle breeds.

Interdigital Hyperplasia in Holstein Cattle Is Associated With a Missense Mutation in the Signal Peptide Region of the Tyrosine-Protein Kinase Transmembrane Receptor Gene

Bovine interdigital hyperplasia (IH) is a typical disease of the foot with varying prevalence depending on age, breed, and environmental factors resulting in different degrees of lameness. In studies based on assessments of claw health status at time of hoof trimming and applying genetic-statistical models to analyze this data, IH consistently exhibits high estimates of heritability in the range of 0.30–0.40. Although some studies have identified chromosomal regions that could possibly harbor causative genes, a clear identification of molecular causes for IH is lacking. While analyzing the large database of claw health status as documented at time of hoof trimming, we identified one herd with extreme prevalence of IH of > 50% of affected Holstein dairy cows. This herd subsequently was chosen as the object of a detailed study. A total of n = 91 cows was assessed and revealed a prevalence of 59.3% and 38.5% for IH cases, documented as “one-sided” or “two-sided”, respectively. Cows were genotyped using the BovineSNP50 BeadChip. A genome wide association study revealed two significantly associated chromosomal positions (-log₁₀P = 5.57) on bovine chromosome 8 (BTA8) located in intron 5 and downstream of the receptor tyrosine kinase-like orphan receptor 2 (ROR2) gene. As ROR2 plays a key role in ossification of the distal limbs and is associated with brachydactylies in humans, it was a reasonable candidate for IH. A comparative sequencing of the ROR2 gene between cases and controls revealed two missense variants in exon 1 (NC_037335.1:g.85,905,534T > A, ARS-UCD1.2) and exon 9 (NC_037335.1:g.86,140,379A > G, ARS-UCD1.2), respectively. Genotyping of both variants in the cohort of 91 cattle showed that the exon 1 variant (rs377953295) remained significantly associated with IH (p < 0.0001) as a risk factor of the disease. This variant resulted in an amino acid exchange (ENSBTAP00000053765.2:p.Trp9Arg) in the N-terminal region of the ROR2 signal peptide which is necessary for proper topology of the polypeptide during translocation. Quantification of ROR2 mRNA and ROR2 protein showed that the variant resulted in a significant suppression of ROR2 expression in homozygous affected compared to wild type and carrier cows.

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

Simple Summary

The Ashidan yak is a new cultivated breed which has polled characteristics and a mild temperament. Improving milk yield is an important aspect of a breeding program for this breed. The mammary gland undergoes dramatic physiological and metabolic changes during the transition from lactation to dry periods, which involves the expression and regulation of a great number of genes. Quantification of gene expression levels by real-time quantitative polymerase chain reaction (RT-qPCR) is important to reveal the molecular mechanisms of mammary gland development and lactation. The accuracy of RT-qPCR is strongly influenced by the expression stability of reference genes, however, a systematic approach for selecting reference genes used for analyzing gene expression of the Ashidan yak has not been developed. In this study, we selected reference genes and analyzed their expression stability at different physiological stages (lactation and dry period). We found the hydroxymethylbilane synthase gene (HMBS) and the tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide gene (YWHAZ) were the most stable genes of the mammary gland of the Ashidan yak. These results help to improve the accuracy of gene expression analysis and provide a basis for future functional studies of target gene expression in the mammary gland of the Ashidan yak.

Abstract

Investigating the critical genes related to milk synthesis is essential for the improvement of the milk yield of the yak. Real-time quantitative polymerase chain reaction (RT-qPCR) is a reliable and widely used method to measure and evaluate gene expression levels. Selection of suitable reference genes is mandatory to acquire accurate normalization of gene expression results from RT-qPCR. To select the most stable reference genes for reliable normalization of mRNA expression by RT-qPCR in the mammary gland of the Ashidan yak, we selected 16 candidate reference genes and analyzed their expression stability at different physiological stages (lactation and dry period). The expression stability of the candidate reference genes was assessed using geNorm, NormFinder, BestKeeper, Delta Ct, and RefFinder methods. The results showed that the hydroxymethylbilane synthase gene (HMBS) and the tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide gene (YWHAZ) were the most stable genes across all treatment samples. The reliability of selected reference genes was validated by normalizing relative expression of the lactation-related 60S ribosomal protein L35 gene (RPL35). The relative expression of RPL35 varied considerably according to the different reference genes. This work provides valuable information to further promote research in the molecular mechanisms involved in lactation and mammary gland development and provides a foundation for the improvement of the milk yield and quality of the Ashidan yak.

Validation of reference genes for gene expression studies in bovine oocytes and cumulus cells derived from in vitro maturation

Quantitative real-time PCR (qPCR) is a valuable tool for gene expression studies and it is necessary to choose an ideal endogenous reference gene for data normalization. This work studied a set of reference genes in oocytes and cumulus cells of COCs (Cumulus-Oocyte Complexes) that are suitable for relative gene expression analyses after in vitro maturation (IVM) in bovine. Immature COCs were collected from ovaries of Nelore cattle (Bos indicus) and submitted to IVM. MII oocytes and cumulus cells were subjected to RNA extraction, reverse transcription and preamplification of cDNA. The expression level of eight reference genes (ACTB, GADPH, B2M, H2AFZ, GUSB, HPRT1, PPIA, and TBP) was measured by real time PCR and analyzed by geNorm software. The gene stability measure (M) was calculated and the ideal number of reference genes (RGs) was determined by the V value (pairwise variation). For oocyte samples, two RGs were the ideal number for relative quantification: HPRT1 and B2M and for bovine cumulus samples four were indicated: HPRT1, PPIA, B2M, and TBP genes. The normalization of a non-reference target gene (SOD1) by these reference genes was shown to be considerably different from normalization by less stable reference genes. Our results strengthen the importance of choosing good normalizing genes in order to analyze gene expression under specific experimental conditions and we suggest the use of these RGs in oocytes and cumulus cells of bovine cattle in in vitro matured COCs.

Reference genes selection for real-time quantitative PCR analysis in mouse germinal vesicle oocytes

Reference gene selection in mouse oocytes is an important task required to perform further adequate analysis of target gene expression levels. In the current work we have analyzed expression stability of the seven most commonly used reference genes (Actb, Eef1e1, Gapdh, H2afz, Ppia, Rpl4 and Ubc) in mouse oocytes at the germinal vesicle (GV) stage. We have performed analysis of expression stability of the above-mentioned reference genes with the three most commonly used software tools: geNorm, BestKeeper and NormFinder. Taking into account the results obtained from all of these programmes Gapdh, Rpl4 and H2afz seem to be suitable candidate reference genes in GV oocytes of mouse.

Reference gene selection for quantitative gene expression analysis in black soldier fly (Hermetia illucens)

Hermetia illucens is an important resource insect for the conversion of organic waste. Quantitative PCR (qPCR) is the primary tool of gene expression analysis and a core technology of molecular biology research. Reference genes are essential for qPCR analysis; however, a stability analysis of H. illucens reference genes has not yet been carried out. To find suitable reference genes for normalizing gene expression data, the stability of eight housekeeping genes (including ATP6V1A, RPL8, EF1, Tubulin, TBP, GAPDH, Actin and RP49) was investigated under both biotic (developmental stages, tissues and sex) and abiotic (heavy metals, food, antibiotics) conditions. Gene expression data were analysed by geNorm, NormFinder, BestKeeper, and ΔCt programs. A set of specific reference genes was recommended for each experimental condition using the results of RefFinder synthesis analysis. This study offers a solid foundation for further studies of the molecular biology of H. illucens.

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 optimal endogenous reference RNAs for RT-qPCR normalization in hindgut of rat models with anorectal malformations

Background

Quantitative real-time polymerase chain reaction (RT-qPCR) is a sensitive method for quantifying mRNA abundance. With relative expression analysis, however, reliable data output is dependent on stably expressed reference genes across the samples being studied. In anorectal malformations (ARMs), there is limited data on the selection of appropriate reference genes.

Purpose

This study was aimed to investigate the optimal reference genes for PCR in ARM rat models.

Methods

We selected 15 commonly used reference genes (Rps18, Actb, B2m, Gapdh, Ppia, Hprt1, Pgk1, Ywhaz, Tbp, Ubc, Rps16, Rpl13a, Rplp1, Sdha, and Hmbs) as candidate reference genes and detected their mRNA expression in ARM samples by RT-qPCR. The expression stability and variability of these transcripts were subsequently evaluated using four methods (geNorm, NormFinder, comparative ΔCt, and BestKeeper).

Results

The abundance of the candidate reference genes was qualified by RT-qPCR and the cycle threshold (Ct) values ranged between 14.07 (Rplp1) and 21.89 (Sdha). In the overall candidate genes, different variations existed across the different algorithms. A comprehensive analysis revealed that Rpl13a ranked first among the relatively stable genes, followed by Ywhaz, Rps18, Sdha, and Hmbs.

Conclusions

The most stable reference genes for RT-qPCR were Rpl13a, Ywhaz, and Rps18 in ETU-induced ARMs in rat fetus. This study provided a foundation for reference gene selection for future gene expression analyses.

Far-off and close-up feeding levels affect immunological performance in grazing dairy cows during the transition period

During the peripartum period, dairy cows often have signs of inflammation. Various stresses, including infectious and metabolic diseases, have been discussed as causative for this inflammation. In this study, expression profiles for 17 immune markers were measured in whole blood preparations from 78 dairy cows over a time frame starting 1 wk before calving to 4 wk after calving. Additionally, the effects of far-off and close-up feeding on immune function of dairy cows during the peripartum period were investigated. Cows were assigned to 1 of 2 feeding levels in late lactation to achieve a low and high BCS at the time of dry-off (approximately 4.25 and 5.0 on a 10-point scale). Following dry-off, both herds were managed to achieve a BCS of 5.0 one month before calving; this involved controlled feeding (i.e., maintenance) and over-feeding of ME during the far-off dry period. Within each far-off feeding-level treatment, cows were offered 65, 90, or 120% of their precalving ME requirements for 3 wk precalving in a 2 × 3 factorial arrangement. Analysis of gene expression profiles from blood cells revealed effects of time indicating that the transition cow's immune system counteracts the peripartum inflammation, whereas later postcalving it becomes activated to provide protection against postpartum infections. Far-off feeding affected (P < 0.05) the expression of 2 of the investigated genes at calving. Interleukin-6 (IL-6) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) expression in unstimulated, peripheral leukocytes were lower (P < 0.05) in animals from the Far-Off_Over-fed group compared with the Far-Off_Control-fed group. Close-up feeding had several effects on gene expression, indicating that immune function in Feed120 animals was distinct from the Feed90 and Feed65. In conclusion, feeding management precalving becomes an important intervention to ensure immunocompetence at and after calving.

Metabolic gene expression and epigenetic effects of the ketone body β-hydroxybutyrate on H3K9ac in bovine cells, oocytes and embryos

The rapid decline in fertility that has been occurring to high-producing dairy cows in the past 50 years seems to be associated with metabolic disturbances such as ketosis, supporting the need for research to improve our understanding of the relations among the diet, metabolism and embryonic development. Recently, the ketone body β-hydroxybutyrate (BOHB) was demonstrated to be a potent inhibitor of histone deacetylases (HDACs). Herein, we performed a series of experiments aiming to investigate the epigenetic effects of BOHB on histone acetylation in somatic cells, cumulus-oocyte complexes (COCs) and somatic cell nuclear transfer (SCNT) embryos. Treatment with BOHB does not increase histone acetylation in cells but stimulates genes associated with ketolysis and master regulators of metabolism. We further demonstrated that maturing COCs with high levels of BOHB does not affect their maturation rate or histone acetylation but increases the expression of PPARA in cumulus cells. Treatment of somatic cell nuclear transfer zygotes with BOHB causes hyperacetylation, which is maintained until the blastocyst stage, causing enhanced FOXO3A expression and blastocyst production. Our data shed light on the epigenetic mechanisms caused by BOHB in bovine cells and embryos and provide a better understanding of the connection between nutrition and reproduction.

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

Selection of reference genes has become an integral step in any real time quantitative PCR (RT-qPCR) based expression studies. The importance of this study stems from the fact that riverine buffaloes are major dairy species of Indian sub-continent and the information generated here will be of great interest to the investigators engaged in functional genomic studies of this important livestock species. In this study, an effort was made to evaluate a panel of 10 candidate reference genes (glyceraldehyde 3-phosphate dehydrogenase (GAPDH), beta- actin (ACTB), ubiquitously expressed transcript (UXT), ribosomal protein S15 (RPS15), ribosomal protein L-4 (RPL4), ribosomal protein S9 (RPS9), ribosomal protein S23 (RPS23), hydroxymethylbilane synthase (HMBS), β2 Microglobulin (β2M) and eukaryotic translation elongation factor 1 alpha 1 (EEF1A1) across 12 tissues (mammary gland, kidney, spleen, liver, heart, intestine, ovary, lung, muscle, brain, subcutaneous fat and testis) of riverine buffaloes. In addition to overall analysis, tissue wise evaluation of expression stability of individual RG was also performed. Three different algorithms provided in geNorm, NormFinder and BestKeeper softwares were used to evaluate the stability of 10 potential reference genes from different functional classes. The M-value given by geNorm ranged from 0.9797 (RPS9 and UXT) to 1.7362 (RPS15). From the most stable to the least stable, genes were ranked as: UXT/RPS9> RPL4> RPS23> EEF1A1> ACTB> HMBS> GAPDH> B2M> RPS15. While NormFinder analysis ranked the genes as: UXT> RPS23> RPL4> RPS9> EEF1A1> HMBS> ACTB> β2M> GAPDH> RPS15. Based on the crossing point SD value and range of fold change expression, BestKeeper analysis ranked the genes as: RPS9> RPS23/UXT> RPL4> GAPDH> EEF1A1> ACTB> HMBS> β2M> RPS15. Overall the study has identified RPS23, RPS9, RPL4 and UXT genes to be the most stable and appropriate RGs that could be utilized for normalization of transcriptional data in various tissues of buffaloes. This manuscript thus provide useful information on panel of reference genes that could be helpful for researchers conducting functional genomic studies in riverine buffaloes.

Reference gene selection for RT-qPCR analysis in Harmonia axyridis, a global invasive lady beetle

Harmonia axyridis is a voracious predator, a biological control agent, and one of the world most invasive insect species. The advent of next-generation sequencing platforms has propelled entomological research into the genomics and post-genomics era. Real-time quantitative PCR (RT-qPCR), a primary tool for gene expression analysis, is a core technique governs the genomic research. The selection of internal reference genes, however, can significantly impact the interpretation of RT-qPCR results. The overall goal of this study is to identify the reference genes in the highly invasive H. axyridis. Our central hypothesis is that the suitable reference genes for RT-qPCR analysis can be selected from housekeeping genes. To test this hypothesis, the stability of nine housekeeping genes, including 18S, 28S, ACTB, ATP1A1, GAPDH, HSP70, HSP90, RP49, and ATP6V1A, were investigated under both biotic (developmental time, tissue and sex), and abiotic (temperature, photoperiod, in vivo RNAi) conditions. Gene expression profiles were analyzed by geNorm, Normfinder, BestKeeper, and the ΔCt method. Our combined results recommend a specific set of reference genes for each experimental condition. With the recent influx of genomic information for H. axyridis, this study lays the foundation for an in-depth omics dissection of biological invasion in this emerging model.

Supplementation with small-extracellular vesicles from ovarian follicular fluid during in vitro production modulates bovine embryo development

Pregnancy success results from the interaction of multiple factors, among them are folliculogenesis and early embryonic development. Failure during these different processes can lead to difficulties in conception. Alternatives to overcome these problems are based on assisted reproductive techniques. Extracellular vesicles are cell-secreted vesicles present in different body fluids and contain bioactive materials, such as messenger RNA, microRNAs (miRNAs), and proteins. Thus, our hypothesis is that extracellular vesicles from follicular fluid from 3-6 mm ovarian follicles can modulate bovine embryo development in vitro. To test our hypothesis follicular fluid from bovine ovaries was aspirated and small-extracellular vesicles (<200 nm) were isolated for further analysis. Additionally, small-extracellular vesicles (EVs) were utilized for functional experiments investigating their role in modulating messenger RNA, microRNA as well as global DNA methylation and hydroxymethylation levels of bovine blastocysts. EVs from 3-6 mm follicles were used for RNA-seq and miRNA analysis. Functional annotation analysis of the EVs transcripts revealed messages related to chromatin remodeling and transcriptional regulation. EVs treatment during oocyte maturation and embryo development causes changes in blastocyst rates, as well as changes in the transcription levels of genes related to embryonic metabolism and development. Supplementation with EVs from 3-6 mm follicles during oocyte maturation and early embryo development (until the 4-cell stage) increased the levels of bta-miR-631 (enriched in EVs from 3-6 mm follicles) in embryos. Interestingly, the addition of EVs from 3-6 mm follicles induced changes in global DNA methylation and hydroxymethylation levels compared to embryos produced by the standard in vitro production system. Our results indicate that the supplementation of culture media with EVs isolated from the follicular fluid of 3-6 mm follicles during oocyte maturation and early embryo development can partially modify metabolic and developmental related genes as well as miRNA and global DNA methylation and hydroxymethylation, suggesting that EVs play an important role during oocyte maturation and early embryo development in vitro.

Fatty Acid Binding Protein 3 And Transzonal Projections Are Involved In Lipid Accumulation During In Vitro Maturation Of Bovine Oocytes

Oocytes that undergo in vitro maturation (IVM) are metabolically abnormal and accumulate excess lipid content. However, the mechanism of lipid accumulation and the role of cumulus cells in this process are unclear. Recently, it was shown that fatty acid binding proteins (FABPs) performed intra- and extracellular fatty acid transport. We postulated that FABP3 might be responsible for fatty acid transport from cumulus cells to the oocytes via transzonal projections (TZPs) during IVM. Transcript and protein levels of FABP3 were analyzed in both in vivo- and in vitro-matured cumulus-oocyte-complexes and were increased in IVM samples. Further analysis showed increased lipid content in oocytes and cumulus cells in IVM samples compared to in vivo-derived. We therefore speculated that altered traffic of fatty acids via FABP3 during IVM was the mechanism leading to the excess of lipids accumulated within IVM oocytes. Furthermore, we demonstrated an increase in FABP3 levels and lipid content during the first 9 h of IVM, further strengthening the possibility of fatty acid transport via FABP3 and TZPs. Additionally, disruptions of TZPs during IVM decreased lipid accumulation in oocytes. Our results shed light on a possible mechanism involving FABP3 and TZPs that causes excess lipid accumulation in oocytes during IVM.

Methods of RNA preparation affect mRNA abundance quantification of reference genes in pig maturing oocytes

To ensure accurate normalization and quantification of target RNA transcripts using reverse transcription quantitative polymerase chain reaction (RT-qPCR), most studies focus on the identification of stably expressed gene(s) as internal reference. However, RNA preparation methods could also be an important factor, especially for test samples of limited quantity (e.g. oocytes). In this study, we aimed to select appropriate reference gene(s), and evaluate the effect of RNA preparation methods on gene expression quantification in porcine oocytes and cumulus cells during in vitro maturation. Expression profiles of seven genes (GAPDH, 18S, YWHAG, BACT, RPL4, HPRT1 and PPIA) were examined, on RNA samples extracted from cumulus cells (RNeasy Kit) and oocytes (RNeasy Kit and Lysis Kit) during in vitro maturation, respectively. Interestingly, different RNA preparation methods were found to potentially affect the quantification of reference gene expression in pig oocytes cultured in vitro. After geNorm analyses, the most suitable genes for normalization were identified, GAPDH/18S for cumulus cells and YWHAG/BACT for oocytes, respectively. Thus, our results provide useful data and information on the selection of better reference genes and RNA preparation method for related functional studies.

Developmental competence and expression profile of genes in buffalo (Bubalus bubalis) oocytes and embryos collected under different environmental stress

The study examined the effects of different environmental stress on developmental competence and the relative abundance (RA) of various gene transcripts in oocytes and embryos of buffalo. Oocytes collected during cold period (CP) and hot period (HP) were matured, fertilized and cultured in vitro to blastocyst hatching stage. The mRNA expression patterns of genes implicated in developmental competence (OCT-4, IGF-2R and GDF-9), heat shock (HSP-70.1), oxidative stress (MnSOD), metabolism (GLUT-1), pro-apoptosis (BAX) and anti-apoptosis (BCL-2) were evaluated in immature and matured oocytes as well as in pre-implantation stage embryos. Oocytes reaching MII stage, cleavage rates, blastocyst yield and hatching rates increased (P < 0.05) during the CP. In MII oocytes and 2-cell embryos, the RA of OCT-4, IGF-2R, GDF-9, MnSOD and GLUT-1 decreased (P < 0.05) during the HP. In 4-cell embryos, the RA of OCT-4, IGF-2R and BCL-2 decreased (P < 0.05) in the HP, whereas GDF-9 increased (P < 0.05). In 8-to 16-cell embryos, the RA of OCT-4 and BCL-2 decreased (P < 0. 05) in the HP, whereas HSP-70.1 and BAX expression increased (P < 0.05). In morula and blastocyst, the RA of OCT-4, IGF-2R and MnSOD decreased (P < 0.05) during the HP, whereas HSP-70.1 was increased (P < 0.05). In conclusion, deleterious seasonal effects induced at the GV-stage carry-over to subsequent embryonic developmental stages and compromise oocyte developmental competence and quality of developed blastocysts.