Immune transcriptomes of highly exposed SARS-CoV-2 asymptomatic seropositive versus seronegative individuals from the Ischgl community

To investigate prevalence of ongoing activation of inflammation following asymptomatic SARS-CoV-2 infection we characterized immune cell transcriptomes from 43 asymptomatic seropositive and 52 highly exposed seronegative individuals with few underlying health issues following a community superspreading event. Four mildly symptomatic seropositive individuals examined three weeks after infection as positive controls demonstrated immunological activation. Approximately four to six weeks following the event, the two asymptomatic groups showed no significant differences. Two seropositive patients with underlying genetic disease impacting immunological activation were included (Cystic Fibrosis (CF), Nuclear factor-kappa B Essential Modulator (NEMO) deficiency). CF, but not NEMO, associated with significant immune transcriptome differences including some associated with severe SARS-CoV-2 infection (IL1B, IL17A, respective receptors). All subjects remained in their usual state of health from event through five-month follow-up. Here, asymptomatic infection resolved without evidence of prolonged immunological activation. Inclusion of subjects with underlying genetic disease illustrated the pathophysiological importance of context on impact of immunological response.

On occasion COVID-19 patients can suffer from longer term sequelae. In one survey, thirty-five percent of patients with mild outpatient-treated disease were reported as having not returned to their usual state of health by two to three weeks following infection 11 . Prolonged myocardial inflammation 12 and subacute thyroiditis 13 post resolution of acute infection have been reported. Although it is increasingly being recognized that up to 96% of infected individuals are asymptomatic 14,15 , their immune responses and the prevalence of unrecognized ongoing inflammation have not been investigated and are not understood.
With the pandemic extending itself into populations with underlying genetic disease, there is an urgent need to understand their immune response to SARS-CoV-2 infections. Towards this end we determined the immune transcriptome of an asymptomatic seropositive patient with Cystic Fibrosis (CF) and one with Nuclear factorkappa B Essential Modulator (NEMO) deficiency within our study population.

