bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
1 A transcriptional regulatory atlas of coronavirus infection of
2 human cells
3 Scott A Ochsner & Neil J McKenna
4 The Signaling Pathways Project and Department of Molecular and
5 Cellular Biology, Baylor College of Medicine, Houston, TX 77030
6 Address Correspondence To:
7 Neil J. McKenna
8 Department of Molecular and Cellular Biology
9 Baylor College of Medicine
10 Houston, TX 77030
11 USA
12 e: [email protected]
13 t: 713-798-8568
14
1
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
15 Abstract
16 Identifying transcriptional responses that are most consistently associated with
17 experimental coronavirus (CoV) infection can help illuminate human cellular signaling
18 pathways impacted by CoV infection. Here, we distilled over 3,000,000 data points from
19 publically archived CoV infection transcriptomic datasets into consensus regulatory
20 signatures, or consensomes, that rank genes based on their transcriptional
21 responsiveness to infection of human cells by MERS, SARS-CoV-1 and SARS-CoV-2
22 subtypes. We computed overlap between genes with elevated rankings in the CoV
23 consensomes against those from transcriptomic and ChIP-Seq consensomes for nearly
24 880 cellular signaling pathway nodes. Validating the CoV infection consensomes, we
25 identified robust overlap between their highly ranked genes and high confidence targets
26 of signaling pathway nodes with known roles in CoV infection. We then developed a
27 series of use cases that illustrate the utility of the CoV consensomes for hypothesis
28 generation around mechanistic aspects of the cellular response to CoV infection. We
29 make the CoV infection datasets and their universe of underlying data points freely
30 accessible through the Signaling Pathways Project web knowledgebase at
31 https://www.signalingpathways.org/datasets/index.jsf.
32
2
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
33 Introduction
34 Infection by coronaviruses (CoV) represents a major current global public health
35 concern. Signaling within and between airway epithelial and immune cells in response
36 to infections by CoV and other viruses is coordinated by a complex network of signaling
37 pathway nodes. These include chemokine and cytokine-activated receptors, signaling
38 enzymes and transcription factors, and the genomic targets encoding their downstream
39 effectors (Takeda et al., 2003; Stark et al., 1998; Darnell et al., 1994)]. We recently
40 described the Signaling Pathways Project (SPP; (Ochsner et al., 2019), an integrated
41 ‘omics knowledgebase designed to assist bench researchers in leveraging publically
42 archived transcriptomic and ChIP-Seq datasets to generate research hypotheses. A
43 unique aspect of SPP is its collection of consensus regulatory signatures, or
44 consensomes, which rank genes based on the frequency of their significant differential
45 expression across transcriptomic experiments mapped to a specific signaling pathway
46 node or node family. By surveying across multiple independent datasets, we have
47 shown that consensomes recapitulate pathway node-genomic target regulatory
48 relationships to a high confidence level (Ochsner et al., 2019).
49 Placing the transcriptional events associated with human CoV infection in context with
50 those associated with other signaling paradigms has the potential to catalyze the
51 development of novel therapeutic approaches. The CoV research community has been
52 active in generating and archiving transcriptomic datasets documenting the
53 transcriptional response of human cells to infection by the three major CoV species,
54 namely, Middle East respiratory syndrome coronavirus (MERS) and severe acute
55 respiratory syndrome coronaviruses 1 (SARS1) and 2 (SARS2) (DeDiego et al., 2011;
3
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
56 Josset et al., 2013; Sims et al., 2013; Yoshikawa et al., 2010). To date however the field
57 has lacked a resource that fully capitalizes on these datasets by, firstly, using them to
58 rank human genes according to their transcriptional response to CoV infection and
59 secondly, contextualizing these transcriptional responses by integrating them with
60 ‘omics data points relevant to host cellular signaling pathways. Here, as a service to the
61 research community to catalyze the development of novel CoV therapeutics, we
62 generated consensomes for infection of human cells by MERS, SARS1 and SARS2
63 CoVs. We then analyzed the extent to which high confidence transcriptional targets for
64 these viruses intersected with genes with elevated rankings in transcriptomic and ChIP-
65 Seq consensomes for cellular signaling pathway nodes. Integration of the CoV
66 consensomes with the existing universes of SPP transcriptomic and ChIP-Seq data
67 points in a series of use cases illuminates previously uncharacterized intersections
68 between CoV infection and human cellular signaling pathways. The CoV infection
69 consensome and its underlying datasets provide researchers with a unique and freely
70 accessible framework within which to generate and pressure test hypotheses around
71 human cellular signaling pathways impacted by CoV infection.
72
73
74
4
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
75 Results
76 Generation of the CoV consensomes
77 We first set out to generate a set of consensomes ranking human genes based on the
78 frequency of their significant differential expression in response to infection by MERS,
79 SARS1 and SARS2 CoVs. To do this we searched the Gene Expression Omnibus
80 (GEO) and ArrayExpress databases to identify datasets involving infection of human
81 cells by these species. From this initial collection of datasets, we next carried out a
82 three step quality control check as previously described (Ochsner et al., 2019), yielding
83 a total of 3,041,047 million data points in 111 experiments from 25 independent CoV
84 infection transcriptomic datasets (Supplementary information, Section 1). Using these
85 curated datasets, we next generated consensomes for each CoV species, as well as
86 one ranking genes across all CoV infection experiments (ALL CoV). As a reference
87 consensome for a virus whose transcriptional impact on human cells has been studied
88 in depth, we also generated a consensome for human influenza A virus (IAV) infection.
89 The Supplementary information files contain the full human ALL CoV (Section 2), MERS
90 (Section 3), SARS1 (Section 4), SARS2 (Section 5) and IAV (Section 6) infection
91 transcriptomic consensomes. To assist researchers in inferring transcriptional regulation
92 of signaling networks, the ALL CoV consensome is annotated to indicate the identity of
93 a gene as encoding a bioactive small molecule, receptor, signaling enzyme or
94 transcription factor (Supplementary information, Section 2, columns W-Z). As an initial
95 benchmark for the ALL CoV consensome, we assembled a list of 20 interferon-
96 stimulated genes (ISGs), which encode many of the key canonical viral response
97 factors (Supplementary information, Section 2, column I) (Schneider et al., 2014). As
5
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
98 shown in the scatterplot representation of the ALL COV consensome in Figure 1, all
99 canonical ISGs were assigned appropriately elevated rankings in the consensome.
100 To gain insight into transcriptional intersections between CoV infection and human
101 cellular signaling pathways, we next computed genes with elevated rankings in the CoV
102 consensomes against genes with high confidence regulatory relationships with cellular
103 signaling pathway nodes. To do this we generated five lists of genes corresponding to
104 the ALL CoV, MERS, SARS1, SARS2 and IAV transcriptomic consensome 95th
105 percentiles. We then retrieved genes in the 95th percentiles of human transcriptomic (n
106 = 30) consensomes and ChIP-Seq (n = 834) consensomes for a collection of cellular
107 signaling pathway nodes, and computed the extent and significance of their overlap with
108 genes in the 95th percentiles of each of the five CoV consensomes. Significant overlap
109 between the 95th percentiles of a node/node family consensome and a CoV
110 consensome would indicate a potential biological relationship between loss or gain of
111 function of that node and the transcriptional response to CoV infection.
112 For clarity and brevity we will refer from here on to the 95th percentile of a ChIP-Seq
113 consensome as “CC95” and the 95th percentile of a transcriptomic consensome as
114 “TC95”. Of the 834 node CC95s evaluated, 377 had significant overlap (p < 0.05) with at
115 least one of the CoV infection TC95s; of the 30 node/node family TC95s evaluated, 25
116 had significant overlap (p< 0.05) with at least one of the CoV infection TC95s. Results of
117 these analyses are shown in Figure 2 (receptor and enzyme transcriptomic analysis),
118 Figure 3 (ChIP-Seq transcription factor analysis), Figure 4 (ChIP-Seq enzyme analysis)
119 and Figure 5 (ChIP-Seq co-node analysis); see also Supplementary information Section
120 7 for the complete numerical data. We next surveyed the significant overlaps to identify
6
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
121 (i) canonical inflammatory signaling pathway nodes with characterized roles in the
122 response to CoV infection, thereby validating the consensome approach in this context;
123 and (ii) evidence for previously uncharacterized transcriptional biology of CoV infection
124 that are consistent with their roles in the response to other viral infections.
125 Receptors Reflecting their well-documented roles in the response to viral infection, we
126 observed appreciable significant overlap between all TC95s and those for the toll-like
127 (TLR (Totura et al., 2015), interferon (IFNR (Hensley et al., 2004)) and tumor necrosis
128 factor (TNFR (W. Wang et al., 2007)) receptor families (Fig. 2). Interestingly, these
129 signatures were particularly highly enriched in the SARS2 TC95 – potentially reflecting a
130 particularly strong antiviral response to infection by SARS2. TC95 overlaps for receptor
131 systems with previously uncharacterized connections to CoV infection, including
132 epidermal growth factor receptors, glucocorticoid receptor and NOTCH receptor
133 signaling, are consistent with the known roles of these receptor systems in the context
134 of other viral infections (Hayward, 2004; Ito et al., 2011; Ng et al., 2013; Ostler et al.,
135 2019; Zheng et al., 2014). The relatively strong enrichment for xenobiotic receptors
136 reflects work describing a role for pregnane X receptor in innate immunity (S. Wang et
137 al., 2014) and points to a potential role for members of this family in the response to
138 CoV infection.
139 Transcription factors In general, ChIP-Seq enrichments for transcription factors and
140 other nodes were more specific for individual CoV infection TC95s (compare Fig. 2 with
141 Figs 3, 4 and 5). This is likely due to the fact that ChIP-Seq consensomes are based on
142 direct promoter binding by a specific node antigen, whereas transcriptomic
143 consensomes encompass both direct and indirect targets of specific receptor and
7
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
144 enzyme node families. Not unexpectedly – and speaking again to validation of the
145 consensomes - the strongest and most significant CoV TC95 overlaps were observed
146 for CC95s for known transcription factor mediators of the transcriptional response to
147 CoV infection, including members of the NFκB (Ludwig & Planz, 2008; Poppe et al.,
148 2017; Ruckle et al., 2012), IRF (Chiang & Liu, 2018) and STAT (Blaszczyk et al., 2016;
149 Frieman et al., 2007; Garcia-Sastre et al., 1998) transcription factor families. Consistent
150 with its known role in the regulation of interferon-stimulated genes (Hasan et al., 2013),
151 we also observed appreciable overlap between the CC95 for TFEB and the ALL CoV
152 TC95. Moreover, the strong overlap between the GTF2B/TFIIB CC95 and all viral
153 TC95s reflects previous evidence identifying GTF2B as a target for orthomyxovirus
154 (Haas et al., 2018), herpes simplex virus (Gelev et al., 2014) and hepatitis B virus
155 (Haviv et al., 1998).
156 Enzymes Compared to the roles of receptors and transcription factors in the response
157 to viral infection, the roles of signaling enzymes are less well illuminated – indeed, in the
158 context of CoV infection, they are entirely unstudied. Through their regulation of cell
159 cycle transitions, cyclin-dependent kinases play important roles in the orchestration of
160 DNA replication and cell division, processes that are critical in the viral life cycle. CDK6,
161 which has been suggested to be a critical G1 phase kinase (Bellail et al., 2014; Grossel
162 & Hinds, 2006), has been shown to be targeted by a number of viral infections, including
163 Kaposi's sarcoma-associated herpesvirus (Kaldis et al., 2001) and HIV-1 (Pauls et al.,
164 2014). Consistent with this common role across distinct viral infections, we observed
165 robust overlap between the CDK95 TC95 (Fig. 2) and the CDK6 CC95 (Fig. 4) and
166 those of all viral TC95s. As with the TLRs, IFNRs and TNFRs, which are known to
8
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
167 signal through CDK6 (Cingoz & Goff, 2018; Handschick et al., 2014; Hennessy et al.,
168 2011), overlap with the CDK6 CC95 was particularly strong in the case of the SARS2
169 TC95 (Fig. 4). CCNT2 is a member of the highly conserved family cyclin family and,
170 along with CDK9, is a member of the viral-targeted p-TEFB complex (Zaborowska et al.,
171 2016). Reflecting a potential general role in viral infection, appreciable overlap was
172 observed between the CCNT2 CC95 and all viral TC95s (Fig. 4). Finally in the context
173 of enzymes, DNA Topoisomerase (TOP1) has been shown to be required for efficient
174 replication of simian virus 40 (Wobbe et al., 1987) and Ebola (Takahashi et al., 2013)
175 viruses. The prominent overlap between its CC95 and those of SARS2 and IAV (Fig. 4)
176 suggest that it may play a similar role in facilitating the replication of these viruses.
177 Co-nodes We have coined the term “co-nodes” as a convenient catch-all for cellular
178 factors that are not members of the three principal signaling pathway node categories
179 (receptors, enzymes and transcription factors; Ochsner et al., 2019). The breadth of
180 functions encompassed by members of this category reflects the diverse mechanisms
181 employed both by viruses to infect and propagate in cells, as well as by hosts to mount
182 an efficient immune response. Consistent with its characterized role in the recruitment
183 of p-TEFB by IV-1 Tat protein (Schulze-Gahmen et al., 2013), we observed consistently
184 strong enrichment of the AFF4 CC95 in all viral TC95s (Fig. 5). The targeting of CNOT3
185 for degradation in response to adenoviral infection (Chalabi Hagkarim et al., 2018) is
186 reflected in the significant overlap between its CC95 and the viral TC95s, particularly in
187 the case of the SARS2 (Fig. 5).
188 By way of additional independent validation of the CoV and IAV consensomes, Gene
189 Set Enrichment Analysis (Subramanian et al., 2005) reflected significant overlap
9
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
190 between the CoV and IAV TC95s and a variety of viral infection and inflammatory
191 transcription factor target gene sets (Supplementary information Section 8).
192
193 Use cases
194 Having established the reliability of the consensome approach in the context of CoV
195 infection, we next wished to use the CoV consensomes to gather evidence to identify
196 previously uncharacterized mediators of the transcriptional response to infection by
197 CoVs and, of particular current interest, SARS2.
198
199 Use case 1: consensome redundancy analysis identifies potential
200 uncharacterized players in the response to CoV infection
201 We previously showed (Figs. 2 and 3) that the overlap of the CoV TC95 genes was
202 most robust among consensomes for members of the IFNR, TLR and TNF receptor
203 families (Fig. 2), and the NFkB, RELA, IRF and STAT transcription factor families (fig.
204 5). To investigate this further, we ranked genes in the ALL CoV consensome by their
205 aggregate 80th percentile rankings across these consensomes (Supplementary
206 information, Section 2, column V). Transcriptomic consensomes for IFNRs (Section 9),
207 TLRs (Section 10) and TNFRs (Section 11) are provided in Supplementary information,
208 as are links to the ChIP-Atlas lists used to generate the ChIP-Seq consensomes
209 (Section 12).
10
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
210 This redundancy ranking elevates genes with known critical roles in the response to
211 viral infection that are acutely responsive to a spectrum of inflammatory signaling
212 nodes, such as NCOA7 (Doyle et al., 2018), STAT (Chapgier et al., 2009) and TAP1
213 (Gruter et al., 1998). Interestingly, genes such as PSMB9, CSRNP1 and MRPL24 have
214 ALL CoV discovery rates that are comparable to or exceed those of many of the classic
215 viral response ISGs (Fig. 1), but are either largely or completely uncharacterized in the
216 context of viral infection. This use case the reflects the potential of consensome-driven
217 data mining to illuminate novel and potentially therapeutically relevant transcription
218 pathway effectors in the response to CoV infection.
219
220 Use case 2: discrimination between genes frequently differentially expressed in
221 response to CoV, but not IAV, infection
222 We next wished to gain insight into signaling paradigms that were selectively
223 characteristic of CoV infection but not IAV infection. To do this we generated a set of
224 genes that were frequently differentially expressed in response to CoV infection (ALL
225 CoV TC99), but less frequently expressed in response to IAV infection (50th percentile)
226 (Table 1). Interestingly, this group contained two genes encoding transcription factors
227 with well characterized roles in the regulation of oscillatory gene expression, PER1 and
228 PER2 (Yamajuku et al., 2010). This led us to speculate that genomic targets of
229 transcriptional regulators of circadian rhythms might be enriched among high
230 confidence CoV-regulated genes, but not IAV-regulated genes. Consistent with this we
231 observed significant overlap between the CC95s of the master circadian transcription
232 factors ARNTL1/BMAL1 and CLOCK (Gaucher et al., 2018) and the CoV TC95, but not
11
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
233 the IAV TC95 (Fig. 2). These data indicate a hitherto uncharacterized intersection
234 between CoV signaling and the cellular oscillatory apparatus. Another gene of note in
235 this group is INHBA, encoding activin A, overexpression of which in murine lung gives
236 rise to a phenotype resembling acute respiratory distress-like syndrome (Apostolou et
237 al., 2012), a condition commonly associated with CoV infection (Totura & Bavari, 2019).
238 Other genes in the group provide evidence for the potential mechanisms of infection
239 and propagation of CoVs. TJAP1, for example, encodes a member of a family of
240 proteins involved in regulation of tight junctions, shown to be route of porcine epidemic
241 diarrhea coronavirus entry into epithelial cells (Luo et al., 2017). In addition, GON7 is a
242 member of the KEOPS complex (Wan et al., 2017), which is involved in
243 threonylcarbamoylation of tRNAs, which represent an important host facet of the
244 retroviral life cycle (Jin & Musier-Forsyth, 2019). Interestingly, another critical gene in
245 this tRNA pathway, YRDC, is the 19th ranked gene in the ALL CoV consensome
246 (Supplementary information Section 2).
247
248 Use case 3: evidence for antagonism between progesterone receptor and
249 interferon receptor signaling in the airway epithelium
250 Although a lack of clinical data has so far prevented a definitive evaluation of the
251 connection between pregnancy and susceptibility to SARS2 infection in COVID-19,
252 pregnancy has been previously associated with the incidence of viral infectious
253 diseases, particularly respiratory infections (Sappenfield et al., 2013; Siston et al.,
254 2010). We were interested therefore to see consistent overlap between the
255 progesterone receptor (PGR) TC95 and all the viral TC95s, with the enrichment being
12
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
256 particularly evident in the case of the SARS2 TC95 (Fig. 2). To investigate the specific
257 nature of the crosstalk implied by this overlap in the context of the airway epithelium, we
258 first identified a set of 16 genes that were in both the ALL CoV and PGR TC95s. We
259 then retrieved two SPP experiments involving treatment of A549 airway epithelial cells
260 with the PGR full antagonist RU486 (RU), alone or in combination with the GR agonist
261 dexamethasone (DEX). As shown in Figure 6 , there was nearly unanimous correlation
262 in the direction of regulation of all 16 genes in response to CoV infection and PGR loss
263 of function. These data indicate that antagonism between PGR and IFNR signaling in
264 the airway epithelium may predispose pregnant women to infection by SARS2.
265
266 Use case 4: evidence for a role for the telomerase catalytic subunit TERT in the
267 interferon response to viral infection
268 Although telomerase activation has been well characterized in the context of cell
269 immortalization by human tumor virus infection (Gewin et al., 2004; Klingelhutz et al.,
270 1996; Yang et al., 2004), no connection has previously been made between CoV or IAV
271 infection and telomerase. We were therefore intrigued to observe robust overlap
272 between all viral TC95s and that of the telomerase catalytic subunit TERT. In support of
273 this finding, NFkB signaling has been shown to induce expression (Yin et al., 2000) and
274 nuclear translocation (Akiyama et al., 2003) of TERT, and direct co-regulation by
275 telomerase of NFkB-dependent transcription has been linked to chronic inflammation
276 (Ghosh et al., 2012). Inspecting the ALL CoV consensome underlying data points (data
277 not shown) we found that the TERT gene was not transcriptionally induced in response
278 to infection by any of the CoVs, indicating that the overlap between its TC95 and those
13
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
279 of the CoVs might occur in response to an upstream regulatory signal. If functional
280 interactions between TERT and inflammatory nodes did indeed take place in response
281 to CoV infection, we anticipated that this would be reflected in close agreement
282 regarding the direction of differential expression of CoV infection-regulated genes in
283 response to perturbation of TERT on the one hand, and on the other, to stimulation of
284 the classic viral response IFNRs. To test this hypothesis, we took the same set of 20
285 ISGs referred to previously (Fig. 1) and compared their direction of regulation across all
286 experiments underlying the ALL CoV, TERT and IFNR consensomes. For reference, the
287 TERT consensome (Section 13) and its underlying data points (Section 14) are
288 provided in the Supplementary information. With respect to the IFNR and TERT data
289 points, we observed a nearly universal alignment in the direction of regulation of all
290 genes tested with those in the CoV infection experiments (Fig. 6), with agreement in the
291 direction of regulation across 99% of the underlying probesets. We should note that of
292 the 1859 p < 0.05 CoV infection ISG data points, we observed repression, rather than
293 induction, in response to CoV infection in 303 (~15%) , which may be attributable to the
294 impact of differences in cell type, cell cycle stage or virus incubation time across the
295 independent experiments. Our results are consistent with a model in which activation of
296 telomerase is a component of the human innate immune response to viral infection.
297
298 Use case 5: SARS2 infection of human cells is specifically associated with an
299 epithelial to mesenchymal transition transcriptional signature
300 Epithelial to mesenchymal transition (EMT) is the process by which epithelial cells lose
301 their polarity and adhesive properties and acquire the migratory and invasive 14
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
302 characteristics of mesenchymal stem cells (Lamouille et al., 2014). EMT is known to
303 contribute to pulmonary fibrosis (Hill et al., 2019) and acute interstitial pneumonia (H. Li
304 et al., 2014), both of which have been reported in connection with SARS2 infection in
305 COVID-19 (Adair & Ledermann, 2020; P. Zhou et al., 2020). Moreover, EMT is widely
306 accepted as a core component of the process by which renal tubular cells transform into
307 mesenchymal cells during the development of fibrosis in kidney disease (Y. Liu, 2006),
308 a signature comorbidity of SARS2 infection (Durvasula et al., 2020). Of the 834 CC95s
309 analyzed, overlap (p < 0.05) was specific to the SARS2 CC95 for only five; SNAI2/Slug,
310 SOX2, GATA6, CTBP1 and PRMT1. Strikingly, a literature search indicated that these
311 five nodes were connected by documented roles in in the promotion of EMT (Avasarala
312 et al., 2015; Herreros-Villanueva et al., 2013; Martinelli et al., 2017; Nieto, 2002; Sahu
313 et al., 2017). In addition to these nodes, whose C95 genes were exclusively enriched (p
314 < 0.05) in the SARS2 C95 genes, we identified several other EMT-linked nodes whose
315 CC95 genes were preferentially enriched (p < 0.05) in the SARS2 TC95 genes (Figs. 3-
316 5), including the homeodomain transcription factor SIX2 (C.-A. Wang et al., 2014),
317 SMAD4 (Siraj et al., 2019), and the co-nodes PYGO2 (Chi et al., 2019), SKI (Tecalco-
318 Cruz et al., 2018), BRD7 (T. Liu et al., 2017) and STAG2 (Nie et al., 2019).
319 To investigate this further, we computed overlap between the individual viral TC99s and
320 a list of 335 genes manually curated from the research literature as signature EMT
321 markers (Supplementary information Section 15; Zhao et al., 2015). Consistent with the
322 node consensome enrichment analysis, we observed significant enrichment of
323 members of this gene set within the SARS2 CC99, but not those of the ALL CoV,
324 SARS1, MERS or IAV consensomes (Supplementary information Section 16). One
15
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
325 possible explanation for this was the fact that the SARS2 consensome was comprised
326 of airway epithelial cell lines, whereas the SARS1 and MERS consensomes included
327 non-epithelial cell biosamples (Supplementary information Section 1). To exclude this
328 possibility therefore, we next calculated airway epithelial cell-specific consensomes for
329 SARS1 and MERS and computed overlap of their TC95s against the 864 pathway
330 node/node family CC95s & TC95s. We found that significant overlap with the CoV
331 TC95s remained specific to SARS2 (data not shown), confirming that significant overlap
332 with the EMT node signature was specific to the SARS2 TC99.
333 We next applied consensome redundancy analysis (see Use case 1) to isolate a set of
334 SARS2 regulated genes (CPV P < 0.05) that were high confidence targets (i.e. in the
335 CC80) for at least two of the EMT nodes (SNAI2, SOX2, GATA6, CTBP1 and PRMT1;
336 Supplementary information, Section 4, column M). A literature survey showed that 13 of
337 these 21 genes had a documented connection to EMT (Supplementary information,
338 Section 5, column N). Figure 9 compares the percentile rankings for these genes across
339 the three CoV infection consensomes and the IAV consensome. Although some EMT
340 genes, such as CXCL2 and IRF9, had elevated rankings across all four consensomes,
341 the collective EMT gene signature had a significantly higher mean percentile value in
342 the SARS2 consensome than in each of the three others (Fig. 9).
343 Although EMT has been associated with infection by transmissible gastroenteritis virus
344 (Xia et al., 2017), this is to our knowledge the first evidence connecting CoV infection,
345 and specifically SARS2 infection, to EMT. Interestingly, several members of the group
346 of SARS2-induced EMT genes have been associated with signature pulmonary
347 comorbidities of CoV infection, including ADAR (Diaz-Pina et al., 2018), CLDN1
16
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
348 (Vukmirovic et al., 2017) and SOD2 (Gao et al., 2008). Of note in the context of these
349 data is the fact that signaling through two SARS2 cellular receptors, ACE2/AT2 and
350 CD147/basigin, has been linked to EMT in the context of organ fibrosis (Kato et al.,
351 2011; Ruster & Wolf, 2011; C. Wang et al., 2018). Moreover, overlap between of the
352 CoV TC95s and the TERT CC95 was particularly robust in the case of SARS2, a finding
353 of potential relevance to the fact that telomerase has been implicated in EMT (Z. Liu et
354 al., 2013). Collectively, our data indicate that EMT warrants further investigation as a
355 SARS2-specific pathological mechanism.
356 Development of a CoV infection cell signaling knowledgebase
357 Having validated the ALL CoV consensome, we next wished to make it freely available
358 to the research community for re-use in the characterization of signaling events
359 associated with CoV infection. Firstly, the viral infection datasets were curated
360 accordingly to our previously described protocol (Ochsner et al., 2019) made available
361 for browsing in the SPP Dataset listing
362 (https://www.signalingpathways.org/datasets/index.jsf). As with other SPP datasets, and
363 per FAIR data best practice, CoV infection datasets were associated with detailed
364 descriptions, assigned a digital object identifier, and linked to the associated article to
365 place the dataset in its original experimental context. Loading the CoV datasets into the
366 SPP also automatically made the underlying data points discoverable by the SPP query
367 tool Ominer (Ochsner et al., 2019). These reports represent a summary of the current
368 state of transcriptomic and ChIP-Seq knowledge on the regulatory relationship of a
369 given gene with upstream regulatory pathway nodes, or in clinical and model
370 experiments. The full value of the integration of the CoV consensome with the existing
17
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
371 SPP framework can perhaps be best appreciated in the context of the one click links to
372 Ominer Regulation Reports from the individual CoV datasets. These Reports provide
373 researchers with a wealth of contextual data on signaling pathways impacted by CoV
374 infection in the context of a specific gene. Table 2 shows links to Regulation Reports for
375 the top twenty ranked genes in the ALL CoV consensome. The order of sections in the
376 Reports is Receptors, Enzymes, Transcription Factors, Co-nodes, Clinical and Models,
377 the last section including data points from the CoV infection model experiments that
378 form the basis of this study.
18
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
379 Discussion
380 An effective research community response to the impact of CoVs on human health
381 demands systematic exploration of the transcriptional interface between viral infection
382 and human cell signaling systems. It also demands routine access to existing datasets
383 that is unhindered either by paywalls or by lack of the informatics training required to
384 manipulate archived datasets in their unprocessed state. Moreover, the substantial
385 logistical obstacles to BSL3 certification only emphasizes the need for fullest possible
386 access to, and re-usability of, existing CoV infection datasets to focus and refine
387 hypotheses prior to carrying out in vivo CoV infection experiments. To this end, we
388 generated a set of CoV infection consensomes that rank human genes by the
389 reproducibility of their significant differential expression in response to infection of
390 human cells by CoVs. We then computed the CoV consensomes against high
391 confidence transcriptional signatures for a broad range of cellular signaling pathway
392 nodes, affording investigators with a broad range of signaling interests an entrez into
393 the study of CoV infection of human cells. The five use cases described here represent
394 illustrative examples of the types of analyses that users are empowered to carry out in
395 the CoV infection knowledgebase.
396 To democratize access to the CoV consensome and its >3,000,000 underlying data
397 points by the broadest possible audience, we have integrated them into the SPP
398 database to create a cell signaling knowledgebase for CoV infection. Incorporation of
399 the CoV data points into SPP in this manner places them in the context of millions more
400 existing SPP data points documenting transcriptional regulatory relationships between
401 pathway nodes and their genomic targets. In doing so, we afford users a unique
19
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
402 appreciation of the cellular signaling pathway nodes whose gain or loss of function in
403 response to CoV infection gives rise to these transcriptional patterns.
404 The human CoV and IAV consensomes and their underlying datasets are “living”
405 resources on SPP that will be updated and versioned with appropriate datasets. This
406 will be particularly important in the case of SARS2, as datasets involving infection of
407 human cells with this virus are necessarily currently limited in number. This will allow
408 for hardening of observations that are intriguing, but whose significance is currently
409 unclear, such as the overlap of the CoV TC95s with the TERT TC95, as well as the
410 enrichment of EMT genes among those with elevated rankings in the SARS2
411 consensome. We welcome feedback and suggestions from the research community for
412 the future development of the SPP CoV infection consensomes.
413
414
415
416
417
20
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
418 Methods
419 Dataset processing and consensome analysis.
420 Differential expression values were calculated for each gene in each experiment using
421 the limma analysis package from Bioconductor then committed to the consensome
422 analysis pipeline as previously described (Ochsner et al., 2019). Briefly, the
423 consensome algorithm surveys each experiment across all datasets and ranks genes
424 according to the frequency with which they are significantly differentially expressed. For
425 each transcript, we counted the number of experiments where the significance for
426 differential expression was ≤0.05, and then generated the binomial probability, referred
427 to as the consensome p-value (CPV), of observing that many or more nominally
428 significant experiments out of the number of experiments in which the transcript was
429 assayed, given a true probability of 0.05. A more detailed description of the
430 transcriptomic consensome algorithm is available (Ochsner et al., 2019). The
431 consensomes and underlying datasets were loaded into an Oracle 13c database and
432 made available on the SPP user interface as previously described (Ochsner et al.,
433 2019).
434 Statistical analysis
435 Gene Overlap analysis was performed using the Bioconductor GeneOverlap analysis
436 package
437 (https://www.rdocumentation.org/packages/GeneOverlap/versions/1.8.0/topics/GeneOv
438 erlap) implemented in R. Briefly, given a whole set I of IDs and two sets A ∈ I and B ∈ I,
439 and S = A ∩ B, GeneOverlap calculates the significance of obtaining S. The problem is
21
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
440 formulated as a hypergeometric distribution or a contingency table, which is solved by
441 Fisher’s exact test. The universe for the overlap was set at a recent estimate of the total
442 number of coding genes in the human genome (21,500; Pertea et al., 2018). Paired
443 Two Sample t-Test for comparing the mean percentile ranking of EMT genes in the
444 MERS, SARS1, SARS2 and IAV consensomes was performed in PRISM at 12 degrees
445 of freedom.
446 Consensome generation
447 The procedure for generating transcriptomic consensomes has been previously
448 described (Ochsner et al., 2019). To generate the ChIP-Seq consensomes, we first
449 retrieved processed gene lists from ChIP-Atlas, which rank human genes based upon
450 their average MACS2 occupancy across all publically archived datasets in which a
451 given transcription factor, enzyme or co-node is the IP antigen. Of the three stringency
452 levels available (10, 5 and 1 kb from the transcription start site), we selected the most
453 stringent (1 kb). According to SPP convention (Ochsner et al., 2019), we then mapped
454 the IP antigen to its pathway node category and class, and the experimental cell line to
455 its appropriate biosample physiological system and organ. We then organized the
456 ranked lists into percentiles to generate the node ChIP-Seq consensome.
457 SPP web application
458 The SPP knowledgebase is a gene-centric Java Enterprise Edition 6, web-based
459 application around which other gene, mRNA, protein and BSM data from external
460 databases such as NCBI are collected. After undergoing semiautomated processing
461 and biocuration as described above, the data and annotations are stored in SPP’s
22
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
462 Oracle 13c database. RESTful web services exposing SPP data, which are served to
463 responsively designed views in the user interface, were created using a Flat UI Toolkit
464 with a combination of JavaScript, D3.JS, AJAX, HTML5, and CSS3. JavaServer Faces
465 and PrimeFaces are the primary technologies behind the user interface. SPP has been
466 optimized for Firefox 24+, Chrome 30+, Safari 5.1.9+, and Internet Explorer 9+, with
467 validations performed in BrowserStack and load testing in LoadUIWeb. XML describing
468 each dataset and experiment is generated and submitted to CrossRef to mint DOIs
469 (Ochsner et al., 2019).
470
23
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
471 Data availability
472 The entire set of data points used to generate the CoV consensome has been uploaded in an R
473 file to Figshare and a link included for reviewer access. The entire set of metadata for these
474 data points is available in Supplementary information Section 1. Consensome data points are in
475 Supplementary information Sections 2-6.
476 SPP is freely accessible at https://www.signalingpathways.org. Programmatic access to all
477 underlying data points and their associated metadata are supported by a RESTful API at
478 https://www.signalingpathways.org/docs/. All SPP datasets are biocurated versions of publically
479 archived datasets, are formatted according to the recommendations of the FORCE11 Joint
480 Declaration on Data Citation Principles 74, and are made available under a Creative Commons
481 CC 3.0 BY license. The original datasets are available are linked to from the corresponding SPP
482 datasets using the original repository accession identifiers. These identifiers are for
483 transcriptomic datasets, the Gene Expression Omnibus (GEO) Series (GSE); and for
484 cistromic/ChIP-Seq datasets, the NCBI Sequence Read Archive (SRA) study identifier (SRP).
485 DOIs for the consensomes and datasets are pending.
486 Code availability
487 The full SPP source code is available in the SPP GitHub account under a Creative Commons
488 CC 3.0 BY license at https://github.com/signaling-pathways-project/ominer/.
24
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
489 Acknowledgements
490 This work was supported by the National Institute of Diabetes, Digestive and Kidney
491 Diseases NIDDK Information Network U24 (DK097748), the National Cancer Institute
492 (CA125123) and by the Brockman Medical Research Foundation. We appreciate the
493 technical assistance of Apollo McOwiti and Shijing Qu of the Biostatistics and
494 Informatics Shared Resource of the Duncan NCI Comprehensive Cancer Center. We
495 thank all investigators who archive their datasets, without whom SPP would not be
496 possible.
497
498 Author contributions
499 Dataset biocuration: SO
500 Data analysis: SO, NM
501 Manuscript drafting: NM
502
503 Competing interests
504 The authors declare no competing interests.
505
506
25
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
507 References
508 Adair, L. B. 2nd, & Ledermann, E. J. (2020). Chest CT Findings of Early and 509 Progressive Phase COVID-19 Infection from a US Patient. In Radiology case 510 reports. https://doi.org/10.1016/j.radcr.2020.04.031
511 Akiyama, M., Hideshima, T., Hayashi, T., Tai, Y.-T., Mitsiades, C. S., Mitsiades, N., 512 Chauhan, D., Richardson, P., Munshi, N. C., & Anderson, K. C. (2003). Nuclear 513 factor-kappaB p65 mediates tumor necrosis factor alpha-induced nuclear 514 translocation of telomerase reverse transcriptase protein. Cancer Research, 63(1), 515 18–21.
516 Apostolou, E., Stavropoulos, A., Sountoulidis, A., Xirakia, C., Giaglis, S., 517 Protopapadakis, E., Ritis, K., Mentzelopoulos, S., Pasternack, A., Foster, M., 518 Ritvos, O., Tzelepis, G. E., Andreakos, E., & Sideras, P. (2012). Activin-A 519 overexpression in the murine lung causes pathology that simulates acute 520 respiratory distress syndrome. American Journal of Respiratory and Critical Care 521 Medicine, 185(4), 382–391. https://doi.org/10.1164/rccm.201105-0784OC
522 Avasarala, S., Van Scoyk, M., Karuppusamy Rathinam, M. K., Zerayesus, S., Zhao, X., 523 Zhang, W., Pergande, M. R., Borgia, J. A., DeGregori, J., Port, J. D., Winn, R. A., & 524 Bikkavilli, R. K. (2015). PRMT1 Is a Novel Regulator of Epithelial-Mesenchymal- 525 Transition in Non-small Cell Lung Cancer. The Journal of Biological Chemistry, 526 290(21), 13479–13489. https://doi.org/10.1074/jbc.M114.636050
527 Barata, J. T., Durum, S. K., & Seddon, B. (2019). Flip the coin: IL-7 and IL-7R in health 528 and disease. Nature Immunology, 20(12), 1584–1593. 529 https://doi.org/10.1038/s41590-019-0479-x
530 Bellail, A. C., Olson, J. J., & Hao, C. (2014). SUMO1 modification stabilizes CDK6 531 protein and drives the cell cycle and glioblastoma progression. Nature 532 Communications, 5, 4234. https://doi.org/10.1038/ncomms5234
533 Blaszczyk, K., Nowicka, H., Kostyrko, K., Antonczyk, A., Wesoly, J., & Bluyssen, H. A. 534 R. (2016). The unique role of STAT2 in constitutive and IFN-induced transcription 535 and antiviral responses. Cytokine & Growth Factor Reviews, 29, 71–81. 536 https://doi.org/10.1016/j.cytogfr.2016.02.010
537 Chalabi Hagkarim, N., Ryan, E. L., Byrd, P. J., Hollingworth, R., Shimwell, N. J., 538 Agathanggelou, A., Vavasseur, M., Kolbe, V., Speiseder, T., Dobner, T., Stewart, 539 G. S., & Grand, R. J. (2018). Degradation of a Novel DNA Damage Response 540 Protein, Tankyrase 1 Binding Protein 1, following Adenovirus Infection. Journal of 541 Virology, 92(12). https://doi.org/10.1128/JVI.02034-17
542 Chang, B., Yang, H., Jiao, Y., Wang, K., Liu, Z., Wu, P., Li, S., & Wang, A. (2016). 543 SOD2 deregulation enhances migration, invasion and has poor prognosis in 544 salivary adenoid cystic carcinoma. Scientific Reports, 6, 25918. 545 https://doi.org/10.1038/srep25918
26
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
546 Chapgier, A., Kong, X.-F., Boisson-Dupuis, S., Jouanguy, E., Averbuch, D., Feinberg, 547 J., Zhang, S.-Y., Bustamante, J., Vogt, G., Lejeune, J., Mayola, E., de 548 Beaucoudrey, L., Abel, L., Engelhard, D., & Casanova, J.-L. (2009). A partial form 549 of recessive STAT1 deficiency in humans. The Journal of Clinical Investigation, 550 119(6), 1502–1514. https://doi.org/10.1172/JCI37083
551 Chi, Y., Wang, F., Zhang, T., Xu, H., Zhang, Y., Shan, Z., Wu, S., Fan, Q., & Sun, Y. 552 (2019). miR-516a-3p inhibits breast cancer cell growth and EMT by blocking the 553 Pygo2/Wnt signalling pathway. Journal of Cellular and Molecular Medicine, 23(9), 554 6295–6307. https://doi.org/10.1111/jcmm.14515
555 Chiang, H.-S., & Liu, H. M. (2018). The Molecular Basis of Viral Inhibition of IRF- and 556 STAT-Dependent Immune Responses. Frontiers in Immunology, 9, 3086. 557 https://doi.org/10.3389/fimmu.2018.03086
558 Cingoz, O., & Goff, S. P. (2018). Cyclin-dependent kinase activity is required for type I 559 interferon production. Proceedings of the National Academy of Sciences of the 560 United States of America, 115(13), E2950–E2959. 561 https://doi.org/10.1073/pnas.1720431115
562 Darnell, J. E. J., Kerr, I. M., & Stark, G. R. (1994). Jak-STAT pathways and 563 transcriptional activation in response to IFNs and other extracellular signaling 564 proteins. Science (New York, N.Y.), 264(5164), 1415–1421. 565 https://doi.org/10.1126/science.8197455
566 DeDiego, M. L., Nieto-Torres, J. L., Jimenez-Guardeno, J. M., Regla-Nava, J. A., 567 Alvarez, E., Oliveros, J. C., Zhao, J., Fett, C., Perlman, S., & Enjuanes, L. (2011). 568 Severe acute respiratory syndrome coronavirus envelope protein regulates cell 569 stress response and apoptosis. PLoS Pathogens, 7(10), e1002315. 570 https://doi.org/10.1371/journal.ppat.1002315
571 Diaz-Pina, G., Ordonez-Razo, R. M., Montes, E., Paramo, I., Becerril, C., Salgado, A., 572 Santibanez-Salgado, J. A., Maldonado, M., & Ruiz, V. (2018). The Role of ADAR1 573 and ADAR2 in the Regulation of miRNA-21 in Idiopathic Pulmonary Fibrosis. Lung, 574 196(4), 393–400. https://doi.org/10.1007/s00408-018-0115-9
575 Doherty, M. R., Cheon, H., Junk, D. J., Vinayak, S., Varadan, V., Telli, M. L., Ford, J. 576 M., Stark, G. R., & Jackson, M. W. (2017). Interferon-beta represses cancer stem 577 cell properties in triple-negative breast cancer. Proceedings of the National 578 Academy of Sciences of the United States of America, 114(52), 13792–13797. 579 https://doi.org/10.1073/pnas.1713728114
580 Doyle, T., Moncorge, O., Bonaventure, B., Pollpeter, D., Lussignol, M., Tauziet, M., 581 Apolonia, L., Catanese, M.-T., Goujon, C., & Malim, M. H. (2018). The interferon- 582 inducible isoform of NCOA7 inhibits endosome-mediated viral entry. Nature 583 Microbiology, 3(12), 1369–1376. https://doi.org/10.1038/s41564-018-0273-9
584 Durvasula, R., Wellington, T., McNamara, E., & Watnick, S. (2020). COVID-19 and 585 Kidney Failure in the Acute Care Setting: Our Experience From Seattle. In
27
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
586 American journal of kidney diseases : the official journal of the National Kidney 587 Foundation. https://doi.org/10.1053/j.ajkd.2020.04.001
588 Ercolano, G., De Cicco, P., Rubino, V., Terrazzano, G., Ruggiero, G., Carriero, R., 589 Kunderfranco, P., & Ianaro, A. (2019). Knockdown of PTGS2 by CRISPR/CAS9 590 System Designates a New Potential Gene Target for Melanoma Treatment. 591 Frontiers in Pharmacology, 10, 1456. https://doi.org/10.3389/fphar.2019.01456
592 Frieman, M., Yount, B., Heise, M., Kopecky-Bromberg, S. A., Palese, P., & Baric, R. S. 593 (2007). Severe acute respiratory syndrome coronavirus ORF6 antagonizes STAT1 594 function by sequestering nuclear import factors on the rough endoplasmic 595 reticulum/Golgi membrane. Journal of Virology, 81(18), 9812–9824. 596 https://doi.org/10.1128/JVI.01012-07
597 Gao, F., Kinnula, V. L., Myllarniemi, M., & Oury, T. D. (2008). Extracellular superoxide 598 dismutase in pulmonary fibrosis. Antioxidants & Redox Signaling, 10(2), 343–354. 599 https://doi.org/10.1089/ars.2007.1908
600 Garcia-Sastre, A., Durbin, R. K., Zheng, H., Palese, P., Gertner, R., Levy, D. E., & 601 Durbin, J. E. (1998). The role of interferon in influenza virus tissue tropism. Journal 602 of Virology, 72(11), 8550–8558.
603 Gaucher, J., Montellier, E., & Sassone-Corsi, P. (2018). Molecular Cogs: Interplay 604 between Circadian Clock and Cell Cycle. Trends in Cell Biology, 28(5), 368–379. 605 https://doi.org/10.1016/j.tcb.2018.01.006
606 Gelev, V., Zabolotny, J. M., Lange, M., Hiromura, M., Yoo, S. W., Orlando, J. S., 607 Kushnir, A., Horikoshi, N., Paquet, E., Bachvarov, D., Schaffer, P. A., & Usheva, A. 608 (2014). A new paradigm for transcription factor TFIIB functionality. Scientific 609 Reports, 4, 3664. https://doi.org/10.1038/srep03664
610 Gewin, L., Myers, H., Kiyono, T., & Galloway, D. A. (2004). Identification of a novel 611 telomerase repressor that interacts with the human papillomavirus type-16 E6/E6- 612 AP complex. Genes & Development, 18(18), 2269–2282. 613 https://doi.org/10.1101/gad.1214704
614 Ghosh, A., Saginc, G., Leow, S. C., Khattar, E., Shin, E. M., Yan, T. D., Wong, M., 615 Zhang, Z., Li, G., Sung, W.-K., Zhou, J., Chng, W. J., Li, S., Liu, E., & Tergaonkar, 616 V. (2012). Telomerase directly regulates NF-kappaB-dependent transcription. 617 Nature Cell Biology, 14(12), 1270–1281. https://doi.org/10.1038/ncb2621
618 Grossel, M. J., & Hinds, P. W. (2006). Beyond the cell cycle: a new role for Cdk6 in 619 differentiation. Journal of Cellular Biochemistry, 97(3), 485–493. 620 https://doi.org/10.1002/jcb.20712
621 Gruter, P., Tabernero, C., von Kobbe, C., Schmitt, C., Saavedra, C., Bachi, A., Wilm, 622 M., Felber, B. K., & Izaurralde, E. (1998). TAP, the human homolog of Mex67p, 623 mediates CTE-dependent RNA export from the nucleus. Molecular Cell, 1(5), 649– 624 659. https://doi.org/10.1016/s1097-2765(00)80065-9
28
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
625 Haas, D. A., Meiler, A., Geiger, K., Vogt, C., Preuss, E., Kochs, G., & Pichlmair, A. 626 (2018). Viral targeting of TFIIB impairs de novo polymerase II recruitment and 627 affects antiviral immunity. PLoS Pathogens, 14(4), e1006980. 628 https://doi.org/10.1371/journal.ppat.1006980
629 Handschick, K., Beuerlein, K., Jurida, L., Bartkuhn, M., Muller, H., Soelch, J., Weber, A., 630 Dittrich-Breiholz, O., Schneider, H., Scharfe, M., Jarek, M., Stellzig, J., Schmitz, M. 631 L., & Kracht, M. (2014). Cyclin-dependent kinase 6 is a chromatin-bound cofactor 632 for NF-kappaB-dependent gene expression. Molecular Cell, 53(2), 193–208. 633 https://doi.org/10.1016/j.molcel.2013.12.002
634 Hasan, M., Koch, J., Rakheja, D., Pattnaik, A. K., Brugarolas, J., Dozmorov, I., Levine, 635 B., Wakeland, E. K., Lee-Kirsch, M. A., & Yan, N. (2013). Trex1 regulates 636 lysosomal biogenesis and interferon-independent activation of antiviral genes. 637 Nature Immunology, 14(1), 61–71. https://doi.org/10.1038/ni.2475
638 Haviv, I., Shamay, M., Doitsh, G., & Shaul, Y. (1998). Hepatitis B virus pX targets TFIIB 639 in transcription coactivation. Molecular and Cellular Biology, 18(3), 1562–1569. 640 https://doi.org/10.1128/mcb.18.3.1562
641 Hayward, S. D. (2004). Viral interactions with the Notch pathway. Seminars in Cancer 642 Biology, 14(5), 387–396. https://doi.org/10.1016/j.semcancer.2004.04.018
643 Hennessy, E. J., Sheedy, F. J., Santamaria, D., Barbacid, M., & O’Neill, L. A. J. (2011). 644 Toll-like receptor-4 (TLR4) down-regulates microRNA-107, increasing macrophage 645 adhesion via cyclin-dependent kinase 6. The Journal of Biological Chemistry, 646 286(29), 25531–25539. https://doi.org/10.1074/jbc.M111.256206
647 Hensley, L. E., Fritz, L. E., Jahrling, P. B., Karp, C. L., Huggins, J. W., & Geisbert, T. W. 648 (2004). Interferon-beta 1a and SARS coronavirus replication. Emerging Infectious 649 Diseases, 10(2), 317–319. https://doi.org/10.3201/eid1002.030482
650 Herreros-Villanueva, M., Zhang, J.-S., Koenig, A., Abel, E. V, Smyrk, T. C., Bamlet, W. 651 R., de Narvajas, A. A.-M., Gomez, T. S., Simeone, D. M., Bujanda, L., & Billadeau, 652 D. D. (2013). SOX2 promotes dedifferentiation and imparts stem cell-like features 653 to pancreatic cancer cells. Oncogenesis, 2, e61. 654 https://doi.org/10.1038/oncsis.2013.23
655 Hill, C., Jones, M. G., Davies, D. E., & Wang, Y. (2019). Epithelial-mesenchymal 656 transition contributes to pulmonary fibrosis via aberrant epithelial/fibroblastic cross- 657 talk. Journal of Lung Health and Diseases, 3(2), 31–35.
658 Ito, T., Allen, R. M., Carson, W. F. 4th, Schaller, M., Cavassani, K. A., Hogaboam, C. 659 M., Lukacs, N. W., Matsukawa, A., & Kunkel, S. L. (2011). The critical role of Notch 660 ligand Delta-like 1 in the pathogenesis of influenza A virus (H1N1) infection. PLoS 661 Pathogens, 7(11), e1002341. https://doi.org/10.1371/journal.ppat.1002341
662 Jiang, R., Li, Y., Xu, Y., Zhou, Y., Pang, Y., Shen, L., Zhao, Y., Zhang, J., Zhou, J., 663 Wang, X., & Liu, Q. (2013). EMT and CSC-like properties mediated by the
29
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
664 IKKbeta/IkappaBalpha/RelA signal pathway via the transcriptional regulator, Snail, 665 are involved in the arsenite-induced neoplastic transformation of human 666 keratinocytes. Archives of Toxicology, 87(6), 991–1000. 667 https://doi.org/10.1007/s00204-012-0933-0
668 Jin, D., & Musier-Forsyth, K. (2019). Role of host tRNAs and aminoacyl-tRNA 669 synthetases in retroviral replication. The Journal of Biological Chemistry, 294(14), 670 5352–5364. https://doi.org/10.1074/jbc.REV118.002957
671 Jordan, N. V., Prat, A., Abell, A. N., Zawistowski, J. S., Sciaky, N., Karginova, O. A., 672 Zhou, B., Golitz, B. T., Perou, C. M., & Johnson, G. L. (2013). SWI/SNF chromatin- 673 remodeling factor Smarcd3/Baf60c controls epithelial-mesenchymal transition by 674 inducing Wnt5a signaling. Molecular and Cellular Biology, 33(15), 3011–3025. 675 https://doi.org/10.1128/MCB.01443-12
676 Josset, L., Menachery, V. D., Gralinski, L. E., Agnihothram, S., Sova, P., Carter, V. S., 677 Yount, B. L., Graham, R. L., Baric, R. S., & Katze, M. G. (2013). Cell host response 678 to infection with novel human coronavirus EMC predicts potential antivirals and 679 important differences with SARS coronavirus. MBio, 4(3), e00165-13. 680 https://doi.org/10.1128/mBio.00165-13
681 Kaldis, P., Ojala, P. M., Tong, L., Makela, T. P., & Solomon, M. J. (2001). CAK- 682 independent activation of CDK6 by a viral cyclin. Molecular Biology of the Cell, 683 12(12), 3987–3999. https://doi.org/10.1091/mbc.12.12.3987
684 Kato, N., Kosugi, T., Sato, W., Ishimoto, T., Kojima, H., Sato, Y., Sakamoto, K., 685 Maruyama, S., Yuzawa, Y., Matsuo, S., & Kadomatsu, K. (2011). Basigin/CD147 686 promotes renal fibrosis after unilateral ureteral obstruction. The American Journal 687 of Pathology, 178(2), 572–579. https://doi.org/10.1016/j.ajpath.2010.10.009
688 Klingelhutz, A. J., Foster, S. A., & McDougall, J. K. (1996). Telomerase activation by the 689 E6 gene product of human papillomavirus type 16. Nature, 380(6569), 79–82. 690 https://doi.org/10.1038/380079a0
691 Kolosionek, E., Savai, R., Ghofrani, H. A., Weissmann, N., Guenther, A., Grimminger, 692 F., Seeger, W., Banat, G. A., Schermuly, R. T., & Pullamsetti, S. S. (2009). 693 Expression and activity of phosphodiesterase isoforms during epithelial 694 mesenchymal transition: the role of phosphodiesterase 4. Molecular Biology of the 695 Cell, 20(22), 4751–4765. https://doi.org/10.1091/mbc.e09-01-0019
696 Kuo, W.-Y., Wu, C.-Y., Hwu, L., Lee, J.-S., Tsai, C.-H., Lin, K.-P., Wang, H.-E., Chou, 697 T.-Y., Tsai, C.-M., Gelovani, J., & Liu, R.-S. (2015). Enhancement of tumor initiation 698 and expression of KCNMA1, MORF4L2 and ASPM genes in the adenocarcinoma 699 of lung xenograft after vorinostat treatment. Oncotarget, 6(11), 8663–8675. 700 https://doi.org/10.18632/oncotarget.3536
701 Lamouille, S., Xu, J., & Derynck, R. (2014). Molecular mechanisms of epithelial- 702 mesenchymal transition. Nature Reviews. Molecular Cell Biology, 15(3), 178–196. 703 https://doi.org/10.1038/nrm3758
30
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
704 Li, H., Zhang, J., Song, X., Wang, T., Li, Z., Hao, D., Wang, X., Zheng, Q., Mao, C., Xu, 705 P., & Lv, C. (2014). Alveolar epithelial cells undergo epithelial-mesenchymal 706 transition in acute interstitial pneumonia: a case report. BMC Pulmonary Medicine, 707 14, 67. https://doi.org/10.1186/1471-2466-14-67
708 Li, N., Gou, J.-H., Xiong, J., You, J.-J., & Li, Z.-Y. (2020). HOXB4 promotes the 709 malignant progression of ovarian cancer via DHDDS. BMC Cancer, 20(1), 222. 710 https://doi.org/10.1186/s12885-020-06725-4
711 Liu, T., Zhao, M., Liu, J., He, Z., Zhang, Y., You, H., Huang, J., Lin, X., & Feng, X.-H. 712 (2017). Tumor suppressor bromodomain-containing protein 7 cooperates with 713 Smads to promote transforming growth factor-beta responses. Oncogene, 36(3), 714 362–372. https://doi.org/10.1038/onc.2016.204
715 Liu, X., Fu, Y., Huang, J., Wu, M., Zhang, Z., Xu, R., Zhang, P., Zhao, S., Liu, L., & 716 Jiang, H. (2019). ADAR1 promotes the epithelial-to-mesenchymal transition and 717 stem-like cell phenotype of oral cancer by facilitating oncogenic microRNA 718 maturation. Journal of Experimental & Clinical Cancer Research : CR, 38(1), 315. 719 https://doi.org/10.1186/s13046-019-1300-2
720 Liu, Y. (2006). Renal fibrosis: new insights into the pathogenesis and therapeutics. 721 Kidney International, 69(2), 213–217. https://doi.org/10.1038/sj.ki.5000054
722 Liu, Z., Li, Q., Li, K., Chen, L., Li, W., Hou, M., Liu, T., Yang, J., Lindvall, C., Bjorkholm, 723 M., Jia, J., & Xu, D. (2013). Telomerase reverse transcriptase promotes epithelial- 724 mesenchymal transition and stem cell-like traits in cancer cells. Oncogene, 32(36), 725 4203–4213. https://doi.org/10.1038/onc.2012.441
726 Ludwig, S., & Planz, O. (2008). Influenza viruses and the NF-kappaB signaling pathway 727 - towards a novel concept of antiviral therapy. Biological Chemistry, 389(10), 1307– 728 1312. https://doi.org/10.1515/BC.2008.148
729 Luo, X., Guo, L., Zhang, J., Xu, Y., Gu, W., Feng, L., & Wang, Y. (2017). Tight Junction 730 Protein Occludin Is a Porcine Epidemic Diarrhea Virus Entry Factor. Journal of 731 Virology, 91(10). https://doi.org/10.1128/JVI.00202-17
732 Martinelli, P., Carrillo-de Santa Pau, E., Cox, T., Sainz, B. J., Dusetti, N., Greenhalf, W., 733 Rinaldi, L., Costello, E., Ghaneh, P., Malats, N., Buchler, M., Pajic, M., Biankin, A. 734 V, Iovanna, J., Neoptolemos, J., & Real, F. X. (2017). GATA6 regulates EMT and 735 tumour dissemination, and is a marker of response to adjuvant chemotherapy in 736 pancreatic cancer. Gut, 66(9), 1665–1676. https://doi.org/10.1136/gutjnl-2015- 737 311256
738 Mendoza-Rodriguez, M., Arevalo Romero, H., Fuentes-Panana, E. M., Ayala-Sumuano, 739 J.-T., & Meza, I. (2017). IL-1beta induces up-regulation of BIRC3, a gene involved 740 in chemoresistance to doxorubicin in breast cancer cells. Cancer Letters, 390, 39– 741 44. https://doi.org/10.1016/j.canlet.2017.01.005
742 Ng, S. S. M., Li, A., Pavlakis, G. N., Ozato, K., & Kino, T. (2013). Viral infection
31
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
743 increases glucocorticoid-induced interleukin-10 production through ERK-mediated 744 phosphorylation of the glucocorticoid receptor in dendritic cells: potential clinical 745 implications. PloS One, 8(5), e63587. https://doi.org/10.1371/journal.pone.0063587
746 Nie, Z., Gao, W., Zhang, Y., Hou, Y., Liu, J., Li, Z., Xue, W., Ye, X., & Jin, A. (2019). 747 STAG2 loss-of-function mutation induces PD-L1 expression in U2OS cells. Annals 748 of Translational Medicine, 7(7), 127. https://doi.org/10.21037/atm.2019.02.23
749 Nieto, M. A. (2002). The snail superfamily of zinc-finger transcription factors. Nature 750 Reviews. Molecular Cell Biology, 3(3), 155–166. https://doi.org/10.1038/nrm757
751 Ochsner, S. A., Abraham, D., Martin, K., Ding, W., McOwiti, A., Kankanamge, W., 752 Wang, Z., Andreano, K., Hamilton, R. A., Chen, Y., Hamilton, A., Gantner, M. L., 753 Dehart, M., Qu, S., Hilsenbeck, S. G., Becnel, L. B., Bridges, D., Ma’ayan, A., 754 Huss, J. M., McKenna, N. J. (2019). The Signaling Pathways Project, an 755 integrated ’omics knowledgebase for mammalian cellular signaling pathways. 756 Scientific Data, 6(1), 252. https://doi.org/10.1038/s41597-019-0193-4
757 Ostler, J. B., Harrison, K. S., Schroeder, K., Thunuguntla, P., & Jones, C. (2019). The 758 Glucocorticoid Receptor (GR) Stimulates Herpes Simplex Virus 1 Productive 759 Infection, in Part Because the Infected Cell Protein 0 (ICP0) Promoter Is 760 Cooperatively Transactivated by the GR and Kruppel-Like Transcription Factor 15. 761 Journal of Virology, 93(6). https://doi.org/10.1128/JVI.02063-18
762 Pauls, E., Ruiz, A., Badia, R., Permanyer, M., Gubern, A., Riveira-Munoz, E., Torres- 763 Torronteras, J., Alvarez, M., Mothe, B., Brander, C., Crespo, M., Menendez-Arias, 764 L., Clotet, B., Keppler, O. T., Marti, R., Posas, F., Ballana, E., & Este, J. A. (2014). 765 Cell cycle control and HIV-1 susceptibility are linked by CDK6-dependent CDK2 766 phosphorylation of SAMHD1 in myeloid and lymphoid cells. Journal of Immunology 767 (Baltimore, Md. : 1950), 193(4), 1988–1997. 768 https://doi.org/10.4049/jimmunol.1400873
769 Pertea, M., Shumate, A., Pertea, G., Varabyou, A., Breitwieser, F. P., Chang, Y.-C., 770 Madugundu, A. K., Pandey, A., & Salzberg, S. L. (2018). CHESS: a new human 771 gene catalog curated from thousands of large-scale RNA sequencing experiments 772 reveals extensive transcriptional noise. Genome Biology, 19(1), 208. 773 https://doi.org/10.1186/s13059-018-1590-2
774 Poppe, M., Wittig, S., Jurida, L., Bartkuhn, M., Wilhelm, J., Muller, H., Beuerlein, K., 775 Karl, N., Bhuju, S., Ziebuhr, J., Schmitz, M. L., & Kracht, M. (2017). The NF- 776 kappaB-dependent and -independent transcriptome and chromatin landscapes of 777 human coronavirus 229E-infected cells. PLoS Pathogens, 13(3), e1006286. 778 https://doi.org/10.1371/journal.ppat.1006286
779 Ruckle, A., Haasbach, E., Julkunen, I., Planz, O., Ehrhardt, C., & Ludwig, S. (2012). 780 The NS1 protein of influenza A virus blocks RIG-I-mediated activation of the 781 noncanonical NF-kappaB pathway and p52/RelB-dependent gene expression in 782 lung epithelial cells. Journal of Virology, 86(18), 10211–10217. 783 https://doi.org/10.1128/JVI.00323-12
32
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
784 Ruster, C., & Wolf, G. (2011). Angiotensin II as a morphogenic cytokine stimulating 785 renal fibrogenesis. Journal of the American Society of Nephrology : JASN, 22(7), 786 1189–1199. https://doi.org/10.1681/ASN.2010040384
787 Sahu, S. K., Tiwari, N., Pataskar, A., Zhuang, Y., Borisova, M., Diken, M., Strand, S., 788 Beli, P., & Tiwari, V. K. (2017). FBXO32 promotes microenvironment underlying 789 epithelial-mesenchymal transition via CtBP1 during tumour metastasis and brain 790 development. Nature Communications, 8(1), 1523. https://doi.org/10.1038/s41467- 791 017-01366-x
792 Sappenfield, E., Jamieson, D. J., & Kourtis, A. P. (2013). Pregnancy and susceptibility 793 to infectious diseases. Infectious Diseases in Obstetrics and Gynecology, 2013, 794 752852. https://doi.org/10.1155/2013/752852
795 Schneider, W. M., Chevillotte, M. D., & Rice, C. M. (2014). Interferon-stimulated genes: 796 a complex web of host defenses. Annual Review of Immunology, 32, 513–545. 797 https://doi.org/10.1146/annurev-immunol-032713-120231
798 Schulze-Gahmen, U., Upton, H., Birnberg, A., Bao, K., Chou, S., Krogan, N. J., Zhou, 799 Q., & Alber, T. (2013). The AFF4 scaffold binds human P-TEFb adjacent to HIV 800 Tat. ELife, 2, e00327. https://doi.org/10.7554/eLife.00327
801 Sims, A. C., Tilton, S. C., Menachery, V. D., Gralinski, L. E., Schafer, A., Matzke, M. M., 802 Webb-Robertson, B.-J. M., Chang, J., Luna, M. L., Long, C. E., Shukla, A. K., 803 Bankhead, A. R. 3rd, Burkett, S. E., Zornetzer, G., Tseng, C.-T. K., Metz, T. O., 804 Pickles, R., McWeeney, S., Smith, R. D., Baric, R. S. (2013). Release of severe 805 acute respiratory syndrome coronavirus nuclear import block enhances host 806 transcription in human lung cells. Journal of Virology, 87(7), 3885–3902. 807 https://doi.org/10.1128/JVI.02520-12
808 Siraj, A. K., Pratheeshkumar, P., Divya, S. P., Parvathareddy, S. K., Bu, R., Masoodi, 809 T., Kong, Y., Thangavel, S., Al-Sanea, N., Ashari, L. H., Abduljabbar, A., Al- 810 Homoud, S., Al-Dayel, F., & Al-Kuraya, K. S. (2019). TGFbeta-induced SMAD4- 811 dependent Apoptosis Proceeded by EMT in CRC. Molecular Cancer Therapeutics, 812 18(7), 1312–1322. https://doi.org/10.1158/1535-7163.MCT-18-1378
813 Siston, A. M., Rasmussen, S. A., Honein, M. A., Fry, A. M., Seib, K., Callaghan, W. M., 814 Louie, J., Doyle, T. J., Crockett, M., Lynfield, R., Moore, Z., Wiedeman, C., Anand, 815 M., Tabony, L., Nielsen, C. F., Waller, K., Page, S., Thompson, J. M., Avery, C., 816 Jamieson, D. J. (2010). Pandemic 2009 influenza A(H1N1) virus illness among 817 pregnant women in the United States. JAMA, 303(15), 1517–1525. 818 https://doi.org/10.1001/jama.2010.479
819 Stark, G. R., Kerr, I. M., Williams, B. R., Silverman, R. H., & Schreiber, R. D. (1998). 820 How cells respond to interferons. Annual Review of Biochemistry, 67, 227–264. 821 https://doi.org/10.1146/annurev.biochem.67.1.227
822 Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., 823 Paulovich, A., Pomeroy, S. L., Golub, T. R., Lander, E. S., & Mesirov, J. P. (2005).
33
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
824 Gene set enrichment analysis: a knowledge-based approach for interpreting 825 genome-wide expression profiles. Proceedings of the National Academy of 826 Sciences of the United States of America, 102(43), 15545–15550. 827 https://doi.org/10.1073/pnas.0506580102
828 Suh, Y., Yoon, C.-H., Kim, R.-K., Lim, E.-J., Oh, Y. S., Hwang, S.-G., An, S., Yoon, G., 829 Gye, M. C., Yi, J.-M., Kim, M.-J., & Lee, S.-J. (2013). Claudin-1 induces epithelial- 830 mesenchymal transition through activation of the c-Abl-ERK signaling pathway in 831 human liver cells. Oncogene, 32(41), 4873–4882. 832 https://doi.org/10.1038/onc.2012.505
833 Takahashi, K., Halfmann, P., Oyama, M., Kozuka-Hata, H., Noda, T., & Kawaoka, Y. 834 (2013). DNA topoisomerase 1 facilitates the transcription and replication of the 835 Ebola virus genome. Journal of Virology, 87(16), 8862–8869. 836 https://doi.org/10.1128/JVI.03544-12
837 Takeda, K., Kaisho, T., & Akira, S. (2003). Toll-like receptors. Annual Review of 838 Immunology, 21, 335–376. 839 https://doi.org/10.1146/annurev.immunol.21.120601.141126
840 Taki, M., Abiko, K., Baba, T., Hamanishi, J., Yamaguchi, K., Murakami, R., Yamanoi, K., 841 Horikawa, N., Hosoe, Y., Nakamura, E., Sugiyama, A., Mandai, M., Konishi, I., & 842 Matsumura, N. (2018). Snail promotes ovarian cancer progression by recruiting 843 myeloid-derived suppressor cells via CXCR2 ligand upregulation. Nature 844 Communications, 9(1), 1685. https://doi.org/10.1038/s41467-018-03966-7
845 Tecalco-Cruz, A. C., Rios-Lopez, D. G., Vazquez-Victorio, G., Rosales-Alvarez, R. E., & 846 Macias-Silva, M. (2018). Transcriptional cofactors Ski and SnoN are major 847 regulators of the TGF-beta/Smad signaling pathway in health and disease. Signal 848 Transduction and Targeted Therapy, 3, 15. https://doi.org/10.1038/s41392-018- 849 0015-8
850 Totura, A. L., & Bavari, S. (2019). Broad-spectrum coronavirus antiviral drug discovery. 851 Expert Opinion on Drug Discovery, 14(4), 397–412. 852 https://doi.org/10.1080/17460441.2019.1581171
853 Totura, A. L., Whitmore, A., Agnihothram, S., Schafer, A., Katze, M. G., Heise, M. T., & 854 Baric, R. S. (2015). Toll-Like Receptor 3 Signaling via TRIF Contributes to a 855 Protective Innate Immune Response to Severe Acute Respiratory Syndrome 856 Coronavirus Infection. MBio, 6(3), e00638-15. https://doi.org/10.1128/mBio.00638- 857 15
858 Vukmirovic, M., Herazo-Maya, J. D., Blackmon, J., Skodric-Trifunovic, V., Jovanovic, D., 859 Pavlovic, S., Stojsic, J., Zeljkovic, V., Yan, X., Homer, R., Stefanovic, B., & 860 Kaminski, N. (2017). Identification and validation of differentially expressed 861 transcripts by RNA-sequencing of formalin-fixed, paraffin-embedded (FFPE) lung 862 tissue from patients with Idiopathic Pulmonary Fibrosis. BMC Pulmonary Medicine, 863 17(1), 15. https://doi.org/10.1186/s12890-016-0356-4
34
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
864 Wan, L. C. K., Maisonneuve, P., Szilard, R. K., Lambert, J.-P., Ng, T. F., Manczyk, N., 865 Huang, H., Laister, R., Caudy, A. A., Gingras, A.-C., Durocher, D., & Sicheri, F. 866 (2017). Proteomic analysis of the human KEOPS complex identifies C14ORF142 867 as a core subunit homologous to yeast Gon7. Nucleic Acids Research, 45(2), 805– 868 817. https://doi.org/10.1093/nar/gkw1181
869 Wang, C.-A., Drasin, D., Pham, C., Jedlicka, P., Zaberezhnyy, V., Guney, M., Li, H., 870 Nemenoff, R., Costello, J. C., Tan, A.-C., & Ford, H. L. (2014). Homeoprotein Six2 871 promotes breast cancer metastasis via transcriptional and epigenetic control of E- 872 cadherin expression. Cancer Research, 74(24), 7357–7370. 873 https://doi.org/10.1158/0008-5472.CAN-14-0666
874 Wang, C., Zhang, J., Fok, K. L., Tsang, L. L., Ye, M., Liu, J., Li, F., Zhao, A. Z., Chan, 875 H. C., & Chen, H. (2018). CD147 Induces Epithelial-to-Mesenchymal Transition by 876 Disassembling Cellular Apoptosis Susceptibility Protein/E-Cadherin/beta-Catenin 877 Complex in Human Endometriosis. The American Journal of Pathology, 188(7), 878 1597–1607. https://doi.org/10.1016/j.ajpath.2018.03.004
879 Wang, K., & Li, J. (2016). Overexpression of ANXA3 is an independent prognostic 880 indicator in gastric cancer and its depletion suppresses cell proliferation and tumor 881 growth. Oncotarget, 7(52), 86972–86984. 882 https://doi.org/10.18632/oncotarget.13493
883 Wang, S., Lei, T., Zhang, K., Zhao, W., Fang, L., Lai, B., Han, J., Xiao, L., & Wang, N. 884 (2014). Xenobiotic pregnane X receptor (PXR) regulates innate immunity via 885 activation of NLRP3 inflammasome in vascular endothelial cells. The Journal of 886 Biological Chemistry, 289(43), 30075–30081. 887 https://doi.org/10.1074/jbc.M114.578781
888 Wang, W., Ye, L., Ye, L., Li, B., Gao, B., Zeng, Y., Kong, L., Fang, X., Zheng, H., Wu, 889 Z., & She, Y. (2007). Up-regulation of IL-6 and TNF-alpha induced by SARS- 890 coronavirus spike protein in murine macrophages via NF-kappaB pathway. Virus 891 Research, 128(1–2), 1–8. https://doi.org/10.1016/j.virusres.2007.02.007
892 Wei, C.-Y., Zhu, M.-X., Yang, Y.-W., Zhang, P.-F., Yang, X., Peng, R., Gao, C., Lu, J.- 893 C., Wang, L., Deng, X.-Y., Lu, N.-H., Qi, F.-Z., & Gu, J.-Y. (2019). Downregulation 894 of RNF128 activates Wnt/beta-catenin signaling to induce cellular EMT and 895 stemness via CD44 and CTTN ubiquitination in melanoma. Journal of Hematology 896 & Oncology, 12(1), 21. https://doi.org/10.1186/s13045-019-0711-z
897 Wobbe, C. R., Dean, F. B., Murakami, Y., Borowiec, J. A., Bullock, P., & Hurwitz, J. 898 (1987). In vitro replication of DNA containing either the SV40 or the polyoma origin. 899 Philosophical Transactions of the Royal Society of London. Series B, Biological 900 Sciences, 317(1187), 439–453. https://doi.org/10.1098/rstb.1987.0071
901 Xia, L., Dai, L., Yu, Q., & Yang, Q. (2017). Persistent Transmissible Gastroenteritis 902 Virus Infection Enhances Enterotoxigenic Escherichia coli K88 Adhesion by 903 Promoting Epithelial-Mesenchymal Transition in Intestinal Epithelial Cells. Journal 904 of Virology, 91(21). https://doi.org/10.1128/JVI.01256-17
35
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
905 Xu, C.-B., Liu, X.-S.-B.-J., Li, J.-Q., Zhao, X., Xin, D., & Yu, D. (2019). microRNA-539 906 functions as a tumor suppressor in papillary thyroid carcinoma via the transforming 907 growth factor beta1/Smads signaling pathway by targeting secretory leukocyte 908 protease inhibitor. Journal of Cellular Biochemistry, 120(6), 10830–10846. 909 https://doi.org/10.1002/jcb.28374
910 Yamajuku, D., Shibata, Y., Kitazawa, M., Katakura, T., Urata, H., Kojima, T., Nakata, O., 911 & Hashimoto, S. (2010). Identification of functional clock-controlled elements 912 involved in differential timing of Per1 and Per2 transcription. Nucleic Acids 913 Research, 38(22), 7964–7973. https://doi.org/10.1093/nar/gkq678
914 Yang, J., Deng, X., Deng, L., Gu, H., Fan, W., & Cao, Y. (2004). Telomerase activation 915 by Epstein-Barr virus latent membrane protein 1 is associated with c-Myc 916 expression in human nasopharyngeal epithelial cells. Journal of Experimental & 917 Clinical Cancer Research : CR, 23(3), 495–506.
918 Yin, L., Hubbard, A. K., & Giardina, C. (2000). NF-kappa B regulates transcription of the 919 mouse telomerase catalytic subunit. The Journal of Biological Chemistry, 275(47), 920 36671–36675. https://doi.org/10.1074/jbc.M007378200
921 Yoshikawa, T., Hill, T. E., Yoshikawa, N., Popov, V. L., Galindo, C. L., Garner, H. R., 922 Peters, C. J., & Tseng, C.-T. K. (2010). Dynamic innate immune responses of 923 human bronchial epithelial cells to severe acute respiratory syndrome-associated 924 coronavirus infection. PloS One, 5(1), e8729. 925 https://doi.org/10.1371/journal.pone.0008729
926 Zaborowska, J., Isa, N. F., & Murphy, S. (2016). P-TEFb goes viral. BioEssays : News 927 and Reviews in Molecular, Cellular and Developmental Biology, 38 Suppl 1, S75- 928 85. https://doi.org/10.1002/bies.201670912
929 Zhang, L., Wang, X., Lai, C., Zhang, H., & Lai, M. (2019). PMEPA1 induces EMT via a 930 non-canonical TGF-beta signalling in colorectal cancer. Journal of Cellular and 931 Molecular Medicine, 23(5), 3603–3615. https://doi.org/10.1111/jcmm.14261
932 Zhao, M., Kong, L., Liu, Y., & Qu, H. (2015). dbEMT: an epithelial-mesenchymal 933 transition associated gene resource. Scientific Reports, 5, 11459. 934 https://doi.org/10.1038/srep11459
935 Zheng, K., Kitazato, K., & Wang, Y. (2014). Viruses exploit the function of epidermal 936 growth factor receptor. Reviews in Medical Virology, 24(4), 274–286. 937 https://doi.org/10.1002/rmv.1796
938 Zhou, F., Geng, J., Xu, S., Meng, Q., Chen, K., Liu, F., Yang, F., Pan, B., & Yu, Y. 939 (2019). FAM83A signaling induces epithelial-mesenchymal transition by the 940 PI3K/AKT/Snail pathway in NSCLC. Aging, 11(16), 6069–6088. 941 https://doi.org/10.18632/aging.102163
942 Zhou, P., Yang, X.-L., Wang, X.-G., Hu, B., Zhang, L., Zhang, W., Si, H.-R., Zhu, Y., Li, 943 B., Huang, C.-L., Chen, H.-D., Chen, J., Luo, Y., Guo, H., Jiang, R.-D., Liu, M.-Q.,
36
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
944 Chen, Y., Shen, X.-R., Wang, X., Shi, Z.-L. (2020). A pneumonia outbreak 945 associated with a new coronavirus of probable bat origin. Nature, 579(7798), 270– 946 273. https://doi.org/10.1038/s41586-020-2012-7
947
37
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
948
38
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
949 Figure Titles Legends
950 Figure 1. Ranking of ISGs in the human ALL CoV infection transcriptomic
951 consensome. See Supplementary information, Section 2, column I.
952 Figure 2. GeneOverlap analysis of CoV and IAV TC95s and SPP receptor or
953 enzyme TC95s. OR: Odds ratio. All ORs are p < 0.05. Numerical data are provided in
954 Supplementary information Section 8.
955 Figure 3. GeneOverlap analysis of CoV and IAV TC95s and ChIP-Atlas
956 transcription factor CC95s. OR: Odds ratio. All ORs are p < 0.05. Numerical data are
957 provided in Supplementary information Section 8.
958 Figure 4. GeneOverlap analysis of CoV and IAV TC95s and ChIP-Atlas enzyme
959 CC95s. OR: Odds ratio. Numerical data are provided in Supplementary information
960 Section 8. All ORs are p < 0.05.
961 Figure 5. GeneOverlap analysis of CoV and IAV TC95s and ChIP-Atlas co-node
962 CC95s. OR: Odds ratio. Numerical data are provided in Supplementary information
963 Section 8. All ORs are p < 0.05.
964 Figure 6. Antagonism between PGR and IFNR signaling in the regulation of viral
965 response genes in the airway epithelium. PGR loss of function experiments were
966 retrieved from GSE17307 (SPP DOI 10.1621/ xigKzGn1se).
967 Figure 7. Conservation of polarity of differential expression in CoV, IFNR and
968 TERT perturbation transcriptomic experiments across diverse canonical
969 interferon-inducible viral response genes. Positive/negative regulatory relationship
970 of TERT with a transcriptional target was inferred from the design of the underlying
39
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
971 experiments: this was loss of function for experiments in GSE77014 (MST132 inhibitor)
972 and GSE60175 (shRNA), and gain of function in E-MEXP-563 (overexpression). All
973 IFNR experiments were gain of function using characterized physiological ligands for
974 members of the family.
975 Figure 8. Identification of a SARS2 epithelial to mesenchymal transcriptional
976 signature. Note that the SARS2 genes are all in the 99th percentile and are therefore
977 superimposed in the scatterplot. Indicated are the results of the Paired Two Sample for
978 Means t-Test comparing the Relative Ranks of the genes in the SARS2 consensome
979 with those in the SARS1, MERS and IAV consensomes.
980
40
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
981 Tables
982 Table 1. Genes in the 99th %ile of the ALL CoV consensome and the 50th
983 percentile of the IAV consensome
984
985
986
41
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
987 Table 2. Links to SPP Regulation Reports for the top 20 ranked genes in the ALL
988 CoV consensome
Symbol Name Link JUN Jun proto-oncogene, AP-1 transcription factor subunit SPP Regulation Report DUSP1 dual specificity phosphatase 1 SPP Regulation Report FOS Fos proto-oncogene, AP-1 transcription factor subunit SPP Regulation Report EGR1 early growth response 1 SPP Regulation Report CXCL2 C-X-C motif chemokine ligand 2 SPP Regulation Report TNFAIP3 TNF alpha induced protein 3 SPP Regulation Report IER2 immediate early response 2 SPP Regulation Report CCNL1 cyclin L1 SPP Regulation Report PIM3 Pim-3 proto-oncogene, serine/threonine kinase SPP Regulation Report CSRNP1 cysteine and serine rich nuclear protein 1 SPP Regulation Report FOSB FosB proto-oncogene, AP-1 transcription factor subunit SPP Regulation Report ZC3H12A zinc finger CCCH-type containing 12A SPP Regulation Report ZFP36 ZFP36 ring finger protein SPP Regulation Report BHLHE40 basic helix-loop-helix family member e40 SPP Regulation Report RELB RELB proto-oncogene, NF-kB subunit SPP Regulation Report EGR2 early growth response 2 SPP Regulation Report NFKBIA NFKB inhibitor alpha SPP Regulation Report BCL3 BCL3 transcription coactivator SPP Regulation Report YRDC yrdC N6-threonylcarbamoyltransferase domain containing SPP Regulation Report IFIT1 interferon induced protein with tetratricopeptide repeats 1 SPP Regulation Report 989
42
bioRxiv preprint was notcertifiedbypeerreview)istheauthor/funder,whohasgrantedbioRxivalicensetodisplaypreprintinperpetuity.Itmade
doi: https://doi.org/10.1101/2020.04.24.059527
available undera ; CC-BY 4.0Internationallicense this versionpostedMay14,2020.
. Figure 1 The copyrightholderforthispreprint(which
bioRxiv preprint was notcertifiedbypeerreview)istheauthor/funder,whohasgrantedbioRxivalicensetodisplaypreprintinperpetuity.Itmade
doi: https://doi.org/10.1101/2020.04.24.059527 Figure 2 available undera ; CC-BY 4.0Internationallicense this versionpostedMay14,2020. . The copyrightholderforthispreprint(which
bioRxiv preprint was notcertifiedbypeerreview)istheauthor/funder,whohasgrantedbioRxivalicensetodisplaypreprintinperpetuity.Itmade
doi: https://doi.org/10.1101/2020.04.24.059527
available undera ; CC-BY 4.0Internationallicense this versionpostedMay14,2020. . The copyrightholderforthispreprint(which
Figure 3 bioRxiv preprint was notcertifiedbypeerreview)istheauthor/funder,whohasgrantedbioRxivalicensetodisplaypreprintinperpetuity.Itmade doi: https://doi.org/10.1101/2020.04.24.059527 available undera ; CC-BY 4.0Internationallicense this versionpostedMay14,2020. . The copyrightholderforthispreprint(which
Figure 3 (contd) bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
Figure 4 bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
Figure 5 bioRxiv preprint was notcertifiedbypeerreview)istheauthor/funder,whohasgrantedbioRxivalicensetodisplaypreprintinperpetuity.Itmade Figure 6
doi: https://doi.org/10.1101/2020.04.24.059527
available undera ; CC-BY 4.0Internationallicense this versionpostedMay14,2020.
. The copyrightholderforthispreprint(which
bioRxiv preprint doi: https://doi.org/10.1101/2020.04.24.059527; this version posted May 14, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
log10 log10 FC
log10 log10 FC
Figure 7 bioRxiv preprint was notcertifiedbypeerreview)istheauthor/funder,whohasgrantedbioRxivalicensetodisplaypreprintinperpetuity.Itmade
Figure 8 doi: https://doi.org/10.1101/2020.04.24.059527
available undera ; CC-BY 4.0Internationallicense this versionpostedMay14,2020.
. The copyrightholderforthispreprint(which