Identification of differential abundances of mRNA transcript in cumulus cells and CCND1 associated with yak oocyte developmental competence.

The development of an accurate and noninvasive preselection process for competent oocytes is essential to achieve a highly efficient in vitro production (IVP) of embryos. Cumulus cells (CCs) have important functions in oocyte growth, development, maturation, and fertilization. It, therefore, is important to know if the quality of oocytes can be ascertained by assessment of gene expression of the surrounding CCs or not. The aim of this study was to identify differentially expressed genes in yak CCs from oocytes with varying developmental competences as possible biomarkers for distinguishing oocyte competence. The isolated CCs were pooled into immature and mature groups in accordance with the maturation outcome of oocytes. A total of 9516 genes were differentially expressed in the two CC categories (P <  0.05). With a minimum change of 2.5-fold, 45 up-regulated and 79 down-regulated genes were observed in CCs belonging to the mature group compared with those in the immature group (P <  0.01). These genes were primarily enriched for the cell cycle, meiosis, cell signaling, metabolism, and apoptosis. The selected candidate genes (CCND1, BMP15, GDF9, H19, KLF4, GPC1, SYCP3, and CTSB) were validated using quantitative real-time polymerase chain reaction (RT-qPCR) and there were expression patterns similar to those detected with transcriptome analysis. The CCs from fertilized oocytes arrested at the 2-cell (2-cell group), or 8-cell (8-cell group) stages or that developed into blastocysts (the blastocyst group) had a 1.5-, 1.8-, and 2.3-fold increase, respectively, in mRNA relative abundance of CCND1 compared with CCs from unfertilized oocytes (P <  0.05). The results with the RT-qPCR analysis confirmed that the relative abundance of CCND1 mRNA in CCs was associated with oocyte developmental competence. In conclusion, RNA-Seq is useful in extracting transcriptomes and selecting markers associated with oocyte developmental competence. Furthermore, the expression of the CCND1 gene in yak CCs can be used to preselect oocytes for IVP efficiency.

Characterization of transcriptional complexity during pre-implantation development of the yak (Bos grunniens) using RNA-Seq.

The objective of this study was to investigate the mechanism that regulates pre-implantation development of the yak (Bos grunniens). We determined the transcriptomes of in vitro-produced yak embryos at two-cell, four-cell, eight-cell stages, and morula and blastocyst using the Illumina RNA-seq for the first time. We obtained 47.36-50.86 million clean reads for each stage, of which, 85.65%-90.02% reads were covered in the reference genome. A total of 17,368 genes were expressed during the two-cell stage to blastocyst of the yak, of which 7,236 genes were co-expressed at all stages, whereas 10,132 genes were stage-specific expression. Transcripts from 9,827 to 14,893 different genes were detected in various developmental stages. When |log2 ratio| ≥ 1 and q-value <0.05 were set as thresholds for identifying differentially expressed genes (DEGs), we detected a total of 6,922-10,555 DEGs between any two consecutive stages. The GO distributions of these DEGs were classified into three categories: biological processes (23 terms), cellular components (22 terms) and molecular functions (22 terms). Pathway analysis revealed 310 pathways of the DEGs that were operative in early pre-implantation yak development, of which 32 were the significantly enriched pathways. In conclusion, this is the first report to investigate the mechanism that regulates yak embryonic development using high-throughput sequencing, which provides a comprehensive framework of transcriptome landscapes of yak pre-implantation embryos.

Comparative analysis of ovarian transcriptomes between prolific and non-prolific goat breeds via high-throughput sequencing.

