Human IRF1 governs macrophagic IFN-γ immunity to mycobacteria.

Inborn errors of human IFN-γ-dependent macrophagic immunity underlie mycobacterial diseases, whereas inborn errors of IFN-α/ß-dependent intrinsic immunity underlie viral diseases. Both types of IFNs induce the transcription factor IRF1. We describe unrelated children with inherited complete IRF1 deficiency and early-onset, multiple, life-threatening diseases caused by weakly virulent mycobacteria and related intramacrophagic pathogens. These children have no history of severe viral disease, despite exposure to many viruses, including SARS-CoV-2, which is life-threatening in individuals with impaired IFN-α/ß immunity. In leukocytes or fibroblasts stimulated in vitro, IRF1-dependent responses to IFN-γ are, both quantitatively and qualitatively, much stronger than those to IFN-α/ß. Moreover, IRF1-deficient mononuclear phagocytes do not control mycobacteria and related pathogens normally when stimulated with IFN-γ. By contrast, IFN-α/ß-dependent intrinsic immunity to nine viruses, including SARS-CoV-2, is almost normal in IRF1-deficient fibroblasts. Human IRF1 is essential for IFN-γ-dependent macrophagic immunity to mycobacteria, but largely redundant for IFN-α/ß-dependent antiviral immunity.

Allele-dependent interaction of LRRK2 and NOD2 in leprosy.

Leprosy, caused by Mycobacterium leprae, rarely affects children younger than 5 years. Here, we studied a multiplex leprosy family that included monozygotic twins aged 22 months suffering from paucibacillary leprosy. Whole genome sequencing identified three amino acid mutations previously associated with Crohn's disease and Parkinson's disease as candidate variants for early onset leprosy: LRRK2 N551K, R1398H and NOD2 R702W. In genome-edited macrophages, we demonstrated that cells expressing the LRRK2 mutations displayed reduced apoptosis activity following mycobacterial challenge independently of NOD2. However, employing co-immunoprecipitation and confocal microscopy we showed that LRRK2 and NOD2 proteins interacted in RAW cells and monocyte-derived macrophages, and that this interaction was substantially reduced for the NOD2 R702W mutation. Moreover, we observed a joint effect of LRRK2 and NOD2 variants on Bacillus Calmette-Guérin (BCG)-induced respiratory burst, NF-κB activation and cytokine/chemokine secretion with a strong impact for the genotypes found in the twins consistent with a role of the identified mutations in the development of early onset leprosy.

The risk of COVID-19 death is much greater and age dependent with type I IFN autoantibodies.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection fatality rate (IFR) doubles with every 5 y of age from childhood onward. Circulating autoantibodies neutralizing IFN-α, IFN-ω, and/or IFN-ß are found in ∼20% of deceased patients across age groups, and in ∼1% of individuals aged <70 y and in >4% of those >70 y old in the general population. With a sample of 1,261 unvaccinated deceased patients and 34,159 individuals of the general population sampled before the pandemic, we estimated both IFR and relative risk of death (RRD) across age groups for individuals carrying autoantibodies neutralizing type I IFNs, relative to noncarriers. The RRD associated with any combination of autoantibodies was higher in subjects under 70 y old. For autoantibodies neutralizing IFN-α2 or IFN-ω, the RRDs were 17.0 (95% CI: 11.7 to 24.7) and 5.8 (4.5 to 7.4) for individuals <70 y and ≥70 y old, respectively, whereas, for autoantibodies neutralizing both molecules, the RRDs were 188.3 (44.8 to 774.4) and 7.2 (5.0 to 10.3), respectively. In contrast, IFRs increased with age, ranging from 0.17% (0.12 to 0.31) for individuals <40 y old to 26.7% (20.3 to 35.2) for those ≥80 y old for autoantibodies neutralizing IFN-α2 or IFN-ω, and from 0.84% (0.31 to 8.28) to 40.5% (27.82 to 61.20) for autoantibodies neutralizing both. Autoantibodies against type I IFNs increase IFRs, and are associated with high RRDs, especially when neutralizing both IFN-α2 and IFN-ω. Remarkably, IFRs increase with age, whereas RRDs decrease with age. Autoimmunity to type I IFNs is a strong and common predictor of COVID-19 death.

Human genetics of Buruli ulcer.

