Skip to main content

Comparison of American mink embryonic stem and induced pluripotent stem cell transcriptomes



Recently fibroblasts of many mammalian species have been reprogrammed to pluripotent state using overexpression of several transcription factors. This technology allows production of induced pluripotent stem (iPS) cells with properties similar to embryonic stem (ES) cells. The completeness of reprogramming process is well studied in such species as mouse and human but there is not enough data on other species. We produced American mink (Neovison vison) ES and iPS cells and compared these cells using transcriptome analysis.


We report the generation of 10 mink ES and 22 iPS cell lines. The majority of the analyzed cell lines had normal diploid chromosome number. The only ES cell line with XX chromosome set had both X-chromosomes in active state that is characteristic of pluripotent cells. The pluripotency of ES and iPS cell lines was confirmed by formation of teratomas with cell types representing all three germ layers. Transcriptome analysis of mink embryonic fibroblasts (EF), two ES and two iPS cell lines allowed us to identify 11831 assembled contigs which were annotated. These led to a number of 6891 unique genes. Of these 3201 were differentially expressed between mink EF and ES cells. We analyzed expression levels of these genes in iPS cell lines. This allowed us to show that 80% of genes were correctly reprogrammed in iPS cells, whereas approximately 6% had an intermediate expression pattern, about 7% were not reprogrammed and about 5% had a "novel" expression pattern. We observed expression of pluripotency marker genes such as Oct4, Sox2 and Rex1 in ES and iPS cell lines with notable exception of Nanog.


We had produced and characterized American mink ES and iPS cells. These cells were pluripotent by a number of criteria and iPS cells exhibited effective reprogramming. Interestingly, we had showed lack of Nanog expression and consider it as a species-specific feature.


Recently mouse and human adult and embryonic fibroblasts (EF) have been reprogrammed into pluripotent state by overexpression of only four transcription factors [13]. At present, induced pluripotent stem (iPS) cells have been derived from somatic cells of mammalian species such as primates [47], rat [8], prairie vole [9], rabbit [10], dog [1116], pig [17], horse [18], sheep [19], cow [20], goat [21] and buffalo [22]. In addition, iPS cells have been produced from endangered species, e.g. rhinoceros [7] and snow leopard [23]. Mouse and human iPS cells are very similar to embryonic stem (ES) cells produced from blastocyst inner cell mass [24]. Apart from extensively studied mouse and human iPS cells, little is known about characteristics of iPS cells of other mammalian species. Additionally, ES cells are produced from a limited number of species and their pluripotent state is rarely well characterized. This complicates comparison between ES and iPS cells. So far, there are no reports of iPS cell derivation from somatic cells of any mustelid species. The aim of the current research is producing mink ES and iPS cells to assess reprogramming completeness by comparing gene expression profiles of mink EF, ES and iPS cells.

We report generation of American mink pluripotent ES and iPS cells. Transcriptome analysis shows efficient EF genome reprogramming. Pluripotent stem cells express key pluripotency markers with notable exception of Nanog.


ES and iPS cell derivation

American mink ES cell line MES12 used in this study was previously described by our group [25]. By using the same protocol we have generated an additional set of mink ES cell lines. Twelve mink morulas and early blastocysts were plated on feeder cells, and in a total of 10 ES cell lines were produced. As in the previous series, they were designated as MES (MES20 - MES29). Unlike mouse ES cells, mink ES cells contain peripheral visible granules and form flat monolayer colonies of epithelial-like morphology [25, 26].

To produce iPS cells from primary mink EF we used the following human reprogramming transcription factors: OCT4, SOX2, C-MYC and KLF4. To facilitate reprogramming process we supplemented culture medium with valproic acid, an epigenetic modifier that inhibits histone deacetylases. Two weeks after first lentiviral transduction of EF we observed many colonies with different morphology. Based on mink ES-like morphology we selected and picked up 25 colonies. From these primary colonies 22 iPS cell lines were successfully produced. We designated them as iNV (iNV1 - iNV22). Morphology of iPS cell line iNV11 colony is shown in Figure 1a.

Figure 1
figure 1

iPS cell morphology, karyotype and examples of cell types in teratomas developed from mink ES and iPS cells. a - morphology of iNV11 iPS cell line colony; b - iNV11 metaphase plate; c - iNV11 karyotype; d - expression of human transgene C-MYC; e - expression of human transgene KLF4; examples of representative cell types in teratomas formed from pluripotent cell lines: f - epidermal epithelium with hair follicle; g - gut-like epithelium; h - bone, adipose and hematopoietic tissues; i - rosettes of neural epithelium; j - gut-like epithelium; k - connective and adipose tissues; l - rosettes of neural epithelium; m - hepatocytes; n - erythropoietic tissue.

We assessed the efficacy of colony formation in a separate experiment with the same conditions as above. We observed 722 colonies (1.2% of initial 6 × 104 transduced cells) on day 11 after transduction of mink EF. Most of them had mink ES-like phenotype.

Cytogenetic analysis

We randomly selected 6 out of 10 ES cell lines and 11 out of 22 iPS cell lines for the analysis. The results of the cytogenetic analysis are shown in Table 1. Five ES cell lines were XY, and MES25 was XX. All analyzed iPS cell lines were XY as they were produced from EF derived from a single male embryo. All analyzed pluripotent cell lines with the exception of MES20 had diploid modal chromosome number. We consider it as an indication of karyotype stability. In MES20 64.3% of cells had tetraploid chromosome number. Tetraploidization of pluripotent cells frequently happens in vitro for human and mouse ES cells and thus some percentage of tetraploid cell lines could be considered as normal [27]. A metaphase plate and a karyotype of iNV11 are presented in Figure 1b-c.

Table 1 Cytogenetic analysis of ES and iPS cell lines.

X-inactivation analysis

Female pluripotent cells have both X-chromosomes in active state. One of the markers of inactive X-linked chromatin modifications is histone H3 trimethylated at lysine 27 (H3K27me3) [28, 29]. Control female mink EF had 82.3% (n = 96) cells positive for H3K27me3 as expected. In the only XX pluripotent cell line, MES25, 87% (n = 285) of the cells were negative for inactive X-linked chromatin marker. Most of the positive cells were situated at the borders of ES cell colonies and might represent differentiation (Additional file 1). Thus, MES25 has both X-chromosomes in active state confirming its pluripotency.

