A single dose of cocaine rewires the 3D genome structure of midbrain dopamine neurons.

Midbrain dopamine neurons (DNs) respond to a first exposure to addictive drugs and play key roles in chronic drug usage 1-3 . As the synaptic and transcriptional changes that follow an acute cocaine exposure are mostly resolved within a few days 4,5 , the molecular changes that encode the long-term cellular memory of the exposure within DNs remain unknown. To investigate whether a single cocaine exposure induces long-term changes in the 3D genome structure of DNs, we applied Genome Architecture Mapping and single nucleus transcriptomic analyses in the mouse midbrain. We found extensive rewiring of 3D genome architecture at 24 hours past exposure which remains or worsens by 14 days, outlasting transcriptional responses. The cocaine-induced chromatin rewiring occurs at all genomic scales and affects genes with major roles in cocaine-induced synaptic changes. A single cocaine exposure triggers extensive long-lasting changes in chromatin condensation in post-synaptic and post-transcriptional regulatory genes, for example the unfolding of Rbfox1 which becomes most prominent 14 days post exposure. Finally, structurally remodeled genes are most expressed in a specific DN sub-type characterized by low expression of the dopamine auto-receptor Drd2 , a key feature of highly cocaine-sensitive cells. These results reveal an important role for long-lasting 3D genome remodelling in the cellular memory of a single cocaine exposure, providing new hypotheses for understanding the inception of drug addiction and 3D genome plasticity.

Pathogenic mutations of human phosphorylation sites affect protein-protein interactions.

Despite their lack of a defined 3D structure, intrinsically disordered regions (IDRs) of proteins play important biological roles. Many IDRs contain short linear motifs (SLiMs) that mediate protein-protein interactions (PPIs), which can be regulated by post-translational modifications like phosphorylation. 20% of pathogenic missense mutations are found in IDRs, and understanding how such mutations affect PPIs is essential for unraveling disease mechanisms. Here, we employ peptide-based interaction proteomics to investigate 36 disease-associated mutations affecting phosphorylation sites. Our results unveil significant differences in interactomes between phosphorylated and non-phosphorylated peptides, often due to disrupted phosphorylation-dependent SLiMs. We focused on a mutation of a serine phosphorylation site in the transcription factor GATAD1, which causes dilated cardiomyopathy. We find that this phosphorylation site mediates interaction with 14-3-3 family proteins. Follow-up experiments reveal the structural basis of this interaction and suggest that 14-3-3 binding affects GATAD1 nucleocytoplasmic transport by masking a nuclear localisation signal. Our results demonstrate that pathogenic mutations of human phosphorylation sites can significantly impact protein-protein interactions, offering insights into potential molecular mechanisms underlying pathogenesis.

The SPOC proteins DIDO3 and PHF3 co-regulate gene expression and neuronal differentiation.

Transcription is regulated by a multitude of activators and repressors, which bind to the RNA polymerase II (Pol II) machinery and modulate its progression. Death-inducer obliterator 3 (DIDO3) and PHD finger protein 3 (PHF3) are paralogue proteins that regulate transcription elongation by docking onto phosphorylated serine-2 in the C-terminal domain (CTD) of Pol II through their SPOC domains. Here, we show that DIDO3 and PHF3 form a complex that bridges the Pol II elongation machinery with chromatin and RNA processing factors and tethers Pol II in a phase-separated microenvironment. Their SPOC domains and C-terminal intrinsically disordered regions are critical for transcription regulation. PHF3 and DIDO exert cooperative and antagonistic effects on the expression of neuronal genes and are both essential for neuronal differentiation. In the absence of PHF3, DIDO3 is upregulated as a compensatory mechanism. In addition to shared gene targets, DIDO specifically regulates genes required for lipid metabolism. Collectively, our work reveals multiple layers of gene expression regulation by the DIDO3 and PHF3 paralogues, which have specific, co-regulatory and redundant functions in transcription.

Transcriptional reprogramming by mutated IRF4 in lymphoma.