Buruli ulcer, the third most common mycobacterial disease worldwide, is caused by Mycobacterium ulcerans and characterized by devastating necrotizing skin lesions. Susceptibility to Buruli ulcer is thought to depend on host genetics, but very few genetic studies have been performed. The identification of a microdeletion on chromosome 8 in a familial form of severe Buruli ulcer suggested a monogenic basis of susceptibility. The role of common host genetic variants in Buruli ulcer development has been investigated in only three candidate-gene studies targeting genes involved in mycobacterial diseases. A recent genome-wide association study suggested a probable role for long non-coding RNAs and strengthened the contribution of autophagy as a major defense mechanism against mycobacteria. In this review, we summarize the history, epidemiological and clinical aspects of Buruli ulcer, focusing particularly on genetic findings relating to susceptibility to this disease. Finally, we discuss exciting new genetic avenues arising, in particular, from studies of mouse models, and the need for different disciplines to work together, to benefit from the extensive work on other mycobacterial diseases, mostly tuberculosis and leprosy. We are convinced that such pooling of effort will lead to the development of efficient novel strategies for combatting Buruli ulcer.

Deciphering the genetic control of gene expression following Mycobacterium leprae antigen stimulation.

Leprosy is a human infectious disease caused by Mycobacterium leprae. A strong host genetic contribution to leprosy susceptibility is well established. However, the modulation of the transcriptional response to infection and the mechanism(s) of disease control are poorly understood. To address this gap in knowledge of leprosy pathogenicity, we conducted a genome-wide search for expression quantitative trait loci (eQTL) that are associated with transcript variation before and after stimulation with M. leprae sonicate in whole blood cells. We show that M. leprae antigen stimulation mainly triggered the upregulation of immune related genes and that a substantial proportion of the differential gene expression is genetically controlled. Indeed, using stringent criteria, we identified 318 genes displaying cis-eQTL at an FDR of 0.01, including 66 genes displaying response-eQTL (reQTL), i.e. cis-eQTL that showed significant evidence for interaction with the M. leprae stimulus. Such reQTL correspond to regulatory variations that affect the interaction between human whole blood cells and M. leprae sonicate and, thus, likely between the human host and M. leprae bacilli. We found that reQTL were significantly enriched among binding sites of transcription factors that are activated in response to infection, and that they were enriched among single nucleotide polymorphisms (SNPs) associated with susceptibility to leprosy per se and Type-I Reaction, and seven of them have been targeted by recent positive selection. Our study suggested that natural selection shaped our genomic diversity to face pathogen exposure including M. leprae infection.

A genome wide association study identifies a lncRna as risk factor for pathological inflammatory responses in leprosy.

Leprosy Type-1 Reactions (T1Rs) are pathological inflammatory responses that afflict a sub-group of leprosy patients and result in peripheral nerve damage. Here, we employed a family-based GWAS in 221 families with 229 T1R-affect offspring with stepwise replication to identify risk factors for T1R. We discovered, replicated and validated T1R-specific associations with SNPs located in chromosome region 10p21.2. Combined analysis across the three independent samples resulted in strong evidence of association of rs1875147 with T1R (p = 4.5x10-8; OR = 1.54, 95% CI = 1.32-1.80). The T1R-risk locus was restricted to a lncRNA-encoding genomic interval with rs1875147 being an eQTL for the lncRNA. Since a genetic overlap between leprosy and inflammatory bowel disease (IBD) has been detected, we evaluated if the shared genetic control could be traced to the T1R endophenotype. Employing the results of a recent IBD GWAS meta-analysis we found that 10.6% of IBD SNPs available in our dataset shared a common risk-allele with T1R (p = 2.4x10-4). This finding points to a substantial overlap in the genetic control of clinically diverse inflammatory disorders.

Reemergence and amplification of tuberculosis in the Canadian arctic.

