miRNA-mediated control of exogenous OCT4 during mesenchymal-epithelial transition increases measles vector reprogramming efficiency

OCT4 is a key mediator of induced pluripotent stem cell (iPSC) reprogramming, but the mechanistic insights into the role of exogenous OCT4 and timelines that initiate pluripotency remain to be resolved. Here, using measles reprogramming vectors, we present microRNA (miRNA) targeting of exogenous OCT4 to shut down its expression during the mesenchymal to the epithelial transition phase of reprogramming. We showed that exogenous OCT4 is required only for the initiation of reprogramming and is dispensable for the maturation stage. However, the continuous expression of SOX2, KLF4, and c-MYC is necessary for the maturation stage of the iPSC. Additionally, we demonstrate a novel application of miRNA targeting in a viral vector to contextually control the vector/transgene, ultimately leading to an improved reprogramming efficiency. This novel approach could be applied to other systems for improving the efficiency of vector-induced processes.


INTRODUCTION
Reprogramming is a multi-dynamic molecular process involving conversion of somatic cells to induced pluripotent stem cells (iPSCs) via overexpression of four reprogramming factors (RFs): Oct4, Sox2, Klf4, and c-Myc (OSKM)), which can be delivered using multiple vector systems. [1][2][3] Based on gene expression profiling, human somatic cell reprogramming is divided into two stages; initiation and maturation, linked by mesenchymal to epithelial transition (MET) phase, which is characterized by the activation of epithelial genes and repression of mesenchymal genes. [4][5][6] MET is orchestrated by exogenous RFs-OSKM. [4][5][6][7][8] Oct4, Sox2, and Klf4 are considered "pioneer factors" for inducing pluripotency. 8 Among the pioneer factors, OCT4 protein, encoded by the Pou5f1 gene, is one of the first to be identified as a master regulator of pluripotency 9 and is found in both Yamanaka and Thomson cocktails to convert skin cells to iPSCs. 1,10 While, recent studies demonstrate that a significant increase in exogenous Oct4 expression over moderate levels negatively affects reprogramming 11,12 and iPSC quality. 13 Others studies have shown that exogenous Oct4 expression can be excluded entirely, inducing low efficiency and slower kinetics, 13 or replaced with other factors in certain conditions. 14 Altogether, these indicate that the role of exogenous OCT4 in reprogramming remains unclear. To address the question if the initial overexpression of OCT4 or its continuous expression throughout the reprogramming is favorable or detrimental to the process, it is essential to develop a flexible vector system in which OCT4 expression can be fine-tuned in a timely manner.
MicroRNA (miRNA) are 20-25-nts-long, small, non-coding RNAs that regulate gene expression by binding to a specific "seed sequence" on the target mRNA and either translationally repress or degrade it, controlling several cellular processes. [15][16][17][18] This regulation system, nicknamed miRNA targeting, is used to restrict (trans)gene expression to a particular cell or tissue type by incorporating miRNA target sequence (MTS) in the 3 0 untranslated region (UTR) of a specific gene. This process has been exploited in viral vectors to 1) control viral replication or propagation, 2) control virulence of oncolytic virus, 3) eliminate viral transgenes, or 4) to prevent it from regulating its natural targets. [19][20][21][22][23][24][25][26] MeV is a negative sense, non-segmented, single-stranded RNA virus that belongs to the Paramyxoviridae family. 27 Recombinant MeVs have been used as an oncolytic in clinical settings, 28 but its application as a viral vector in reprogramming has only started recently. [29][30][31] MiRNA-controlled oncolytic MeVs have previously been developed to increase tumor specificity and protect the surrounding organs. Specific MTSs were introduced either in the 3 0 UTR of the F gene to abolish fusion and propagation 32,33 or in the 3 0 UTR of the P, N, or L gene to control replication. 33,34 In 2019, our lab developed a single-cycle MeV, by substituting the hemagglutinin gene with GFP or RFs. 29,30 The MeV-derived iPSC could re-differentiate into the three germ lineages, and vector-free iPSCs were achieved in 3-5 clonal passages, indicating an excellent future for the single-cycle MeV reprogramming platform. 29,30 However, the first generation of MeV reprogrammed with low efficiency and additional effort in vector design was required to bring MeV reprogramming to a competitive level.
Here, we present a more efficient reprogramming MeV vector and a novel way to use miRNA targeting. In this study, we take advantage of the modular gene arrangement in MeV to separate the OCT4 gene from the SOX2-KLF4 bicistron to increase reprogramming efficiency. We then selected a miRNA, which is specifically upregulated during MET, 35-37 miR-375, and showed that we can use it to control the expression of OCT4 by MET. This direct control of OCT4 increases both the kinetics and efficiency of reprogramming. Through indirect silencing of OCT4, using P gene silencing, we demonstrate that OCT4 is dispensable for the maturation phase of reprogramming, while all three other factors, SOX2, KLF4, and c-MYC, are not. Finally, our study demonstrates a novel application for miRNA targeting in a viral vector in an autoregulatory loop to improve the overall efficiency of the vector-induced process.