Generation and analysis of teratomas

To analyze the pluripotency of ES and iPS cells we used the in vivo teratoma formation assay. The summary of teratoma analysis data is presented in Table 2. Cell types in the individual teratomas are presented in Additional file 2. Teratomas generated from ES cell lines MES22, MES24 and MES29, and iPS cell lines iNV7, iNV11, iNV19 and iNV20 had most or all analyzed cell types representing three germ layers (ectoderm, endoderm and mesoderm). Thus, we consider these cell lines pluripotent. Representative cell types in MES29, iNV7 and iNV11 derived teratomas are shown in Figure 1f-n. These three cell lines and the previously characterized MES12 [25] were chosen for the transcriptome analysis.

Table 2 Summary of histological analysis of teratomas formed after injection of ES and iPS cells into immunodeficient mice.

Transgene silencing analysis

Transgenes inserted into pluripotent stem cells with viral vectors are usually silenced [28, 30, 31]. We performed analysis of the human transgene silencing in iPS cell lines iNV7 and iNV11 by RT-PCR using human specific primers (Additional file 3) for OCT4, C-MYC, KLF4 and SOX2. Human transgenes KLF4 and C-MYC were silenced, as no product was visible in the agarose gel with the exception of a slight band in the iNV7 line for C-MYC (Figure 1d-e). In addition, we showed C-MYC and KLF4 transgene silencing for iPS cell lines iNV5, iNV15 and iNV20 (data not shown). In the case of SOX2 and OCT4 genes, due to their high sequence conservation between human and mink in the coding region, we were not able to design human specific primers. Therefore, we have no information about the silencing of the SOX2 and OCT4 transgenes. We can thus conclude that the human transgenes are mostly silenced in the mink iPS cell lines.

Quantitative real-time PCR (qPCR) results

We verified expression of several mink genes by qPCR (Additional file 4 Additional file 5). We found significant differences between the five cell lines for the gene expression: Oct4 (P = 2.00E-08), Nanog (P = 3.00E-08), Gdf3 (P = 8.00E-08), Sox2 (P = 1.63E-05) and Nestin (P = 0.0351). Interestingly, looking at the pairwise post comparisons (Additional file 6 Figure 2a), we can see that gene expression varies between pluripotent cell lines. As expected, the mink EF show a pattern of expression very different from the pluripotent cell lines.

Figure 2
figure 2

Expression levels of several genes of interest in mink EF, ES and iPS cells based on transcriptome analysis and qPCR. a - results of qPCR for Oct4, Sox2, Gdf3, Nestin (Nes) and Nanog; b- expression levels of Oct4, Sox2, Rex1, Nanog, Cer1, Otx1, Gata6 and Gapdh genes in mink pluripotent cells and EF. The vertical axis represents counts determined for each sample by transcriptome analysis.

Transcriptome analysis

To assess completeness of reprogramming in mink iPS cells we compared transcriptomes of ES cells, iPS cells and EF. We selected cell lines for transcriptome analysis based on cytogenetic analysis and teratoma formation test. For transcriptome analysis all cell lines were cultured for two passages in three independent replicates. Mink EF were collected on passage 7, ES cells lines MES12 and MES29 on passages 21 and 25, respectively, iPS cell lines iNV7 and iNV11 on passages 10 and 11, respectively. To eliminate contamination with mink EF transcripts ES and iPS cells were grown on mouse strain CD-1 feeder cells. We produced between 4.0 and 23.9 million reads for each cell line replicate (see Additional file 7).

Hierarchical clustering, as well as principal component analysis (PCA) of transcriptome data, shows a sharp contrast between EF and pluripotent cells (Figure 3a-b). 80% of all differences in expression levels observed between cell lines can be explained by difference between differentiated and pluripotent cells, whereas only 13% of dissimilarities reflect difference between ES and iPS cells (Figure 3b). This result clearly indicates that the reprogramming process results in cells that lose most of EF signature genes and are very similar to the pluripotent ES cells.

Figure 3
figure 3

Transcriptome gene expression analysis. a - hierarchical cluster dendrogram based on the expression levels of 100 genes with a highest variance between samples; b - principal component analysis of expression based on the same set of genes as in a. The horizontal axis represents the first principal component of variation - a projection of the array data representing the maximal variance among them (80% of total variance between samples). The vertical axis represents the second principal component of variation (13% of total variance); c - example of heat-map presentation of gene expression profiles of correctly reprogrammed genes in mink EF, ES and iPS cells. Each row of the heat-map corresponds to one gene, and expression levels are shown using color-code (blue color represents the lowest expression level, red color - the highest); d - example of heat-map presentation of gene expression profiles of not reprogrammed genes in mink EF, ES and iPS cells.

The observed slight difference between ES and iPS cells might reflect incomplete reprogramming, a phenomenon well known for mouse and human iPS cells. A common source of incomplete reprogramming is so-called "somatic memory" determined as expression of certain genes in iPS cells on a level of original somatic cells [32]. To study "somatic memory" in iPS cells we decided to perform detailed analysis of genes that are differently expressed in EF and ES cells. We annotated 11831 contigs out of 30490, generated by Trinity. Among them we identified 6891 unique genes and found 3201 showing significant difference in expression between EF and either MES12 or MES29 cells (Additional file 8). Of these, the majority of genes were expressed at the same level in iPS and ES cells, thus being correctly reprogrammed. List and examples of such genes are presented in Additional file 9, Figure 2b (Oct4, Sox2 and Rex1) and Figure 3c. Some genes show intermediate expression pattern, i.e. characterized by expression level between EF and ES cells. A small number of genes were not reprogrammed, displaying EF-like expression pattern in iPS cells with expression level significantly different from both MES12 and MES29 cell lines (Figure 3)d. Finally, few genes do not fall into correctly reprogrammed, intermediate or not reprogrammed categories. We designed these genes as "novel", since they show a novel expression pattern in iPS cells, different from both EF and ES cells (Additional file 10). In contrast to genes from intermediate category, the expression level of these "novel" genes does not lie between levels of EF and ES cells, suggesting that unexpected expression of these genes is not due to incomplete activation or repression of respective promoters during reprogramming process, but rather is a result of some reprogramming mistakes. The lists and number of genes with different patterns of reprogramming are presented in Additional file 8 and Table 3.