To increase the current understanding of the gene expression in the pre-ovulatory ovary and identify the key genes involved in the regulation of ovulation rate, we compared the transcriptomes of ovaries from the prolific Jintang black goat (JTG) and the non-prolific Tibetan goat (TBG) during the follicular phase using the Illumina RNA-Seq method. Three ovarian libraries were constructed for each breed. On average, we obtained approximately 49.2 and 45.9 million reads for each individual ovary of TBGs and JTGs, respectively, of which 79.76% and 78.67% reads were covered in the genome database. A total of 407 differentially expressed genes (DEG) were detected between these two breeds, in which 316 were upregulated, and 91 were downregulated in the ovaries of JTGs versus TBGs. Based on the results of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment, some of these DEGs potentially play an important role in controlling the development of ovarian follicles. SRD5A2, MSMB, STAR and 3BHSD, etc. were the most significantly differentially expressed between these two distinct breeds. In addition, each ovary expressed 1,066 versus 989 novel transcripts, and 171,829 versus 140,529 putative SNPs in TBGs versus JTBs, respectively. All data sets (GEO and dbSNP) were available via public repositories. Our study provides insight into the transcriptional regulation of the ovaries of two distinct breeds of goats that might serve as a key resource for understanding goat fecundity. SRD5A2, MSMB, STAR and 3BHSD may be associated with the high fecundity of JTGs.

Identification of candidate miRNAs and expression profile of yak oocytes before and after in vitro maturation by high-throughput sequencing.

Small RNA represents several unique non-coding RNA classes that have important function in a wide range of biological processes including development of germ cells and early embryonic, cell differentiation, cell proliferation and apoptosis in diverse organisms. However, little is known about their expression profiles and effects in yak oocytes maturation and early development. To investigate the function of small RNAs in the maturation process of yak oocyte and early development, two small RNA libraries of oocytes were constructed from germinal vesicle stage (GV) and maturation in vitro to metaphase II-arrested stage (M II) and then sequenced using small RNA high-throughput sequencing technology. A total of 9,742,592 and 12,168,523 clean reads were obtained from GV and M II oocytes, respectively. In total, 801 and 1,018 known miRNAs were acquired from GV and M II oocytes, and 75 miRNAs were found to be significantly differentially expressed: 47 miRNAs were upregulated and 28 miRNAs were downregulated in the M II oocytes compared to the GV stage. Among the upregulated miRNAs, miR-342 has the largest fold change (9.25-fold). Six highly expressed miRNAs (let-7i, miR-10b, miR-10c, miR-143, miR-146b and miR-148) were validated by real-time quantitative PCR (RT-qPCR) and consistent with the sequencing results. Furthermore, the expression patterns of two miRNAs and their potential targets were analysed in different developmental stages of oocytes and early embryos. This study provides the first miRNA profile in the mature process of yak oocyte. Seventy-five miRNAs are expressed differentially in GV and M II oocytes as well as among different development stages of early embryos, suggesting miRNAs involved in regulating oocyte maturation and early development of yak. These results showed specific miRNAs in yak oocytes had dynamic changes during meiosis. Further functional and mechanistic studies on the miRNAs during meiosis may beneficial to understanding the role of miRNAs on meiotic division.

Oocyte extract improves epigenetic reprogramming of yak fibroblast cells and cloned embryo development.

The objective was to investigate the effects of bovine oocyte extract (BOE) on epigenetic reprogramming of yak fibroblast cells, based on their cell cycle status, histone acetylation, DNA methylation, gene expression, and cloned blastocyst formation. Permeabilization of yak fibroblasts after treatment with 10 or 50 µL of BOE (treated-S and treated-L groups, respectively) for 24 hours increased (P < 0.05) the cell population at the G(0)/G(1) phase (85.2 ± 2.3% and 89.6 ± 1.5%, respectively) compared with controls (75.4 ± 1.1%). Acetylation at lysine 9 of histone H3 was also higher (26.1 ± 1.4 and 33.5 ± 2.1) than in the control group (15.3 ± 1.6; P < 0.05). Moreover, BOE reduced methylation of the promoter regions of Oct-4 and Nanog (76.4% and 72.2%; and 35.6% and 30.0%, respectively) compared with the control group (92.1% and 47.8%; P < 0.05). In addition, the relative expression levels of HDAC-1, HADC-2, Dnmt-1, and Dnmt-3a were downregulated (P < 0.05) after yak fibroblasts were treated with BOE. Furthermore, when yak fibroblasts were used for interspecies somatic cell nuclear transfer after BOE treatment, 8-cell and blastocyst formation rates significantly exceeded those of the control. In conclusion, BOE induced epigenetic reprogramming of yak fibroblasts, making them suitable donors for yak interspecies somatic cell nuclear transfer.