Disease-causing mutations in genes encoding transcription factors (TFs) can affect TF interactions with their cognate DNA-binding motifs. Whether and how TF mutations impact upon the binding to TF composite elements (CE) and the interaction with other TFs is unclear. Here, we report a distinct mechanism of TF alteration in human lymphomas with perturbed B cell identity, in particular classic Hodgkin lymphoma. It is caused by a recurrent somatic missense mutation c.295 T > C (p.Cys99Arg; p.C99R) targeting the center of the DNA-binding domain of Interferon Regulatory Factor 4 (IRF4), a key TF in immune cells. IRF4-C99R fundamentally alters IRF4 DNA-binding, with loss-of-binding to canonical IRF motifs and neomorphic gain-of-binding to canonical and non-canonical IRF CEs. IRF4-C99R thoroughly modifies IRF4 function by blocking IRF4-dependent plasma cell induction, and up-regulates disease-specific genes in a non-canonical Activator Protein-1 (AP-1)-IRF-CE (AICE)-dependent manner. Our data explain how a single mutation causes a complex switch of TF specificity and gene regulation and open the perspective to specifically block the neomorphic DNA-binding activities of a mutant TF.

Cardiovascular disease biomarkers derived from circulating cell-free DNA methylation.

Acute coronary syndrome (ACS) remains a major cause of worldwide mortality. The syndrome occurs when blood flow to the heart muscle is decreased or blocked, causing muscle tissues to die or malfunction. There are three main types of ACS: Non-ST-elevation myocardial infarction, ST-elevation myocardial infarction, and unstable angina. The treatment depends on the type of ACS, and this is decided by a combination of clinical findings, such as electrocardiogram and plasma biomarkers. Circulating cell-free DNA (ccfDNA) is proposed as an additional marker for ACS since the damaged tissues can release DNA to the bloodstream. We used ccfDNA methylation profiles for differentiating between the ACS types and provided computational tools to repeat similar analysis for other diseases. We leveraged cell type specificity of DNA methylation to deconvolute the ccfDNA cell types of origin and to find methylation-based biomarkers that stratify patients. We identified hundreds of methylation markers associated with ACS types and validated them in an independent cohort. Many such markers were associated with genes involved in cardiovascular conditions and inflammation. ccfDNA methylation showed promise as a non-invasive diagnostic for acute coronary events. These methods are not limited to acute events, and may be used for chronic cardiovascular diseases as well.

Disintegration of the NuRD Complex in Primary Human Muscle Stem Cells in Critical Illness Myopathy.

Critical illness myopathy (CIM) is an acquired, devastating, multifactorial muscle-wasting disease with incomplete recovery. The impact on hospital costs and permanent loss of quality of life is enormous. Incomplete recovery might imply that the function of muscle stem cells (MuSC) is impaired. We tested whether epigenetic alterations could be in part responsible. We characterized human muscle stem cells (MuSC) isolated from early CIM and analyzed epigenetic alterations (CIM n = 15, controls n = 21) by RNA-Seq, immunofluorescence, analysis of DNA repair, and ATAC-Seq. CIM-MuSC were transplanted into immunodeficient NOG mice to assess their regenerative potential. CIM-MuSC exhibited significant growth deficits, reduced ability to differentiate into myotubes, and impaired DNA repair. The chromatin structure was damaged, as characterized by alterations in mRNA of histone 1, depletion or dislocation of core proteins of nucleosome remodeling and deacetylase complex, and loosening of multiple nucleosome-spanning sites. Functionally, CIM-MuSC had a defect in building new muscle fibers. Further, MuSC obtained from the electrically stimulated muscle of CIM patients was very similar to control MuSC, indicating the impact of muscle contraction in the onset of CIM. CIM not only affects working skeletal muscle but has a lasting and severe epigenetic impact on MuSC.

A multimorphic mutation in IRF4 causes human autosomal dominant combined immunodeficiency.

Interferon regulatory factor 4 (IRF4) is a transcription factor (TF) and key regulator of immune cell development and function. We report a recurrent heterozygous mutation in IRF4, p.T95R, causing an autosomal dominant combined immunodeficiency (CID) in seven patients from six unrelated families. The patients exhibited profound susceptibility to opportunistic infections, notably Pneumocystis jirovecii, and presented with agammaglobulinemia. Patients' B cells showed impaired maturation, decreased immunoglobulin isotype switching, and defective plasma cell differentiation, whereas their T cells contained reduced TH17 and TFH populations and exhibited decreased cytokine production. A knock-in mouse model of heterozygous T95R showed a severe defect in antibody production both at the steady state and after immunization with different types of antigens, consistent with the CID observed in these patients. The IRF4T95R variant maps to the TF's DNA binding domain, alters its canonical DNA binding specificities, and results in a simultaneous multimorphic combination of loss, gain, and new functions for IRF4. IRF4T95R behaved as a gain-of-function hypermorph by binding to DNA with higher affinity than IRF4WT. Despite this increased affinity for DNA, the transcriptional activity on IRF4 canonical genes was reduced, showcasing a hypomorphic activity of IRF4T95R. Simultaneously, IRF4T95R functions as a neomorph by binding to noncanonical DNA sites to alter the gene expression profile, including the transcription of genes exclusively induced by IRF4T95R but not by IRF4WT. This previously undescribed multimorphic IRF4 pathophysiology disrupts normal lymphocyte biology, causing human disease.

