Cattle In Vitro Induced Pluripotent Stem Cells Generated and Maintained in 5 or 20% Oxygen and Different Supplementation

The event of cellular reprogramming into pluripotency is influenced by several factors, such as in vitro culture conditions (e.g., culture medium and oxygen concentration). Herein, bovine iPSCs (biPSCs) were generated in different levels of oxygen tension (5% or 20% of oxygen) and supplementation (bFGF or bFGF + LIF + 2i—bFL2i) to evaluate the efficiency of pluripotency induction and maintenance in vitro. Initial reprogramming was observed in all groups and bFL2i supplementation initially resulted in a superior number of colonies. However, bFL2i supplementation in low oxygen led to a loss of self-renewal and pluripotency maintenance. All clonal lines were positive for alkaline phosphatase; they expressed endogenous pluripotency-related genes SOX2, OCT4 and STELLA. However, expression was decreased throughout the passages without the influence of oxygen tension. GLUT1 and GLUT3 were upregulated by low oxygen. The biPSCs were immunofluorescence-positive stained for OCT4 and SOX2 and they formed embryoid bodies which differentiated in ectoderm and mesoderm (all groups), as well as endoderm (one line from bFL2i in high oxygen). Our study is the first to compare high and low oxygen environments during and after induced reprogramming in cattle. In our conditions, a low oxygen environment did not favor the pluripotency maintenance of biPSCs.