Table 3 Percentage of different categories of genes in mink iPS cells.


We have produced American mink ES and iPS cell lines. Their pluripotency was confirmed by teratoma formation test. We have observed silencing of two of four transgenes in analyzed iPS cell lines. That is in line with reports on transgene silencing in mouse and human iPS cells [28, 30, 31]. American mink represents one of a few species for which both ES and iPS cells were obtained. This had allowed us to address several important questions. One of them is completeness of somatic cell reprogramming, whether genes of EF origin are effectively reprogrammed to pluripotent state in iPS cells. Another question refers to the nature of the pluripotent state in the American mink ES and iPS cells, naïve or primed [33, 34].

We had shown that the majority of genes were correctly reprogrammed (Table 3). About 400 genes were not reprogrammed or have "novel" expression pattern in iPS cells. That is comparable with data on human iPS cells where approximately 600 genes were not reprogrammed in iPS cells compared to ES cells [32]. In another study, it was shown that in 12 human iPS cell lines an average of 550 genes were not reprogrammed compared to ES cell reference. It should be noted that ES cell lines differed in gene expression. In the same study an analysis of 20 ES cell lines revealed on average 387 differentially expressed genes in each line compared to ES cell reference [24]. Thus, such number of differentially expressed genes could be attributed to "normal" gene expression fluctuations or differences. In general, we conclude that the majority of the genes in both analyzed iPS cell lines were effectively reprogrammed.

The second question concerns the status of pluripotency in mink cells. Recently, pluripotent cells were divided into two categories: naïve and primed, distinguished by colony morphology, pattern of gene expression and other features [33, 34]. Mouse ES cells are considered naïve while human ES cells resemble primed epiblast stem cells. The appliance of 2i inhibitors (MEK inhibitor PD 0325901 and GSK3 inhibitor CHIR 99021) was shown to switch the pluripotent state from primed to naïve. Among species with reported iPS cell derivation the closest to N. vison is dog, Canis familiaris, order Carnivora. Interestingly, in canine ES [3539] and iPS cells [1116] different research groups were able to observe both morphological cell types, although their pluripotent state was not identified. However, it should be noted that only three research groups were able to produce teratomas with all three germ layer derivatives from canine ES cells [38] and iPS cells [15, 16]. Lack of teratoma formation raise a question about pluripotency of these cell lines. In our study, ES and iPS cells are pluripotent by this criterion. Both types of the cells form morphologically monolayer colonies. In mouse and human pluripotent stem cells, this morphology is characteristic of primed pluripotent state. However, the morphology of mink ES and iPS cells is very similar to that of mink blastocyst inner cell mass [26]. This characteristic could be species-specific and does not point out a pluripotent state.

Expression of pluripotency-associated genes was assessed based on transcriptome data showing that both mink ES and iPS cells represent one cell type. Moreover, we observed expression of several known key pluripotency genes, including Rex1 (Figure 2b), which is considered a marker of naïve pluripotent state [40]. One notable exception is the lack of Nanog expression as only several copies of its transcript were found in solitary replicates of pluripotent cell lines. This finding was also verified by qPCR. Taking into consideration the proven pluripotency of analyzed cell lines, we presume that the low level of Nanog expression is a species-specific feature. Interestingly, Gata6 is expressed in all four pluripotent cell lines. Gata6 is Nanog's antagonist in mouse early embryo development, Nanog-positive cells tend to form epiblast while Gata6-positive - primitive endoderm [41]. Thus, coexpression of Rex1 and Nanog could be considered contradictory.

We have also analyzed the expression of several genes such as Cer1 [42] and Otx2 [43] that are characteristic of mouse primed pluripotent cells (Figure 2b). We can observe the differences in these genes expression between the different ES cell lines as well as when compared to the iPS cell lines. Interestingly, the comparative levels of Oct4, Sox2 and Rex1 expression in MES12 and MES29 are reciprocal to Cer1 and Otx2. It might point out to different pluripotent states of these ES cell lines. Nevertheless, based on these gene expression levels we cannot assess the pluripotent state of the analyzed cell lines. It was shown that in mouse these genes are expressed both in ES and epiblast stem cells but on different levels [42, 44]. Due to the fact that we do not have a control with known pluripotency status the expression itself is not an indicator.

As it was shown for various mouse pluripotent cell lines, addition of 2i could shift primed cells into naïve [33, 34]. Interestingly, to produce and culture canine pluripotent cells investigators used supplementation with substantially different factors, e.g. LIF as used for mouse ES cells with bFGF as for human ES cells [12, 14, 15, 38]. In addition, some groups were able to obtain pluripotent cells using 2i + LIF + bFGF [16] and LIF + bFGF + 2i + valproic acid + TGH-β antagonist A83-01 [11]. Some researchers used mix of all mentioned factors for iPS cell production but cultured iPS cells with LIF only [13]. To test whether the change of culture condition could change morphology of mink iPS colonies we applied various combinations: 2i, (2i + LIF), (2i + bFGF) and (2i + LIF + bFGF) respectively to iNV11 cells for two weeks. The morphology of the colonies remained unchanged. If mink iPS cells are in primed pluripotent state, it maybe that additional factors are needed to shift it to naïve. Alternatively, they could already be in naïve state as indicated by Rex1 expression.


We produced and characterized American mink ES and iPS cells. These cell lines have diploid chromosome number, and are pluripotent based on teratoma formation test. The transcriptome analysis shows efficient reprogramming of the mink EF genome to the pluripotent state in iPS cells. Colony morphology and expression of several marker genes are not enough to conclude whether the cells are in naïve or primed pluripotent state. We have found that Nanog is nearly absent in these pluripotent stem cells and consider it as species-specific feature.


Production of mink embryonic fibroblasts

