Human inherited complete STAT2 deficiency underlies inflammatory viral diseases

STAT2 is a transcription factor activated by type I and III IFNs. We report 23 patients with loss-of-function variants causing autosomal recessive (AR) complete STAT2 deficiency. Both cells transfected with mutant STAT2 alleles and the patients' cells displayed impaired expression of IFN-stimulated genes and impaired control of in vitro viral infections. Clinical manifestations from early childhood onward included severe adverse reaction to live attenuated viral vaccines (LAV) and severe viral infections, particularly critical influenza pneumonia, critical COVID-19 pneumonia, and herpes simplex virus type 1 (HSV-1) encephalitis. The patients displayed various types of hyperinflammation, often triggered by viral infection or after LAV administration, which probably attested to unresolved viral infection in the absence of STAT2-dependent types I and III IFN immunity. Transcriptomic analysis revealed that circulating monocytes, neutrophils, and CD8+ memory T cells contributed to this inflammation. Several patients died from viral infection or heart failure during a febrile illness with no identified etiology. Notably, the highest mortality occurred during early childhood. These findings show that AR complete STAT2 deficiency underlay severe viral diseases and substantially impacts survival.

Encephalitis and poor neuronal death-mediated control of herpes simplex virus in human inherited RIPK3 deficiency

Inborn errors of TLR3-dependent type I IFN immunity in cortical neurons underlie forebrain herpes simplex virus-1 (HSV-1) encephalitis (HSE) due to uncontrolled viral growth and subsequent cell death. We report an otherwise healthy patient with HSE who was compound heterozygous for nonsense (R422*) and frameshift (P493fs9*) RIPK3 variants. Receptor-interacting protein kinase 3 (RIPK3) is a ubiquitous cytoplasmic kinase regulating cell death outcomes, including apoptosis and necroptosis. In vitro, the R422* and P493fs9* RIPK3 proteins impaired cellular apoptosis and necroptosis upon TLR3, TLR4, or TNFR1 stimulation and ZBP1/DAI-mediated necroptotic cell death after HSV-1 infection. The patient's fibroblasts displayed no detectable RIPK3 expression. After TNFR1 or TLR3 stimulation, the patient's cells did not undergo apoptosis or necroptosis. After HSV-1 infection, the cells supported excessive viral growth despite normal induction of antiviral IFN-β and IFN-stimulated genes (ISGs). This phenotype was, nevertheless, rescued by application of exogenous type I IFN. The patient's human pluripotent stem cell (hPSC)-derived cortical neurons displayed impaired cell death and enhanced viral growth after HSV-1 infection, as did isogenic RIPK3-knockout hPSC-derived cortical neurons. Inherited RIPK3 deficiency therefore confers a predisposition to HSE by impairing the cell death-dependent control of HSV-1 in cortical neurons but not their production of or response to type I IFNs.

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

BACKGROUND

We previously reported that impaired type I IFN activity, due to inborn errors of TLR3- and TLR7-dependent type I interferon (IFN) immunity or to autoantibodies against type I IFN, account for 15-20% of cases of life-threatening COVID-19 in unvaccinated patients. Therefore, the determinants of life-threatening COVID-19 remain to be identified in ~ 80% of cases.

METHODS

We report here a genome-wide rare variant burden association analysis in 3269 unvaccinated patients with life-threatening COVID-19, and 1373 unvaccinated SARS-CoV-2-infected individuals without pneumonia. Among the 928 patients tested for autoantibodies against type I IFN, a quarter (234) were positive and were excluded.

RESULTS

No 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.1 × 10-4) for 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.1 × 10-4). This enrichment was further strengthened by (1) adding the recently reported TYK2 and TLR7 COVID-19 loci, particularly under a recessive model (OR = 19.65[95%CI 2.1-2635.4], P = 3.4 × 10-3), and (2) considering as pLOF branchpoint variants with potentially strong impacts on splicing among the 15 loci (OR = 4.40[9%CI 2.3-8.4], P = 7.7 × 10-8). Finally, the patients with pLOF/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.68 × 10-5).

CONCLUSIONS

Rare 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.

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

Background

We 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.

Methods

We 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.

Results

No 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.1×10 -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.1×10 -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.4×10 -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.7×10 -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.68×10 -5 ).

Conclusions

Rare 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.

Recessive inborn errors of type I IFN immunity in children with COVID-19 pneumonia

Recessive or dominant inborn errors of type I interferon (IFN) immunity can underlie critical COVID-19 pneumonia in unvaccinated adults. The risk of COVID-19 pneumonia in unvaccinated children, which is much lower than in unvaccinated adults, remains unexplained. In an international cohort of 112 children (<16 yr old) hospitalized for COVID-19 pneumonia, we report 12 children (10.7%) aged 1.5-13 yr with critical (7 children), severe (3), and moderate (2) pneumonia and 4 of the 15 known clinically recessive and biochemically complete inborn errors of type I IFN immunity: X-linked recessive TLR7 deficiency (7 children) and autosomal recessive IFNAR1 (1), STAT2 (1), or TYK2 (3) deficiencies. Fibroblasts deficient for IFNAR1, STAT2, or TYK2 are highly vulnerable to SARS-CoV-2. These 15 deficiencies were not found in 1,224 children and adults with benign SARS-CoV-2 infection without pneumonia (P = 1.2 × 10-11) and with overlapping age, sex, consanguinity, and ethnicity characteristics. Recessive complete deficiencies of type I IFN immunity may underlie ∼10% of hospitalizations for COVID-19 pneumonia in children.