Introduction
The success of reprogramming somatic cells into pluripotent cells, or iPSCs, has long been reported and studied, particularly in humans and mice. This technology has provided researchers with a valuable source of different cell types [1][2][3] to use in regenerative medicine, as well as for biotechnology and reproductive purposes in domestic and wild animals [4].
Pluripotent stem cells are mainly classified into two different states, naïve and primed. Naïve ESCs and iPSCs, mostly observed in murine models, rely on leukemia inhibitory factor (LIF) supplementation to retain pluripotency, whereas primed mouse EpiSC and human ESCs and iPSCs rely on fibroblast growth factor 2 (FGF2) stimulation [5]. Additionally, two small molecules are widely used to improve the pluripotent state by blocking the mitogen-activated protein kinase (MEK) and Glycogen synthase kinase 3-beta (GSKB3) activity pathways; these molecules are known as 2i (2 inhibitors) [6].
iPSCs from large animals are important models for both biomedical and agricultural applications due to their physiological and morphological similarity to humans and, in particular, cattle genetically manipulated by different genome editing tools would provide a powerful platform for the production of modified herds with superior genetic traits for production of reproduction. Indeed, multiple genetic alterations in livestock could be used to improve traits connected with animal production and create large animal models for in vitro and in vivo disease modeling [4,[7][8][9].
The generation of livestock iPSCs has been reported, including in sheep [10,11], horses [12,13], pigs [14][15][16] and cattle [17][18][19]. These livestock iPSCs revealed pluripotent characteristics, such as a high proliferation capacity, the expression of pluripotency marks and the ability to generate three germ layers tissues on in vitro and in vivo tests; however, most of the characteristics and protocols are still divergent and often describe either an intermediate or "pluripotency-like" state [20].
Different methodologies have been explored to improve efficiency during the reprogramming process. The in vitro culture environment is one promising field that enables the use of different molecules as inhibitors and supplements and allows the control of oxygen tension. In mouse ESCs, hypoxia conditions (5% oxygen) increased the number of colonies with alkaline phosphatase activity, enhanced cell proliferation and prevented differentiation [21]. Additionally, for mouse and human iPSCs, hypoxia increased reprogramming efficiency by promoting an increased number of colonies [22]. However, livestock iPSCs have not yet been described on low oxygen.
This study aimed to reprogram bovine fibroblasts under different culture conditions and oxygen tensions to generate biPSCs. These were then evaluated through their pluripotency-related transcripts at different passages, immunostaining and in vitro differentiation capacity. As such, we aimed to obtain a better understanding of the failures in the establishment of the pluripotency state in biPSCs that are mainly associated with mechanisms of self-renewal and related to in vitro culture conditions.

Materials and Methods
All experiments were approved by the Ethics Committee on Animal Experimentation at the Faculty of Animal Science and Food Engineering of the University of São Paulo (CEUA 2192250918 and 3526250717).

Feeder Cells
MEFs were isolated from 13.5d mouse fetuses and cultured on complete IMDM. Cells were treated with mitomycin (M4287, Sigma-Aldrich, St. Louis, MO, USA) and used as a feeder monolayer at 1 × 10 5 cells in 6-well dishes. The dishes were treated with 0.1% gelatin (G1393, Sigma-Aldrich, St. Louis, MO, USA) prior to seeding.

Bovine Induced Pluripotent Stem Cells Generation
To reprogram fetal fibroblasts, the cells were seeded 24 h before lentivirus transduction at 2 × 10 4 per well on a 6-well plate. A fresh medium with lentivirus stock was added to the culture medium supplemented with 8 µg/mL hexadimethrine bromide (polybrene, Sigma-Aldrich, St. Louis, MO, USA. After overnight transduction, the culture medium was replaced by fresh medium. On day 5 after transduction, the cells were harvested with TrypLE (12604-021, Thermo Fisher Scientific, Waltham, MA, USA) and split into plates with feeder layers and iPSC medium supplemented with bFGF or bFGF + LIF + 2i (bFL2i). The culture was performed at two different oxygen levels (5%-low O 2 tension-and 20%-high O 2 ) at 38.5 • C and 5% CO 2 in a humidified atmosphere. Three clonal lineages were analyzed for each experimental group.

In Vitro Embryo Production
One pool of in vitro produced blastocysts was used as positive control for pluripotency analyses. For that, bovine ovaries were collected from a local abattoir and cumulusoocyte complexes (COCs) were aspirated from selected follicles (2-6 mm in diameter) using a 21-gauge needle. COCs were incubated in tissue culture medium 199 (TCM-199) supplemented with 10% FBS, 0.2 mM pyruvate, 50 µg/mL gentamicin, 0.5 mg/mL FSH and hCG, 5 IU/mL (Vetecor, Ceva Hertape, Brazil) for 22 h. The in vitro fertilization was performed with frozen-thawed semen from a single Holstein bull. Capacitated sperm were obtained after Percoll gradient (45% and 90%) separation. The fertilization medium was composed of Tyrode's lactate stock, 50 µg/mL gentamicin, 22 µg/mL sodium pyruvate, 40 µL/mL PHE (2 mM D-penicillamine, 1 mM hypotaurine and 245 µM epinephrine), 5.5 IU/mL heparin and 6 mg/mL bovine serum albumin (BSA). Presumptive zygotes were partially denuded after 18 h and cultured in SOF medium in an incubator with a humidified atmosphere of 5% CO 2 and 5% O 2 in air at 38.5 • C for 7 days (insemination is considered as day 0).

RNA Extraction and Reverse Transcription
All samples were submitted to total RNA extraction using the Trizol ® protocol (Thermo Fisher Scientific, Waltham, MA, USA). After extraction, total RNA samples were quantified with a NanoDrop ® (Thermo Fisher Scientific, Waltham, MA, USA). For DNA digestion and reverse transcription, we adjusted the concentration to 1.000 ng RNA per sample. All the samples were subjected to DNA digestion using a DNAse I-Amplification Grade ® (18068015, Thermo Fisher Scientific, Waltham, MA, USA). For reverse transcription (RT) of the samples, a high-capacity cDNA reverse transcription kit was utilized (4368814, Thermo Fisher Scientific, Waltham, MA, USA).

Quantitative PCR
For gene expression analysis on biPSCs (generated in 5% and 20% O 2 conditions and different supplements), samples were collected at five different moments during the reprogramming process (passages 5, 10, 15, 20 and 25). The targets were pluripotent and metabolism-related genes and the exogenous OSKM expression. The reference genes were ACTB and PPIA [24] and normalization was performed using the geometric means from both. The primers were designed using the Primer-BLAST (NCBI) software based upon sequences available in GenBank (Table 1). Relative quantification of transcribed genes was analyzed by qPCR using the PowerUP SYBR Green ® PCR Master Mix reagent (A25742, Thermo Fisher Scientific, Waltham, MA, USA) in Applied Biosystems QuantStudio 6 equipment. Then, qPCR reactions were run in a volume of 10 µL containing 100 nM of each primer, 1X Power UP SYBR Green, 2.5 µL H 2 O and 1 µL template (four-fold diluted cDNA; 25 ng). Cycling conditions for amplification were as follows: 95 • C for 10 min, followed by 45 cycles at 95 • C for 15 s, 57 • C for 20 s and 60 • C for 40 s. Each sample was analyzed in duplicate for each of the genes. Additionally, ultrapure DNAse and RNAse free water was used as a negative control of the reaction. The relative gene expression was determined by 2 −∆Cq .

Undirected Differentiation in Embryoid Bodies
Undirected differentiation of bovine iPSCs at p25 through embryoid bodies (EBs) formation was performed in plates treated with agarose 0.1% and iPSCs media without bFGF, 2i, or LIF. After 5 days in suspension culture, EBs were collected for qPCR analysis. EBs were placed on gelatin-coated culture plates and cultured in IMDM medium; EBs were documented 7 days after plating.

Statistical Analysis
Statistical analysis was performed using GraphPad Prism 8 (GraphPad Software, San Diego, CA, USA). Data were tested for normality of residuals and homogeneity of variances using the Shapiro-Wilk test and analyzed. Gene expression on biPSCs lines data were analyzed by two-way ANOVA followed by Tukey's post hoc test. Gene expression on EBs data were analyzed by one-way ANOVA followed by Tukey's post hoc test. Differences with probabilities of p < 0.05 were considered significant. In the text, values are presented as the mean ± the mean standard error (SEM). All experiments were repeated at least three times unless otherwise stated.

Influence of Oxygen Tension and Different Pluripotency Inducers during the Reprogramming Process of Bovine Fibroblasts
Colonies were observed on day 14 after transduction and the number of colonies that appeared in one well of each treatment was counted for each experimental group. A higher number of colonies was observed with bFL2i supplementation in both oxygen tensions (0.075 and 0.007%, respectively, for hypoxia and normoxia, meaning that 15 and 14 colonies were observed at D14 after plating 2 × 10 4 transduced cells per well; Table 2), whereas 0.02% and 0.005% were observed for the bFGF group (4 and 1 colonies observed at D14, respectively, for hypoxia and normoxia). Three colonies from each experimental group were selected at day 21 post-transduction for expansion and clonal culture. All lines were cultured for 25 passages (±200 days), with the exception of the clonal lines derived from the bFL2i group in hypoxia conditions, which lost its iPSCs characteristics after 8 passages, suggesting a loss of its self-renewal ability and the maintenance of pluripotency.
All biPSCs lines generated in low oxygen or high oxygen tensions were small and compact colonies, with clearly defined borders ( Figure 1) and were positive for alkaline phosphatase (Figure 1). Additionally, we did not observe morphological differences between the biPSCs groups.
All biPSCs lines generated in low oxygen or high oxygen tensions were small and compact colonies, with clearly defined borders ( Figure 1) and were positive for alkaline phosphatase (Figure 1). Additionally, we did not observe morphological differences between the biPSCs groups. Immunofluorescence revealed that all lineages were positive for SOX2 and OCT4 ( Figure S1); however, both genes are also known to be expressed from the exogenous cassette.

Molecular Profile of biPSCs Generated in Different Oxygen Tensions and Different Supplementations
The generated biPSCs were analyzed at five different time points during culture: p5, p10, p15, p20 and p25. At initial passages, (p5) no differences were observed. The SOX2 relative gene expression varied throughout all passages (p = 0.0074); however, experimental groups did not differ. Relative expression of OCT4 and STELLA increased in the bFGF group that was cultured in high O2 compared to low O2; in normoxia, bFGF supplementation showed higher expression than bFL2i for these genes. However, both genes decreased after passaging (p = 0.0287 and p = 0.0367 for OCT4 and STELLA, respectively) ( Figure 2). Immunofluorescence revealed that all lineages were positive for SOX2 and OCT4 ( Figure S1); however, both genes are also known to be expressed from the exogenous cassette.

Molecular Profile of biPSCs Generated in Different Oxygen Tensions and Different Supplementations
The generated biPSCs were analyzed at five different time points during culture: p5, p10, p15, p20 and p25. At initial passages, (p5) no differences were observed. The SOX2 relative gene expression varied throughout all passages (p = 0.0074); however, experimental groups did not differ. Relative expression of OCT4 and STELLA increased in the bFGF group that was cultured in high O 2 compared to low O 2 ; in normoxia, bFGF supplementation showed higher expression than bFL2i for these genes. However, both genes decreased after passaging (p = 0.0287 and p = 0.0367 for OCT4 and STELLA, respectively) ( Figure 2). SLC2A1 and SLC2A3 (also called GLUT1 and GLUT3) presented higher relative expression in hypoxic conditions (p = 0.0002 and p = 0.0014, respectively). In normoxia, no differences in bFGF x bFL2i supplementation were observed ( Figure 2). The exogenous expression of mOKSM was reduced throughout passaging (p = 0.0002); however, groups were not different (p = 0.6745) (Figure 2).

Undirected Differentiation of biPSCs Generated in Hypoxia and Normoxia Conditions and Different Pluripotent Supplements
In vitro differentiation of cells was performed via the undirected differentiation of embryoid body (EB) formation. All biPSCs lines generated in hypoxia and normoxia conditions and different pluripotent supplements were cultured in suspension for 5 days without supplementation. All lineages were able to produce EBs ( Figure 3A). Additionally, these EBs were evaluated by qPCR regarding ectoderm, mesoderm and endoderm markers. SLC2A1 and SLC2A3 (also called GLUT1 and GLUT3) presented higher relative expression in hypoxic conditions (p = 0.0002 and p = 0.0014, respectively). In normoxia, no differences in bFGF x bFL2i supplementation were observed (Figure 2). The exogenous expression of mOKSM was reduced throughout passaging (p = 0.0002); however, groups were not different (p = 0.6745) (Figure 2).

Undirected Differentiation of biPSCs Generated in Hypoxia and Normoxia Conditions and Different Pluripotent Supplements
In vitro differentiation of cells was performed via the undirected differentiation of embryoid body (EB) formation. All biPSCs lines generated in hypoxia and normoxia conditions and different pluripotent supplements were cultured in suspension for 5 days without supplementation. All lineages were able to produce EBs ( Figure 3A). Additionally, these EBs were evaluated by qPCR regarding ectoderm, mesoderm and endoderm markers.  The relative expression of the ectoderm marker (Vimentin-VIM) and mesoderm markers (BMP4 and PECAM1) was observed for all lines. However, we did not observe the relative expression of TUBB3 and FOXA2, both ectoderm and endoderm markers, respectively ( Figure 3B).
Interestingly, we observed the relative gene expression of AFP, an endoderm marker, in EBs derived from biPSCs culture in bFL2i in 20% O2 conditions only. Oxygen tension did not seem to affect cell differentiation in our experimental conditions; nevertheless, we observed a higher expression of PECAM1 and VIM on EBs produced from the bFL2i high The relative expression of the ectoderm marker (Vimentin-VIM) and mesoderm markers (BMP4 and PECAM1) was observed for all lines. However, we did not observe the relative expression of TUBB3 and FOXA2, both ectoderm and endoderm markers, respectively ( Figure 3B).
Interestingly, we observed the relative gene expression of AFP, an endoderm marker, in EBs derived from biPSCs culture in bFL2i in 20% O 2 conditions only. Oxygen tension did not seem to affect cell differentiation in our experimental conditions; nevertheless, we observed a higher expression of PECAM1 and VIM on EBs produced from the bFL2i high oxygen tension group ( Figure 3B).
Additionally, EBs were transferred to gelatin-coated dishes and cultured for another 7 days with IMDM supplemented with FBS. During the culture, it was possible to identify several different types of cells presenting varying morphologies ( Figure 3C).

Discussion
The generation of bovine in vitro reprogrammed cells has been already reported [17,18,25]. Indeed, it may contribute to the generation of genetically modified herds for optimized milk and meat production, as well as to the production of therapeutic proteins, the generation of animals resistant to several syndromes or diseases and even to biomedical models for basic and applied medicine due to cattle's more similar characteristics to humans in embryology, size and longevity when compared to rodents [26,27].
Herein, we describe the generation of bovine iPSCs generated in 20% and 5% O 2 . We performed a factorial experimental design of different oxygen tensions (5 and 20%) with different supplements (bFGF and bFL2i) to maintain in vitro bovine reprogrammed cells. In our conditions, we observed the formation of biPSCs colonies in both hypoxia and normoxia conditions after 14 days of transduction, with a higher percentage of colonies formed per transduced cell when using bFL2i supplementation independently of oxygen tension; however, further investigation and replicates must be carried out to confirm such initial observations.
In an elegant review, Nichols and Smith described two different states of pluripotency, naïve and primed; each state was described as having typical colony morphology. Pluripotent stem cells presenting the naïve profile grow as dome-shaped and compact colonies with well-defined borders, whereas primed profiles grow as flat colonies [5]. Indeed, pluripotency maintenance using different supplementation relies on several pluripotency pathways, leading to transcriptional and morphological changes during the reprogramming process. The fibroblast growth factors (FGFs) phosphorylate and activate FGF receptors (FGFRs), which leads to protein kinase C (PKC) signaling, regulating cell self-renewal, metabolism, survival, proliferation and differentiation [27]. LIF inhibits differentiation activity in iPSCs and stimulates self-renewal, directing them to pluripotency by acting in three ways: STAT3, PI (3)-kinase and MAPK [28]. Moreover, 2i differentiation inhibitors (MEKi: PD0325901 and GSK3i: CHIR99021) inhibit both MAPK and GSK-3 signaling, regulating the maintenance of naïve pluripotency by inhibiting the fibroblast growth factor receptor (FGFR) and repressing DNA methyltransferase (Dnmt3) expression [29].
The first reports of bovine iPSCs (in 2011-2012) differed regarding colony morphology; when biPSCs showed dome-shaped [17,30] or flat colonies [25,31], authors did not classify those lines as belonging to naïve or primed profiles. More recently, Kawaguchi et al. [32] demonstrated the differences in the morphology of colonies with different supplementations; bFGF + LIF lineages presented a flattened shape, whereas lineages cultured with PD0325901 (MEK pathway inhibitor) changed to a dome-shaped morphology. Additionally, a higher efficiency in colony generation was reported by using a combination of bFGF + LIF + 2i + forskolin.
Furthermore, oxygen tension can modulate different pathways, including those related to the maintenance of pluripotency. During the maintenance of human embryonic stem cells in culture, Ezashi et al. [33] suggested that the hypoxia environment was able to preserve the pluripotency of these cells and diminish the differentiated cells in culture; initial studies in iPSCs described the improvement of reprogramming efficiency in hypoxia [22,34].
Conversely, in this study, hypoxia tension could not maintain the pluripotent state of biPSCs when they were supplemented with bFL2i and these cells differentiated before 10 passages.
De Castro et al. [35] demonstrated that hypoxic conditions did not improve the reprogramming efficiency in equine iPSCs, as has previously been described in humans and mice; however, they also described different profiles of eiPSCs generated in hypoxic compared to normoxia. Indeed, the reprogramming of iPSC cells in a low oxygen atmosphere has not yet been studied in many species, including cattle. Under our experimental conditions, significant differences in the biPSCs colonies morphology were not observed, including when they were derived in either low or high O 2 or when they were derived in different supplementations.
The relative gene expression of pluripotency targets was evaluated and, in our conditions, no differences were observed in initial passages (p5). We observed the endogenous expression of SOX2, OCT4 and STELLA (also known as Dppa3); for these genes, we observed relative expression decreasing in early to late passages, which could be influenced by a freeze-thaw cycle (not shown) or differentiation due to exogenous vector silencing at 20% O 2 .
Lengner et al. [36] described the inefficiency of different oxygen tensions for promoting significant changes in the expression of core pluripotency genes OCT4, SOX2 and NANOG in human ESCs. In this work, we did not observe SOX2 and OCT4 expression changes when cultured in hypoxia or normoxia conditions. Another possible explanation for the pluripotency core transcription inefficiency is related to oxygen levels-lower oxygen levels (2% or less of oxygen) could potentially promote pluripotency core transcription regulation.
Great advances have been made in understanding pluripotency acquisition during the reprogramming process, which is itself considered to be a highly heterogenous process. For example, the mesenchymal-to-epithelial transition has long been reported as an important step for pluripotency acquisition [37]. More recently, Guo et al. 2019 demonstrated two other important pathways important to determine the reprogrammed fate through mouse reprogramming. Using the single-cell resolution analysis, they showed that Klf4 transcription factor contributes to the non-reprogrammed Cd34+/Fxyd5+/Psca+ keratinocyte-like fate and that IFN-γ impedes the final transition to chimera-competent pluripotency along the reprogrammed cells [38]. Nonetheless, a detailed description of a time-lapse gene expression analysis during reprogramming and pluripotency maintenance (comparing passages) in other species is still lacking in the literature. Several speciesspecific differences could influence the response of oxygen tension and the reprogramming process. Under our conditions, we maintained pluripotent-like cells for at least 25 passages in both oxygen conditions.
In low O 2 conditions, GLUT1 and GLUT3 expression was upregulated in our biP-SCs. GLUT3 demonstrated a higher affinity for glucose than GLUT1 and its high turnover makes it an efficient transporter [39]. Additionally, GLUT3 is the functional transporter for maternal glucose in mouse embryos and GLUT3 expression is required for blastocyst formation [40]. Somewhat similar to our result, Harvey et al. [41] described the same upregulation in GLUT1 expression in bovine blastocysts cultured under 2% oxygen and blastocysts cultured under 7% and 20% oxygen. Additionally, in in vitro hypoxia conditions (2% oxygen), GLUT1 expression was similar to that of in vivo-derived embryos. Controversially, hypoxia conditions (5% oxygen) decreased GLUT3 expression in equine iPSCs and human embryonic stem cells [35,42]. These controversial results demonstrate the different effects of hypoxia on metabolism pathways between species.
Under our conditions, we hypothesize, based on literature data, the molecular pathways ( Figure 4) being modulated and discuss the possible action of CHIR99021 in its inhibition of the GSK3 pathway and rising beta-catenin nuclear concentration, therefore leading to OCT4, SOX2 and NANOG expression. We suggest that the resultant C-MYC expression increase, an essential gene for increasing proliferation and work in the reprogramming process, together with low oxygen tension, may inhibit the mTORC1 pathway, thus decreasing biogenesis and mitochondrial function and potentially explaining the increase in GLUT3 once glycolysis becomes the only option for cells with low oxygen tension. ferent effects of hypoxia on metabolism pathways between species.
Under our conditions, we hypothesize, based on literature data, the molecular pathways ( Figure 4) being modulated and discuss the possible action of CHIR99021 in its inhibition of the GSK3 pathway and rising beta-catenin nuclear concentration, therefore leading to OCT4, SOX2 and NANOG expression. We suggest that the resultant C-MYC expression increase, an essential gene for increasing proliferation and work in the reprogramming process, together with low oxygen tension, may inhibit the mTORC1 pathway, thus decreasing biogenesis and mitochondrial function and potentially explaining the increase in GLUT3 once glycolysis becomes the only option for cells with low oxygen tension. Figure 4. Graphical hypothesis of molecular pathways involved on the cell signaling process during the reprogramming process using different inhibitors of differentiation and low oxygen tension, based on our results, as well as literature data [28,[43][44][45][46][47].
It is important to stress that most of the iPSC lines established in non-rodent species, especially livestock species, are dependent on continuous transgene expression to maintain their pluripotency [20]. In bovine iPSCs, Kawaguchi et al. [32] described the same characteristics of naïve-type iPSCs under the control of transgenes; after Dox removal, biPSCs readily changed their morphology and no longer expressed both exogenous and endogenous pluripotent genes. More studies are necessary to optimize culture conditions for the establishment and maintenance of ESCs and iPSCs, while naïve, intermediate and primed states seem to vary among species, which may prove a challenge for solving these issues in livestock species.
In conclusion, one biPSCs line was positive for differentiation in all three embryonic tissues. This lineage was cultured at a high O 2 tension and with bFL2i. These undirected differentiation results may present an incomplete ability of our EBs to express ectoderm, mesoderm and endoderm markers after the reprogramming process in our conditions, which is probably due to an insufficient period of in vitro differentiation or transgene expression. Indeed, we underline that, however the hypoxia condition was described in mice cells as an enhancer to promote differentiation in the endoderm lineage [45], the same effects were not found in bovine EBs derived from biPSCs that were produced and cultured in low O 2 conditions.

Conclusions
The acquisition and maintenance of pluripotency in vitro offers an important tool for improving knowledge on early development. For the first time, our findings show biPSCs characteristics when they are generated under low oxygen tension. In our conditions, oxygen tension may lead to less pronounced differences in in vitro pluripotency maintenance in biPSCs when compared to previous reports on humans and mice. Herein it was evident that bFGF combined with LIF and 2i did not favor pluripotency maintenance at 5% O 2 . In addition, our findings may contribute for demonstrating the ways in which a combination of oxygen and different supplementation sets impact bovine iPSCs.

Informed Consent Statement: Not applicable.
Data Availability Statement: The authors confirm that all data and materials support the published claims and comply with field standards.