Primary EF of American mink were obtained from individual 29-day embryos by standard protocol [45]. Mink of wild type genotype were used as donors of embryos. The EF culture medium consisted of DMEM (Invitrogen, USA) supplemented with 10% fetal bovine serum (Invitrogen, USA), and 1x penicillin and streptomycin (Invitrogen, USA).

Production of mink ES cell lines

To produce mink ES cells, the previously published protocol was followed [25]. Embryos were obtained from Public Center "Fur-bearing and farm animals" of Federal State Budget Scientific Institution "The Federal Research Center Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences" (ICG SB RAS), Novosibirsk, Russia. Briefly, embryos of standard (wild type) genotypes at morula and early blastocyst stage were plated on plastic dishes coated with 0.1% gelatin on mitomycin C inactivated mink EF. Zona pellucidae of embryos was previously removed by treatment in 0.5% pronase solution. Within a few days the embryos attached to the feeder layer of EF and formed colonies of morphologically homogeneous cells similar to the ICM cells. These primary colonies were passaged by trypsinization with 0.25% Trypsin-EDTA (Invitrogen, USA) on the fresh feeder. ES cell culture medium contained α-MEM (Invitrogen, USA) with 20% ES cell qualified FBS (Invitrogen, USA), 1x NEAA (Invitrogen, USA), 1x GlutaMAX (Invitrogen, USA), 0.1 mM β-mercaptoethanol (Sigma, USA) and 1x Penicillin-Streptomycin (Invitrogen, USA). For subsequent culture we used 15% ES cell qualified FBS.

Production of mink iPS cell lines

