Interface Gain-of-Function Mutations in TLR7 Cause Systemic and Neuro-inflammatory Disease

TLR7 recognizes pathogen-derived single-stranded RNA (ssRNA), a function integral to the innate immune response to viral infection. Notably, TLR7 can also recognize self-derived ssRNA, with gain-of-function mutations in human TLR7 recently identified to cause both early-onset systemic lupus erythematosus (SLE) and neuromyelitis optica. Here, we describe two novel mutations in TLR7, F507S and L528I. While the L528I substitution arose de novo, the F507S mutation was present in three individuals from the same family, including a severely affected male, notably given that the TLR7 gene is situated on the X chromosome and that all other cases so far described have been female. The observation of mutations at residues 507 and 528 of TLR7 indicates the importance of the TLR7 dimerization interface in maintaining immune homeostasis, where we predict that altered homo-dimerization enhances TLR7 signaling. Finally, while mutations in TLR7 can result in SLE-like disease, our data suggest a broader phenotypic spectrum associated with TLR7 gain-of-function, including significant neurological involvement. Supplementary Information The online version contains supplementary material available at 10.1007/s10875-024-01660-6.


Introduction
Systemic lupus erythematosus (SLE) describes a heterogeneous set of clinical phenotypes associated with type I interferon (IFN) signaling upregulation and the presence of autoantibodies targeting nuclear autoantigens [1].Familial aggregation and higher concordance rates between monozygotic versus dizygotic twins suggest a major hereditary component, with rare monogenic forms of SLE providing important insights into disease pathogenesis [2,3].Expressed in sentinel cells such as macrophages and dendritic cells, Toll-like receptors (TLRs) are a family of single-pass membrane-spanning proteins that engage structurally conserved microbial features as part of a coordinated innate and adaptive immune response to pathogens.Endosomal TLR3, TLR7/8, and TLR9 traffic from the ER to the endosomal compartment where they sense nucleic acids, respectively, dsRNA, ssRNA, and CpG DNA.Inherited TLR3 deficiency underlies herpes simplex encephalitis [4] and severe influenza or COVID-19 pneumonia [5], and inherited TLR7 deficiency predisposes to critical COVID-19 pneumonia [6].Importantly, TLR7/8 and TLR9 can signal upon sensing both viral-and self-derived nucleic acid [7].Reflective of the latter situation, disease in murine models of SLE is attenuated in animals deficient for TLR7 [8], while SLE-like pathology is driven by TLR7 overexpression [9].Further, gain-of-function (GOF) mutations in human TLR7 were recently identified to cause severe autoimmune phenotypes, specifically, SLE and neuromyelitis optica [10].These observations indicate the importance of endosomal TLR signaling in immunological homeostasis.Here, we describe two novel mutations in TLR7, causing a spectrum of systemic and neuro-inflammatory diseases.

Samples Obtained from Patients
Samples were obtained from the probands and parents with written informed consent.The study was approved by the Comité de Protection des Personnes (ID-RCB/EUDRACT: 2014-A01017-40, revalidated in 2022).