BACKGROUND: Between November 2011 and November 2012, a Canadian village of 933 persons had 50 culture-positive cases of tuberculosis, with 49 sharing the same genotype. METHODS: We performed Illumina-based whole-genome sequencing on Mycobacterium tuberculosis isolates from this village, during and before the outbreak. Phylogenetic trees were generated using the maximum likelihood method. RESULTS: Three distinct genotypes were identified. Strain I (n = 7) was isolated in 1991-1996. Strain II (n = 8) was isolated in 1996-2004. Strain III (n = 62) first appeared in 2007 and did not arise from strain I or II. Within strain III, there were 3 related but distinct clusters: IIIA, IIIB, and IIIC. Between 2007 and 2010, cluster IIIA predominated (11 of 22 vs 2 of 40; P < .001), whereas in 2011-2012 clusters IIIB (n = 18) and IIIC (n = 20) predominated over cluster IIIA (n = 11). Combined evolutionary and epidemiologic analysis of strain III cases revealed that the outbreak in 2011-2012 was the result of ≥6 temporally staggered events, spanning from 1 reactivation case to a point-source outbreak of 20 cases. CONCLUSIONS: After the disappearance of 2 strains of M. tuberculosis in this village, its reemergence in 2007 was followed by an epidemiologic amplification, affecting >5% of the population.

Association of TNFSF8 regulatory variants with excessive inflammatory responses but not leprosy per se.

BACKGROUND: Type 1 reactions (T1R) affect a considerable proportion of patients with leprosy. In those with T1R, the host immune response pathologically overcompensates for the actual infectious threat, resulting in nerve damage and permanent disability. Based on the results of a genome-wide association study of leprosy per se, we investigated the TNFSF15 chromosomal region for a possible contribution to susceptibility to T1R. METHODS: We performed a high-resolution association scan of the TNFSF15 locus to evaluate the association with T1R in 2 geographically and ethnically distinct populations: a family-based sample from Vietnam and a case-control sample from Brazil, comprising a total of 1768 subjects. RESULTS: In the Vietnamese sample, 47 single-nucleotide polymorphisms (SNPs) overlapping TNFSF15 and the adjacent TNFSF8 gene were associated with T1R but not with leprosy. Of the 47 SNPs, 39 were cis-expression quantitative trait loci (cis-eQTL) for TNFSF8 including SNPs located within the TNFSF15 gene. In the Brazilian sample, 18 of these cis-eQTL SNPs overlapping the TNFSF8 gene were validated for association with T1R. CONCLUSIONS: Taken together, these results indicate TNFSF8 and not TNFSF15 as an important T1R susceptibility gene. Our data support the need for infection genetics to go beyond genes for pathogen control to explore genes involved in a commensurate host response.

Different selective pressures shape the evolution of Toll-like receptors in human and African great ape populations.

The study of the genetic and selective landscape of immunity genes across primates can provide insight into the existing differences in susceptibility to infection observed between human and non-human primates. Here, we explored how selection has driven the evolution of a key family of innate immunity receptors, the Toll-like receptors (TLRs), in African great ape species. We sequenced the 10 TLRs in various populations of chimpanzees and gorillas, and analysed these data jointly with a human data set. We found that purifying selection has been more pervasive in great apes than in humans. Furthermore, in chimpanzees and gorillas, purifying selection has targeted TLRs irrespectively of whether they are endosomal or cell surface, in contrast to humans where strong selective constraints are restricted to endosomal TLRs. These observations suggest important differences in the relative importance of TLR-mediated pathogen sensing, such as that of recognition of flagellated bacteria by TLR5, between humans and great apes. Lastly, we used a population genetics-phylogenetics method that jointly analyses polymorphism and divergence data to detect fine-scale variation in selection pressures at specific codons within TLR genes. We identified different codons at different TLRs as being under positive selection in each species, highlighting that functional variation at these genes has conferred a selective advantage in immunity to infection to specific primate species. Overall, this study showed that the degree of selection driving the evolution of TLRs has largely differed between human and non-human primates, increasing our knowledge on their respective biological contribution to host defence in the natural setting.

The evolutionary landscape of cytosolic microbial sensors in humans.