The SPOC domain is a phosphoserine binding module that bridges transcription machinery with co- and post-transcriptional regulators.

The heptad repeats of the C-terminal domain (CTD) of RNA polymerase II (Pol II) are extensively modified throughout the transcription cycle. The CTD coordinates RNA synthesis and processing by recruiting transcription regulators as well as RNA capping, splicing and 3'end processing factors. The SPOC domain of PHF3 was recently identified as a CTD reader domain specifically binding to phosphorylated serine-2 residues in adjacent CTD repeats. Here, we establish the SPOC domains of the human proteins DIDO, SHARP (also known as SPEN) and RBM15 as phosphoserine binding modules that can act as CTD readers but also recognize other phosphorylated binding partners. We report the crystal structure of SHARP SPOC in complex with CTD and identify the molecular determinants for its specific binding to phosphorylated serine-5. PHF3 and DIDO SPOC domains preferentially interact with the Pol II elongation complex, while RBM15 and SHARP SPOC domains engage with writers and readers of m6A, the most abundant RNA modification. RBM15 positively regulates m6A levels and mRNA stability in a SPOC-dependent manner, while SHARP SPOC is essential for its localization to inactive X-chromosomes. Our findings suggest that the SPOC domain is a major interface between the transcription machinery and regulators of transcription and co-transcriptional processes.

Herpesviruses mimic zygotic genome activation to promote viral replication.

DUX4 is a germline transcription factor and a master regulator of zygotic genome activation. During early embryogenesis, DUX4 is crucial for maternal to zygotic transition at the 2-8-cell stage in order to overcome silencing of genes and enable transcription from the zygotic genome. In adult somatic cells, DUX4 expression is silenced and its activation in adult muscle cells causes the genetic disorder Facioscapulohumeral Muscular Dystrophy (FSHD). Here we show that herpesviruses from alpha-, beta- and gamma-herpesvirus subfamilies as well as papillomaviruses actively induce DUX4 expression to promote viral transcription and replication. We demonstrate that HSV-1 immediate early proteins directly induce expression of DUX4 and its target genes including endogenous retroelements, which mimics zygotic genome activation. We further show that DUX4 directly binds to the viral genome and promotes viral transcription. DUX4 is functionally required for herpesvirus infection, since genetic depletion of DUX4 by CRISPR/Cas9 abrogates viral replication. Our results show that herpesviruses induce DUX4 expression and its downstream germline-specific genes and retroelements, thus mimicking an early embryonic-like transcriptional program that prevents epigenetic silencing of the viral genome and facilitates herpesviral gene expression.

Multi-Omics Alleviates the Limitations of Panel Sequencing for Cancer Drug Response Prediction.

Comprehensive genomic profiling using cancer gene panels has been shown to improve treatment options for a variety of cancer types. However, genomic aberrations detected via such gene panels do not necessarily serve as strong predictors of drug sensitivity. In this study, using pharmacogenomics datasets of cell lines, patient-derived xenografts, and ex vivo treated fresh tumor specimens, we demonstrate that utilizing the transcriptome on top of gene panel features substantially improves drug response prediction performance in cancer.

SARS-CoV-2 infection dynamics revealed by wastewater sequencing analysis and deconvolution.

