Search

During mammalian pre-implantation embryonic development, dramatic and orchestrated changes occur in gene transcription. Pregnancy rates were low when yak females were crossbred with cattle breeds, but few studies exist to describe the unique molecular network regulation behind the pre-implantation development of these embryos. We determined the transcriptomes of crossbred embryos derived from yak oocytes in vitro fertilized with Jersey sperm using Illumina RNA-seq for the first time in this study. Embryos were sampled at the 2-, 4-, and 8-cell, morula and blastocyst stages. The results showed that in total, 291.9 million short reads were generated from the five libraries of 2-, 4-, and 8-cell, morula and blastocyst stages, with 276.2 million high-quality reads selected for further analysis. Eighty to 91% of the clean reads were aligned against the yak reference genome. A total of 19,072 transcripts were identified in five libraries, of which 7,785 transcripts were co-expressed in each stage and 2,013 transcripts were stage-specific. When a |log 2 ratio| ≥1 and q-value ≤ 0.05 were set as thresholds for identifying differentially expressed genes (DEGs), we detected a total of 3,690 to 10,298 DEGs between any two consecutive stages. Based on the results of GO and KEGG enrichment, some of these DEGs potentially play an important role in regulating pre-implantation development, but they are most likely stage-specific. There were 2,960, 7,287, 6,420, 7,724 and 10,417 DEGs in 2-, 4-, 8-cell, morula and blastocyst stages between the crossbred embryos and purebred embryos of the yak, respectively, leading to a large difference in GO terms and pathways. In conclusion, we sequenced transcriptomes of in vitro-produced crossbred embryos of yak and cattle during pre-implantation and provided comprehensive examinations of gene activities. These will be helpful for development of assisted reproductive technology and better understanding the early maternal-fetal or maternal-embryonic dialog in inter-species crossbreeding.

The kidding rate is one of the most important economic traits for goat production, but the genetic mechanism that is associated with ovulation rate is poorly understood. Recently, increasing evidence has suggested that microRNAs (miRNAs) influence ovarian biological processes. The present study provides the first comparison of the ovarian miRNAs of prolific Jintang black goats (JTGs) and non-prolific Tibetan goats (TBGs) during the follicular phase using RNA-Seq technology. We generated 11.19 million (M) and 11.34 M clean reads from the TBG and JTG libraries, respectively, from which a total of 389 known miRNAs were identified and 142 novel miRNAs were predicted. A total of 191 miRNAs were differentially expressed between the two breeds. Among the 10 most abundant miRNAs, miR-21-5p was defined as differentially expressed miRNA with a higher level in the JTG library than in the TBG library, but the other miRNAs were not different between the breeds. The predicted miRNA-targeted genes were further analyzed by Gene Ontology and KEGG pathway analyses. The results revealed that miR-21, miR-99a, miRNA-143, let-7f, miR-493 and miR-200b may affect follicular development. These findings will increase the current understanding of the role of ovarian miRNAs in the regulation of ovulation rate in goats.