Host-pathogen interactions are generally initiated by host recognition of microbial components or danger signals triggered by microbial invasion. This recognition involves germline-encoded microbial sensors or pattern-recognition receptors (PRRs). By studying the way in which natural selection has driven the evolution of these microbial sensors in humans, we can identify genes playing an essential role and distinguish them from other, more redundant genes. We characterized the sequence diversity of the NOD-like receptor family, including the NALP and NOD/IPAF subfamilies, in various populations worldwide and compared this diversity with that of other PRR families, such as Toll-like receptors (TLRs) and RIG-I-like receptors (RLRs). We found that most NALPs had evolved under strong selective constraints, suggesting that their functions are essential and possibly much broader than previously thought. Conversely, most NOD/IPAF subfamily members were subject to more relaxed selective constraints, suggesting greater redundancy. Furthermore, some NALP genes, including NLRP1, NLRP14, and CIITA, were found to have evolved adaptively. We identified those variants conferring a selective advantage on some human populations as the most likely targets of positive selection. More generally, the strength of selection differed considerably between the major families of microbial sensors. Endosomal TLRs and most NALPs were found to evolve under stronger purifying selection than most NOD/IPAF subfamily members and cell-surface TLRs and RLRs, suggesting some degree of redundancy in the signaling pathways triggered by these molecules. This study provides novel perspectives and experimentally testable hypotheses concerning the relative biological relevance of the various families of microbial sensors in humans.

The Basque paradigm: genetic evidence of a maternal continuity in the Franco-Cantabrian region since pre-Neolithic times.

Different lines of evidence point to the resettlement of much of western and central Europe by populations from the Franco-Cantabrian region during the Late Glacial and Postglacial periods. In this context, the study of the genetic diversity of contemporary Basques, a population located at the epicenter of the Franco-Cantabrian region, is particularly useful because they speak a non-Indo-European language that is considered to be a linguistic isolate. In contrast with genome-wide analysis and Y chromosome data, where the problem of poor time estimates remains, a new timescale has been established for the human mtDNA and makes this genome the most informative marker for studying European prehistory. Here, we aim to increase knowledge of the origins of the Basque people and, more generally, of the role of the Franco-Cantabrian refuge in the postglacial repopulation of Europe. We thus characterize the maternal ancestry of 908 Basque and non-Basque individuals from the Basque Country and immediate adjacent regions and, by sequencing 420 complete mtDNA genomes, we focused on haplogroup H. We identified six mtDNA haplogroups, H1j1, H1t1, H2a5a1, H1av1, H3c2a, and H1e1a1, which are autochthonous to the Franco-Cantabrian region and, more specifically, to Basque-speaking populations. We detected signals of the expansion of these haplogroups at ∼4,000 years before present (YBP) and estimated their separation from the pan-European gene pool at ∼8,000 YBP, antedating the Indo-European arrival to the region. Our results clearly support the hypothesis of a partial genetic continuity of contemporary Basques with the preceding Paleolithic/Mesolithic settlers of their homeland.

Evidence of pre-Roman tribal genetic structure in Basques from uniparentally inherited markers.

Basque people have received considerable attention from anthropologists, geneticists, and linguists during the last century due to the singularity of their language and to other cultural and biological characteristics. Despite the multidisciplinary efforts performed to address the questions of the origin, uniqueness, and heterogeneity of Basques, the genetic studies performed up to now have suffered from a weak study design where populations are not analyzed in an adequate geographic and population context. To address the former questions and to overcome these design limitations, we have analyzed the uniparentally inherited markers (Y chromosome and mitochondrial DNA) of ~900 individuals from 18 populations, including those where Basque is currently spoken and populations from adjacent regions where Basque might have been spoken in historical times. Our results indicate that Basque-speaking populations fall within the genetic Western European gene pool, that they are similar to geographically surrounding non-Basque populations, and also that their genetic uniqueness is based on a lower amount of external influences compared with other Iberians and French populations. Our data suggest that the genetic heterogeneity and structure observed in the Basque region result from pre-Roman tribal structure related to geography and might be linked to the increased complexity of emerging societies during the Bronze Age. The rough overlap of the pre-Roman tribe location and the current dialect limits support the notion that the environmental diversity in the region has played a recurrent role in cultural differentiation and ethnogenesis at different time periods.

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.

Evolution of the TIR domain-containing adaptors in humans: swinging between constraint and adaptation.