To produce iPS cells from the mink EF we used lentiviral vectors LeGO ( with GFP and human reprogramming transcription factors: OCT4, SOX2, C-MYC and KLF4, courtesy of Dr. Sergei L. Kiselev. Lentiviruses were produced in Phoenix cell line using Lipofectamine LTX (Invitrogen, USA) according to manufacturer's recommendations. Multiplicity of infection was estimated as 4.8 using GFP lentiviral vector. Mink EF (3 × 105 cells, 15 × 103 cells/cm2) plated the day before were transduced with viruses containing four reprogramming transcription factors and 4 µg/ml Polybrene. Transduction was performed three consecutive days, on second and third day with halved amount of C-MYC lentivirus. Until day 10 the medium was changed daily with addition of 1 mM valproic acid. On day 7 transduced cells were seeded onto 10 cm culture dishes on feeder in ES cell culture medium. From days 13 to 20 colonies with mink ES cell-like morphology were picked up and expanded. All cell cultures were maintained at 37°C and 5% CO2.

All animal studies were undertaken with prior approval from Interinstitutional Bioethical Committee of ICG SB RAS.

Cytogenetic analysis

Cytogenetic analysis for ES was carried out on passages 7-23 and for iPS cell lines on passages 4-6. Preparation of metaphase chromosomes from ES cells was performed as previously described with minor modifications: cells were treated with hypotonic solution (0.56% KCl for 15 min) and after that fixated with methanol/acetic acid (3:1) solution [46]. For each cell line an average of 30 metaphase plates were analyzed on Carl Zeiss Axioscop 2 imaging microscope with CoolCube 1 CCD-camera (Meta Sistems). Digital images were analyzed using ISIS (In Situ Imaging System, MetaSystems GmbH) software in collective Microscopic Center of ICG SB RAS, Novosibirsk, Russia.

Immunofluorescence analysis of X-chromosome inactivation

Immunofluorescence analysis was carried out according to previously published protocol [47]. In short, cells cultured on glass coverslips were fixed with 3% formaldehyde (Fluka, Germany) and permeabilized with 0.1% Triton X-100 (Fluka, Germany). Non-specific binding was blocked with 2% bovine serum albumin (Sigma, USA) and cells were incubated with rabbit polyclonal or anti-H3K27me3 antibodies (1:500) (Molecular Probes, USA). The primary antibodies were visualized with the secondary goat anti-rabbit IgG antibodies conjugated with Alexa Fluor 488 (1:400) (Molecular Probes, USA). The cells were stained with DAPI, mounted in a glycerol solution containing 1,4 diazobicyclo-[2.2.3] octane (DABCO) (Sigma, USA) and visualized on LSM 780 NLO (Zeiss) based on AxioObserver Z1 (Zeiss) using ZEN software in collective Microscopic Center of ICG SB RAS, Novosibirsk.

Teratoma formation analysis

For teratoma formation, we used SCID hairless outbred mice (Crl:SHO-PrkdcscidHrhr) of SPF status and BALB/c-nu. Experiments with SCID mice were performed in the Center for Genetic Resources of Laboratory Animals (RFMEFI61914X0005) at ICG SB RAS; BALB/c-nu mice were kept in the Vivarium for conventional animals at ICG SB RAS. Teratomas were produced using standard protocol [45]. Between 3 and 17.7 × 106 ES cells on passages 7-23 and 1.5 to 7 × 106 iPS cells on passages 5-14 were injected subcutaneously into immunodeficient mice. Teratomas generated by ES or iPS cells were dissected after 3-12 weeks and were fixed in Bouin solution. Paraffin sections were prepared according to standard protocol and were stained with histological dyes Picro-Mallory trichromica (04-021822), Masson trichromica (04-011802), P.T.A.H.-hematoxyline (04-060802), Luxol fast blue Krever Barrera (04-200812), Azan trichromica (04-001802), Picrofuchsin Van Gizon (04-030802) (Bio-Optica Milano S.P.A., Italy) and with hematoxylin-eosin. Images were analyzed on Carl Zeiss Axioscop 2+ microscope with AxioCam HRc CCD-camera. Digital images were taken using AxioVision software in collective Microscopic Center of ICG SB RAS.

All animal studies were undertaken with prior approval from Interinstitutional Bioethical Committee of ICG SB RAS.

RNA isolation and cDNA synthesis

RNA was isolated from each sample using TriReagent® (MRC Inc., USA) according to manufacturer's recommendations. We had three technical replicates for each cell line (a, b and c, resulting in 15 samples in total). Genomic DNA was removed using RNeasy MinElute Cleanup kit (Qiagen, Germany). Quantity and quality were assessed in a Nanodrop ND-1000 machine (Thermo Scientific). The RNA integrity was assessed by gel electrophoresis and by an Experion system (BioRad Laboratories) using the ExperionTM RNA StdSens analysis kit (BioRad, Sweden). Average RNA quality indicator values (RQI) were 9.7 for the American mink ES cells, 9.6 for iPS cells and 9.3 for EF. This RNA was used for RNAseq library preparation and for cDNA synthesis, gene silencing analysis and qPCR.

One microgram of DNase I treated total RNA was used for cDNA synthesis [48]. cDNA synthesis was done in duplicates for each RNA sample, thus resulting in 30 samples. All samples were diluted eight times before using in qPCR.

Primer design

We used Primer3Plus software to design primers for the following mink genes (Additional file 3): Sox2, Oct4, Gdf3, Nanog and Nestin. For the reference genes Gapdh and Hprt1 mink-specific primers were from Rouvinen-Watt et al. [49].

Human-specific primers for C-MYC and KLF4 were from Mathew et al. [50]. For OCT4 and SOX2 primer design was based on the maximum number of mismatches between human and mink sequences.


qPCR was performed as described by Cirera et al. [51]. Briefly, PCR efficiency was calculated from the log-linear portion of the standard curve for each assay, which was done with dilution serials of a pool of cDNA from all the samples to test. Efficiencies between 80 and 110% were accepted. GeneEx Professional (MultiD) software was used for pre-processing the qPCR data. Normalisation was undertaken using the two stable reference genes according to GeNorm and NormFinder software [52, 53], namely GAPDH and HPRT1. Subsequently, technical replicates were averaged and relative quantities (fold changes, FC) were calculated based in the less expressed sample for each assay. Prior to statistical analysis, FC were log2 transformed to ensure normal distribution. Transcriptional differences between cell lines were evaluated using analysis of variance (ANOVA). A post hoc test (Tukey-Kramer) was performed in order to see which pairwise comparisons were significant. P values below 0.05 were considered significant. Results are reported as mean of the Log2 of the FC. The qPCR experiments, as well as the data analysis, were all compliant with MIQE guidelines [54].

RNA-Seq data processing, reference assembly and alignment

Prior to library preparation the RNA quality and integrity was assessed according to Illumina guidelines. Library preparation was done using the TruSeq® Stranded mRNA sample preparation 96 rxn kit (Illumina™) following the low sample protocol according to manufacturer's recommendations. Briefly, approximately 2.5 µg of total RNA was diluted and purified using RNA purification beads targeting the poly-A tail of the mRNA and subsequently was fragmented by means of the enzymes provided in the kit. After the cDNA synthesis adenylation of 3' ends and ligation of the adaptors were performed. Adaptors were ligated in 12-plex formations, allowing the pooling of 12 samples after the PCR enrichment of the library. Subsequently, the library was quantified using PicoGreen® dye (Life Technologies™) as described in the manufacturer's protocol. Thereafter 12 samples were pooled at equal concentrations to create the eight pools. In order to accurately quantify the concentration in nM of our pools, the Kapa SYBR® FAST universal qPCR kit (Kapa Biosystems™) for Illunima™ sequencing was used to quantify the number of the amplifiable molecules in the pools and the Bioanalyzer® machine (Agilent Technologies™) to determine the average fragment size of the pools. These measurements allowed optimizing the flow cell clustering and proceed with the Run. The samples were sequenced in five lanes of an Illumina HiSeq2000 sequencer, for 50 bp pair-end read.

Since mink pluripotent cells were maintained on mouse feeder cells raw reads produced by next generation sequencing might contain transcripts of mouse feeder origin. Considering this, we aligned all data produced from iPS and ES cells sequencing to mm10-based mouse transcriptome (obtained from Illumina iGenomes project, using bowtie2 with parameters "--no-unal --no-sq --no-hd --np 0 --rdg 5,1 --rfg 5,1 --score-min L,-6.1,0 -X 500 --mp 3,2". All concordantly aligned read-pairs were considered as of mouse origin and eliminated from the analysis.

We used filtered data from mink iPS and ES cells as well as raw data from mink fibroblasts to perform de novo transcriptome assembly. Transcriptome was assembled using Trinity [55] with default parameters. We blasted Trinity-produced data against mouse transcripts and manually checked all highly similar (similarity >= 97%) sequences to ensure absence of contamination.

Annotation of de novo assembled contigs was done by alignment coding sequences against the human RefSeq mRNA database (Homo Sapiens GRCh37.66.cdna; The contigs that did not aligned against the human RefSeq were then blasted against the Mustela putorius furo protein sequences (Ensemble Genes 81; The annotation was made only for sequences translating to minimum of 50 amino acids in length. For efficient assembly, analysis and annotation of the data, scripts developed in-house and described earlier were used [56].

To perform differential expression we used mouse-filtered mink iPS and ES cells data. To make it comparable with fibroblasts data we also performed filtering of fibroblasts reads aligning it to mouse transcriptome as described above.

Resulting filtered sequencing data was aligned to the transcriptome assembly using bowtie [57] with options "-aS -X 800 --offrate 1". Produced alignment results were passed to the eXpress software [58] to obtain bias-corrected counts. To analyze expression on gene-level counts of all contigs representing the same gene were summed using self-made python script. Obtained counts were analyzed using DESeq2 package, that included library size normalization, dispersions estimation, differential expression tests etc. [59].

Availability of supporting data

The raw sequenced data sets supporting the results of this article are available in the NCBI BioProject repository,


  1. Takahashi K, Yamanaka S: Induction of Pluripotent Stem Cells from Mouse Embryonic and Adult Fibroblast Cultures by Defined Factors. Cell. 2006, 126 (4): 663-676.

    Article  CAS  Google Scholar 

  2. Takahashi K, Tanabe K, Ohnuki M, Narita M, Ichisaka T, Tomoda K, Yamanaka S: Induction of Pluripotent Stem Cells from Adult Human Fibroblasts by Defined Factors. Cell. 2007, 131 (5): 861-872.

    Article  CAS  Google Scholar 

  3. Yu J, Vodyanik MA, Smuga-Otto K, Antosiewicz-Bourget J, Frane JL, Tian S, et al: Induced pluripotent stem cell lines derived from human somatic cells. Science. 2007, 318 (5858): 1917-1920.

    Article  CAS  Google Scholar 

  4. Marchetto MC, Narvaiza I, Denli AM, Benner C, Lazzarini TA, Nathanson JL, et al: Differential L1 regulation in pluripotent stem cells of humans and apes. Nature. 2013, 503 (7477): 525-529.

    Article  CAS  Google Scholar 

  5. Wu Y, Zhang Y, Mishra A, Tardif SD, Hornsby PJ: Generation of induced pluripotent stem cells from newborn marmoset skin fibroblasts. Stem Cell Res. 2010, 4 (3): 180-188.

    Article  CAS  Google Scholar 

  6. Wunderlich S, Kircher M, Vieth B, Haase A, Merkert S, Beier J, et al: Primate iPS cells as tools for evolutionary analyses. Stem Cell Res. 2014, 12 (3): 622-629.

    Article  Google Scholar 

  7. Ben-Nun IF, Montague SC, Houck ML, Tran HT, Garataonandia I, Leonardo TR: Induced pluripotent stem cells from highly endangered species. Nature Methods. 2011, 8 (10): 2-6.

    Google Scholar 

  8. Li W, Wei W, Zhu S, Zhu J, Shi Y, Lin T, et al: Generation of Rat and Human Induced Pluripotent Stem Cells by Combining Genetic Reprogramming and Chemical Inhibitors. Cell Stem Cell. 2009, 4 (1): 370-

    Article  CAS  Google Scholar 

  9. Manoli DS, Subramanyam D, Carey C, Sudin E, van Westerhuyzen JA, Bales KL, et al: Generation of induced pluripotent stem cells from the prairie vole. PLoS One. 2012, 7 (5): e38119-

    Article  CAS  Google Scholar 

  10. Honda A, Hirose M, Hatori M, Matoba S, Miyoshi H, Inoue K, Ogura A: Generation of induced pluripotent stem cells in rabbits: potential experimental models for human regenerative medicine. J Biol Chem. 2010, 285 (41): 31362-31369.

    Article  CAS  Google Scholar 

  11. Shimada H, Nakada A, Hashimoto Y, Shigeno K, Shionoya Y, Nakamura T: Generation of canine-induced pluripotent stem cells by retroviral transduction and chemical inhibitors. Mol Reprod Dev. 2010, 77 (1): 2-

    Article  CAS  Google Scholar 

  12. Luo J, Suhr ST, Chang EA, Wang K, Ross PJ, Nelson LL, et al: Generation of Leukemia Inhibitory Factor and Basic Fibroblast Growth Factor-Dependent Induced Pluripotent Stem Cells from Canine Adult Somatic Cells. Stem Cells Dev. 2011, 20 (10): 1669-1678.

    Article  CAS  Google Scholar 

  13. Whitworth DJ, Ovchinnikov DA, Wolvetang EJ: Generation and Characterization of LIF-dependent Canine Induced Pluripotent Stem Cells from Adult Dermal Fibroblasts. Stem Cells Dev. 2012, 21 (12): 2288-2297.

    Article  CAS  Google Scholar 

  14. Baird AEG, Barsby T, Guest DJ: Derivation of Canine Induced Pluripotent Stem Cells. Reprod Domest Anim. 2015, 50 (4): 669-676.

    Article  CAS  Google Scholar 

  15. Lee AS, Xu D, Plews JR, Nguyen PK, Nag D, Lyons JK, et al: Preclinical derivation and imaging of autologously transplanted canine induced pluripotent stem cells. J Biol Chem. 2011, 286 (37): 32697-32704.

    Article  CAS  Google Scholar 

  16. Koh S, Thomas R, Tsai S, Bischoff S, Lim J-H, Breen M, et al: Growth requirements and chromosomal instability of induced pluripotent stem cells (iPSC) generated from adult canine fibroblasts. Stem Cells Dev. 2012, 22 (6): 951-963.

    Article  Google Scholar 

  17. Wu Z, Chen J, Ren J, Bao L, Liao J, Cui C, et al: Generation of pig induced pluripotent stem cells with a drug-inducible system. J Mol Cell Biol. 2009, 1 (1): 46-54.

    Article  CAS  Google Scholar 

  18. Nagy K, Sung HK, Zhang P, Laflamme S, Vincent P, Agha-Mohammadi S, et al: Induced pluripotent stem cell lines derived from equine fibroblasts. Stem Cell Rev. 2011, 7 (3): 693-702.

    Article  Google Scholar 

  19. Li Y, Cang M, Lee AS, Zhang K, Liu D: Reprogramming of sheep fibroblasts into pluripotency under a drug-inducible expression of mouse-derived defined factors. PLoS One. 2011, 6 (1): e15947-

    Article  CAS  Google Scholar 

  20. Han X, Han J, Ding F, Cao S, Lim SS, Dai Y, et al: Generation of induced pluripotent stem cells from bovine embryonic fibroblast cells. Cell Res. 2011, 21 (10): 1509-1512.

    Article  CAS  Google Scholar 

  21. Ren J, Pak Y, He L, Qian L, Gu Y, Li H, et al: Generation of hircine-induced pluripotent stem cells by somatic cell reprogramming. Cell Res. 2011, 21 (5): 849-853.

    Article  CAS  Google Scholar 

  22. Deng Y, Liu Q, Luo C, Chen S, Li X, Wang C, et al: Generation of induced pluripotent stem cells from buffalo (Bubalus bubalis) fetal fibroblasts with buffalo defined factors. Stem Cells Dev. 2012, 21 (13): 2485-2494.

    Article  CAS  Google Scholar 

  23. Verma R, Holland MK, Temple-Smith P, Verma PJ: Inducing pluripotency in somatic cells from the snow leopard (Panthera uncia), an endangered felid. Theriogenology. 2012, 77 (1): 220-228. e1-228.e2

    Article  CAS  Google Scholar 

  24. Bock C, Kiskinis E, Verstappen G, Gu H, Boulting G, Smith ZD, et al: Reference maps of human es and ips cell variation enable high-throughput characterization of pluripotent cell lines. Cell. 2011, 144 (3): 439-452.

    Article  CAS  Google Scholar 

  25. Sukoyan MA, Vatolin SY, Golubitsa AN, Zhelezova AI, Semenova LA, Serov OL: Embryonic stem cells derived from morulae, inner cell mass, and blastocysts of mink: comparisons of their pluripotencies. Mol Reprod Dev. 1993, 36 (2): 148-158.

    Article  CAS  Google Scholar 

  26. Sukoyan MA, Golubitsa AN, Zhelezova AI, Shilov AG, Vatolin SY, Maximovsky LP, et al: Isolation and cultivation of blastocyst-derived stem cell lines from American mink (Mustela vison). Mol Reprod Dev. 1992, 33 (4): 418-431.

    Article  CAS  Google Scholar 

  27. Mantel C, Guo Y, Lee MR, Kim MK, Han MK, Shibayama H, et al: Checkpoint-apoptosis uncoupling in human and mouse embryonic stem cells: a source of karyotpic instability. Blood. 2007, 109 (10): 4518-4527.

    Article  CAS  Google Scholar 

  28. Maherali N, Sridharan R, Xie W, Utikal J, Eminli S, Arnold K, et al: Directly Reprogrammed Fibroblasts Show Global Epigenetic Remodeling and Widespread Tissue Contribution. Cell Stem Cell. 2007, 1 (1): 55-70.

    Article  CAS  Google Scholar 

  29. Silva SS, Rowntree RK, Mekhoubad S, Lee JT: X-chromosome inactivation and epigenetic fluidity in human embryonic stem cells. Proc Natl Acad Sci USA. 2008, 105 (12): 4820-4825.

    Article  CAS  Google Scholar 

  30. Wernig M, Meissner A, Foreman R, Brambrink T, Ku M, Hochedlinger K, et al: In vitro reprogramming of fibroblasts into a pluripotent ES-cell-like state. Nature. 2007, 448: 318-324.

    Article  CAS  Google Scholar 

  31. Okita K, Ichisaka T, Yamanaka S: Generation of germline-competent induced pluripotent stem cells. Nature. 2007, 448: 313-317.

    Article  CAS  Google Scholar 

  32. Ohi Y, Qin H, Hong C, Blouin L, Polo JM, Guo T, et al: Incomplete DNA methylation underlies a transcriptional memory of somatic cells in human iPS cells. Nat Cell Biol. 2011, 13 (5): 541-549.

    Article  CAS  Google Scholar 

  33. Nichols J, Smith A: Naive and Primed Pluripotent States. Cell Stem Cell. 2009, 4 (6): 487-492.

    Article  CAS  Google Scholar 

  34. Hanna J, Saha K, Jaenisch R: Somatic cell reprogramming and transitions between pluripotent states: facts, hypotheses, unresolved issues. Cell. 2010, 143 (4): 508-525.

    Article  CAS  Google Scholar 

  35. Hatoya S, Torii R, Kondo Y, Okuno T, Kobayashi K, Wijewardana V, et al: Isolation and characterization of embryonic stem-like cells from canine blastocysts. Mol Reprod Dev. 2006, 73 (3): 298-305.

    Article  CAS  Google Scholar 

  36. Schneider MR, Wolf E, Braun J, Kolb HJ, Adler H: Canine embryonic stem cells: State of the art. Theriogenology. 2010, 74 (4): 492-497.

    Article  CAS  Google Scholar 

  37. Hayes B, Fagerlie SR, Ramakrishnan A, Baran S, Harkey M, Graf L, et al: Derivation, characterization, and in vitro differentiation of canine embryonic stem cells. Stem Cells. 2008, 26 (2): 465-473.

    Article  CAS  Google Scholar 

  38. Vaags AK, Rosic-Kablar S, Gartley CJ, Zheng YZ, Chesney A, Villagomez DA, et al: Derivation and characterization of canine embryonic stem cell lines with in vitro and in vivo differentiation potential. Stem Cells. 2009, 27 (2): 329-340.

    Article  CAS  Google Scholar 

  39. Wilcox JT, Semple E, Gartley C, Brisson BA, Perrault SD, Villagómez DA, et al: Characterization of canine embryonic stem cell lines derived from different niche microenvironments. Stem Cells Dev. 2009, 18 (8): 1167-1178.

    Article  CAS  Google Scholar 

  40. Pelton TA, Sharma S, Schulz TC, Rathjen J, Rathjen PD: Transient pluripotent cell populations during primitive ectoderm formation: correlation of in vivo and in vitro pluripotent cell development. J Cell Sci. 2002, 115 (Pt 2): 329-339.

    CAS  Google Scholar 

  41. Cockburn K, Rossant J: Making the blastocyst: Lessons from the mouse. J Clin Invest. 2010, 120 (4): 995-1003.

    Article  CAS  Google Scholar 

  42. Factor DC, Najm FJ, Tesar PJ: Generation and characterization of epiblast stem cells from blastocyst-stage mouse embryos. Methods Mol Biol. 2013, 1074: 1-13.

    Article  CAS  Google Scholar 

  43. Kalkan T, Smith A: Mapping the route from naive pluripotency to lineage specification. Philos Trans R Soc Lond B Biol Sci. 2014, 369 (1657): pii: 20130540

    Google Scholar 

  44. Acampora D, Di Giovannantonio LG, Simeone A: Otx2 is an intrinsic determinant of the embryonic stem cell state and is required for transition to a stable epiblast stem cell condition. Development. 2013, 140 (1): 43-55.

    Article  CAS  Google Scholar 

  45. Hogan B, Beddington R, Costantini F, Lacy L: Manipulating the Mouse Embryo. 1994, Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press, 2nd

    Google Scholar 

  46. Pristyazhnyuk IE, Matveeva NM, Graphodatsky AS, Serdukova NA, Serov OL: Analysis of chromosome composition in interspecific embryonic stem hybrid cells of mice. Tsitologiya. 2010, 52 (2): 136-143.

    Google Scholar 

  47. Gridina MM, Serov OL: Bidirectional reprogramming of mouse embryonic stem cell/fibroblast hybrid cells is initiated at the heterokaryon stage. Cell Tissue Res. 2010, 342 (3): 377-389.

    Article  Google Scholar 

  48. Nygard A-B, Jørgensen CB, Cirera S, Fredholm M: Selection of reference genes for gene expression studies in pig tissues using SYBR green qPCR. BMC Mol Biol. 2007, 8: 67-

    Article  Google Scholar 

  49. Rouvinen-Watt K, Harris L, Dick M, Pal C, Lei S, Mustonen AM, Nieminen P: Role of hepatic de novo lipogenesis in the development of fasting-induced fatty liver in the American mink (Neovison vison). British Journal of Nutrition. 2012, 108 (8): 1360-1370.

    Article  CAS  Google Scholar 

  50. Mathew R, Jia W, Sharma A, Zhao Y, Clarke LE, Cheng X, Wang H, et al: Robust activation of the human but not mouse telomerase gene during the induction of pluripotency. FASEB J. 2010, 24 (8): 2702-2715.

    Article  CAS  Google Scholar 

  51. Cirera S, Jensen MS, Elbrønd VS, Moesgaard SG, Christoffersen B, Kadarmideen HN, et al: Expression studies of six human obesity-related genes in seven tissues from divergent pig breeds. Anim Genet. 2014, 45 (1): 59-66.

    Article  CAS  Google Scholar 

  52. Andersen CL, Jensen JL, Ørntoft TF: Normalization of real-time quantitative reverse transcription-PCR data: A model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004, 64 (15): 5245-5250.

    Article  CAS  Google Scholar 

  53. Mestdagh P, Van Vlierberghe P, De Weer A, Muth D, Westermann F, Speleman F, Vandesompele J: A novel and universal method for microRNA RT-qPCR data normalization. Genome Biol. 2009, 10 (6): R64-

    Article  Google Scholar 

  54. Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al: The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin. 2009, 55 (4): 611-622.

    CAS  Google Scholar 

  55. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29: 644-652.

    Article  CAS  Google Scholar 

  56. Christensen K, Anistoroaei R: An American mink (Neovison vison) transcriptome. Anim Genet. 2014, 45 (2): 301-303.

    Article  CAS  Google Scholar 

  57. Langmead B, Salzberg SL: Fast gapped-read alignment with Bowtie 2. Nature Methods. 2012, 9 (4): 357-359.

    Article  CAS  Google Scholar 

  58. Roberts A, Pachter L: Streaming fragment assignment for real-time analysis of sequencing experiments. Nat Methods. 2013, 10: 71-73.

    Article  CAS  Google Scholar 

  59. Love MI, Huber W, Anders S: Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. 2014

    Book  Google Scholar 

Download references


We acknowledge Dr. Sergei L. Kiselev (VIGG RAS, Moscow) for providing lentiviral vectors, Dr. Oleg V. Trapezov (ICG SB RAS, Novosibirsk) for providing American mink embryos and Antonina I. Zhelezova (ICG SB RAS, Novosibirsk) for work with mice. We also thank Karla K. Freude from University of Copenhagen (Denmark) for providing human iPS samples and for the good advice. We appreciate the help of Novosibirsk State University IT-center administration who provided us server cluster at Publication of this article was supported by state project No. VI.60.1.3 of Russian Academy of Sciences.

This article has been published as part of BMC Genomics Volume 16 Supplement 13, 2015: Selected articles from the 7th International Young Scientists School "Systems Biology and Bioinformatics" (SBB'2015): Genomics. The full contents of the supplement are available online at

Author information

Authors and Affiliations


Corresponding author

Correspondence to Aleksei G Menzorov.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

AGM conceived the study, participated in cell culture, provided analysis and interpretations of the data and is the principal investigator of the project; NMM participated in cell culture and produced all mink ES cells; MNM carried out Illumina libraries preparation, sequence processing and analysis; VSF carried out sequence processing and transcriptome-based gene expression analysis; KC carried out bioinformatics tasks, sequence analysis and gene annotation; AAK participated in cell culture; IEP carried out cytogenetic and immunocytochemical analyses, participated in cell culture; EAK carried out work with animals and teratoma histochemical analysis; SC carried out RNA preparation, qPCR analysis and interpretation; RA carried out bioinformatics tasks, sequence analysis and gene annotation, project coordination in Western Europe; OLS carried out interpretations of the data and project coordination in Russia. AGM did most of the writing with contributions from all authors. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Mink EF and MES25 stained with H3K27me3 antibodies and visualized with the secondary antibodies conjugated with Alexa Fluor 488, counterstained with DAPI. a - mink EF; b - MES25. (Additional file 1.jpg) (JPG 700 KB)


Additional file 2: Histological analysis of individual teratomas formed after injection of ES and iPS cells into immunodeficient mice. (XLSX 14 KB)

Additional file 3: Primers for qPCR and RT-PCR. (DOCX 16 KB)

Additional file 4: qPCR raw data of selected gene expression. (DOCX 14 KB)


Additional file 5: qPCR data of selected gene expression corrected by PCR efficiency and normalized to reference genes Gapdh and Hprt1. (DOCX 14 KB)


Additional file 6: P values of the pairwise comparison (Tukey-Kramer post-test) of selected gene expression. (DOCX 14 KB)

Additional file 7: Number of pair-end reads and base in mink cell lines. (DOCX 13 KB)

Additional file 8: Lists of differentially expressed genes between mink EF, ES and iPS cells. (ZIP 126 KB)


Additional file 9: Expression levels of selected pluripotency-associated genes in mink EF, ES and iPS cells. Vertical axis represents counts determined for each sample by transcriptome analysis. (PDF 15 KB)


Additional file 10: Expression levels of shared between iNV7 and iNV11 genes with "novel" expression pattern in EF, ES and iPS cells. Vertical axis represents counts determined for each sample by transcriptome analysis. (PDF 69 KB)

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Menzorov, A.G., Matveeva, N.M., Markakis, M.N. et al. Comparison of American mink embryonic stem and induced pluripotent stem cell transcriptomes. BMC Genomics 16 (Suppl 13), S6 (2015).

Download citation

  • Published:

  • DOI: