Hippo signaling pathway activation during SARS-CoV-2 infection contributes to host antiviral response

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), responsible for the Coronavirus Disease 2019 (COVID-19) pandemic, causes respiratory failure and damage to multiple organ systems. The emergence of viral variants poses a risk of vaccine failures and prolongation of the pandemic. However, our understanding of the molecular basis of SARS-CoV-2 infection and subsequent COVID-19 pathophysiology is limited. In this study, we have uncovered a critical role for the evolutionarily conserved Hippo signaling pathway in COVID-19 pathogenesis. Given the complexity of COVID-19-associated cell injury and immunopathogenesis processes, we investigated Hippo pathway dynamics in SARS-CoV-2 infection by utilizing COVID-19 lung samples and human cell models based on pluripotent stem cell-derived cardiomyocytes (PSC-CMs) and human primary lung air–liquid interface (ALI) cultures. SARS-CoV-2 infection caused activation of the Hippo signaling pathway in COVID-19 lung and in vitro cultures. Both parental and Delta variant of concern (VOC) strains induced Hippo pathway. The chemical inhibition and gene knockdown of upstream kinases MST1/2 and LATS1 resulted in significantly enhanced SARS-CoV-2 replication, indicating antiviral roles. Verteporfin, a pharmacological inhibitor of the Hippo pathway downstream transactivator, YAP, significantly reduced virus replication. These results delineate a direct antiviral role for Hippo signaling in SARS-CoV-2 infection and the potential for this pathway to be pharmacologically targeted to treat COVID-19.

To replicate, viruses hijack cellular machinery and signaling pathways. Thus, understanding the molecular basis of dysregulated signaling circuitry can provide potential therapeutic targets. Utilizing a library of kinase inhibitors, we previously demonstrated critical pathways that are important for SARS-CoV-2 virus replication such as the DNA-damage response (DDR), mTOR-PI3K-AKT, and ABL-BCR/MAPK pathways [19]. COVID-19 pathogenesis is a complex process involving crosstalk with multiple cellular pathways. To understand the fundamental principles underlying RNA viral mediated cell injury processes, we hypothesize that the conserved cellular circuitry of the Hippo signaling pathway, which is involved in tissue growth, immune response, and inflammation, can contribute to regulating SARS-CoV-2 replication and COVID-19 pathogenesis. Previously, we observed that the Zika virus (ZIKV), a member of the Flaviviridae family, activates this pathway [20]. The core components of the Hippo pathway in mammals are STE20-like kinases STK4/MST1 and STK3/MST2 (Hpo in Drosophila), large tumor suppressors (LATS1/2, Wts in Drosophila), MOB kinase activator 1A and 1B (Mats in Drosophila), Neurofibromatosis 2 (NF2), and Salvador 1 (SAV1) [21]. MST1/2 and LATS1/2 are serine/threonine kinases, whereas NF2 and SAV1 are adaptors or co-factors. Upon activation by stimuli, MST1/2 phosphorylate and activates the downstream serine/threonine kinase LATS1/2. LATS1/2 can be phosphorylated by MAP4Ks independent of MST1/2 axis. The phospho-LATS1/2, along with SAV1, then phosphorylates YAP (Yes-associated protein), a transcription co-activator and oncogene. Phospho-YAP, bound to a 14-3-3 protein and retained in the cytoplasm, is subsequently ubiquitinated and degraded by proteasomes. While the Hippo pathway is inactive, the unphosphorylated YAP and its homologous partner TAZ translocate to the nucleus, where they bind to the TEAD family of transcription factors and mediate expression of target genes involved in cell survival and proliferation (Birc5, Id, CTCF, Axin2, Myc, CycD, etc.) [22]. As the Hippo signaling pathway controls progenitor cell proliferation and differentiation, and the size of organs [23][24][25][26][27], the pathway is closely linked to mitochondrial energy metabolism and biogenesis. Studies suggest that the Hippo pathway modulates host antiviral immune responses [28][29][30][31][32][33], where YAP inhibits antiviral defense mechanisms by antagonizing the function of pro-innate immune factor TBK1 [34], which is a key signal transducer of cytosolic RNA-sensing retinoic acid-inducible gene I (RIG-I)-like receptor (RLR) pathway. Inherited autosomal recessive mutations of MST1 (STK4) in human have been attributed to primary immunodeficiency with T-and B cell lymphopenia, neutropenia, and defective regulatory T cells [30,[35][36][37]. Thus, we postulate that Hippo signaling pathway plays a critical role in SARS-CoV-2 antiviral response and COVID-19 pathogenesis.

Results/Discussion
To understand the pathophysiological effect of SARS-CoV-2 on the Hippo signaling pathway, first we performed transcriptomic analysis of previously published RNA sequencing (RNAseq) datasets of lung samples from COVID-19 patients (n = 5) [38]. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicates that many genes (45 genes) in the Hippo pathway were differentially regulated (p < 0.01 and log2FC>1) (Fig 1A and S1 Data), and 33 (73%) of these genes are up-regulated, likely due to an antiviral host response. WWTR1/TAZ, a YAP homologue and downstream transcriptional coactivator of the Hippo signaling pathway, was found to be significantly down-regulated, suggesting the activation of Hippo signaling pathway. Phosphorylation of YAP (phospho-YAP) at Serine (Ser or S) 127 is a hallmark of Hippo pathway activation and results in cytoplasmic retention and inactivation of its transactivation function. We investigated YAP/TAZ phosphorylation status by immunohistochemical staining and analysis of 4 autopsy lung tissue samples from COVID-19 patients and control lungs ( Fig 1B). Three of the samples were positive for the phospho-YAP (Ser127) phenotype. Using RNA-fluorescence in situ hybridization (RNA-FISH), SARS-CoV-2 viral genomes were detected in the COVID-19 lung tissue ( Fig 1B). In addition, immunohistochemistry demonstrated a higher level of phospho-YAP (Ser127) in the COVID-19 lung compared to control lung ( Fig 1B). Furthermore, we observed that the COVID-19 lung had heavy infiltration of CD68-positive (+) mononuclear inflammatory cells (S1A Fig). CD68 marks macrophages in the inflammatory milieu. The influx of CD68+ cells in the COVID-19 lung could be due to tissue damage caused by virus replication. The Hippo pathway activation in COVID-19 lung can induce Type I IFN signaling pathway that can contribute to host antiviral inflammatory response.
The host response to virus infection is a dynamic process. The temporal differences in the time of virus infection and the progression to severe disease in COVID-19 patients can contribute to the expression and phosphorylation levels of YAP. Thus, to investigate SARS-CoV-2 viral replication and the Hippo signaling pathway in a well-controlled experimental system, we utilized a biologically relevant primary human proximal airway cell culture model, which consisted of mucociliary epithelial cells grown at an air-liquid interface (ALI). Human lung airway basal stem cells (ABSC) harvested from the trachea and bronchi of healthy lungs were differentiated into a ciliated pseudostratified columnar epithelium in a collagen-coated transwell membrane culture system [19,39,40]. The ciliated and secretory airway epithelial cell types in the ALI culture were verified using immunostaining for the cell type-specific markers, acetylated tubulin and Mucin5AC, respectively ( Fig 1C). The apical surface of the ALI cultures was exposed to SARS-CoV-2 infection to understand how the virus affects the Hippo pathway at a post-translational level. Parental SARS-CoV-2 (Isolate USA-WA1/2020 at MOI of 1) from BEI Resources was used for infection studies in the UCLA BSL3 high-containment facility. In all in vitro experiments, we used 48 hours post-infection as a common endpoint. Since primary ALI is presented. Bar graph shows Log2FC (FC >1) of 45 Hippo signaling pathway genes that were differentially expressed in COVID-19 lung samples (p < 0.01). The adjacent heatmap depicts Z scores as expression levels of these 45 DEGs. Blue and red colors represent 12 down-regulated and 33 up-regulated genes, respectively. The KEGG pathway database was used to examine the Hippo signaling pathway genes in the differentially expressed gene list. Individual quantitative observations that underlie the data summarized can be located under the Supporting Information File as S1 Data. (B) IHC of COVID-19 lung autopsy tissue shows high level of phospho-YAP Ser127 protein (dark brown). SARS-CoV-2 RNA (red) presence in COVID-19 lung was confirmed by RNA-FISH (RNAscope) (scale bar: 100 μm). (C) Western blot analyses of lung ALI cells show activation of Hippo-STAT1 pathways during SARS-CoV-2 infection (MOI 1). Density of bands compared to that of 2 dpi mock is provided for pYAP127 and YAP. Proximal lung air interface cultures are susceptible to SARS-CoV-2 infection. Antibody targeting Spike protein (green) was used for probing infection by immunohistochemistry. Cell-specific markers, such as Ac-tubulin (ciliated cells) and Mucin5 (mucus cells) were detected by antibody probes (red) to define infected cell types (6 dpi). Scale bar: 10 μm. (D) Human airway epithelial cells (Calu-3) were analyzed with IHC for pYAP127 and SARS-CoV-2 Spike protein (red) in uninfected (Mock) is less susceptible to parental SARS-CoV-2 infection [39], we followed the viral infection for up to 6 days. Western blot analysis of protein samples collected at 2 and 6 days post-infection (dpi), showed an increased level of phosphorylated YAP (Ser127) suggesting activation of the Hippo signaling pathway during SARS-CoV-2 infection (Fig 1C). Interestingly, we have noticed an increase in total YAP1 protein level in the infected samples. In addition, we observed a concurrent activation of the innate immune pathway as evidenced by an increase in phosphorylated STAT1 (Tyr701). SARS-CoV-2 infection was confirmed in the lung ALI culture by western blot and immunohistochemistry that revealed both ciliated and mucus secreted cells were infected at 6 dpi ( Fig 1C). Next, we examined the Hippo signaling pathway in the Calu-3 human airway epithelial cell line during infection with SARS-CoV-2 parental as well as the Delta variant of concern (VOC) (Fig 1D). We observed that both the parental and Delta variant strains of SARS-CoV-2 increased phospho-YAP (Ser127) level in the infected Calu-3 cell cultures (MOI 0.1). Interestingly, the phospho-YAP (Ser127) protein aggregates formed punctate structures in the cytoplasm, likely for potential degradation in autophagosomes [41]. Taken together, the Hippo signaling pathway is activated in COVID-19 infected lungs and SARS-CoV-2 infected in vitro lung culture systems.
Although SARS-CoV-2 is a major respiratory pathogen, COVID-19 manifestations in the cardiovascular system are also well documented [17,18,[42][43][44]. To further understand how the Hippo signaling pathway may be modulated by SARS-CoV-2 and the host response, we studied this pathway in cardiomyocytes. Human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) have been shown to be efficient at recapitulating cardiovascular diseases at a cellular level [45,46], have demonstrated susceptibility to SARS-CoV-2 replication [47], and have been used as a platform to study potential therapeutics against SARS-CoV-2 infection [19]. Transcriptomic analysis of SARS-CoV-2-infected hPSC-CM system [48] showed that 25 genes in the Hippo pathway are significantly differentially expressed (p < 0.01), of which 19 (76%) of the genes are up-regulated relative to uninfected mock hPSC-CM control (S1B Fig and S4 Data). While comparing the Hippo pathway DEGs of COVID-19 lung samples (Fig 1A and S1 Data) and of hPSC-CMs (S1B Fig and S4 Data), we observed that genes such as BMPR1B and GDF6 are commonly up-regulated. Subsequently, we used these highly susceptible hPSC-CM to further investigate the kinetics of the Hippo signaling pathway during early and late phases of viral infection. SARS-CoV-2-infected cardiomyocytes were confirmed by specific cell-staining of cardiac troponin T (cTnT) and viral Spike proteins (Fig 2A). SARS-CoV-2 induced pronounced cytopathic effect (CPE), characterized by cell fusion and rounding of cells, and fragmentation of cTnT filaments, as well as apoptosis by activation of cleavage of caspase-3 (CC3) indicated the high susceptibility of hPSC-CM to cell death by infection. Western blot revealed that YAP1 is induced at 2 and 24 hpi but degraded by 48 hpi (Fig 2B). Simultaneously, there is an induction of phospho-YAP1 (S127) in the infected cells, which indicate a strong activation of the Hippo signaling pathway. Also, SARS-CoV-2 infection results in increased levels of phospho-TBK1 and phospho-IRF3 indicating activation of the RLR innate immune signaling pathway as early as 2 hpi. In addition, the antiviral STAT1 Type I interferon (IFN) pathway is induced in hPSC-CM following infection ( Fig 2B). Incoming viral RNA and envelope glycoprotein can be detected by TLR and RLR-IFN pathway sensors in minutes [49][50][51], which can allow the cells to mount effective antiviral defense. Furthermore, and infected cells at 48 hours post-infection. Parental and Delta strains were used for infection studies (MOI 0.1). Scale bar: 50 μm. Insets highlight phospho-YAP S127 in infected cells, and white asterisks denote uninfected cells. Representative data from 3 independent studies are provided. ALI, air-liquid interface; COVID- 19

PLOS BIOLOGY
immunohistochemical analysis showed that YAP/TAZ has been depleted in the infected cells with a cytoplasmic localization ( Fig 2C) and the level of pYAP127 has increased at 48 hpi ( Fig  2D). Phosphorylation of YAP at Serine 127 residue by LATS1/2 leads to cytoplasmic distribution of the YAP protein and prevents its transactivation function [52]. Collectively, our data concludes that the Hippo signaling cascade is also active during SARS-CoV-2 infection of cardiomyocyte systems; therefore, we next sought to determine the functional relevance of this pathway.

PLOS BIOLOGY
In order to uncover the functional significance of the Hippo pathway, we utilized both genetic and pharmacological ablation approaches. The core Hippo signaling pathway is a kinase cascade where MST1/2 kinases and Sav1 form a complex to phosphorylate and activate LATS1/2 kinases. After activation, LATS1/2 can phosphorylate and inhibit YAP and TAZ transcription co-activators. YAP/TAZ, when dephosphorylated, interacts with TEAD1-4 and additional transcription factors to stimulate gene expression for cell proliferation and to inhibit apoptosis. Therefore, the Hippo pathway is regulated at multiple levels. MST1/2 and LATS1/2 are regulated upstream by other molecules like Merlin and KIBRA, and YAP/TAZ is regulated downstream by protein ubiquitination. Thus, we postulated that inhibition of MST1/2 and LATS1/2 would lead to an increase in YAP/TAZ activity, and consequently an increase in SARS-CoV-2 replication. To test this mechanism, we performed shRNA-mediated gene silencing experiments in hPSC-CM. For knockdown experiments, we transduced the cells with lentiviral vectors expressing nonspecific control shRNA or shRNAs targeting YAP/ TAZ or LATS1. Approximately 72 hours post-transduction, the cells were infected with SARS-CoV-2 at MOI of 0.01. At 48 hpi, the cells were fixed for IHC or harvested in RIPA protein lysis buffer for western blot analysis. Immunohistochemical analysis revealed that the shRNA-mediated partial knockdown of YAP1 gene in hPSC-CM cells resulted in significantly reduced SARS-CoV-2 infection, whereas LATS1 knockdown, increased SARS-CoV-2 infection (Fig 3A and 3D). The quantification of spike antigen positive infected cells is provided in Fig 3D and S2 Data. These observations were also verified at the protein level ( Fig 3B). Our observation that the induction of YAP protein at early stages of SARS-CoV-2 infection ( Fig  2B) could be due to viral-mediated up-regulation of pro-viral YAP protein as a result of a complex viral-host interaction response and signaling feedback regulation. Taken together, our data provide evidence that the YAP/TAZ is a pro-viral factor, whereas LATS1 has antiviral function.
In addition, we investigated well-characterized direct acting chemical inhibitors on the Hippo signaling pathway. XMU-MP-1 blocks MST1/2 kinase activities and thereby enhances YAP/TAZ activity. MST1 and MST2 are core components of the Hippo pathway and critical targets for tissue repair and regeneration [53]. Therefore, XMU-MP-1, a selective and reversible inhibitor, increases tissue repair and regeneration properties following tissue injuries [53,54]. XMU-MP-1 increases cardiomyocyte survival and reduces apoptosis following oxidative stress [54]. In mice, XMU-MP-1 has been shown to increase intestinal repair, liver repair and regeneration [53], and preserve cardiac function and inhibit cardiomyocyte apoptosis in mice with transverse aortic constriction (TAC) [54]. Infection of cardiomyocytes and Calu-3 lung cells pre-treated with the MST1/2 inhibitor, XMU-MP-1 (10 μm), led to an increase in SARS-CoV-2 replication (Fig 3C-3E). XMU-MP-1 compound is non-toxic to the Calu-3 cells (S2A Fig). Interestingly, XMU-MP-1 treatment induced a coarse granular morphology to the nuclear localized YAP/TAZ (Fig 3E), and in infected cells the YAP/TAZ protein was depleted. Furthermore, XMU-MP1 treatment stimulated an increase in SARS-CoV-2 replication. Similar to the activity of XMU-MP-1, MST1/2 loss of function mutations in humans result in immune deficiency disorders and the affected individuals are more susceptible to human papilloma infections [33]. Consistent with the clinical observation, our study indicates that pharmacological inhibition of MST1/2 promotes the susceptibility of cardiac and lung cells to SARS-CoV-2 infection. Thus, MST1/2 or LATS1/2 loss of function mutations can potentially increase susceptibility to SARS-CoV-2 infections in humans.
Since YAP/TAZ is a proviral factor, pharmacological inhibition can provide additional therapeutic options for COVID-19 treatments. Thus, we focused on pharmacological modulation of YAP/TAZ. Verteporfin, a small molecule derivative of porphyrin, is a YAP-TEAD inhibitor. Verteporfin is a US Food and Drug Administration (FDA) approved drug used in patients with predominantly classic, subfoveal choroidal neovascularization (CNV) caused by age-related macular degeneration (AMD) (VAM) [55]. In addition, because YAP is a protumorigenic factor expressed in multiple cancers, Verteporfin has been used to inhibit YAP as an anti-cancer chemotherapeutic and adjuvant drug [56], which increased phosphorylation of YAP1 (S127) in Calu-3 cells (S2B Fig). Therefore, modulating YAP/TAZ with Verteporfin can provide more insight into YAP's proviral properties. Calu-3 cells were pretreated with 1 μm Verteporfin (non-toxic concentration; S2A Fig) for 24 hours and subsequently the cells were infected with SARS-CoV-2. At 48 hpi, the cells were harvested for protein analysis and the culture supernatants were collected for measuring viral titer. Verteporfin treatment resulted in a reduction of YAP/TAZ protein level compared to vehicle-treated cells (Fig 3F and S3 Data) as well as a decrease in SARS-CoV-2 replication at 48 hpi. Viral titer (TCID50/ml) measurement of Verteporfin-treated infected cell culture media revealed there was a significant decrease of viral production in Calu-3 cells (Fig 3F and S3 Data). Our observation supports the results of an independent study that used an FDA-approved drug library of 3,200 small molecules to discover anti-SARS-CoV-2 compounds and identified Verteporfin and verified its antiviral properties in a mouse model [57]. Together, our results indicate a direct antiviral role of Hippo signaling on SARS-CoV-2 infection and disease pathogenesis processes, which can be therapeutically targeted. A schematic diagram of our working hypothetical model is shown in Fig 3G. In conclusion, our results show that SARS-CoV-2 infection caused activation of the Hippo signaling pathway in hPSC-CMs, human primary lung ALI cultures, and human airway epithelial cells (Calu-3). These in vitro cultures are efficiently infected by SARS-CoV-2, consequently causing activation of immune, and inflammatory responses, and altered Hippo signaling cascade. shRNA-mediated partial knockdown of LATS1 and pharmacological modulation of core upstream MST1/2 kinases resulted in enhanced SARS-CoV-2 replication. Loss of function mutations in these antiviral kinases could enhance susceptibility to SARS-CoV-2 infections in human. While Verteporfin, a pharmacological inhibitor targeting downstream transactivator, YAP, significantly reduced viral replication and production. Moving forward, additional mechanistic investigations are required to elucidate the underpinnings of the Hippo pathway role in antiviral responses across pandemic potential RNA viruses and additional preclinical small and large animal safety and efficacy studies of Verteporfin therapy for COVID-19. Taken together, our results indicate a direct antiviral role for Hippo signaling in SARS--CoV-2 infection, thus creating potentially new avenues this pathway can be pharmacologically targeted to treat COVID-19.  Hippo on). c-GAS, cyclic GMP-AMP synthase; IKKe, inhibitor of nuclear factor kappa B kinase subunit epsilon; MAVS, mitochondrial antiviral-signaling protein; RIG-I, retinoic acid inducible gene I protein; STING, stimulator of interferon response cGAMP interactor 1; TEAD, TEA domain transcription factors. Individual quantitative observations that underlie the data summarized in Fig 3D and 3F

PLOS BIOLOGY
Hippo signaling antiviral response to SARS-CoV-2

Ethics statement
Our study does not involve human subjects, but used deidentified human tissues. For human tissue procurement, large airways and bronchial tissues were acquired after lung transplantations at the Ronald Reagan UCLA Medical Center from deidentified normal human donors following Institutional Review Board [58] exemption. Pluripotent stem cell-related work received UCLA ESCRO-IRB approval. All the SARS-CoV-2 live virus experiments were performed at the UCLA BSL3 High containment facility.

Cell lines, hPSC-CM cultures, and ALI cultures
Vero E6 cells were obtained from ATCC [VERO C1008 (ATCC CRL-1586)] or DSMZ (Braunschweig, Germany). Cells were cultured in EMEM growth media containing 10% fetal bovine serum (FBS) and penicillin (100 units/ml). Human lung adenocarcinoma epithelial cell line (Calu-3) was purchased from ATCC (ATCCHTB-55) and cultured in Dulbecco's Modified Eagles Medium (DMEM), supplemented with 20% FBS, 1% L-glutamine (L-glu), and 1% penicillin/streptomycin (P/S). Cells were incubated at 37˚C with 5% CO 2 . hPSC-CMs were derived from hESC line H9 using the previously described method [59]. Temporarily, hPSCs were maintained in mTeSR1 (STEMCELL Technology). Shortly after, RPMI1640 was supplemented with B27 minus insulin (Invitrogen) and used as differentiation medium. From Day 0 to 1, 6 μm CHIR99021 (Selleckchem) was added into differentiation medium. Next, Day 3 to 5, 5 μm IWR1 (Sigma-Aldrich) was added to differentiation medium. After Day 7, RPMI 1640 plus B27 maintenance medium was added. From Day 10 to 11, RPMI 1640 without D-glucose supplemented with B27 was transiently used for metabolic purification of CMs [60]. After selection, hPSC-CMs were replated for viral infection. ALI cultures derived from primary human proximal ABSCs were used as described previously [40]. The 24-well 6.5-mm transwells with 0.4-mm pore polyester membrane inserts were used for culturing ALI cells. A total of 500 μl ALI media (PneumaCult-ALI Medium, STEMCELL Technologies) was used in the basal chamber for ALI cultures, and cells were cultured at 37˚C with 5% CO 2 . Use of these established cell lines for the study was approved by the Institutional Biosafety Committee at UCLA.

Virus
SARS-Related Coronavirus 2 (SARS-CoV-2), Isolate USA-WA1/2020, and Isolate hCoV-19/ USA/MD-HP05647/2021 (Lineage B.1.617.2; Delta variant) were obtained from BEI Resources of National Institute of Allergy and Infectious Diseases (NIAID). All the studies involving live virus was conducted in UCLA BSL3 high-containment facility. SARS-CoV-2 was passaged once in Vero E6 cells and viral stocks were aliquoted and stored at −80˚C. Virus titer was measured in Vero E6 cells by established plaque assay or TCID50 assay.

SARS-CoV-2 infection
Calu-3 cells were seeded at 30 × 10 3 cells per well in 0.2 ml volumes using a 96-well plate and hPSC-CMs were replated at 1 × 10 5 cells per well in a 48-well plate. Viral inoculum (MOI of 0.01 or 0.1; 100 μl/well) was added using serum-free base media. The conditioned media from each well and condition was removed and 100 μl of prepared inoculum was added onto cells. After 1 hour incubation at 37˚C with 5% CO2, inoculum was replaced for Calu-3 cells with serum-supplemented media (200 μl per well) and for hPSC-CM, cell culture medium was replaced with RPMI 1640 + B27 supplement with insulin. In all our experiments, we included 48 hpi as a standard endpoint unless otherwise stated. For ALI cultures, 100 μl of viral inoculum (MOI of 1) prepared in PneumaCult media was added to the apical chamber of ALI culture insert and incubated for 1 hour at 37˚C with 5% CO2. For mock infection, PneumaCult media (100 μl/well) alone was added. The inoculum was dispersed throughout the hour every 15 minutes by gently tilting the plate sideways. At the end of incubation, the inoculum was removed from the apical side for ALI cultures to maintain the ALI [19,40]. At selected time points, cells were fixed with 4% PFA, collected by 1xRIPA for protein analysis, and/or supernatant collected for viral titer. Viral infection was examined by immunostaining or western blot analysis using SARS-CoV-2 antibodies. BEI Resources NR-616 monoclonal anti-SARS-CoV Spike protein (Similar to 240C) SARS coronavirus for immunostaining and GeneTex SARS--CoV-2 (COVID-19) Nucleocapsid antibody [HL344] (Cat#GTX635679) for western blot.

Viral titer by TCID50 (median tissue culture infectious dose) assay
Viral production by infected cells was measured by quantifying TCID50 as previously described [61] methodology for quantifying TCID50; 10 × 10 3 cells/well of Vero E6 cells were seeded in 96-well plates. The next day, culture media samples of supernatant collected from Calu-3 cells were at 48-hour time point was subjected to 10-fold serial dilutions (10 1 to 10 8 ) and inoculated onto Vero E6 cells. The cells were incubated at 37˚C with 5% CO2 for 3 to 4 days to evaluate for presence or absence of viral CPE. After measuring the percent infected dilutions immediately above and immediately below 50%, the TCID50 was calculated based on the method of Reed and Muench.

Histopathology
COVID-19 patient autopsy lung samples and normal lung samples (n = 4 each) were procured from the UCLA Translational Pathology Core Lab for research use. Samples were processed and sectioned, after a pathologist confirmed the section quality by H&E staining, and subsequent confirmation of COVID-19 positivity by RNAscope V-nCoV2019-S probe (ACD, Cat#: 848568, ready to use). Immunohistochemistry stainings were also performed on this lung tissue: Paraffin-embedded sections were cut at 4-μm thickness and paraffin was removed with xylene and the sections were rehydrated through graded ethanol. Endogenous peroxidase activity was blocked with 3% hydrogen peroxide in methanol for 10 minutes. Heat-induced antigen retrieval [9] was carried out for all sections in AR9 buffer (AR9001KT Akoya) using a Biocare decloaker at 95˚C for 25 minutes. The slides were then stained with YAP (S127) antibody (Cell Signaling, 13008, 1-100) and CD68 antibody (Dako, m0876, 1-200) at 4 degree overnight; the signal was detected using Bond Polymer Refine Detection Kit (Leica Microsystems, catalogue #DS9800) with a diaminobenzidine reaction to detect antibody labeling and hematoxylin counterstaining.

Immunohistochemistry
Cells were fixed with 4% PFA for 30 to 60 minutes or methanol (incubated in −20˚C freezer until washed with 1XPBS). Cells were washed 3 times with 1× PBS before permeabilizing with blocking buffer (0.3% Triton X-100, 2% BSA, 5% Goat Serum, 5% Donkey Serum in 1× PBS) for 1 hour at room temperature. After adding specific primary antibodies (refer to S1 Table), cells were incubated overnight at 4˚C. The next day, cells were washed with 1× PBS 3 times and incubated with respective secondary antibody (refer to S1 Table) at room temperature for 1 hour. DAPI (4 0 ,6-diamidino-2-phenylindole, dihydrochloride) (Life Technologies) was used to stain nuclei at a dilution of 1:5,000 in 1× PBS. Image acquisition was done using Leica DMi1 fluorescent microscopes and using the Leica Application Suite X (LAS X). The LSM 700 confocal microscopes and Zeiss Software programs available at the UCLA Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research Microscopy Core at Center for Health Sciences Building was used for capturing images as well.

Image analysis/quantification
The confocal images were obtained using the Leica Application Suite X (LAS X) and/or the Zeiss LSM 700 Confocal Microscopy by Zeiss Software Program with maximum intensity projection feature. Using a double-blinded approach to count the positively stained cells, Image J's plugin Cell Counter program was used. Approximately 3 to 4 independent images per condition of the hPSC-CM were analyzed. For all the analysis, we used confocal images acquired at 63× objective and fluorescent microscope images acquired at 20 and 40×. Total number of cells in each image was obtained by manually counting all the DAPI-stained nuclei using Ima-geJ program (Version 1.8.0; https://imagej.nih.gov/ij/index.html). Total cell count was then used for normalization and to calculate the percentage of individual marker positive cell populations in respective images. The mean percentage of positively stained cells from 3 to 8 independent images (about 150 to 1,000 cells total) were quantified and presented as graph format.

RNA sequencing data analysis
We used a previously published [38] RNA sequencing dataset of 5 normal and autopsy lung samples from COVID-19 patients were obtained from the NCBI Gene Expression Omnibus (GEO) under the accession number GSE150316 [38]. The SARS-CoV-2 infected pluripotent stem cell-derived cardiomyocytes gene expression data used in this study were also retrieved from GEO under the accession number GSE150392 [47]. The raw read counts per gene were used as inputs for differential expression gene analysis using DESeq2 v1.28.1 in R v4.0.3 [62]. Median of ratios method was used to normalize expression counts for each gene in all experimental samples. Each gene in the samples was fitted into a negative binomial generalized linear model. Genes were considered as differentially expressed (DEG) if they were supported by a false discovery rate (FDR) p < 0.01. The heatmaps were prepared using shinyheatmap web interface [63], which uses normalized read counts as inputs for Z-score normalization, which is computed on a gene-by-gene basis by subtracting the mean and then dividing by the standard deviation. The KEGG pathway database was used to examine the hippo signaling pathway in the DEG datasets [64].

Data analysis
All statistical testing was performed at the 2-sided alpha level of 0.05. To test statistical significance, unpaired Student t test was used to compare 2 groups (uninfected versus infected). GraphPad Prism software, version 8.1.2 (GraphPad Software, US) was used.
Supporting information S1 Fig. (A) Immunohistochemistry of COVID-19 lung autopsy tissue shows heavy infiltration of CD68-positive inflammatory cells (red) (scale bar: 100 μm). (B) Heatmap depicting Z scores as expression levels of the 25 differentially expressed genes (p < 0.01) involved in Hippo signaling pathway of SARS-CoV-2-infected hPSC-cardiomyocytes relative to uninfected cells at 3 dpi. Blue and red colors represent 6 down-regulated and 19 up-regulated genes, respectively. These differentially expressed genes in the Hippo signaling pathway were identified using the KEGG pathway database. The gene expression data was retrieved at Gene Expression Omnibus with accession number GSE150392. Individual quantitative observations that underlie the data summarized can be located under the Supporting Information File as S4 Data.