Natural selection is expected to act strongly on immune system genes as hosts adapt to novel, diverse, and coevolving pathogens. Population genetic studies of host defense genes with parallel functions in model organisms have revealed distinct evolutionary histories among the different components-receptors, adaptors, and effectors-of the innate immune system. In humans, however, detailed evolutionary studies have been mainly confined to the receptors and in particular to Toll-like receptors (TLRs). By virtue of a toll/interleukin-1 receptor (TIR) domain, TLRs activate distinct signaling pathways, which are mediated by the five TIR-containing adaptors: myeloid differentiation factor-88 (MyD88), myeloid differentiation factor-88 adaptor-like protein (MAL), toll/interleukin-1 receptor domain-containing adaptor protein inducing interferon (IFN)ß (TRIF), toll/interleukin-1 receptor domain-containing adaptor protein inducing IFNß-related adaptor molecule (TRAM), and sterile α- and armadillo motif-containing protein (SARM). Here, we have examined the extent to which natural selection has affected immune adaptors in humans, using as a paradigm the TIR-containing adaptors. To do so, we characterized their levels of naturally occurring genetic variation in various human populations. We found that MyD88 and TRIF have mainly evolved under purifying selection, suggesting that their role in the early stages of signal transduction is essential and nonredundant for host survival. In addition, the adaptors have been targeted by multiple episodes of positive selection, differing in timing and spatial location. MyD88 and SARM display signatures of a selective sweep that has occurred in all humans, whereas for the other three adaptors, we detected signatures of adaptive evolution that are restricted to specific populations. Our study provides evidence that the contemporary diversity of the five TIR-containing adaptors results from the intermingling of different selective events, swinging between constraint and adaptation.

[Population genetics and human immunity: the interferon paradigm]. / Génétique des populations et immunité chez l'homme - Le cas des interférons.

The evolutionary genetics dissection of immunity-related genes provides insights into immunological defence mechanisms and highlight host pathways playing an important role in pathogen resistance. Recent population genetic data have increased knowledge of the biological relevance of human interferons (IFN), cytokines released by host cells in response to pathogen presence or tumour cells. Some IFN-α subtypes as well as IFN-γ are strongly evolutionarily constrained, suggesting that the functions they fulfil are essential and non redundant. Other IFN, the most extreme cases being IFN-α10 and IFN-ε, can accumulate missense or nonsense mutations at high population frequencies, suggesting higher redundancy. Furthermore, genetic variation at some IFN genes can be advantageous for the host and increase in frequency by positive selection. This has been shown for type III IFN, where mutations at IL28A, IL28B and IL29 have been positively selected in Europeans and Asians, most likely by increasing resistance to viral infection. This review uses the IFN paradigm to illustrate the value of the evolutionary approach in highlighting important determinants of host immune responsiveness in the natural setting.

The risk of COVID-19 death is much greater and age-dependent with type I IFN autoantibodies.

SARS-CoV-2 infection fatality rate (IFR) doubles with every five years of age from childhood onward. Circulating autoantibodies neutralizing IFN-α, IFN-ω, and/or IFN-ß are found in ~20% of deceased patients across age groups. In the general population, they are found in ~1% of individuals aged 20-70 years and in >4% of those >70 years old. With a sample of 1,261 deceased patients and 34,159 uninfected individuals, we estimated both IFR and relative risk of death (RRD) across age groups for individuals carrying autoantibodies neutralizing type I IFNs, relative to non-carriers. For autoantibodies neutralizing IFN-α2 or IFN-ω, the RRD was 17.0[95% CI:11.7-24.7] for individuals under 70 years old and 5.8[4.5-7.4] for individuals aged 70 and over, whereas, for autoantibodies neutralizing both molecules, the RRD was 188.3[44.8-774.4] and 7.2[5.0-10.3], respectively. IFRs increased with age, from 0.17%[0.12-0.31] for individuals <40 years old to 26.7%[20.3-35.2] for those ≥80 years old for autoantibodies neutralizing IFN-α2 or IFN-ω, and from 0.84%[0.31-8.28] to 40.5%[27.82-61.20] for the same two age groups, for autoantibodies neutralizing both molecules. Autoantibodies against type I IFNs increase IFRs, and are associated with high RRDs, particularly those neutralizing both IFN-α2 and -ω. Remarkably, IFR increases with age, whereas RRD decreases with age. Autoimmunity to type I IFNs appears to be second only to age among common predictors of COVID-19 death.

Autoantibodies against type I IFNs in patients with critical influenza pneumonia.