The use of RNA sequencing from wastewater samples is a valuable way for estimating infection dynamics and circulating lineages of SARS-CoV-2. This approach is independent from testing individuals and can therefore become the key tool to monitor this and potentially other viruses. However, it is equally important to develop easily accessible and scalable tools which can highlight critical changes in infection rates and dynamics over time across different locations given sequencing data from wastewater. Here, we provide an analysis of lineage dynamics in Berlin and New York City using wastewater sequencing and present PiGx SARS-CoV-2, a highly reproducible computational analysis pipeline with comprehensive reports. This end-to-end pipeline includes all steps from raw data to shareable reports, additional taxonomic analysis, deconvolution and geospatial time series analyses. Using simulated datasets (in silico generated and spiked-in samples) we could demonstrate the accuracy of our pipeline calculating proportions of Variants of Concern (VOC) from environmental as well as pre-mixed samples (spiked-in). By applying our pipeline on a dataset of wastewater samples from Berlin between February 2021 and January 2022, we could reconstruct the emergence of B.1.1.7(alpha) in February/March 2021 and the replacement dynamics from B.1.617.2 (delta) to BA.1 and BA.2 (omicron) during the winter of 2021/2022. Using data from very-short-reads generated in an industrial scale setting, we could see even higher accuracy in our deconvolution. Lastly, using a targeted sequencing dataset from New York City (receptor-binding-domain (RBD) only), we could reproduce the results recovering the proportions of the so-called cryptic lineages shown in the original study. Overall our study provides an in-depth analysis reconstructing virus lineage dynamics from wastewater. While applying our tool on a wide range of different datasets (from different types of wastewater sample locations and sequenced with different methods), we show that PiGx SARS-CoV-2 can be used to identify new mutations and detect any emerging new lineages in a highly automated and scalable way. Our approach can support efforts to establish continuous monitoring and early-warning projects for detecting SARS-CoV-2 or any other pathogen.

Rapid Inflammasome Activation Is Attenuated in Post-Myocardial Infarction Monocytes.

Inflammasomes are crucial gatekeepers of the immune response, but their maladaptive activation associates with inflammatory pathologies. Besides canonical activation, monocytes can trigger non-transcriptional or rapid inflammasome activation that has not been well defined in the context of acute myocardial infarction (AMI). Rapid transcription-independent inflammasome activation induced by simultaneous TLR priming and triggering stimulus was measured by caspase-1 (CASP1) activity and interleukin release. Both classical and intermediate monocytes from healthy donors exhibited robust CASP1 activation, but only classical monocytes produced high mature interleukin-18 (IL18) release. We also recruited a limited number of coronary artery disease (CAD, n=31) and AMI (n=29) patients to evaluate their inflammasome function and expression profiles. Surprisingly, monocyte subpopulations isolated from blood collected during percutaneous coronary intervention (PCI) from AMI patients presented diminished CASP1 activity and abrogated IL18 release despite increased NLRP3 gene expression. This unexpected attenuated rapid inflammasome activation was accompanied by a significant increase of TNFAIP3 and IRAKM expression. Moreover, TNFAIP3 protein levels of circulating monocytes showed positive correlation with high sensitive troponin T (hsTnT), implying an association between TNFAIP3 upregulation and the severity of tissue injury. We suggest this monocyte attenuation to be a protective phenotype aftermath following a very early inflammatory wave in the ischemic area. Damage-associated molecular patterns (DAMPs) or other signals trigger a transitory negative feedback loop within newly recruited circulating monocytes as a mechanism to reduce post-injury tissue damage.

Identifying tumor cells at the single-cell level using machine learning.

Tumors are complex tissues of cancerous cells surrounded by a heterogeneous cellular microenvironment with which they interact. Single-cell sequencing enables molecular characterization of single cells within the tumor. However, cell annotation-the assignment of cell type or cell state to each sequenced cell-is a challenge, especially identifying tumor cells within single-cell or spatial sequencing experiments. Here, we propose ikarus, a machine learning pipeline aimed at distinguishing tumor cells from normal cells at the single-cell level. We test ikarus on multiple single-cell datasets, showing that it achieves high sensitivity and specificity in multiple experimental contexts.

Single-cell-resolved dynamics of chromatin architecture delineate cell and regulatory states in zebrafish embryos.

