Rare predicted loss-of-function variants of type I IFN immunity genes are associated with critical COVID-19

BackgroundWe previously reported inborn errors of TLR3- and TLR7-dependent type I interferon (IFN) immunity in 1-5% of unvaccinated patients with life-threatening COVID-19, and auto-antibodies against type I IFN in another 15-20% of cases. MethodsWe report here a genome-wide rare variant burden association analysis in 3,269 unvaccinated patients with life-threatening COVID-19 (1,301 previously reported and 1,968 new patients), and 1,373 unvaccinated SARS-CoV-2-infected individuals without pneumonia. A quarter of the patients tested had antibodies against type I IFN (234 of 928) and were excluded from the analysis. ResultsNo gene reached genome-wide significance. Under a recessive model, the most significant gene with at-risk variants was TLR7, with an OR of 27.68 (95%CI:1.5-528.7, P=1.1x10-4), in analyses restricted to biochemically loss-of-function (bLOF) variants. We replicated the enrichment in rare predicted LOF (pLOF) variants at 13 influenza susceptibility loci involved in TLR3-dependent type I IFN immunity (OR=3.70 [95%CI:1.3-8.2], P=2.1x10-4). Adding the recently reported TYK2 COVID-19 locus strengthened this enrichment, particularly under a recessive model (OR=19.65 [95%CI:2.1-2635.4]; P=3.4x10-3). When these 14 loci and TLR7 were considered, all individuals hemizygous (n=20) or homozygous (n=5) for pLOF or bLOF variants were patients (OR=39.19 [95%CI:5.2-5037.0], P=4.7x10-7), who also showed an enrichment in heterozygous variants (OR=2.36 [95%CI:1.0-5.9], P=0.02). Finally, the patients with pLOF or bLOF variants at these 15 loci were significantly younger (mean age [SD]=43.3 [20.3] years) than the other patients (56.0 [17.3] years; P=1.68x10-5). ConclusionsRare variants of TLR3- and TLR7-dependent type I IFN immunity genes can underlie life-threatening COVID-19, particularly with recessive inheritance, in patients under 60 years old.

Multicenter analysis of neutrophil extracellular trap dysregulation in adult and pediatric COVID-19

Dysregulation in neutrophil extracellular trap (NET) formation and degradation may play a role in the pathogenesis and severity of COVID-19; however, its role in the pediatric manifestations of this disease including MIS-C and chilblain-like lesions (CLL), otherwise known as "COVID toes", remains unclear. Studying multinational cohorts, we found that, in CLL, NETs were significantly increased in serum and skin. There was geographic variability in the prevalence of increased NETs in MIS-C, in association with disease severity. MIS-C and CLL serum samples displayed decreased NET degradation ability, in association with C1q and G-actin or anti-NET antibodies, respectively, but not with genetic variants of DNases. In adult COVID-19, persistent elevations in NETs post-disease diagnosis were detected but did not occur in asymptomatic infection. COVID-19-affected adults displayed significant prevalence of impaired NET degradation, in association with anti-DNase1L3, G-actin, and specific disease manifestations, but not with genetic variants of DNases. NETs were detected in many organs of adult patients who died from COVID-19 complications. Infection with the Omicron variant was associated with decreased levels of NETs when compared to other SARS-CoV-2 strains. These data support a role for NETs in the pathogenesis and severity of COVID-19 in pediatric and adult patients. SummaryNET formation and degradation are dysregulated in pediatric and symptomatic adult patients with various complications of COVID-19, in association with disease severity. NET degradation impairments are multifactorial and associated with natural inhibitors of DNase 1, G-actin and anti-DNase1L3 and anti-NET antibodies. Infection with the Omicron variant is associated with decreased levels of NETs when compared to other SARS-CoV-2 strains.

ImmunoTyper-SR: A Novel Computational Approach for Genotyping Immunoglobulin Heavy Chain Variable Genes using Short Read Data