Production of a MeV vector expressing OCT4 under the control of MTS
We have previously produced a single-cycle MeV expressing the four RFs, MV4F N , to produce iPSCs. 29 MV4F N expresses tricistron OSK instead of MeV hemagglutinin (H) and c-MYC in an additional transcription unit (ATU) after the H gene ( Figure 1A, top genome). We modified this vector by isolating the OCT4 from a SOX2 and KLF4 bicistron (SK) in an individual ATU to make the MV(O)(SK)(M) vector ( Figure 1A, middle genome). In MV(O) MTS375 (SK) (M), we introduced three copies of an MTS for miR-375 (MTS375) in the 3 0 UTR of the OCT4 gene ( Figure 1A, bottom genome). For both, MV(O) (SK)(M) and MV(O) MTS375 (SK)(M) vectors, GFP was inserted in the ATU after the P gene to track the vector. To confirm the expression and correct processing of the RFs from the three individual ATUs, Western blot analysis of neonatal human fibroblasts (NHFs) infected with both vectors was performed ( Figure 1B). Expression of OCT4, SOX2, and c-MYC with appropriate molecular weight was observed, indicating that the addition of the MTS375 in the 3 0 UTR region of OCT4 did not affect its expression. Confocal analysis confirmed the nuclear localization of the RFs ( Figure 1C). To address the possible effect of an additional ATU on the MeV vector propagation, a one-step growth curve was performed. All vectors showed comparable growth kinetics with the replication-competent MV(GFP) virus and replicated to maximum titers of greater than 10 6 by 48 h, confirming that vector propagation remains unaltered by the insertion of an additional ATU and the MTS375 in the 3 0 UTR of OCT4 ( Figure 1D).
Subsequently, we verified that the MTS375 control of OCT4 by miR-375 was functional. Cells (293T) were initially transfected with either miR-375, scramble miR precursors, or water, followed by transduc-  (Figure 2). From day 0 to day 6, Vimentin continued to be assembled into a network of filaments in the cytoplasm in both NHF and AHF (Figure 2A, second row). While MeVs can be detected by day 3 (GFP+ cells), it is only by day 9 that the GFP+ cells gained more epithelial morphology, with downregulation of Vimentin. The loss of Vimentin was synchronized with a weak expression of Occludin, indicating the initial formation of tight junctions. By days 12-15, the loss of Vimentin was evident, and the expression of Occludin stronger, suggesting loss of mesenchymal characteristics and gain of epithelial junctions. Subcellular localization of b-catenin is a marker of MET. By day 9, as GFP+ cells started losing their long protrusions and decreasing in size, a transition of b-catenin from the cytoplasm toward the nucleus was observed (Figure 2A, fourth row). By day 12, some GFP+ cells displayed both nuclear and cellular membrane expression of b-catenin. b-catenin's cortical distribution was reported to help in the formation of adherens junctions in emerging iPSC colonies. 38 The cortical expression of b-catenin and E-cadherin was observed by day 15 in most of the GFP+ cells, confirming completion of MET. Finally, expression of pluripotency markers NANOG and TRA-1-60 were observed starting at day 12 and onwards for both NHFs and AHFs (Figure 2A, eighth and nine rows). Confocal approach was validated by showing establishment of MET between days 6 and 9 in NHF lentivirus (LV)-mediated reprogramming (Figure S1), as previous reported. [39][40][41] The establishment of MET in MeV-mediated reprogramming was also verified using RT-qPCR in NHFs and AHFs ( Figures 2B and 2C). E-cadherin, Occludin, and Error bars indicate mean ± SD. A two-way ANOVA was used followed by Tukey's multiple comparison test (ns, not significant) (E and F) qPCR analysis of OCT4 (E) and KLF4 (F) expression from indicated vectors in cells transfected with miR-375, scrambled miRNA, or water. Data were normalized to GAPDH and represent the average ±SD of the mean from three independent replicates. A two-way ANOVA was used followed by Tukey's multiple comparison test (ns, not significant, ***p % 0.001 and ****p % 0.0001). (G and H) Immunofluorescence analysis of OCT4 (G) and KLF4 (H) expression from indicated vectors in 293T cells transfected with miR-375, scrambled miRNA, or water. Scale bars represent 100 mm.
EPCAM's relative expression increased between days 9 and 12 and onwards, supporting that MET was initiated by day 9 (Figures 2B,  2C, and 3G). Expression of pluripotency markers, NANOG or DPPA2, started by days 12-15 ( Figures 2B and 2C). On the other hand, there was no significant decrease in the expression of Vimentin ( Figure S2), and this is mostly due to the presence of a high number of untransduced fibroblasts. 6 Taken together, these results indicate that MET occurs between days 9 and 15 and is also cell maturity independent, as it occurs at a similar time frame in both NHFs and AHFs.
Previous reports demonstrated that miR-375 enhances MET in various cancer models and b cell-derived cells. [35][36][37] However, there is no existing report on the expression or involvement of miR-375 during reprogramming, making it an ideal candidate to regulate exogenous OCT4 without affecting reprogramming process. We determined the kinetics of miR-375 expression during MeV reprogramming of NHFs and AHFs using RT-qPCR analysis. While more modest in NHFs, a significant upregulation of miR-375 was observed in both NHFs and AHFs starting between days 9 and 12 and onwards. In both cells, the upregulation was directly correlated with the timing of the MET.  Figure 3B and not shown), indicating reprogramming was also increased by the control of OCT4 by miR-375.
We next evaluated the influence of the presence of the MTS375 in the 3 0 UTR of OCT4 on the MET timelines ( Figures 3C and S3A). Starting at day 6, a weak signal for Occludin and a loss of Vimentin were observed, and this became more significant as time progresses by days 9 and 12 for AHF reprogramming ( Figure 3C, second and third rows, and Figure S3D). Similar results were observed on NHFs (Figures S3A and S3C). While the subcellular localization of b-catenin followed the same pattern as previously described for MV(O)(SK)(M) reprogramming; once again, these processes started at day 6 instead of day 9 ( Figures 3C and S3A, fifth row). However, there was no change in the timing of the expression of E-cadherin between both vectors at day 12, and the subsequent increase in its expression is comparable ( Figure 3C We next determined the status of the exogenous and endogenous OCT4 during reprogramming. Specific primers for either the codon-optimized exogenous or endogenous OCT4 were used to differentiate between both transcript populations. While the expression kinetics pattern for exogenous KLF4 mRNA followed the MeV nucleoprotein (N) mRNA, there was a gradual loss in exogenous OCT4 expression during MV(O) MTS375 (SK)(M) reprogramming ( Figure 3D, right half). This gradual loss was directly correlated with the increase in miR-375 ( Figure 3E). In contrast, during MV(O)(SK)(M) reprogramming, all three mRNA levels, for exogenous OCT4, KLF4, and N, followed a similar expression kinetics (Figure 3D, left half). Expression of N, OCT4, or KLF4 was undetectable in both iPSCs (data not shown). This is attributed to both the addition of the antiviral at D20 and to the known elimination of our single cycle measles vector from iPSC reported in Wang et al. 29 Additionally, these vectors do not produce permanently infected iPSC, in contrast to the infectious measles virus that can lead to persistent infection. 42 Interestingly, the overall kinetics of the activation of endogenous and AHF (C) reprogramming. All values are relative to day 0 and normalized to GAPDH. Error bars represent mean ± SD with n = 3 independent experiments. A two-way ANOVA was used followed by Tukey's multiple comparison test (*p % 0.05, **p % 0.01, ***p % 0.001, ****p % 0.0001). (D and E) qPCR analysis of miR-375 during NHF (D) and AHF (E) reprogramming. Error bars represent mean ± SD. Data represent values from three independent experiments. All values are relative to day 0 and normalized to RNU6B. A one-way ANOVA was used followed by Tukey's multiple comparison test (ns, not significant, **p % 0.01, ***p % 0.001, **** p % 0.0001). OCT4 remained unaffected by the early elimination of exogenous OCT4 throughout reprogramming ( Figure 3F).
Recent studies showed that Claudin 7 (CLDN7)-mediated modulation of EPCAM, a transmembrane glycoprotein, results in the nuclear translocation of its intracellular domains (EpICD) by forming a complex with b-catenin and FHL2 to regulate gene expression of endogenous OCT4 and c-MYC. [43][44][45][46][47] Keeping this in mind, we next investigated the expression patterns of EPCAM and CLDN7 using qPCR. While EPCAM was significantly upregulated during early stages in MV(O) MTS375 (SK)(M) reprogramming compared with MV(O)(SK)(M) ( Figure 3G, days 3 and 6), there was no increase in CLDN7 ( Figure 3H).
The capability of iPSC-like clones to self-renew, proliferate, and differentiate was tested using confocal microscopy and RT-PCR analyses on isolated individual clones. All pluripotency markers were expressed in all clones tested, and there was no difference in the expression between the clones issued from both vectors ( Figures S4A and  S4B). The multi-lineages propensity of all clones was confirmed by the formation of the embryoid body and spontaneous differentiation into mesoderm (CD31) endoderm (FOXA2) or ectoderm (b-III tubulin) ( Figure S4C, left, middle, and right panel, respectively). The embryonic stem-cell-specific (ESCC) miRNAs are a group of miRNAs highly expressed in embryonic stem cells and promote cell cycle progression. 48  We next explored the effect of eliminating SOX2, KLF4, and c-MYC by MET during MeV reprogramming. Two additional vectors MV4F P and MV4F PMTS375 were produced. MV4F P expresses triciston OSK instead of H and c-MYC in the ATU after the tricistron (Figure 5A, top genome). MV4F PMTS375 was modified by inserting three repeats of the MTS375 in the 3 0 UTR of P gene ( Figure 5A, bottom genome). P protein of MeV serves as a cofactor to the virus polymerase L and plays a major role in replicating and transcribing the vector. 34,[49][50][51] Hence, targeting the MeV P using miRNA will control viral replication and transcription and, with it, the expression of all four RFs.
The expression of the RFs, OCT4, SOX2, and c-MYC, as well as the expression of the viral proteins, were not affected by the presence of the MTS375 in the 3 0 UTR of the P gene ( Figure 5B), indicating that the presence of the MTS375 did not affect the cofactor function of the P. Confocal microscopy confirmed the nuclear localization of the RFs in human fibroblasts transduced with both vectors (Figure 5C). To determine the effect of MTS375 insertion on vector propagation, we performed a one-step growth curve analysis on the MV4F P and MV4F PMTS375 and compared it with a replication competent MV(GFP) virus. All vectors replicated to comparable titers, indicating that the insertion of the MTS375 in the 3 0 UTR of the P gene does not affect vector replication and propagation ( Figure 5D). We next determined the functionality of the MTS375 insertion on the 3 0 UTR of P gene ( Figure 5E). Cells (293-H) expressing the MeV H gene were transfected with miR-375 or scrambled miRNA precursors or water and later transduced with either MV4F P or MV4F PMTS375 vectors. Propagation of MV4F PMTS375 was significantly decreased in the presence of miR-375, as shown by the reduction in titers and syncytia formation, compared with cells treated with control miRNA or water ( Figures 5E and S5). Unlike MV4F PMTS375 , MV4F P replicated to high titers in all conditions ( Figures 5E and S5), showing that the targeting of P by miR-375 is functional.
We next analyzed the ability of MV4F P and MV4F PMTS375 vectors to reprogram NHFs. While the MV4F P vector reprogrammed NHFs, the efficiency of MV4F PMTS375 was drastically decreased ( Figure 5F). To determine if this effect was due to the early elimination of all exogenous RFs by MET, reprogramming of MV4F PMTS375 was performed in the presence of individual LV-expressing SOX2, KLF4, from three independent experiments. Comparison made using unpaired two-tailed T tests (*p % 0.05). (C) Immunofluorescent labeling of specified markers at days 0, 3, 6, 9, 12, and 15 of AHF reprogramming. Merge is presented as overlays of indicated antibody staining, GFP, and Dapi staining. GFP+ cells within areas of interest are highlighted with yellow dashed lines for better visualization. The scale bars represent 100 mm. (D) qPCR analysis of relative expression of exogenous OCT4 (pink), KLF4 (blue), and Measles N (green) during AHF reprogramming with indicated vectors. All values are relative to day 0 and normalized to GAPDH. Error bars represent mean ± SD with n = 3, independent experiments. A two-way ANOVA was used followed by Sidak post hoc multiple comparisons test to compare the relative expression levels of OCT4 and KLF4 (**p % 0.01, ***p % 0.001, ****p % 0.0001). (E) qRT-PCR analysis of endogenous miR-375 expression in MV(O) MTS375 (SK) (M) reprogramming. All values are relative to day 0 and normalized to RNU6B. Error bars represent mean ± SD. Data represent values from three independent experiments. A one-way ANOVA was used followed by Tukey's multiple comparison test (***p % 0.001, ****p % 0.0001). (F-H) qPCR analysis of indicated endogenous OCT4 (F), EPCAM (G), and CLDN7 (H) during AHF reprogramming. All values are relative to day 0 and normalized to GAPDH. A two-way ANOVA was used followed by Tukey's multiple comparison test. Error bars represent mean ± SD, n = 3, independent experiments (*p % 0.05, **p % 0.01, ****p % 0.0001). and c-MYC. This LV co-transduction ensured the expression of SKM after MET, while OCT4 expression was suppressed through MTS375 targeting. Reprogramming with MV4F PMTS375 was rescued by adding the three LVs and reached similar efficiencies than MV4F P . However, transduction with a combination of one or two LVs could not rescue MV4F PMTS375 reprogramming ( Figure 5F), indicating that exogenous OCT4 can be shut down by MET but not SOX2, KLF4, and c-MYC ( Figure 5F).
The self-renewal, proliferation, and differentiation ability of the iPSClike clones derived from MV4F P and MV4F PMTS375 + LV were performed using immunofluorescence or RT-PCR on isolated individual clones. All clones tested expressed the pluripotency markers, and there was no difference in the expression between the clones issued from both vectors (Figures S6A and S6B). Multi-lineage propensity was established by the formation of embryoid body and spontaneously differentiating them into mesoderm, endoderm, or ectoderm ( Figure S6C), indicating that all clones produce are iPSCs.

DISCUSSION
This study demonstrates an autoregulatory system that amplifies the efficiency of the process initiated by the vector. We validate a new application for miRNA targeting in a viral vector that activates a vector-induced miRNA, which in turn controls the vector itself or one of the transgenes in an autoregulatory loop, resulting in increased efficiency of the process started initially by the vector. In this case, the MeV vector expressing the four RFs starts the reprogramming process and induces miR-375. Mir-375 can then bind to the MTS375 in the 3 0 UTR of the OCT4 gene expressed from the vector, causing its elimination during the MET phase, increasing reprogramming efficiency. In this work, we further demonstrate that only OCT4, and not SKM, is dispensable for the maturation stage of reprogramming.
Previous studies using partially and intermediately reprogrammed iPSC have facilitated context-dependent studies of reprogramming. [52][53][54] However, their production is sporadic and unpredictable, since reprogramming is transient and asynchronous. 55,56 Unlike the  previous methods, the miRNA controllable MeV vector system allows for context-dependent expression of RFs at a single cell level, catering to the asynchronous nature of reprogramming. Additionally, as MeV vector replicates in the cytoplasm, it is free of positional effects, unlike lentivectors. 57 Another hurdle in cell reprogramming is the low efficiency and slow kinetics of iPSC production. Many modulators (miRNA mimics and small molecules) have been identified to enhance iPSC generation. 58 In comparison with other strategies to enhance efficiency, the main advantage of this technology is simplicity. We showed that the expression of OCT4 can be fine-tuned in MeV by exploiting endogenous miRNA to ensure a precise and optimized level of transgene to improve efficiency. Hence, in this study, the new and improved MeV vector can reprogram AHFs with an average efficiency of 0.16% compared with 0.046%. Further, our novel approach could be incorporated into existing reprogramming systems.
During the initiation phase, the RFs, OCT4, SOX2, and KLF4, act as "pioneer factors" and bind to the inaccessible chromatin regions, leading to subsequent activation or repression of genes. 8 However, their role during the maturation phase is not well understood. We demonstrate that with the exception of OCT4, the other RFs (SOX2, KLF4, and c-MYC) are important for iPSC production, even in the maturation phase. This is in accordance with the previous reports that show the dependency of reprogrammed cells on transgene expression during early maturation for survival and growth. 59 However, silencing of the RFs during later stages of reprogramming is crucial to establish pluripotency. 59,60 In MeV reprogramming systems, the RF expression remains high, even during the maturation phase. It is most likely that this high expression of the RFs, especially OCT4, is detrimental for the transduced cell to complete reprogramming. By silencing the exogenous OCT4 promptly, by MET, we have allowed more cells to complete the process and adopt iPSC phenotype successfully.
We attribute the increase in reprogramming efficiency and the kinetics of MV(O) MTS375 (SK)(M) to an early expression of EPCAM and ESCC miRNA-miR-302a, miR-372, and miR-373. Indeed, an increase or decrease in EPCAM and ESCC miRNA expression has been shown to improve or decrease the reprogramming efficiency individually. 43,45,48,61,62 We believe that early EPCAM expression could be contributing to the faster kinetic of reprogramming by recruiting CLDN7 to the cell-cell junctions. 63 Further, CLDN7-mediated modulation of EpCAM could result in EpICD (intracellular domain) release, which complexes in the nucleus with proteins of the Wnt signal pathway (b-catenin, Lef-1), initiating gene transcription, leading to an early MET. Additionally, ESCC miRNA, by either targeting nuclear receptor subfamily 2 through indirect positive regulation of OCT4 or inhibiting the TGF-b signaling pathway, synergistically enhance MET and reprogramming of MV(O) MTS375 (SK) (M). 61 Our work demonstrates a novel application of miRNA targeting and that shutting down the expression of exogenous OCT4 by MET increases the reprogramming efficiency and kinetics. Studying the ef-fects of the temporal expression of the other exogenous RFs on the reprogramming process using the MeV vector systems and miRNAs upregulated during reprogramming could yield insights that can further enhance the efficiency and/or quality of reprogramming.

MATERIALS AND METHODS
De-identified human cells in this study were either obtained from the American Type Culture Collection (ATCC) or approved through Mayo Institutional Review Board. The viral vectors, viruses, and experiments associated with them were permitted by the Mayo Clinic Institutional Biosafety Committee.

Full-length measles virus cDNA plasmid production
Full-length cDNA vector p(+) MVvac2 DH(O)(SK)(cM)H was produced by splitting the one ATU containing the OSK polycistron in two ATUs containing the (O) and (SK) bicistron. Addition of the new ATU was performed using PCR and addition of the unique restriction site between the two ATUs to facilitate the cloning. All cloning steps were performed in accord with the "rule of six," using the intermediate vector pCG containing a PacI-SpeI fragment from full-length p(+) MVvac2DH(OSK) (cM)H30 (MV, Figure 1A). Then, the PacI-SpeI fragment was cloned back into the MV fulllength genome containing a GFP in an ATU between the P and M gene. The resulting full-length vector was called p(+) MVvac2

Viral vector production
Rescue of recombinant MeV vectors from transfected plasmids was carried out, as previously described, using Vero-H2 and helper 293-3-46-H2 cells. 30 In brief, the MV genome and MV polymerase were cotransfected into the helper 293-3-46 H2 cells and were transferred onto Vero-H2 cells 3 days later. The rescue was monitored for appearance and spread of GFP, and the virus was further expanded on Vero-H2 cells. All vectors were propagated on Vero-H2 cells, and stocks from the second or third passage were used for all experiments. Virus titers were determined by titration on Vero-H2 cells, by 50% endpoint dilution (TCID50) with individual infection events counted by GFP fluorescence at 72 h post infection using the Spearman-Kärber method. 64 Reprogramming of human fibroblasts cells NHF (7 Â 10 4 ) and AHF cells (5 Â 10 4 ) were seeded on Matrigel (Corning, Corning, NY, USA)-coated 12-well plates. Cells were either transduced with MeV vectors alone (MOI of 0.5) or in combination with LV vectors (at best optimal volume). 30 The virus and cells were spinoculated at 1,100 rpm for one hour at 25 C, after which the inoculum was left O/N at 37 C. The next day, infected cells were washed, and media 1 was added. Subsequently, the media were changed every other day until day 8, after which it was switched to media 2, with daily replacement until iPSC-like clones appear around days 20-25. For AHF, the cells were split onto two 6-well Matrigel-coated plates in media 1. On day 7, and for the next 7 days, cells were switched to media 2 containing small molecules (sm); SB431542 (5mM), PD0325901 (0.2mM), and Thiazovivin (0.5mM) (All Stemgent, Cambridge, MA, USA). After day 14, daily media change was performed with media 2 without sm, until iPSC-like clones were ready to be picked or fixed. At day 20, anti-MeV inhibitor, AS-136A (Sigma Aldrich, St. Louis, MO, USA), was added to eliminate the viral vector.
Reprogramming efficiency was calculated as the percent of TRA-1-60-positive iPSC colonies generated divided by the number of input cells. TRA-1-60-positive colonies were visualized with NovaRED HRP substrate (Vector Laboratories, Burlingame, California, USA), according to the manufacture.

One-step growth curves
Vero-H2 cells (4 Â 10 5 ) were infected with MeV vectors or MeV control virus with an MOI of 0.05 in OptiMEM for two hours at 37 C. Following viral adsorption, the inoculum was aspirated, cells were washed, and DMEM-10 was added. Both cells and supernatant were collected (24,48, and 72 h post infection), and viral progeny was tittered, as described above.

Immunostaining and confocal microscopy
For immunostaining, cells on Thermo Fisher Scientific Lab-Tek chamber slides (Sigma-Aldrich, St. Louis, MO, USA) or a 12-well tissue culture plate (MatTek corporation, Ashland, MA, USA) were fixed with 4% PFA, permeabilized, and stained for appropriate primary antibodies in 5% FBS/PBS O/N at 4 C and corresponding secondary antibodies for 1 h at room temperature (RT). Once stained, they were mounted with Prolong Gold Antifade (Life Technologies, Carlsbad, CA, USA). Further analyses were performed using a Zeiss LSM 780 confocal microscope followed by image processing with Zen black software (Zeiss). Primary and secondary antibodies are listed in Table S1.

Western blot
BJs (2.1 Â 10 5 ) were transduced with MeV at MOI 0.5. Cells were processed according to previously described procedures, following 36 h. 30 Separation of protein samples was performed on SDS page gels (Bio-Rad Laboratories, Hercules, CA, USA), followed by a transfer to polyvinylidene difluoride membranes (Immobilon-P, Bio-Rad Laboratories, Hercules, CA, USA). After blocking the membranes, they were incubated with primary antibodies. Following washes, membranes were subjected to peroxidase-conjugated secondary antibodies for 2 h at room temperature. Three washes later with TBS-Tween 0.1%, membrane was subjected to ECL2 substrate (Thermo Pierce, Waltham, MA, USA). Table S2 lists the primary and secondary antibodies.

Cellular and viral gene transcription by qRT-PCR and RT-PCR
Cells were subjected to Trizol reagent (Life Technologies, Carlsbad, CA, USA) to extract total RNA. EcoDry TM Premix Oligo dT kits (Takara Bio, Shiga, Japan) were used for cDNA synthesis. The transcript level of the target gene was determined using TaqMan PreAmp Master Mix (Thermo Fisher Scientific, Waltham, MA, USA). The gene expression was normalized to GAPDH and was calculated using the delta-delta-Ct algorithm. Platinum Taq DNA polymerase (Life Technologies, Carlsbad, CA, USA) was used to amplify the cDNA, subsequently used for PCR. Primers and probes are listed in Table  S3. RT-PCR to verify pluripotency markers in iPSC clones was carried out using Platinum Taq DNA polymerase (Life Technologies, Carlsbad, CA, USA). Primers for pluripotency markers were described in previous literature. 30 miRNA targeting assay HEK293LV-H cells or HEK293T cells (2 Â 10 5 ) were seeded in 12well Matrigel plates. The next day, we transfected miRNA mimics (mirVana miRNA Mimics, Life Technologies, Carlsbad, CA, USA) at a final concentration of 40 nM using Lipofectamine 2000 transfection reagent (Life Technologies, Carlsbad, CA, USA). Cells were transduced with MeV vectors at an MOI of 0.05, 4 h post-transfection. Viral vectors inoculum was incubated with the cells for 2 h at 37C, later removed, followed by a wash, and media addition. To determine production of virus, transduced cells were scraped, and samples were tittered or qPCR analyzed, as described above. For confocal analysis, HEK293T cells were fixed and immunostained, as described above.

Spontaneous differentiation assay
Briefly, the iPSCs were detached using EZ-LiFT Stem Cell Passaging Reagent (Sigma-Aldrich, St. Louis, MO, USA). After which, the cells were cultured to form EBs in a non-adherent 6-well plate and subsequently differentiated on Matrigel-coated chamber slides (LAB-TECKR -II, Thermo Fisher Scientific, Waltham, MA, USA). Subsequently, the cells were subjected to immunostaining and confocal microscopy (as described above) to verify that the iPSCs can differentiate in endoderm (FOXA2), ectoderm (b-III tubulin), and mesoderm (CD-31) lineages. Primary and secondary antibodies are listed in Table S1.

Statistics
Data were processed in Microsoft Excel. GraphPad Prism 9 was used to graph as well as perform statistical analysis. Single comparisons were made using unpaired two-tailed T tests. Several comparisons were analyzed by one-way or two-way ANOVA, following which Sidak or Tukey's multiple comparison test was used. For reprogramming, growth curve, and transfection/transduction, all experiments are presented as the average of three independent experiments. For qRT-PCR, results are presented as technical triplicate of three independent experiments. Data are graphed as group mean ± SD. Statistical significance cut off at p % 0.05 and ns > 0.05 for all experiments.

ACKNOWLEDGMENTS
We thank Drs. Navaratnarajah and Beyene for reading the manuscript. We thank Dr. Yasuhiro Ikeda for helpful discussion. We thank Megan Rasmussen and Amber Lopez for technical assistance. This work was supported by Regenerative Medicine Minnesota (RMM 102516 002 and RMM 092319 DS 005 to P.D.), Mayo Graduate School of Biomedical Sciences, and the National Institutes of Health (R21AI105233, R56HL147852-01, and R01HL147852-01A to P.D.). R.R. was supported by the Mayo Clinic Center for Regenerative Medicine and Michael S. and Mary Sue Shannon Foundation. The salary of C.B.D. was supported in part by grant T32 (AI132165) from the National Institutes of Health. The graphical abstract was created using Biorender.com.

DECLARATION OF INTERESTS
P.D., R.R., Q.W., and C.B.D. are inventors of an IP that was licensed by National Resilience Inc, and for which they may receive payments. The other authors declared no conflict of interest.