Identification of differentially expressed placental transcripts during multiple gestations in the Eurasian beaver (Castor fiber L.).

The Eurasian beaver is one of the largest rodents that, despite its high impact on the environment, is a non-model species that lacks a reference genome. Characterising genes critical for pregnancy outcome can serve as a basis for identifying mechanisms underlying effective reproduction, which is required for the success of endangered species conservation programs. In the present study, high-throughput RNA sequencing (RNA-seq) was used to analyse global changes in the Castor fiber subplacenta transcriptome during multiple pregnancy. De novo reconstruction of the C. fiber subplacenta transcriptome was used to identify genes that were differentially expressed in placentas (n=5) from two females (in advanced twin and triple pregnancy). Analyses of the expression values revealed 124 contigs with significantly different expression; of these, 55 genes were identified using MegaBLAST. Within this group of differentially expressed genes (DEGs), 18 were upregulated and 37 were downregulated in twins. Most DEGs were associated with the following gene ontology terms: cellular process, single organism process, response to stimulus, metabolic process and biological regulation. Some genes were also assigned to the developmental process, the reproductive process or reproduction. Among this group, four genes (namely keratin 19 (Krt19) and wingless-type MMTV integration site family - member 2 (Wnt2), which were downregulated in twins, and Nik-related kinase (Nrk) and gap junction protein β2 (Gjb2), which were upregulated in twins) were assigned to placental development and nine (Krt19, Wnt2 and integrin α7 (Itga7), downregulated in twins, and Nrk, gap junction protein β6 (Gjb6), GATA binding protein 6 (Gata6), apolipoprotein A-I (ApoA1), apolipoprotein B (ApoB) and haemoglobin subunit α1 (HbA1), upregulated in twins) were assigned to embryo development. The results of the present study indicate that the number of fetuses affects the expression profile in the C. fiber subplacental transcriptome. Enhancement of transcriptomic resources for C. fiber will improve understanding of the pathways relevant to proper placental development and successful reproduction.