Human immunoglobulin heavy chain (IGH) locus on chromosome 14 includes more than 40 functional copies of the variable gene (IGHV), which, together with the joining genes (IGHJ), diversity genes (IGHD), constant genes (IGHC) and immunoglobulin light chains, code for antibodies that identify and neutralize pathogenic invaders as a part of the adaptive immune system. Because of its highly repetitive sequence composition, the IGH locus has been particularly difficult to assemble or genotype through the use of standard short read sequencing technologies. Here we introduce ImmunoTyper-SR, an algorithmic method for genotype and CNV analysis of the germline IGHV genes using Illumina whole genome sequencing (WGS) data. ImmunoTyper-SR is based on a novel combinatorial optimization formulation that aims to minimize the total edit distance between reads and their assigned IGHV alleles from a given database, with constraints on the number and distribution of reads across each called allele. We have validated ImmunoTyper-SR on 12 individuals with Illumina WGS data from the 1000 Genomes Project, whose IGHV allele composition have been studied extensively through the use of long read and targeted sequencing platforms, as well as nine individuals from the NIAID COVID Consortium who have been subjected to WGS twice. We have then applied ImmunoTyper-SR on 585 samples from the NIAID COVID Consortium to investigate associations between distinct IGHV alleles and anti-type I IFN autoantibodies which have been linked to COVID-19 severity.

Genetic regulation of OAS1 nonsense-mediated decay underlies association with risk of severe COVID-19

Genomic regions have been associated with COVID-19 susceptibility and outcomes, including the chr12q24.13 locus encoding antiviral proteins OAS1-3. Here, we report genetic, functional, and clinical insights into genetic associations within this locus. In Europeans, the risk of hospitalized vs. non-hospitalized COVID-19 was associated with a single 19Kb-haplotype comprised of 76 OAS1 variants included in a 95% credible set within a large genomic fragment introgressed from Neandertals. The risk haplotype was also associated with impaired spontaneous but not treatment-induced SARS-CoV-2 clearance in a clinical trial with pegIFN-{lambda}1. We demonstrate that two exonic variants, rs10774671 and rs1131454, affect splicing and nonsense-mediated decay of OAS1. We suggest that genetically-regulated loss of OAS1 expression contributes to impaired spontaneous clearance of SARS-CoV-2 and elevated risk of hospitalization for COVID-19. Our results provide the rationale for further clinical studies using interferons to compensate for impaired spontaneous SARS-CoV-2 clearance, particularly in carriers of the OAS1 risk haplotypes.

Magnitude and Dynamics of the T-Cell Response to SARS-CoV-2 Infection at Both Individual and Population Levels

T cells are involved in the early identification and clearance of viral infections and also support the development of antibodies by B cells. This central role for T cells makes them a desirable target for assessing the immune response to SARS-CoV-2 infection. Here, we combined two high-throughput immune profiling methods to create a quantitative picture of the T-cell response to SARS-CoV-2. First, at the individual level, we deeply characterized 3 acutely infected and 58 recovered COVID-19 subjects by experimentally mapping their CD8 T-cell response through antigen stimulation to 545 Human Leukocyte Antigen (HLA) class I presented viral peptides (class II data in a forthcoming study). Then, at the population level, we performed T-cell repertoire sequencing on 1,815 samples (from 1,521 COVID-19 subjects) as well as 3,500 controls to identify shared "public" T-cell receptors (TCRs) associated with SARS-CoV-2 infection from both CD8 and CD4 T cells. Collectively, our data reveal that CD8 T-cell responses are often driven by a few immunodominant, HLA-restricted epitopes. As expected, the T-cell response to SARS-CoV-2 peaks about one to two weeks after infection and is detectable for at least several months after recovery. As an application of these data, we trained a classifier to diagnose SARSCoV-2 infection based solely on TCR sequencing from blood samples, and observed, at 99.8% specificity, high early sensitivity soon after diagnosis (Day 3-7 = 85.1% [95% CI = 79.9-89.7]; Day 8-14 = 94.8% [90.7-98.4]) as well as lasting sensitivity after recovery (Day 29+/convalescent = 95.4% [92.1-98.3]). These results demonstrate an approach to reliably assess the adaptive immune response both soon after viral antigenic exposure (before antibodies are typically detectable) as well as at later time points. This blood-based molecular approach to characterizing the cellular immune response has applications in clinical diagnostics as well as in vaccine development and monitoring.