RNA sequencing and transcriptome analysis of buffalo (Bubalus bubalis) blastocysts produced by somatic cell nuclear transfer and in vitro fertilization.

Across farm animal species, the live birth rate obtained with somatic cell nuclear transfer (SCNT) embryos is only <2% compared with >40% obtained with in vitro fertilization (IVF) embryos, primarily due to incomplete nuclear reprogramming which results in aberrant embryonic gene expression. We used RNA sequencing to compare the global transcriptome profile of SCNT and IVF buffalo blastocysts. SCNT blastocysts expressed 17,061 transcripts, of which 941 were unique whereas, IVF blastocysts expressed 17,303 transcripts, of which 1,183 were unique. At ≥2-folds change (p < .05), 331 transcripts were differentially expressed in the two groups among which, 19 were unique, 188 were downregulated and 143 were upregulated in SCNT compared with IVF blastocysts. Many genes affecting pluripotency, trophectoderm development, developmental regulation, and epigenetic modifications were upregulated in SCNT compared with IVF blastocysts. Among the four functional categories analyzed, epigenetic regulators were the most affected. Most of the WNT signaling pathway genes were upregulated whereas, the inhibitors of this pathway, such as DKK1, were downregulated in SCNT blastocysts, suggesting that this pathway is overexpressed in SCNT embryos. Gene Ontology analysis revealed that 25 biological processes, 20 molecular functions, and 24 cellular compartment categories were enriched in SCNT blastocysts. This data can help identify reprogramming errors for improving cloning efficiency.

Selection of Reference miRNAs for Relative Quantification in Buffalo (Bubalus bubalis) Blastocysts Produced by Hand-Made Cloning and In Vitro Fertilization.

Very low birth rate and a high incidence of abnormalities in offspring born from cloned embryos, which have limited the application of cloning technology on a wide scale, are believed to be because of incomplete or aberrant nuclear reprogramming. MicroRNAs (miRNAs) are involved in regulating a wide range of biological processes including reprogramming and embryonic development. Selection of suitable reference miRNAs is critical for normalization of data for accurate relative quantification of miRNAs by quantitative real-time polymerase chain reaction (qRT-PCR), which is currently the most widely used technique for quantifying miRNAs. This study was aimed at identification of reference miRNAs suitable for normalization of qRT-PCR data from blastocyst-stage buffalo embryos produced by handmade cloning and in vitro fertilization (IVF). RNA isolated from cloned and IVF blastocysts was subjected to next-generation sequencing based on which, 12 highly and most consistently expressed miRNAs, which included miR-92a, miR-423, miR-151, Let-7a, miR-103a, miR-93, miR-16b, miR-25, miR-30e, miR-101, miR-127, and miR-197, were selected as candidates for identification of suitable reference miRNAs using three statistical algorithms namely geNorm, NormFinder, and BestKeeper. Based on consensus of the three algorithms, the combination of miRNAs found to be suitable as reference miRNAs were miR-127 and miR-103 for IVF blastocysts; miR-92a and miR-103 for cloned blastocysts, and miR-103, miR-423, and miR-93 across both IVF and cloned blastocysts. The data of this study can be very useful in miRNA expression analysis of blastocyst-stage cloned and IVF embryos.

Buffalo SCNT embryos exhibit abnormal gene expression of ERK/MAPK pathway and DNA methylation.

Inhibition of ERK/MAPK pathway has been shown to decrease DNA methylation via down-regulation of DNA methyltransferases (DNMTs) in several studies suggesting that this pathway plays an important role in regulation of DNA methylation. We examined the relative expression level of seven important genes related to ERK/MAPK pathway and DNMTs (DNMT1, DNMT3a and DNMT3b) by quantitative real-time PCR in buffalo blastocysts produced by Hand-made cloning and compared it with that in blastocyst-stage embryos produced by in vitro fertilization (IVF). The expression level of six of seven genes related to ERK/MAPK pathway examined i.e., p21RAS, RAF1, AKT1, ERK2, PIK3R2 and c-Myc was significantly higher (p < 0.05) in cloned than in IVF embryos. However, the expression level of FOS was lower (p < 0.005) in cloned than in IVF embryos. The relative expression level of DNMT3a and DNMT3b but not that of DNMT1 was significantly higher (p < 0.05) in cloned than in IVF embryos. These results indicate that the cloned embryos exhibit an abnormal expression of several important genes related to ERK/MAPK pathway and DNMTs. Although a direct link between ERK/MAPK pathway and DNMTs was not examined in the present study, it can be speculated that ERK/MAPK pathway may have a role in regulating the expression of DNMTs in embryos, as also observed in other tissues.

Selection of Suitable Internal Control Genes for Accurate Normalization of Real-Time Quantitative PCR Data of Buffalo (Bubalus bubalis) Blastocysts Produced by SCNT and IVF.

We evaluated the suitability of 10 candidate internal control genes (ICGs), belonging to different functional classes, namely ACTB, EEF1A1, GAPDH, HPRT1, HMBS, RPS15, RPS18, RPS23, SDHA, and UBC for normalizing the real-time quantitative polymerase chain reaction (qPCR) data of blastocyst-stage buffalo embryos produced by hand-made cloning and in vitro fertilization (IVF). Total RNA was isolated from three pools, each of cloned and IVF blastocysts (n = 50/pool) for cDNA synthesis. Two different statistical algorithms geNorm and NormFinder were used for evaluating the stability of these genes. Based on gene stability measure (M value) and pairwise variation (V value), calculated by geNorm analysis, the most stable ICGs were RPS15, HPRT1, and ACTB for cloned blastocysts, HMBS, UBC, and HPRT1 for IVF blastocysts and RPS15, GAPDH, and HPRT1 for both the embryo types analyzed together. RPS18 was the least stable gene for both cloned and IVF blastocysts. Following NormFinder analysis, the order of stability was RPS15 = HPRT1>GAPDH for cloned blastocysts, HMBS = UBC>RPS23 for IVF blastocysts, and HPRT1>GAPDH>RPS15 for cloned and IVF blastocysts together. These results suggest that despite overlapping of the three most stable ICGs between cloned and IVF blastocysts, the panel of ICGs selected for normalization of qPCR data of cloned and IVF blastocyst-stage embryos should be different.