Transcriptome profiling of rabbit parthenogenetic blastocysts developed under in vivo conditions

Analyses of histological and transcriptome differences in the skin of short-hair and long-hair rabbits

BACKGROUND

Hair fibre length is an important economic trait of rabbits in fur production. However, molecular mechanisms regulating rabbit hair growth have remained elusive.

RESULTS

Here we aimed to characterise the skin traits and gene expression profiles of short-hair and long-hair rabbits by histological and transcriptome analyses. Haematoxylin-eosin staining was performed to observe the histological structure of the skin of short-hair and long-hair rabbits. Compared to that in short-hair rabbits, a significantly longer anagen phase was observed in long-hair rabbits. In addition, by RNA sequencing, we identified 951 genes that were expressed at significantly different levels in the skin of short-hair and long-hair rabbits. Nine significantly differentially expressed genes were validated by quantitative real-time polymerase chain reaction. A gene ontology analysis revealed that epidermis development, hair follicle development, and lipid metabolic process were significantly enriched. Further, we identified potential functional genes regulating follicle development, lipid metabolic, and apoptosis as well as important pathways including extracellular matrix-receptor interaction and basal cell carcinoma pathway.

CONCLUSIONS

The present study provides transcriptome evidence for the differences in hair growth between short-hair and long-hair rabbits and reveals that lipid metabolism and apoptosis might constitute major factors contributing to hair length.

Disruption of NNAT, NAP1L5 and MKRN3 DNA methylation and transcription in rabbit parthenogenetic fetuses

Parthenogenetically activated oocytes cannot develop to term in mammals due to lack of paternal gene expression. Disruption of imprinted gene expression and DNA methylation status in parthenogenetic fetuses has been reported in mice and pigs, but not in rabbits. In this study, the genomic imprinting status of the paternally expressed genes Neuronatin (NNAT), Nucleosome assembly protein 1-like 5 (NAP1L5), and Makorin ring finger protein 3 (MKRN3) was compared between rabbit parthenogenetic (PA) and normally fertilized fetuses (Con) using quantitative real-time PCR (qRT-PCR) and bisulfite sequencing PCR (BSP). The results revealed a significantly reduced expression of NNAT, NAP1L5, and MKRN3 in rabbit PA fetuses compared with Con fetuses (p<0.05). In addition, the BSP results demonstrated hypermethylation in the differentially methylated regions (DMRs) of NNAT, NAP1L5, and MKRN3 in rabbit PA fetuses. Taken together, these results suggest that hypermethylation of DMRs is associated with decreased NNAT, NAP1L5, and MKRN3 expression, which may be responsible for developmental failure of rabbit PA fetuses.

Transcriptome Analysis of Thermal Parthenogenesis of the Domesticated Silkworm

Thermal induction of parthenogenesis (also known as thermal parthenogenesis) in silkworms is an important technique that has been used in artificial insemination, expansion of hybridization, transgenesis and sericultural production; however, the exact mechanisms of this induction remain unclear. This study aimed to investigate the gene expression profile in silkworms undergoing thermal parthenogenesis using RNA-seq analysis. The transcriptome profiles indicated that in non-induced and induced eggs, the numbers of differentially expressed genes (DEGs) for the parthenogenetic line (PL) and amphigenetic line (AL) were 538 and 545, respectively, as determined by fold-change ≥ 2. Gene ontology (GO) analysis showed that DEGs between two lines were mainly involved in reproduction, formation of chorion, female gamete generation and cell development pathways. Upregulation of many chorion genes in AL suggests that the maturation rate of AL eggs was slower than PL eggs. Some DEGs related to reactive oxygen species removal, DNA repair and heat shock response were differentially expressed between the two lines, such as MPV-17, REV1 and HSP68. These results supported the view that a large fraction of genes are differentially expressed between PL and AL, which offers a new approach to identifying the molecular mechanism of silkworm thermal parthenogenesis.

Effect of in vitro and in vivo conditions on development of parthenogenetic rabbit embryos after vitrification