Results
for use under a CC0 license. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted September 2, 2020. . https://doi.org/10.1101/2020.09.01.20185884 doi: medRxiv preprint In an attempt to define the immune response in asymptomatic SARS-CoV-2 seropositive and highly exposed seronegative individuals within an isolated population, we conducted a study on residents from the ski resort of Ischgl that experienced a superspreading event in early March of 2020. This explosive local outbreak led to the spread of the virus throughout Austria and to many other European countries and worldwide 16 . Ischgl and the Paznaun valley were quarantined on March 13, 2020 and remained under lockdown for six weeks. An epidemiologic study targeting 79% (n=1473) of the population of Ischgl Adults in the Austrian Ski Resort Ischgl, submitted) with approximately 17% of these being asymptomatic. Our study encompassed 43 seropositive asymptomatic individuals (Group A) and 52 highly exposed seronegative individuals (Group B) with an equal gender and age distribution ( Fig. 1a and Table 1). Six households had both seropositive asymptomatic and highly exposed seronegative members (Supplementary Table 1). Only a few asymptomatic seropositive individuals had conditions that increased their risk of severe illness from COVID-19 17 , one patient with Cystic Fibrosis (CFTR G551D mutation) and one with Nuclear factor-kappa B Essential Modulator (NEMO) deficiency (Incontinentia pigmenti, IKBKG exon4_10del mutation) ( Table 1). All infected individuals in the study remained asymptomatic throughout their infection and were confirmed in in person interviews as having their usual state of health at the time of the study (4-6 weeks after the infection) by telephone and remained so at the second phone call five months following the outbreak.
To validate that RNA sequencing (RNA-seq) performed on RNA extracted from peripheral blood mononuclear cells (PBMCs) could be used to identify gene expression changes in patients following SARS-CoV-2 infection, we conducted an unbiased RNAseq analysis on PBMCs from four patients with mild symptoms (Group D), about three weeks after they had tested PCR-positive. These were compared to four highly exposed seronegative individuals (Group E) ( Fig. 1b and  This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted September 2, 2020. . https://doi.org/10.1101/2020.09.01.20185884 doi: medRxiv preprint roles in immune regulation: TNF-a/NFkB, mTORC1 signaling, IL2-STAT5 signaling, TGFb signaling and inflammatory response ( Fig. 2 Table 2). This demonstrates that the RNA-seq approach utilizing buffy-coat-isolated PBMCs identified statistically significant differences in inflammatory gene expression between infected and non-infected individuals. Notably, we observed significantly elevated expression of IL10 in the symptomatic patients, a cytokine whose elevated expression has been associated with disease severity 18 .

and Supplementary
We then conducted an unbiased RNA-seq transcriptome analysis (Supplementary Table 3) and plasma cytokine profiling ( Fig. 3 and Supplementary Table 4) from the asymptomatic seropositive (Group A) and highly exposed seronegative (Group B) Ischgl resident cohorts. Very few statistically significant changes in gene expression were found (11 induced, 7 down-regulated) (Supplementary Table 3). Quantification of inflammatory cytokines and chemokines revealed no significant differences in protein levels between Group A and B ( Fig. 3 and Supplementary Table 4).
The PBMC transcriptome data from the CF (Supplementary Table 5) and NEMO (Supplementary Table 6) patients in Group A were analyzed separately. While the NEMO patient showed no statistically significant differences as compared to Group A, the CF patient demonstrated statistically significant differences in expression for approximately 4670 genes in comparison with Group A. Overall, expression levels of 3020 genes were significantly higher and expression levels of 1648 genes were significantly lower (Supplementary Table 5). Genes were enriched in 32 Hallmark gene sets, 11 of which have defined roles in immune regulation (Fig. 4a). Expression of key immune signature genes, including interferon response genes, IL1B, IL17A and their respective receptors, and JAK-STAT pathway genes, were significantly induced ( Fig. 4b-d).

Discussion
The ability to assess the immune response of both symptomatic and asymptomatic SARS-CoV-2 infected individuals might provide critical information on dysregulated immune-response signatures that could foretell disease trajectories. While longitudinal studies on hospitalized patients demonstrate elevated levels of pro-inflammatory cytokines and signatures associated with ongoing inflammation 3 , there is a parallel need for use under a CC0 license.
This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted September 2, 2020. . https://doi.org/10.1101/2020.09.01.20185884 doi: medRxiv preprint to pinpoint the immune response of infected, yet asymptomatic individuals. With the exception of a CF patient, we did not detect an aberrant immune transcriptome in 41 asymptomatic seropositive individuals that were infected during a super spreading event as compared to 52 highly exposed seronegative individuals from the same community.
Several households housed both seropositive and seronegative individuals. In contrast to the seropositive asymptomatic individuals, a proinflammatory immune signature was detected in patients with mild symptoms. These results demonstrate that development of an antibody response to COVID-19 following viral exposure and seroconversion in asymptomatic cases is not necessarily associated with sustained alterations in the immune system transcriptome.
While other studies have focused on single cells RNA-seq (scRNA-seq) profiling of small numbers of individuals with COVID-19 disease 4-7,9,10 , the use of RNA-seq from mononuclear cells permitted us to analyze larger cohorts at a great depth. This approach was validated through the detection of inflammatory immune signatures in COVID-19 patients exhibiting mild symptoms and one asymptomatic seropositive CF patient.
Patients with CF manifest cytokine dysfunction and hyperinflammation that overlaps with the pathophysiology of COVID-19 19 . While there are limited data on the immune response of CF patients to COVID-19 infection, preliminary information suggests that the course of disease may be milder than expected 20,21 . The immune transcriptome of the asymptomatic seropositive CF patient provided evidence of highly activated cytokine signaling pathways. Expression of key components of interferon, interleukin and JAK-STAT pathways is highly elevated. Many of these genes, such as IL1B and its receptors are also highly activated in COVID-19 patients 9 . In contrast, the immune transcriptome of a SRAS-CoV-2 infected asymptomatic patient with a NEMO deficiency syndrome, a rare primary immunodeficiency, was indistinguishable from asymptomatic seropositive controls. It remains to be understood how elevated cytokine signaling, documented here in an asymptomatic CF patient, contributes to disease progression in non-CF patients and why CF patients, with chronic high expression, do not invariably experience disease progression. A critical question for development of targeted therapies is to discern between direct pathogenic immune determinants of severe disease and correlates of inflammation.
for use under a CC0 license.
This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted September 2, 2020. . https://doi.org/10.1101/2020.09.01.20185884 doi: medRxiv preprint Limitations of our study include the translatability of our findings to other populations. There were few underlying health issues in this rural alpine population living at an altitude of 1,400 meters. For example, obesity, which is associated with an inflammatory state and is recognized as risk factor for severe COVID-19 disease 22,23 , was less than 10% in our study population, differing greatly from higher prevalence rate in other infected populations. Similarly, other defined risk factors for severe disease such as diabetes and chronic kidney disease, were also comparatively low. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Methods
The copyright holder for this preprint this version posted September 2, 2020. Extraction of the buffy coat and purification of RNA. Extraction of the buffy coat and subsequent RNA purification will be performed as described 10 . In short, the drawn blood is centrifuged at 1,600g for 10 min at 4°C. After vacuming off the plasma layer, the buffy coat layer is carefully collected. The obtained buffy coat is mixed with 1 mL RBC lysis buffer and incubated for 10 min at room temperature. After a centrifugation step, supernatants were discarded, and the pellets mixed with 1 mL RBC lysis buffer. The pellet is washed with PBS buffer and then mixed with 1 mL TRIzol ® . For extraction of the RNA the TRIzol reagent single-step method will be used 11 . After addition of 2 M sodium acetate (pH 4), the tube is mixed thoroughly by invertion before adding chloroform/isoamyl alcohol This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted September 2, 2020. . https://doi.org/10.1101/2020.09.01.20185884 doi: medRxiv preprint mRNA-seq read quality control was done using Trimmomatic 24 (version 0.36) and STAR RNA-seq 25 (version STAR 2.5.4a) using 150bp paired-end mode was used to align the reads (hg19). HTSeq 26 was to retrieve the raw counts and subsequently, R (https://www.R-project.org/), Bioconductor 27 and DESeq2 28 were used. Additionally, the RUVSeq 29 package was applied to remove confounding factors. The data were prefiltered keeping only those genes, which have at least ten reads in total. Genes were categorized as significantly differentially expressed with an adjusted p-value (pAdj) below 0.05 and a fold change > 2 for up-regulated genes and a fold change of < -2 for downregulated ones. The visualization was done using dplyr (https://CRAN.Rproject.org/package=dplyr) and ggplot2 30 . The genes were cutoff in the standard, less than 5 value and less than 1 log2 fold change and then conducted gene enrichment analysis (https://www.gsea-msigdb.org/gsea/msigdb). for use under a CC0 license.
This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted September 2, 2020.   This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted September 2, 2020.  for use under a CC0 license. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted September 2, 2020. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

X IQR (whiskers).
for use under a CC0 license. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted September 2, 2020. . https://doi.org/10.1101/2020.09.01.20185884 doi: medRxiv preprint This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.