DNA accessibility of cis-regulatory elements (CREs) dictates transcriptional activity and drives cell differentiation during development. While many genes regulating embryonic development have been identified, the underlying CRE dynamics controlling their expression remain largely uncharacterized. To address this, we produced a multimodal resource and genomic regulatory map for the zebrafish community, which integrates single-cell combinatorial indexing assay for transposase-accessible chromatin with high-throughput sequencing (sci-ATAC-seq) with bulk histone PTMs and Hi-C data to achieve a genome-wide classification of the regulatory architecture determining transcriptional activity in the 24-h post-fertilization (hpf) embryo. We characterized the genome-wide chromatin architecture at bulk and single-cell resolution, applying sci-ATAC-seq on whole 24-hpf stage zebrafish embryos, generating accessibility profiles for ∼23,000 single nuclei. We developed a genome segmentation method, ScregSeg (single-cell regulatory landscape segmentation), for defining regulatory programs, and candidate CREs, specific to one or more cell types. We integrated the ScregSeg output with bulk measurements for histone post-translational modifications and 3D genome organization and identified new regulatory principles between chromatin modalities prevalent during zebrafish development. Sci-ATAC-seq profiling of npas4l/cloche mutant embryos identified novel cellular roles for this hematovascular transcriptional master regulator and suggests an intricate mechanism regulating its expression. Our work defines regulatory architecture and principles in the zebrafish embryo and establishes a resource of cell-type-specific genome-wide regulatory annotations and candidate CREs, providing a valuable open resource for genomics, developmental, molecular, and computational biology.

Cell-type specialization is encoded by specific chromatin topologies.

The three-dimensional (3D) structure of chromatin is intrinsically associated with gene regulation and cell function1-3. Methods based on chromatin conformation capture have mapped chromatin structures in neuronal systems such as in vitro differentiated neurons, neurons isolated through fluorescence-activated cell sorting from cortical tissues pooled from different animals and from dissociated whole hippocampi4-6. However, changes in chromatin organization captured by imaging, such as the relocation of Bdnf away from the nuclear periphery after activation7, are invisible with such approaches8. Here we developed immunoGAM, an extension of genome architecture mapping (GAM)2,9, to map 3D chromatin topology genome-wide in specific brain cell types, without tissue disruption, from single animals. GAM is a ligation-free technology that maps genome topology by sequencing the DNA content from thin (about 220 nm) nuclear cryosections. Chromatin interactions are identified from the increased probability of co-segregation of contacting loci across a collection of nuclear slices. ImmunoGAM expands the scope of GAM to enable the selection of specific cell types using low cell numbers (approximately 1,000 cells) within a complex tissue and avoids tissue dissociation2,10. We report cell-type specialized 3D chromatin structures at multiple genomic scales that relate to patterns of gene expression. We discover extensive 'melting' of long genes when they are highly expressed and/or have high chromatin accessibility. The contacts most specific of neuron subtypes contain genes associated with specialized processes, such as addiction and synaptic plasticity, which harbour putative binding sites for neuronal transcription factors within accessible chromatin regions. Moreover, sensory receptor genes are preferentially found in heterochromatic compartments in brain cells, which establish strong contacts across tens of megabases. Our results demonstrate that highly specific chromatin conformations in brain cells are tightly related to gene regulation mechanisms and specialized functions.

Lymphocyte access to lymphoma is impaired by high endothelial venule regression.

Blood endothelial cells display remarkable plasticity depending on the demands of a malignant microenvironment. While studies in solid tumors focus on their role in metabolic adaptations, formation of high endothelial venules (HEVs) in lymph nodes extends their role to the organization of immune cell interactions. As a response to lymphoma growth, blood vessel density increases; however, the fate of HEVs remains elusive. Here, we report that lymphoma causes severe HEV regression in mouse models that phenocopies aggressive human B cell lymphomas. HEV dedifferentiation occurrs as a consequence of a disrupted lymph-carrying conduit system. Mechanosensitive fibroblastic reticular cells then deregulate CCL21 migration paths, followed by deterioration of dendritic cell proximity to HEVs. Loss of this crosstalk deprives HEVs of lymphotoxin-ß-receptor (LTßR) signaling, which is indispensable for their differentiation and lymphocyte transmigration. Collectively, this study reveals a remodeling cascade of the lymph node microenvironment that is detrimental for immune cell trafficking in lymphoma.

PHF3 regulates neuronal gene expression through the Pol II CTD reader domain SPOC.