Parthenote embryos offer multiple opportunities in biotechnological research, so it is important to analyse the possibilities for their cryopreservation in order to establish a biobank. The aim of this experiment was to determine the effect of culture conditions and vitrification on rabbit parthenogenetic embryos. Parthenotes were cultured under in vivo and in vitro conditions until day 3 (late morula/early blastocyst), when they were vitrified. Immediately after warming, they were newly cultured under in vivo and in vitro conditions till day 6 (blastocyst stage). Both culture conditions showed similar late morula/early blastocyst (0.39±0.056 vs. 0.46±0.043, for in vivo and in vitro, respectively) and blastocyst rates (0.12±0.068 vs. 0.13±0.070, for in vivo and in vitro, respectively). However, no parthenote was recovered when a combination of culture conditions was performed. To our best knowledge, this is the first demonstration of the ability of rabbit parthenogenetic embryos to develop after vitrification, with similar embryo development after in vivo or in vitro culture. Nevertheless, our results highlight the importance of culture conditions on the morphology of parthenote embryos. Therefore, we have described that special attention should be paid on culture conditions to generate parthenote embryos, with a view to their subsequent use, for example in embryonic stem cell production.

Characterization of the altered gene expression profile in early porcine embryos generated from parthenogenesis and somatic cell chromatin transfer

The in vitro production of early porcine embryos is of particular scientific and economic interest. In general, embryos produced from in vitro Assisted Reproductive Technologies (ART) manipulations, such as somatic cell chromatin transfer (CT) and parthenogenetic activation (PA), are less developmentally competent than in vivo–derived embryos. The mechanisms underlying the deficiencies of embryos generated from PA and CT have not been completely understood. To characterize the altered genes and gene networks in embryos generated from CT and PA, comparative transcriptomic analyses of in vivo (IVV) expanded blastocysts (XB), IVV hatched blastocyst (HB), PA XB, PA HB, and CT HB were performed using a custom microarray platform enriched for genes expressed during early embryonic development. Differential expressions of 1492 and 103 genes were identified in PA and CT HB, respectively, in comparison with IVV HB. The “eIF2 signalling”, “mitochondrial dysfunction”, “regulation of eIF4 and p70S6K signalling”, “protein ubiquitination”, and “mTOR signalling” pathways were down-regulated in PA HB. Dysregulation of notch signalling–associated genes were observed in both PA and CT HB. TP53 was predicted to be activated in both PA and CT HB, as 136 and 23 regulation targets of TP53 showed significant differential expression in PA and CT HB, respectively, in comparison with IVV HB. In addition, dysregulations of several critical pluripotency, trophoblast development, and implantation-associated genes (NANOG, GATA2, KRT8, LGMN, and DPP4) were observed in PA HB during the blastocyst hatching process. The critical genes that were observed to be dysregulated in CT and PA embryos could be indicative of underlying developmental deficiencies of embryos produced from these technologies.

Microarray analysis of gene expression in parthenotes and in vitro-derived goat embryos

The present work was carried out to investigate the global gene expression profile to search differentially expressed candidate transcripts between parthenogenetic and in vitro-fertilized (IVF) caprine morula. For this study, total RNA was isolated from diploid parthenogenetic and IVF embryos, and complementary DNA was synthesized. Microarray and relative real-time polymerase chain reaction analysis were performed to check global gene expression profile and validation, respectively. According to the microarray analysis, the total number of upregulated (UR) and downregulated (DR) genes was 613 and 220, respectively in diploid parthenogenetic morula as compared with IVF morula. The number of genes showing about two-, two- to five-, five- to 10-, 10- to 20-, and above 20-fold UR and DR genes was 147, 229, 122, 59, and 56 and 94, 73, 18, 13, and 22, respectively. Five UR genes validated (PTEN, PHF3, CTNNB1, SELK, and NPDC1) and all of them were significantly higher in parthenotes, which was in accordance with microarray results, whereas the expression of DR (AURKC and KLF15) genes were downregulated in parthenotes as observed in microarray results but the difference was not significant (P < 0.05). In conclusion, our findings demonstrate differential expression of a large number of genes in parthenotes compared with IVF embryos, which may be the reason for aberrant parthenogenetic embryo development in caprine species.