Genetic Studies
DNA was extracted from whole blood using standard methods.In family AGS571, a Twist Custom Panel Kit (Twist Bioscience), comprising around 8000 disease-causing genes, was used for library preparation, and a NovaSeq6000 platform (Illumina) for sequencing.In family AGS3740, exome sequencing was performed on genomic DNA using a Sure-Select Human All Exon Kit (Agilent Technologies) for targeted enrichment and Illumina HiSeq2000 for sequencing.Variants were assessed using the in silico programs SIFT (http:// sift.jcvi.org) and Polyphen2 (http:// genet ics.bwh.harva rd.edu/ pph2/).Population allele frequencies were obtained from the gnomAD database (http:// gnomad.broad insti tute.org).Sanger sequencing was performed to confirm the TLR7 variant identified in family AGS3740.The reference sequence used for primer design and nucleotide numbering was TLR7: ENST00000380659.4/NM_016562.4;NP_057646.1.

Protein Structure Modeling
We searched the Protein Data Bank [11] on 16/11/2023 for protein structures that have > 90% homology with TLR7 (UniProt ID: Q9NYK1) over a stretch of at least 50 amino acids.Twenty-three structures were analyzed with FoldX 5.0 [12] to estimate the change in Gibbs free energy (ΔΔG) of the missense variants as the average of 3 replicates.The FoldX "RepairPDB" command was run before modeling the mutations.The buried surface area was defined as the solvent-accessible surface area between the monomer and the complex and was calculated at the residue level using AREAIMOL from the CCP4 software suite [13].Structure visualization was performed with PyMOL [14].

Plasmids
pCMV6 vector encoding wild-type (WT) TLR7 was used as the parental vector for mutagenesis [6].Mutant plasmids of TLR7 were generated via site-directed mutagenesis using the Q5 kit (E0554S; New England Biolabs) according to the manufacturer's instructions and using oligonucleotides listed in the table below.NEB 5-alpha competent Escherichia coli were transformed with the ligated product, and colonies were screened by sequencing for the presence of the desired variants.pcDNA 3.1™ vector encoding V5-tagged human WT UNC93B1 has been described previously [5].

Luciferase Reporter Assay
For each condition, 2.10 4 HEK293T cells per well were plated in a 96-well plate in duplicate and co-transfected with a plasmid containing the Firefly luciferase gene under the control of the human NF-κB promoter (pGL4.32,Promega), a plasmid constitutively expressing the Renilla luciferase gene for normalization (pRL-SV40, Promega), as well as a plasmid encoding WT, empty vector (EV) or variant TLR7 and a plasmid encoding UNC93B1, using the TransIT®-293 Transfection Reagent (Mirus, # MIR2700) according to the manufacturer's instructions.After incubation for 24 h, cells were left either unstimulated or stimulated with the TLR7 agonist R848 0.01, 0.1, and 1 µg/mL for 24 h.Cells were then lysed, and luciferase levels were measured with the Dual-Luciferase® Reporter assay system (Promega, #E1980) according to the manufacturer's protocol.Luminescence intensity was acquired on an Infinite® F200 Pro microplate reader (TECAN).
Firefly luciferase activity values were normalized against Renilla luciferase activity before further data processing (see figure legend).

Family AGS571
The proband (II:1) was a female born to non-consanguineous parents of Italian ancestry (Fig. 1A).She was well until age 4 years, with normal physical and neurological development (Table 1, Table S1).At that age, she presented with fever, malar rash, and cytopenia (lymphopenia and anemia).
One year later, she began to experience severe headaches and hallucinations.She was subsequently diagnosed with extensive cerebral vasculitis and suffered a cerebrovascular infarct, considered secondary to high positivity for antiphospholipid antibodies, resulting in a major neurological deficit.Widespread micro-calcifications were noted at that time on cranial computed tomography.Titers of antinuclear antibodies (ANA) and anti-dsDNA were elevated, and there was C3 and C4 hypocomplementemia.She subsequently developed a class III/IV lupus nephritis and intestinal ischemia.She was treated with steroids, hydroxychloroquine, mycophenolate mofetil, and fondaparinux (the latter as secondary thromboprophylaxis).During follow-up, she continued to show low levels of complement, but autoantibody screening was negative.She died of an acute myocardial infarction at the age of 17 years.After a term pregnancy and normal delivery (birth weight 2.63 kg, head circumference 49 cm), the proband's brother  The patient has no dysphagia or respiratory problems.He exhibits a dystonic movement disorder.His current neurological state appears to be stable.The mother of II:1 and II:2 had a history of malar rash and ANA positivity identified at age 12 years.At 24 years of age, she was diagnosed with immune thrombocytopenic purpura and treated with steroids for 2 years.Aged 44 years she experienced a major arterial thrombotic event in the right lower limb, at which time proteinuria was noted.Immunological work-up showed low levels of complement in the absence of ANA, ENA, and antiphospholipid autoantibodies.
Interferon signaling in whole blood tested in the proband at age 7 years was elevated, with an interferon score of 14.416 (normal < 2.446) (Fig. S1A).Further testing in all three affected individuals was also grossly elevated (Fig. S1B).

Family AGS3740
This now 11-year-old female (II:1) was born to non-consanguineous parents of Italian ancestry with no family history of note (Fig. 1B).She was delivered vaginally at 36 weeks gestation weighing 3.3 kg.During infancy, she exhibited mild motor developmental delay, achieving complete head control at 5 months, independent sitting at 8 months, and independent walking at 18 months.At 1 year of age, she developed immune thrombocytopenia refractory to multiple therapies (steroids, intravenous immunoglobulins, mycophenolate mofetil, rituximab, eltrombopag) except sirolimus and, later, autoimmune hemolytic anemia (negative Coombs reaction, positive anti-SSA), leading to a diagnosis of Evanslike/SLE-like syndrome.
At age 8 years, she was referred to pediatric neurology because of a stereotypic movement disorder (i.e., head nodding), present since age 18 months.Cranial MRI revealed multiple hyperintense foci in cortico-subcortical areas, mainly parasagittal in the frontal and parietal lobes, also involving the temporal lobes and periventricular white matter (Fig. 1D).Serial brain MRIs from 8 years of age, including after transplantation (Fig. 1E), demonstrated a progressive, diffuse leukoencephalopathy of the deep white matter, associated with brain atrophy and both supra-and infra-tentorial scattered calcifications (stone-like in the basal ganglia), in the absence of a clear worsening of neurological status.
By the age of 10 years, both hematological conditions had become refractory to all previous therapies (including the addition of bortezomib), and life-threatening, with platelet counts consistently < 10 × 10 9 /L and hemoglobin < 7 g/dl, prompting peripheral blood hematopoietic stem cell transplantation (HSCT) from an HLA-matched unrelated donor.The conditioning regimen included thiotepa, treosulfan, fludarabine, and rabbit anti-T lymphocyte globulin and graft-versus-host disease (GVHD) prophylaxis comprising cyclosporine A and mycophenolate mofetil.HSCT was successful, with prompt achievement of 100% donor chimerism in peripheral blood mononuclear cells as early as day + 10 after HSCT.Neutrophil engraftment (first of 3 consecutive days with neutrophils > 0.5 × 10 9 /L) was achieved on day + 11, while platelet engraftment (first of 3 consecutive days with platelets > 50 × 10 9 /L) was achieved on day + 105 post-HSCT.Grade I acute cutaneous GVHD developed on day + 29 but was completely resolved with topical steroids.Six months after transplantation, the patient developed skin maculopathy and panniculitis, both thought to be due to limited chronic GVHD.Histological confirmation was obtained 10 months after HSCT.Initial treatment included imatinib and rituximab, subsequently changed to ruxolitinib, started 16 months after HSCT.At the last hematological examination, 22 months after HSCT, the patient had 100% donor chimerism and a completely normal blood count.After 6 months of treatment with ruxolitinib, an initial improvement in chronic cutaneous GVHD was observed.
On the last neurological examination, at the age of 11 years and 10 months, there was a slight reduction in muscle strength and processing speed and some hypokinesia.
The same residue was reported to be mutated to a leucine (F507L) in Brown et al. [10] (again, maternally inherited).

In Vitro Data
To investigate the effect of these two missense substitutions on TLR signaling, we expressed WT and variant (c.1520 T > C (F507S) and c.1582C > A (L528I)) TLR7 in HEK293T cells, which lack TLRs, using an NF-κB luciferase reporter as the read-out.We measured reporter activity following stimulation with the TLR7 agonist R848.Both the F507S and the L528I variants were found to confer a gainof-function in TLR7 activity (Fig. 3, Fig. S3), similar to two recently described gain-of-function mutations (F507L and Y264H) in TLR7 [10].

Protein Structure Modeling
The F507, L528, and Y264 residues all lie at the endosomal leucine-rich repeat (LRR) domain of TLR7 (Fig. 4A).We used a meta-structural analysis approach to explore the likely structural impacts of missense variants by modeling their energetic effects (ΔΔG) across 23 different experimentally determined TLR7 structures, considering both the stability of the monomeric protein and their intermolecular effects on the TLR7 homodimer (Fig. 4B).All variants were predicted to have minimal effects on protein stability, based on previous observations suggesting ~ 1.5 kcal/mol as an approximate ΔΔG threshold above which mutations are less likely to be tolerated in monomers [15,16].In contrast, when considering intermolecular effects, L528I is predicted to be highly disruptive to the TLR7 dimerization interface, with much smaller intermolecular effects predicted for the other variants.However, despite the small predicted intermolecular effects of these variants, it is notable that the F507, L528, and Y264 residues are situated directly at the homodimeric interface, compared to < 10% (35/365) of missense variants recorded in gnomAD v4 (P = 0.0001, Fisher's exact test).Moreover, the residues associated with pathogenic variants tend to bury significantly more intermolecular surface area  3 In vitro assessment of TLR7 signaling.NF-κB luciferase reporter activity following co-transfection of HEK293T cells with empty vector (EV), wild-type (WT) TLR7 or variants of TLR7, and UNC93B1 plasmid and stimulation with R848 0.1 µg/mL.Data are normalized to WT response to R848.Mean ± SEM of n = 3 experiments.Two-way ANOVA with Dunnett's post hoc test: **p < 0.01, *p < 0.5.NS: non-stimulated in the dimer interface than the sites of the gnomAD variants (Fig. 4C).Given that dimerization of TLR7 is necessary for its activation [17], this suggests that intermolecular effects are crucial for pathogenicity.

Discussion
TLRs are type I transmembrane proteins with an extracellular LRR domain, a single transmembrane helix, and a cytoplasmic Toll/interleukin-1 receptor (TIR) domain [18].On recognition of ligand by the LRR domain, inactive monomeric or preformed dimeric TLRs are induced to form a face-to-face activated dimer [19].In this way, TIR domains are brought into proximity with each other, thereby recruiting downstream adaptor proteins to initiate transcriptionally-driven immune responses.TLR7 contains two distinct ligand-binding sites, recognizing ssRNAs in the form of degradation products, nucleosides, and oligoribonucleotides [20].Site 1, which is highly conserved between TLR7 and TLR8, recognizes both nucleosides (guanosine for TLR7, uridine for TLR8) and nucleoside analogs and is essential for receptor dimerization.In contrast, site 2 is not conserved, is spatially distinct between TLR7 and TLR8, and appears to play an auxiliary role in receptor dimerization by enhancing the binding affinities of site 1 ligands.
Here, we describe two novel mutations in TLR7, F507S and L528I.While the L528I mutation is predicted to impact intermolecular contacts at the TLR7 dimer interface, this does not necessarily mean that dimerization will be completely disrupted.Notably, there is considerable variability in the precise orientation of the TLR7 dimer in the available experimental structures, suggesting that it might be robust to a single amino acid change.Comparable to the F507L mutation described by Brown et al. [10], the F507S substitution is predicted to have a minimal effect on monomeric and dimeric TLR7 structures.Neither of the L528 and F507 residues are in direct contact with the guanosine (or other RNA ligand) binding site in any structure.Brown et al. described a Y264H substitution to map to ligand-binding site 1 [10], suggesting an enhanced affinity of the Y264H for guanosine.However, since Y264 lies at the TLR7 dimer interface, it remains possible that this substitution might also play a role in homo-dimerization.We highlight this point by noting that four (F507L, Y264H, F507S, L528I) of the five mutations in human TLR7 described to date, not including R28G (possibly a signal peptide-related variant: Brown et al. [10]), occur directly at the homodimer interface in multiple structures.Of possible further interest, although TLR7-TLR7 molecules are in contact with one another, no noticeable interactions (hydrogen bond, intimate contact, etc.), that might be important in dimer stabilization, were observed in the structure reported by Ishida et al. [17], indicating that dimerization is loosely achieved (and perhaps partly explaining why the isolated endosomal domain of TLR7 mainly exists as a monomer in solution [21]).All told, these observations suggest the importance of the TLR7 interface in maintaining immune homeostasis, where we predict that altered homo-dimerization enhances TLR7 signaling.It will be interesting to map the position of further mutations in human TLR7 as they are reported.
While GOF mutations in TLR7 can result in a phenotype consistent with SLE, our data emphasize a broader spectrum of disease, including significant neurological involvement.Thus, the child with a de novo Y264H mutation described by Brown et al. [10] demonstrated recurrent hemichorea requiring treatment with haloperidol (Table S2), and in a second case reported by the same authors, a child with a maternally inherited F507L mutation, the phenotype was exclusively neurological, consistent with neuromyelitis optica and positive for AQP4 autoantibodies (a phenomenon recognized in the context of enhanced type I interferon signaling [22]).In our series of four affected patients, we observed cerebral vasculitis and cerebral calcification in the female proband in family AGS571 and developmental delay, epilepsy, and marked intracranial calcification in her brother.Further, white matter disease, atrophy, and calcification were seen in the proband in family AGS3740, which was progressive despite an apparently excellent response of life-threatening immune cytopenias to HSCT.
Finally, given that the TLR7 gene is situated on the X chromosome, the observation of a severely affected male in our series is notable, with all other cases so far described being female (and where variable levels of X-inactivation might affect the clinical phenotype).
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Fig. 1
Fig. 1 Genetic and clinical data.Family pedigrees (A AGS571; B AGS3740) where an affected individual harbors a rare non-synonymous missense substitution in TLR7, in either the hemizygous (male) or heterozygous (female) state.Circles and squares indicate female and male family members, respectively.Filled shapes indicate affected status.The line across II:1 in AGS571 indicates deceased

Fig. 4
Fig. 4 Modeling of TLR7 and assessment of identified variants.A Top view and side views of TLR7 structure based upon PDB ID 7cyn, with the location of indicated amino acids (variants at F507 and L528 found in this study) highlighted.B Predicted change in the Gibbs free energy of folding (ΔΔG) of TLR7 variants measured in monomeric TLR7 structures (left).Predicted intermolecular ΔΔG defined as the difference between the ΔΔG of the TLR7 monomer

Table 1
Patients harboring heterozygous TLR7 missense substitutions identified in this study AAB autoantibody, CVA cerebrovascular accident, DD developmental delay, HSCT hematopoietic stem cell transplantation, LN lupus nephritis, MI myocardial infarction, SLE systemic lupus erythematosus, WMD white matter disease a The denominator (112,190) is derived from a p.Phe507LeufsTer16 recorded on gnomAD.Further, note that a p.Ile505Thr and a p.Lys509Glu are each seen once on gnomAD out of 457,192 and 456,952 alleles respectively b The denominator (640,595) is derived from a p.Leu528Leu recorded on gnomAD Currently, the boy is 3 years of age and cannot sit, walk or talk.Visual evoked potentials showed a bilateral increase in cortical latency with a normal auditory brainstem response.