The C-terminal domain (CTD) of the largest subunit of RNA polymerase II (Pol II) is a regulatory hub for transcription and RNA processing. Here, we identify PHD-finger protein 3 (PHF3) as a regulator of transcription and mRNA stability that docks onto Pol II CTD through its SPOC domain. We characterize SPOC as a CTD reader domain that preferentially binds two phosphorylated Serine-2 marks in adjacent CTD repeats. PHF3 drives liquid-liquid phase separation of phosphorylated Pol II, colocalizes with Pol II clusters and tracks with Pol II across the length of genes. PHF3 knock-out or SPOC deletion in human cells results in increased Pol II stalling, reduced elongation rate and an increase in mRNA stability, with marked derepression of neuronal genes. Key neuronal genes are aberrantly expressed in Phf3 knock-out mouse embryonic stem cells, resulting in impaired neuronal differentiation. Our data suggest that PHF3 acts as a prominent effector of neuronal gene regulation by bridging transcription with mRNA decay.

Predicting lethal courses in critically ill COVID-19 patients using a machine learning model trained on patients with non-COVID-19 viral pneumonia.

In a pandemic with a novel disease, disease-specific prognosis models are available only with a delay. To bridge the critical early phase, models built for similar diseases might be applied. To test the accuracy of such a knowledge transfer, we investigated how precise lethal courses in critically ill COVID-19 patients can be predicted by a model trained on critically ill non-COVID-19 viral pneumonia patients. We trained gradient boosted decision tree models on 718 (245 deceased) non-COVID-19 viral pneumonia patients to predict individual ICU mortality and applied it to 1054 (369 deceased) COVID-19 patients. Our model showed a significantly better predictive performance (AUROC 0.86 [95% CI 0.86-0.87]) than the clinical scores APACHE2 (0.63 [95% CI 0.61-0.65]), SAPS2 (0.72 [95% CI 0.71-0.74]) and SOFA (0.76 [95% CI 0.75-0.77]), the COVID-19-specific mortality prediction models of Zhou (0.76 [95% CI 0.73-0.78]) and Wang (laboratory: 0.62 [95% CI 0.59-0.65]; clinical: 0.56 [95% CI 0.55-0.58]) and the 4C COVID-19 Mortality score (0.71 [95% CI 0.70-0.72]). We conclude that lethal courses in critically ill COVID-19 patients can be predicted by a machine learning model trained on non-COVID-19 patients. Our results suggest that in a pandemic with a novel disease, prognosis models built for similar diseases can be applied, even when the diseases differ in time courses and in rates of critical and lethal courses.

Parallel genetics of regulatory sequences using scalable genome editing in vivo.

How regulatory sequences control gene expression is fundamental for explaining phenotypes in health and disease. Regulatory elements must ultimately be understood within their genomic environment and development- or tissue-specific contexts. Because this is technically challenging, few regulatory elements have been characterized in vivo. Here, we use inducible Cas9 and multiplexed guide RNAs to create hundreds of mutations in enhancers/promoters and 3' UTRs of 16 genes in C. elegans. Our software crispr-DART analyzes indel mutations in targeted DNA sequencing. We quantify the impact of mutations on expression and fitness by targeted RNA sequencing and DNA sampling. When applying our approach to the lin-41 3' UTR, generating hundreds of mutants, we find that the two adjacent binding sites for the miRNA let-7 can regulate lin-41 expression independently of each other. Finally, we map regulatory genotypes to phenotypic traits for several genes. Our approach enables parallel analysis of regulatory sequences directly in animals.

Transcriptomic profiling of SARS-CoV-2 infected human cell lines identifies HSP90 as target for COVID-19 therapy.

Detailed knowledge of the molecular biology of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is crucial for understanding of viral replication, host responses, and disease progression. Here, we report gene expression profiles of three SARS-CoV- and SARS-CoV-2-infected human cell lines. SARS-CoV-2 elicited an approximately two-fold higher stimulation of the innate immune response compared to SARS-CoV in the human epithelial cell line Calu-3, including induction of miRNA-155. Single-cell RNA sequencing of infected cells showed that genes induced by virus infections were broadly upregulated, whereas interferon beta/lambda genes, a pro-inflammatory cytokines such as IL-6, were expressed only in small subsets of infected cells. Temporal analysis suggested that transcriptional activities of interferon regulatory factors precede those of nuclear factor κB. Lastly, we identified heat shock protein 90 (HSP90) as a protein relevant for the infection. Inhibition of the HSP90 activity resulted in a reduction of viral replication and pro-inflammatory cytokine expression in primary human airway epithelial cells.