Autoantibodies neutralizing type I interferons (IFNs) can underlie critical COVID-19 pneumonia and yellow fever vaccine disease. We report here on 13 patients harboring autoantibodies neutralizing IFN-α2 alone (five patients) or with IFN-ω (eight patients) from a cohort of 279 patients (4.7%) aged 6-73 yr with critical influenza pneumonia. Nine and four patients had antibodies neutralizing high and low concentrations, respectively, of IFN-α2, and six and two patients had antibodies neutralizing high and low concentrations, respectively, of IFN-ω. The patients' autoantibodies increased influenza A virus replication in both A549 cells and reconstituted human airway epithelia. The prevalence of these antibodies was significantly higher than that in the general population for patients <70 yr of age (5.7 vs. 1.1%, P = 2.2 × 10-5), but not >70 yr of age (3.1 vs. 4.4%, P = 0.68). The risk of critical influenza was highest in patients with antibodies neutralizing high concentrations of both IFN-α2 and IFN-ω (OR = 11.7, P = 1.3 × 10-5), especially those <70 yr old (OR = 139.9, P = 3.1 × 10-10). We also identified 10 patients in additional influenza patient cohorts. Autoantibodies neutralizing type I IFNs account for ∼5% of cases of life-threatening influenza pneumonia in patients <70 yr 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.

Vaccine breakthrough hypoxemic COVID-19 pneumonia in patients with auto-Abs neutralizing type I IFNs.

Life-threatening 'breakthrough' cases of critical COVID-19 are attributed to poor or waning antibody response to the SARS-CoV-2 vaccine in individuals already at risk. Pre-existing autoantibodies (auto-Abs) neutralizing type I IFNs underlie at least 15% of critical COVID-19 pneumonia cases in unvaccinated individuals; however, their contribution to hypoxemic breakthrough cases in vaccinated people remains unknown. Here, we studied a cohort of 48 individuals (age 20-86 years) who received 2 doses of an mRNA vaccine and developed a breakthrough infection with hypoxemic COVID-19 pneumonia 2 weeks to 4 months later. Antibody levels to the vaccine, neutralization of the virus, and auto-Abs to type I IFNs were measured in the plasma. Forty-two individuals had no known deficiency of B cell immunity and a normal antibody response to the vaccine. Among them, ten (24%) had auto-Abs neutralizing type I IFNs (aged 43-86 years). Eight of these ten patients had auto-Abs neutralizing both IFN-α2 and IFN-ω, while two neutralized IFN-ω only. No patient neutralized IFN-ß. Seven neutralized 10 ng/mL of type I IFNs, and three 100 pg/mL only. Seven patients neutralized SARS-CoV-2 D614G and the Delta variant (B.1.617.2) efficiently, while one patient neutralized Delta slightly less efficiently. Two of the three patients neutralizing only 100 pg/mL of type I IFNs neutralized both D61G and Delta less efficiently. Despite two mRNA vaccine inoculations and the presence of circulating antibodies capable of neutralizing SARS-CoV-2, auto-Abs neutralizing type I IFNs may underlie a significant proportion of hypoxemic COVID-19 pneumonia cases, highlighting the importance of this particularly vulnerable population.

Evolutionary genetics evidence of an essential, nonredundant role of the IFN-γ pathway in protective immunity.

Identifying how natural selection has affected immunity-related genes can provide insights into the mechanisms that have been crucial for our survival against infection. Rare disorders of either chain of the IFN-γ receptor, but not of IFN-γ itself, have been shown to confer predisposition to mycobacterial disease in patients otherwise normally resistant to most viruses. Here, we defined the levels of naturally occurring variation in the three specific genes controlling the IFN-γ pathway (IFNG, IFNGR1, IFNGR2) and assessed whether and how natural selection has acted on them. To this end, we resequenced the three genes in 186 individuals from sub-Saharan Africa, Europe, and East-Asia. Our results show that IFNG is subject to strong purifying selection against nonsynonymous variants. Conversely, IFNGR1 and IFNGR2 evolve under more relaxed selective constraints, although they are not completely free to accumulate amino acid variation having a major impact on protein function. In addition, we have identified signatures of population-specific positive selection, including at one intronic variant known to be associated with higher production of IFN-γ. The integration of our population genetic data into a clinical framework demonstrates that the IFN-γ pathway is essential and nonredundant in host defense, probably because of its role in protective immunity against mycobacteria.