Lowering storage temperature during ovary transport is beneficial to the developmental competence of bovine oocytes used for somatic cell nuclear transfer.

The objective of this study was to determine the effect of storage temperature during ovary transport on the developmental competence of bovine oocytes for use in somatic cell nuclear transfer (SCNT). Ovaries obtained from a slaughterhouse were stored in physiological saline for 3-4h at one of the three temperatures: 15 °C, 25 °C, or 35 °C. The developmental competence of oocytes used for SCNT was ascertained by cleavage and blastocyst formation rate, total cell number, apoptosis index, and the relative abundance of Bax and Hsp70.1 in day 7 blastocysts. Ovaries stored at 35 °C for 3-4h reduced the recovery rate of grade I and II oocytes compared with those stored at 25 °C or 15 °C (45.1±0.7% vs. 76.7±1.2% or 74.8±2.0%, P<0.05). The proportion of oocytes matured to the MII stage (maturation rate) for oocytes stored at 35 °C was significantly lower than those stored at 25 °C or 15 °C (51.3±0.9% vs. 75.1±1.4% or 71.7±1.3%, P<0.05). Cleavage rate (77.7±2.1%, 77.9±1.1% and 72.1±0.7% for 15 °C, 25 °C and 35 °C groups, respectively) and blastocyst formation rate (39.1±0.5%, 36.8±1.4% and 32.2±0.9% for 15 °C, 25 °C and 35 °C groups, respectively) following SCNT were not significantly different between treatments. Oocytes from ovaries stored at 15 °C, however, produced blastocysts with higher cell numbers (97.3±8.6 vs. 80.2±10.8 or 77.4±11.7; P<0.05) and lower apoptotic index (5.1±1.3 vs. 13.5±1.6 or 18.6±1.1, P<0.05) than those stored at 25 °C or 35 °C. The relative abundance of Bax and Hsp70.1 in day 7 blastocysts produced from oocytes derived from ovaries stored at 15 °C was lower than those stored at 25 °C or 35 °C (P<0.05). It was concluded that a storage temperature of 15 °C for a 3-4h period had a significant beneficial effect on the quality and developmental competence of oocytes used for SCNT due to the alleviation of stresses on the oocytes compared with those subjected to storage temperatures of 25 °C or 35 °C.

Production of cloned calves by combination treatment of both donor cells and early cloned embryos with 5-aza-2/-deoxycytidine and trichostatin A.

We previously reported that treatment of both donor cells and early cloned embryos with a combination of 0.01 µM 5-aza-2(/)-Deoxycytidine (5-aza-dC) and 0.05 µM trichostatin A (TSA) significantly improved development of cloned bovine embryos in vitro. In the present study, we investigated the effect of this combination treatment on the in vivo development potency and postnatal survivability of cloned calves. Blastocysts (77 and 82 blastocysts derived from untreated (control) and treated groups, respectively) were individually transferred to recipient cows. Relative to the control group, the combination treatment of both donor cells and early embryos with 5-aza-dC and TSA dramatically increased the cleavage rate (49.2 vs 63.6%, P < 0.05) at 24 h of culture, and blastocyst development rate on Days 6 and 7 of culture (18.8 vs 33.9% and 27.1 vs 38.5% respectively, P < 0.05). Although pregnancy rate did not differ 40 d after transfer, it was lower in the treated than control group 90 d after transfer (7.8 vs 29.3%, P < 0.05). In the control group, there were three calves born to 77 recipients (only two survived beyond 60 d), whereas in the treated group, 17 calves were born to 82 recipients, and 11 survived beyond 60 d. In conclusion, a combination treatment of donor cells and early cloned embryos with 5-aza-dC and TSA significantly enhanced development of somatic cell cloned bovine embryos in vivo; cloning efficiency (number of surviving calves at 60 d of birth/number of recipient cows) was increased from 2.6 to 13.4%.