Single-Cell and Bulk RNA-Sequencing Reveal Differences in Monocyte Susceptibility to Influenza A Virus Infection Between Africans and Europeans

There is considerable inter-individual and inter-population variability in response to viruses. The potential of monocytes to elicit type-I interferon responses has attracted attention to their role in viral infections. Here, we use single-cell RNA-sequencing to characterize the role of cellular heterogeneity in human variation of monocyte responses to influenza A virus (IAV) exposure. We show widespread inter-individual variability in the percentage of IAV-infected monocytes. Notably, individuals with high cellular susceptibility to IAV are characterized by a lower activation at basal state of an IRF/STAT-induced transcriptional network, which includes antiviral genes such as IFITM3, MX1 and OAS3. Upon IAV challenge, we find that cells escaping viral infection display increased mRNA expression of type-I interferon stimulated genes and decreased expression of ribosomal genes, relative to both infected cells and those never exposed to IAV. We also uncover a stronger resistance of CD16⁺ monocytes to IAV infection, together with CD16⁺ -specific mRNA expression of IL6 and TNF in response to IAV. Finally, using flow cytometry and bulk RNA-sequencing across 200 individuals of African and European ancestry, we observe a higher number of CD16 ⁺ monocytes and lower susceptibility to IAV infection among monocytes from individuals of African-descent. Based on these data, we hypothesize that higher basal monocyte activation, driven by environmental factors and/or weak-effect genetic variants, underlies the lower cellular susceptibility to IAV infection of individuals of African ancestry relative to those of European ancestry. Further studies are now required to investigate how such cellular differences in IAV susceptibility translate into population differences in clinical outcomes and susceptibility to severe influenza.

X-linked recessive TLR7 deficiency in ~1% of men under 60 years old with life-threatening COVID-19

Autosomal inborn errors of type I IFN immunity and autoantibodies against these cytokines underlie at least 10% of critical COVID-19 pneumonia cases. We report very rare, biochemically deleterious X-linked TLR7 variants in 16 unrelated male individuals aged 7 to 71 years (mean: 36.7 years) from a cohort of 1,202 male patients aged 0.5 to 99 years (mean: 52.9 years) with unexplained critical COVID-19 pneumonia. None of the 331 asymptomatically or mildly infected male individuals aged 1.3 to 102 years (mean: 38.7 years) tested carry such TLR7 variants (p = 3.5 × 10-5). The phenotypes of five hemizygous relatives of index cases infected with SARS-CoV-2 include asymptomatic or mild infection (n=2, 5 and 38 years), or moderate (n=1, 5 years), severe (n=1, 27 years), or critical (n=1, 29 years) pneumonia. Two boys (aged 7 and 12 years) from a cohort of 262 male patients with severe COVID-19 pneumonia (mean: 51.0 years) are hemizygous for a deleterious TLR7 variant. The cumulative allele frequency for deleterious TLR7 variants in the male general population is < 6.5x10-4 We also show that blood B cell lines and myeloid cell subsets from the patients do not respond to TLR7 stimulation, a phenotype rescued by wild-type TLR7 The patients' blood plasmacytoid dendritic cells (pDCs) produce low levels of type I IFNs in response to SARS-CoV-2. Overall, X-linked recessive TLR7 deficiency is a highly penetrant genetic etiology of critical COVID-19 pneumonia, in about 1.8% of male patients below the age of 60 years. Human TLR7 and pDCs are essential for protective type I IFN immunity against SARS-CoV-2 in the respiratory tract.

A computational approach for detecting physiological homogeneity in the midst of genetic heterogeneity

The human genetic dissection of clinical phenotypes is complicated by genetic heterogeneity. Gene burden approaches that detect genetic signals in case-control studies are underpowered in genetically heterogeneous cohorts. We therefore developed a genome-wide computational method, network-based heterogeneity clustering (NHC), to detect physiological homogeneity in the midst of genetic heterogeneity. Simulation studies showed our method to be capable of systematically converging genes in biological proximity on the background biological interaction network, and capturing gene clusters harboring presumably deleterious variants, in an efficient and unbiased manner. We applied NHC to whole-exome sequencing data from a cohort of 122 individuals with herpes simplex encephalitis (HSE), including 13 individuals with previously published monogenic inborn errors of TLR3-dependent IFN-α/β immunity. The top gene cluster identified by our approach successfully detected and prioritized all causal variants of five TLR3 pathway genes in the 13 previously reported individuals. This approach also suggested candidate variants of three reported genes and four candidate genes from the same pathway in another ten previously unstudied individuals. TLR3 responsiveness was impaired in dermal fibroblasts from four of the five individuals tested, suggesting that the variants detected were causal for HSE. NHC is, therefore, an effective and unbiased approach for unraveling genetic heterogeneity by detecting physiological homogeneity.