Human Germline Heterozygous Gain-of-Function STAT6 Variants Cause Severe Allergic Disease

Mehul Sharma, Henry Y. Lu, Maryam Vaseghi-Shanjani, Kate L. Del Bel, Oriol Fornes, Robin van der Lee, Phillip A. Richmond, Susan Lin, Joshua Dalmann, Jessica J. Lee, Allison Matthews, Géraldine Blanchard-Rohner, Clara D M van Karnebeek, H. Melanie Bedford, Wyeth W. Wasserman, Michael Seear, Margaret L. McKinnon, Hanan Ahmed, Stuart E. Turvey
doi: https://doi.org/10.1101/2022.04.25.22274265

ABSTRACT

STAT6 (Signal transducer and activator of transcription 6) is a transcription factor that plays a central role in the pathophysiology of allergic inflammation. STAT6 mediates the biological effects of IL-4, a cytokine necessary for type 2 differentiation of T cells and B cell survival, proliferation and class switching to IgE. We have identificated two unrelated patients with a phenotype notable for their early-life onset of profound allergic immune dysregulation, widespread treatment-resistant atopic dermatitis, hypereosinophilia with esosinophilic esophagitis, elevated serum IgE, IgE-mediated food allergies, and vascular anomalies of the brain. Both patients harbored heterozygous de novo missense variants in the DNA binding domain of STAT6 (c.1144G>C, p.E382Q; and c.1256A>G, p.D419G). Functional studies established that both variants caused a gain-of-function (GOF) phenotype associated with enhanced phosphorylation and transcriptional activity of STAT6, in addition to increased transcript abundance of known STAT6 target genes and other genes implicated in allergic disease. JAK inhibitors decreased the enhanced STAT6 responses associated with both these STAT6 GOF variants. This study identifies heterozygous GOF variants in STAT6 as a novel autosomal dominant allergic disorder. We anticipate that our discovery of the first humans with germline STAT6 GOF variants will facilitate the recognition of more affected individuals and the full definition of this new primary atopic disorder.

INTRODUCTION

Signal transducer and activator of transcription (STAT) proteins are a family of seven mammalian transcription factors (STAT1, STAT2, STAT3, STAT4, STAT5A/B, and STAT6) responsible for regulating the expression of genes involved in cellular immunity, survival, proliferation, and differentiation1,2. Although they each have unique functions, they share sequence homology and possess similar activation mechanisms3. In general, upon extracellular cytokine/growth factor-receptor engagement, intracellular membrane receptor-associated Janus kinases (JAK) undergo transphosphorylation and activation, leading to receptor tyrosine phosphorylation4. Phosphorylated tyrosine residues on the receptor then create binding sites for the normally cytosolic-located Src homology 2 (SH2) domain-containing STATs to bind and become phosphorylated and activated. Phosphorylated STATs, dimerize and translocate into the nucleus, where they bind to specific DNA sequences and regulate transcription4-7. As such, STATs are critical for transducing extracellular cytokine-induced responses to transcriptional changes in the nucleus.

Signal transducer and activator of transcription 6 (STAT6) is the main transcription factor that mediates the biological effects of IL-4, a cytokine necessary for type 2 differentiation of T cells, as well as B cell survival, proliferation and class switching to IgE8-11. IL-4 is predominantly secreted by T cells, mast cells, basophils, and type 2 innate lymphoid cells12. To mediate its effects, IL-4 binds the widely expressed IL-4Rα and this complex then binds secondary receptors chains13. In hematopoietic cells, the secondary receptor chain is the common γc chain also used by IL-2, IL-7, IL-9 and IL-15 and this forms the type 1 IL-4R system14. In non-hematopoietic cells such as epithelial cells, the secondary receptor chain is IL-13Rα, which forms the type 2 IL-4R system13. Activation of these two receptor systems lead to diverse functions. In B cells, IL-4-STAT6 signalling supports the growth and differentiation of B cells as well as secretion of pro-allergic immunoglobulin E (IgE)9. In T cells, the IL-4-STAT6 axis is critical for the differentiation of T helper 2 (TH2) cells15, a subset of CD4+ helper T cells that is a major contributor to the pathogenesis of allergic disease. In epithelial cells, the IL-4/IL-13-STAT6 pathway has been implicated in barrier dysfunction16,17 and chronic pruritus18. Collectively, the IL-4-STAT6 axis drives the development and progression of allergic disease and asthma. Indeed, a monoclonal antibody directed against IL-4Rα (dupilumab) is clinically approved for treating allergic inflammation19,20. Mouse studies further emphasize the central role that STAT6 plays in the allergic diathesis. Mice lacking STAT6 fail to develop many of the features of allergic inflammation21-23.

In 2018, the term primary atopic disorders (PADs) was coined to describe heritable genetic disorders presenting with dysregulated pathogenic allergic effector responses24-26. PADs are a subgroup of inborn errors of immunity (IEIs) that manifest with prominent allergic inflammation, in the presence or absence of other clinical features associated with IEIs, such as enhanced susceptibility to infections, autoimmunity, and malignancy. Of the ∼40 monogenic causes of PADs described to date, 12 impact the JAK-STAT signaling cascade that is critical for cytokine responsiveness. For example, both loss-of-function (LOF) and gain-of-function (GOF) variants in STAT327,28 and STAT5B29,30, and GOF variants in JAK131,32 and STAT133 cause PADs. Furthermore, defects in the genes encoding upstream and downstream proteins in the STAT pathway (e.g. ZNF341, DOCK8, IL6ST, IL6R, and ERBIN) also result in severe allergic diseases26. The recent increased recognition of PADs is valuable because: (i) the identification of affected pathways provide important insights into the pathogenesis of allergic inflammation and may uncover new treatment targets; (ii) these conditions are likely under-diagnosed since allergies are traditionally believed to be complex polygenic disorders; and (iii) diagnosis of PADs can be transformative for affected individuals and their families by guiding treatment options, and informing relevant considerations tied to genetic counselling including family planning, long term prognosis, and connection to other families with the same genetic disease.

In this study, we describe a novel human PAD caused by germline heterozygous GOF variants in the DNA binding domain of STAT6 found in two unrelated individuals with severe allergic disease and vascular anomalies of the brain.

METHODS

Ethical Considerations

All study participants provided written informed consent. Research study protocols were approved by The University of British Columbia Clinical Research Ethics Board.

Identification of STAT6 variant via whole exome sequencing

Trio whole exome sequencing (WES) was performed as previously described34. Briefly, WES was performed on genomic DNA extracted from the first patient (P1) and both his parents and a novel de novo predicted damaging STAT6 variant (c.1144G>C p.E382Q) was selected for further analysis. WES was also performed on genomic DNA extracted from the second patient (P2), her mother, and one of her brothers. A novel de novo predicted damaging variant for STAT6 (c.1256A>G, p.D419G) was identified. Sanger sequencing on saliva samples collected from the remaining family members of P2 revealed that no other individuals carried the variant and that it segregates with disease.

Generation of STAT6 variant plasmids

Plasmids used for transfection studies contained full-length STAT6 in a pCMV6 entry vector with a C-terminal GFP tag (Cat#: RG210065, OriGene Technologies; Rockville, Maryland, USA). To generate p.E382Q and p.D419G STAT6, a Q5 site-directed mutagenesis kit (Cat# E0554S, New England Biolabs, Ipswich, MA, USA) was used according to manufacturer’s recommendations along with the following primer pairs: 5’-TGTCACAGAGCAGAAGTGCGC-3’, 5’-GACTCAGTGCCCTTCCGC-3’, and 5’-GGCAACCAAGGCAACAATGCC-3’, 5’-ATGGACGATGACCACCAG-3’, respectively.

To generate lentivirus vectors, wild type (WT), p.E382Q, and p.D419G STAT6 from the above plasmids were cloned into a GFP-tagged Lenti vector (Cat#: PS100071, OriGene Technologies, Rockville, MD, United States) using EcoRI-HF (Cat#: R3101) and NotI-HF (Cat#: R3189) both from New England BioLabs. The three STAT6 plasmids were packaged using 3rd generation packaging plasmids and transfected into HEK293T cells. Culture media was collected, centrifuged, filtered, concentrated, and stored at - 80°C before use.

All new variant plasmids were confirmed by Sanger sequencing and purified from 10-beta competent E. coli using a QIAprep Spin Miniprep Kit (Qiagen, ON, Canada).

Transient and stable expression of STAT6 variants

Transient expression of STAT6 variants in HEK293 cells were accomplished using a Lipofectamine 3000 kit (Thermo Fisher Scientific, Waltham, MA, USA) according to manufacturer’s recommendations. Briefly, HEK293 cells were seeded at 8.0 × 105 cells/well in a 6-well plate in 1.5 mL of Dulbecco modified Eagle medium (DMEM) with 10% FBS (Gibco, Life Technologies; Rockville, MD, USA) and incubated for 24 hours at 37°C. Cells were transfected with 2.5μg of plasmid DNA using the P3000 and Lipofectamine 3000 reagents and harvested after 24h.

Stable expression of STAT6 in Jurkat T cells was accomplished using the Lentivirus approach as previously described35,36. Briefly, Jurkat T cells were infected with lentiviral particles in the presence of 5μg/ml polybrene (Sigma-Aldrich, St Louis, MO, USA) and spinoculated at 800 x g for 30mins at 32°C, cultured, and expanded in complete RPMI-1640 (GE Healthcare, Chicago, IL, United States) supplemented with 10% FBS. Expanded cells were sorted on GFP expression using a BD FACS Aria (BD Biosciences) cell sorter.

Determining STAT6 variant activity by phospho-flow cytometry

Phospho-STAT6 activity was quantified by flow cytometry in transfected cells stimulated with 100ng/mL of IL-4 (Cat# 204-IL-020, R & D Systems, Minneapolis, MN, USA) for 1h or transfected cells pre-treated with either 10μM ruxolitinib (Cat code: tlrl-rux, Invivogen, San Diego, CA, USA) or 4μM tofacitinib (Cat# S50001, Selleckchem, Houston, TX, USA) for 2h before stimulation as previously described37. Briefly, transfected and stimulated cells were fixed using BD Cytofix (Cat# 554655, BD Biosciences, Ontario, Canada) for 20 min at 4°C and permeabilized using Perm III for 30 min on ice (Cat# 558050 BD Biosciences). The cells were then stained with STAT6 PE-CF594 (Cat# 564148, BD Biosciences) and p-STAT6 AF647 (Cat# 612601, BD Biosciences). pSTAT6 expression was measured in GFP+STAT6+ cells from WT and STAT6 variants samples on an LSRII flow cytometer (BD Biosciences) and analyzed using FlowJo software (BD Biosciences).

Luciferase reporter assays

A luciferase reporter plasmid encoding a 4x STAT6 binding site [TTCCCAAGAA] was used to assess WT and variant STAT6 promoter activity38. The p4xSTAT6-Luc2P plasmid was a gift from Axel Nohturfft (Addgene plasmid #35554; http://n2t.net/addgene:35554; RRID:Addgene_35554). Briefly, HEK293 cells were seeded in 24-well plates overnight at a density of 150k cells and transfected with 250ng of p4xSTAT6-Luc2P, 250ng of a plasmid encoding GFP-tagged WT, p.E382Q, or p.D419G STAT6, and 10ng Renilla luciferase (R-Luc) using lipofectamine 3000 as described above. Transfected cells were subsequently stimulated with 100ng/mL of IL-4 (Cat# 6507-IL-010, R&D Systems, Minneapolis, MN, USA) or left untreated. Cell lysates were prepared and processed using a Dual-Glo Luciferase Assay Kit (Cat#: E2920, Promega, Madison, WI, USA) according to manufacturer’s recommendations. Luciferase activity was measured on the Infinite M200 plate reader (Tecan; Männedorf, Switzerland).

RNA sequencing

To investigate the global transcriptome, STAT6 variant-transduced Jurkat T cells were left unstimulated, or stimulated with 100ng/mL IL-4 for 4h. RNA was extracted in triplicate as previously described37 using a RNeasy Mini Plus Kit (Qiagen) according to manufacturer’s recommendations. RNA was prepared following the standard protocol for the NEBNext Ultra II Stranded mRNA (New England Biolabs) and sequenced on the Illumina NextSeq 500 with Paired End 42 bp × 42 bp reads. De-multiplexed read sequences were aligned to a reference sequence using RNA-Seq Alignment app (v1.1.1) on Illumina Basespace, using Spliced Transcripts Alignment to a Reference (STAR) aligner and Cufflinks 2 for assembly and estimation of gene expression.

Expression data were normalized to reads between samples using the edgeR package in R (R Foundation, Vienna, Austria). Normalized counts were filtered to remove low counts using the filterByExpr function in edgeR39. Principal component analysis (PCA) was done on log2(normalized counts+0.25) in R using the PCA function. Differential expression between unstimulated and stimulated samples for all three STAT6-transduced Jurkat T cells was accomplished using Limma40. Differentially expressed genes were defined as those with fold change (FC) greater than 1.25 and adjusted p-value less than 0.05.

Pathway analysis was done by first performing Gene Set Enrichment Analysis (GSEA) with 1000 permutations using the Molecular Signatures Database Hallmark module. Signal-to-noise ratio was used for gene ranking and the obtained NESs and P-values were further adjusted using the Benjamini-Hochberg method. Pathways with an adjusted P-value < 0.05 were considered significant. Leading edge genes from significant pathways between STAT6 WT and STAT6 c.1256A>G Jurkat were identified (no significantly upregulated pathways were identified between STAT6 Wild Type and STAT6 c.1144G>C). Expression levels of these genes were then determined in each of the three groups (Wild Type, c.1144G>C, c.1256A>G) under both stimulated and unstimulated conditions. Sample level enrichment analyses (SLEA) scores were computed as previously described41. Briefly, z-scores were computed for gene sets of interest for each sample. The mean expression levels of significant genes were compared to the expression of 1000 random gene sets of the same size. The difference between observed and expected mean expression was then calculated and represented on heatmaps.

Gene set enrichment analysis was also done on two other gene sets: i) STAT6 targets; and ii) IL-4 TH2 targets. STAT6 targets were defined as the set of genes that were significantly upregulated upon IL-4 stimulation in WT-STAT6 transduced Jurkat T cells. Enrichment of this gene set was determined at both baseline and after IL-4 stimulation between WT and p.E382Q STAT6, and WT and p.D419G STAT6 Jurkat T cells. IL-4-STAT6 TH2 skewing targets were previously reported42 STAT6 target genes that are significantly upregulated in response to IL-4 treatment and that lead to TH2 skewing.

Raw data will be deposited on Gene Expression Omnibus.

Structural modeling

The effects of the STAT6 variants on the protein function and structure were analyzed using three-dimensional models. We used SWISS-MODEL43 to model the variants based on a template structure of the human STAT6 transcription factor solved as a homodimer and in complex with DNA (PDB: 4Y5W, resolution: 3.1 Å, chains A, C, M and N)38. Structures were visualized with UCSF Chimera44.

RESULTS

Clinical features of two patients with severe atopic disease

Patient A-II-1 (P1) was a male patient who had life-long severe allergic disease. The patient was not dysmorphic and was born to unrelated parents. Soon after birth, P1 developed widespread, treatment-resistant atopic dermatitis (eczema), gastroesophageal reflux with episodes of aspiration pneumonia, as well as chronic diarrhea. He underwent Nissen fundoplication a few months later (with a surgical revision in early childhood, age 6-10 years), following which gastrostomy feeding was commenced. Over time, he was also diagnosed with multiple food and drug allergies with anaphylactic reactions (combined with positive skin prick testing or RAST testing) documented in response to peanuts, cows milk protein, bananas, sesame seeds, fish, and eggs. He was fed with extensively hydrolyzed and elemental formulas via gastrostomy with very limited oral feeding for much of his life. He was also diagnosed with asthma. He experienced episodes of secondary bacterial skin infections from which Staphylococcus aureus, group A streptococci and Candida spp. were isolated. Widespread molluscum contagiosum was another notable skin finding. In his second and third decade of life, he had persisting eczema, asthma, and perennial rhinitis which worsened during pollen seasons. Repeated upper GI endoscopies confirmed the diagnosis of eosinophilic esophagitis. During this time, most of his feeding consisted of elemental formula via gastrostomy supplemented with limited oral intake (mainly rice, potatoes, and chicken). He was treated with oral steroids for many years and experienced symptom flares whenever the dose of prednisone was dropped below 10mg/day. He developed side-effects related to corticosteroids including cataracts, osteoporosis, and pathogenic fractures. He had an episode of thromboembolic ischemia of the right first and third toes with pain and cyanosis, ultimately requiring a lumbar sympathectomy to relieve the pain. In his mid-30s (age 30-35 years), he was hospitalized for non-bloody diarrhea, profound weight loss, malnutrition, weakness, and vertigo. He was treated with high dose of prednisone, antibiotics, and elemental feeding. During the hospitalization, he died after developing a spontaneous subarachnoid hemorrhage with obstructive hydrocephalus, presumed secondary to an aneurysm (Table 1). Blood testing over the years confirmed eosinophilia and high serum IgE levels (Fig 1). Lymphocyte phenotyping and proliferation to a panel of mitogens were found to be unremarkable when measured in early adulthood (ages 20-25 years). An initial genetic assessment revealed normal sequences for STAT3 and DOCK8, prompting a more extensive genetic evaluation by WES.

View this table:
Table 1: Major Clinical Features.

Tabulation and comparison of the clinical phenotype of patient 1 and 2.

Figure 1:
Figure 1: Clinical phenotype of two patients with severe atopic disease.

A) Family pedigree of the two patients. Filled symbols=affected individual, arrow=proband, unfilled symbols=unaffected individual. P1 has the variant p.E382Q, P2 has the variant p.D419G. B) Eosinophil frequency in whole blood of P1 and P2. C) IgE concentration in whole blood of P1 and P2. D) Schematic illustrating the protein domains of STAT6. Amino acid location of the variants shown in red. Affected region was aligned to other species. Asterisks indicate full conservation. E) Structural model of the DNA-STAT6 homodimer complex. Green=p.E382, Purple=p.D419. Residues marked with an asterisk (*) were demonstrated to affect the binding of STAT6 to DNA in a previous study38. F) STAT6 expression in WT-, p.E382Q-, and p.D419G-STAT6-transfected HEK293 cells. G) Quantification of relative protein expression in F), n=3. Significance determined through a one-way ANOVA and Tukey’s post-hoc test.

Patient B-II-1 (P2) is a female who had similar life-long struggles with severe allergic disease. She presented with severe chronic, treatment-resistant atopic dermatitis since birth. She was diagnosed with IgE-mediated food allergy to egg, milk, peanut and tree nuts early in life (age 0-5 years). She was diagnosed with asthma as an infant. In her early life (age 0-5 years) she also had a pneumothorax requiring chest tube insertion, required intubation for an asthma exacerbation, and was diagnosed with giant papillary conjunctivitis which has remained an ongoing challenge. Other eye manifestations included unilateral retinal detachment in her adolescence (age 15-20 years) and the development of cataracts requiring surgery in early adulthood (age 20-25 years). Through her second and third decades of life she continued to suffer from eczema, asthma, recurrent bronchitis, rhinoconjunctivitis, and secondary staphylococcal and candida skin infections. She had multiple dental abscess in the second and third decades of her life. In her 20s she began experiencing gastrointestinal symptoms culminating in biopsy diagnosed eosinophilic esophagitis and eosinophilic gastroenteritis. Joint examination confirmed the presence of hypermobility. Serial blood testing over the years confirmed eosinophilia and high serum IgE levels (Fig 1). Brain MRI in her mid-30s (ages 30-35 years) revealed multiple anatomical variants including anomalous vasculature in the circle of Willis with hypoplastic vertebrobasilar arteries, a persistent left-sided congenital trigeminal artery with tortuosity, and a hypoplastic A1 segment of the right anterior carotid artery (Table 1). Multiple therapies were trialed over the years with corticosteroids offering the most obvious benefits. To seek a diagnosis for this spectrum of severe allergic manifestations, WES was performed on the patient, her mother and one of her brothers.

Whole-exome sequencing reveals novel de novo heterozygous variants in the DNA binding domain of STAT6

The patients were found to carry heterozygous variants in the gene STAT6 (NM_001178079.2): c.1144G>C, p.E382Q for P1 and c.1256A>G, p.D419G for P2. The variants were de novo as they were not found in any other family members (Fig 1A). Both variants have not been reported in population databases (i.e. gnomAD), but p.D419G can be found in the Catalogue of Somatic Mutations in Cancer (COSMIC) database as a recurrent somatic variant in lymphoma with experimental evidence of a GOF phenotype for this variant45,46. Pathogenicity prediction models predict both variants to be pathogenic, including CADD (Version 1.6, p.D419G: 29.3, p.E382Q: 27.7)47, SIFT (p.D419G: 0.001, p.E382Q: 0.026)48, and Polyphen-2 (p.D419G: 1.0, p.E382Q: 1.0)49 (Supp Table 1). Clinical laboratory tests showed that both patients had eosinophilia and markedly elevated IgE (Fig 1B-C). Both variants localized to the highly conserved DNA binding domain of STAT6 (Fig 1D). Glu382 and Asp419 lie close to the regions responsible for protein-DNA interaction38. Both of these variants decrease the electro-negativity of the protein near the DNA-binding interface (due to loss of glutamic acid (E) and aspartic acid (D)), and are predicted to enhance STAT6 binding to DNA (Fig 1E). These variants did not affect STAT6 protein expression as measured by flow cytometry in transfected HEK293 cells (Fig 1F-G).

p.E382Q and p.D419G STAT6 variants lead to increased phosphorylation and transcription activity of STAT6 in HEK293 cells

To assess the functional impact of the p.E382Q and p.D419G STAT6 variants, we selected HEK293 cells as our system, as these cells lack endogenous STAT6 but express other components of the IL-4R pathway (Fig 2A)50. HEK293 cells were transfected with wild-type (WT), p.E382Q, or p.D419G STAT6 and the phosphorylation status of STAT6 was quantified by flow cytometry after 24h (Fig 2B). At baseline, we observed significantly increased phosphorylated STAT6 in cells expressing the p.E382Q and p.D419G variants when compared to WT. IL-4 stimulation led to increased phospho-STAT6 in all conditions, with the variant expressing cells maintaining significantly higher phospho-STAT6 than WT (Fig 2C). To investigate whether increased phosphorylation of STAT6 in both variants lead to increased STAT6 transcription factor activity, we conducted luciferase assays with a reporter plasmid containing a canonical 4x STAT6 binding site, TTCCCAAGAA in HEK293 cells38. Both p.E382Q and p.D419G STAT6 variants resulted in significantly higher promoter activity at baseline (Fig 2D). However, there was no significant difference in promoter activity after 4h IL-4 stimulation (Fig 2D). Together these data confirm a GOF phenotype caused by both p.E382Q and p.D419G STAT6 variants.

Figure 2:
Figure 2: p.E382Q and p.D419G STAT6 variants lead to increased STAT6 activity in HEK293 cells.

A) Schematic illustrating classical IL-4-mediated STAT6 activation, dimerization, and phosphorylation. B) Phospho-STAT6 expression in WT-, p.E382Q-, and p.D419G-transfected HEK293 cells before and after treatment with IL-4 (100ng/mL for 1h). C) Time course of phospho-STAT6-positive cells after 100ng/mL IL-4 stimulation in transfected HEK293 cells, n=3. D) Luciferase assay of STAT6 activity on a plasmid containing a 4x STAT6 binding site [TTCCCAAGAA] for WT-, p.E382Q-, and p.D419G STAT6-transfected HEK293 cells before and after stimulation with IL-4 (100ng/mL for 4h), n=3. One-way ANOVA and Tukey’s post-hoc test. *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001.

p.E382Q and p.D419G STAT6 variants exert a gain-of-function phenotype in Jurkat T cells

Given the increased promoter activity of the p.E382Q and p.D419G STAT6 variants, we were interested in evaluating how this may impact the global transcriptome. As transcriptomic studies on HEK293 cells after IL-4 stimulation have been challenging to interpret45, we stably expressed p.E382Q and p.D419G STAT6 by lentivirus transduction in Jurkat T leukemia cells, which express endogenous STAT6 and serve as a model of heterozygosity51. Here, we observed significantly higher phospho-STAT6-positivity at baseline for the p.E382Q-cells but not p.D419G-cells, when compared to WT (Fig 3A). However, there was no significant difference in phospho-STAT6 positivity after IL-4 stimulation (Fig 3A).

Figure 3:
Figure 3: p.E382Q and p.D419G STAT6 variants lead to increased STAT6 activity in Jurkat T cells.

A) Quantification of phospho-STAT6-positive WT-, p.E382Q-, or p.D419G-transduced Jurkat T cells before and after treatment with IL-4 (100ng/mL for 1h), n=3. B) Principal component analysis (PCA) comparing unstimulated and stimulated (100 ng/mL IL-4 for 4 h) WT (green), p.E382Q (blue), p.D419G (purple) STAT6 transduced Jurkat T cells. PC1 and PC2 contribution is shown in brackets. C-H) Normalized counts comparing stimulated WT (green) vs. p.E382Q (blue) or p.D419G (purple), for C) XBP1, D) TGFBR2, E) EPAS1, F) BATF, G) IL4R, H) CISH. I) Heatmap representation of normalized counts of a transcription set defined as IL-4 targets in each STAT6 Jurkat T cells. *adjusted P-val<0.05. J) GSEA plots for curated STAT6 and IL-4-TH2 targets genes, comparing WT vs. either p.E382Q (blue) or p.D419G (purple) at both baseline and after stimulation with IL-4. Normalized enrichment score (NES) and adjusted p-value are shown. K) Sample level enrichment analyses (SLEA) of significantly enriched immune pathways from MSigDB Hallmark in unstimulated and IL-4-stimulated samples, comparing WT vs either p.E382Q or p.D419G. Heatmap is normalized across the rows and shown as relative expression.

Transcriptomically, WT-, p.E382Q-, and p.D419G-transduced cells clustered separately from each other both at baseline and after stimulation with IL-4 (Fig 3B). This confirmed that all variant cells were responsive to IL-4, and identified a unique transcriptomic signature for each. Differential gene expression analysis revealed significantly increased transcript abundance of known STAT6 target genes, including IL4R52, CISH45 and EPAS42, and genes that have been implicated in atopic disease, including XBP153 and TGFBR254 in p.E382Q and p.D419G transduced cells when compared to WT transduced controls (Fig 3C-H).

In response to IL-4 stimulation, we observed 54, 145, and 159 significantly upregulated genes in WT-, p.E382Q-, and p.D419G-transduced Jurkat T cells, respectively (Fig 3I, Supp Fig 2A). Interestingly, p.E382Q and p.D419G had 67 and 80 uniquely increased hits, which did not overlap with WT, nor with each other. On the other hand, p.E382Q and p.D419G had 108 and 54 uniquely downregulated genes (Fig 3I, Supp Fig 2A). This suggests that the altered activity of both p.E382Q and p.D419G is not restricted to enhanced activity on known STAT6 targets alone.

To further investigate the differential gene expression patterns, we performed gene set enrichment analyses (GSEA). We first defined IL-4-STAT6 targets as genes that were significantly upregulated in stimulated WT STAT6-transduced Jurkat T cells (Supp Table 2). We found enrichment of these targets in both p.E382Q and p.D419G Jurkats only at baseline (p.E382Q NES 1.34, adj-pvalue 0.084; p.D419G NES 1.83, adj-pvalue <0.001), but not after IL-4 stimulation. Notably, after IL-4 stimulation both p.E382Q and p.D419G STAT6 Jurkats showed significantly increased enrichment of IL-4 responsive TH2 skewing targets42 (p.E382Q NES 1.56, adj-pvalue <0.001; p.D419G NES 1.52, adj-pvalue 0.007)(Fig 3J).

GSEA using the MSigDB Hallmark gene sets showed no significant difference in enrichment in immunological pathways between WT and p.E382Q STAT6 both at baseline and in response to stimulation (Fig 3K, Supp Table 3). However, p.D419G STAT6 exhibited significantly increased enrichment in proliferation (Myc, G2M, DNA repair, E2F) pathways, consistent with the fact that p.D419G (but not p.E382Q) is reported as a recurrent somatic mutation in lymphoma45,55,56, establishing this variant as oncogenic.

JAK inhibitors reduce enhanced STAT6 phosphorylation induced by p.E382Q and p.D419G STAT6 in HEK293 cells

Having demonstrated that both p.E382Q and p.D419G STAT6 are GOF variants, we tested whether small molecule inhibitors that target the STAT6 pathway could potentially be clinically beneficial. We selected JAK inhibitors, ruxolitinib and tofacitinib, as both of these drugs are already used clinically for treating atopic disease57,58. Both inhibitors effectively decreased the enhanced STAT6 phosphorylation observed at baseline and after IL-4 stimulation in cells expressing the p.E382Q or p.D419G STAT6 variants (Fig 4A-B). This suggests that JAK inhibitors may be promising treatment options for patients with GOF STAT6 variants.

Figure 4:
Figure 4: JAK inhibitor treatment decreases gain-of-function STAT6 activity caused by p.E382Q and p.D419G STAT6.

A) phospho-STAT6 expression in transfected HEK293 cells pre-treated with ruxolitinib (10μM 2hr) or tofacitinib (4μM 2hr), before and after stimulation with IL-4 (100ng/mL 1h). B) Quantification of A), n=3. One-way ANOVA and Tukey’s post-hoc test. *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001.

DISCUSSION

In this study, we present a combination of clinical, molecular, and transcriptional evidence of a new human disorder caused by germline autosomal dominant GOF STAT6 variants in two patients with life-long severe allergic disease. These variants led to increased STAT6 phosphorylation, increased STAT6 target gene expression, and TH2 skewed transcriptional profile. Although the full phenotype of GOF STAT6 variants will only be uncovered through the identification of additional affected individuals, we propose to classify human germline autosomal dominant GOF STAT6 syndrome as a novel primary atopic disorder24-26. Based on our study, possible clinical ‘red flags’ for this new disorder include: (i) early life onset; (ii) peripheral blood eosinophilia, (iii) elevated serum IgE, (iii) widespread, treatment-resistant atopic dermatitis, (iv) multiple food and drug allergies; (v) recurrent skin and respiratory infections; (vi) eosinophilic gastrointestinal disorder, including eosinophilic esophagitis; (vii) allergic rhinoconjunctivitis; and (viii) vascular malformations of the brain.

STAT6 is intimately linked to the biology of allergic inflammation. The central and most studied role of STAT6 is in mediating the biological effects of IL-4, a cytokine necessary for TH2 differentiation, B cell survival, proliferation and class switching to IgE9,42, and driving M2 macrophage polarization59. In T cells, STAT6 activation induces the expression of GATA3, the master regulator of TH2 differentiation, which in turn enhances expression of IL-4, IL-5 and IL-13, cytokines necessary for promoting allergic responses by activating mast cells and eosinophils60. Presence of greater TH2 cell populations, or TH2 cells producing copious amounts of IL-4/IL-5/IL-13, could be a driver of the observed allergic phenotype presented in our patients. Elevated IgE in partnership with mast cells is important for both acute and chronic manifestations of allergic disorders and can be an additional driver of the allergic diathesis61. STAT6 hyperactivation in airway epithelial cells and resident dendritic cells can further create an environment favouring asthma and chronic lung disease, as this would induce production of chemokines that promote TH2 cells and eosinophil recruitment62,63. Population genetics provide further support for the central role that STAT6 plays in the development of human allergic disease. Multiple independent genome wide association studies (GWAS) have found that polymorphisms in STAT6 associate with many allergic conditions (Table 2). Our study further expands this appreciation of the role of STAT6 in human disease by establishing that heterozygous GOF variants cause a monogenic form of severe allergic disease.

View this table:
Table 2. Number of published genome-wide association (GWAS) studies linking polymorphisms (SNPs) in STAT6 to common allergic diseases in the population.

Significant genome-wide associations (P < 5×10−8) between STAT6 SNPs and all relevant allergic traits were obtained through NHGRI-EBI GWAS Catalog (ebi.ac.uk/gwas/).

The fatal cerebral aneurysm in P1 (p.E382Q) was not clinically anticipated, but it is possible that the STAT6 GOF variant also increased the risk of developing cerebral aneurysms. Intracranial aneurysms have been reported in patients with both STAT3 LOF and STAT1 GOF64-66. Increased activation of other STAT family members, including STAT2, STAT3 and STAT5 have also been observed in human abdominal aortic aneurysms (STAT6 was not evaluated), although it is not clear whether enhanced STAT phosphorylation contributes to aneurysms or is the result of inflammation caused by aneurysms67. As more individuals with STAT6 GOF variants are identified, the possibility of cerebral vascular anomalies warrants investigation.

A GOF STAT6 model (designated STAT6VT) has previously been described in vitro68 and has been used to study chronic atopic dermatitis in mouse models69,70. STAT6VT has the substitution of two amino acid residues, at positions 547 and 548, in the SH2 domain resulting in a STAT6 mutant that is constitutively active in an IL-4 independent manner and is unresponsive to IL-4 stimulation68. STAT6VT differs from the variants we describe as the latter are in the DNA binding domain (Fig 1) and do respond to IL-4 stimulation (Fig 2-4). Nevertheless, the humans we identified with STAT6 GOF variants and STAT6VT mice share a number of key features of the allergic diathesis, including elevated serum IgE and chronic atopic dermatitis.

There is now a growing list of human single gene defects that cause the classic hyper-IgE phenotypic triad of eczema, recurrent skin and lung infections, and elevated serum IgE26,71. Autosomal dominant hyper-IgE syndrome caused by dominant negative variants in STAT3 (aka. Job’s syndrome or STAT3 LOF) is the best characterized of these conditions, but this list of disorders does include both autosomal dominant (STAT3, CARD11, ERBIN) and autosomal recessive (DOCK8, PGM3, IL6ST) conditions72,73. Notably, the patients we identified with STAT6 GOF variants did have some of the extra-immunological features typical of STAT3 LOF (i.e. hyperextensible joints, fractures, vascular anomalies73).

Based on the findings reported in this study, we suggest that heterozygous GOF variants in STAT6 be added to the list of autosomal dominant causes of the hyper-IgE phenotype. While each of the conditions known to cause a hyper-IgE phenotype has some specific clinical features (e.g. viral skin infections are a defining feature of DOCK8 deficiency74,75), there is considerable clinical overlap and clinically-approved testing of these pathways is rarely available. Therefore, we recommend genetic testing as the most efficient initial diagnostic approach to patients who experience severe allergic disease beginning very early in life. Finally, in common with other primary atopic disorders31, we demonstrate that JAK inhibition may be an effective treatment for patients with GOF STAT6 variants.

While this study has many strengths, notably the extreme allergic phenotype of the two patients combined with in-depth functional assessment of their de novo STAT6 variants, a limitation is our lack of validation in primary patient cells due to the death of P1 and challenges with phlebotomy in P2 due to the severity of her atopic dermatitis. Despite this limitation, our study does identify GOF variants in STAT6 as a novel monogenic allergic disorder. We anticipate that this discovery will facilitate the recognition of more affected individuals and, with time, a full definition of the human phenotype caused by germline human STAT6 GOF variants will emerge.

Data Availability

All data produced in the present study are available upon reasonable request to the authors.

Contribution

Me.S., G.B.R, H.M.B, Mi.S, M.L.M, H.A, K.L.D., and S.E.T enrolled patients and analyzed clinical data; J.J.L, A.M, W.W.W., C.D.M.v.K conducted genetic screening and identified variants of significance. M.S., O.F, R.v.d.L, and P.A.R performed, structural modeling, bioinformatic and statistical analysis; M.S., H.Y.L., M.V.S, K.L.D, S.L., and J.D. performed the laboratory experiments; M.S., H.Y.L., M.V.S, and S.E.T. wrote and edited the manuscript. All authors reviewed the final manuscript.

Conflict of Interest Disclosures

Authors report no conflicts of interest in relation to this manuscript.

SUPPLEMENTAL FIGURES

Supplementary Figure 1:
Supplementary Figure 1: Extended data on STAT6 variant-transduced Jurkat T cells.

A) Significantly upregulated (i) and downregulated (ii) genes upon IL-4 treatment in WT (green), p.E382Q (blue) and p.D419G (purple) in Jurkat cells as shown via Venn diagram.

SUPPLEMENTAL TABLES

View this table:
Supplementary Table 1:

Variant annotation and pathogenicity prediction of the variants reported in the two patients.

View this table:
Supplementary Table 2:

IL-4/STAT6 target genes in Jurkat. IL-4/STAT6 target genes are shown that were either upregulated (green) or downregulated (red) in WT-STAT6 transduced Jurkat after stimulation with IL-4 (100ng/mL) for 4 hours.

View this table:
Supplementary Table 3:

Statistics for enrichment analysis done on Jurkat T cells, comparing WT-STAT6 with either p.E382Q or p.D419G at baseline or after stimulation.

Acknowledgements

This work was supported by grants from the Canadian Institutes of Health Research (PJT 178054) (S.E.T.), Genome British Columbia (SIP007) (S.E.T.), and BC Children’s Hospital Foundation. S.E.T. holds a Tier 1 Canada Research Chair in Pediatric Precision Health and the Aubrey J. Tingle Professor of Pediatric Immunology. M.S. is supported by a CIHR Frederick Banting and Charles Best Canada Graduate Scholarships Doctoral Award (CGS-D) and University of British Columbia Four Year Doctoral Fellowship (4YF). H.Y.L. is supported by a CGS-D, 4YF, Killam Doctoral Scholarship, Friedman Award for Scholars in Health, and a BC Children’s Hospital Research Institute Graduate Studentship. M.V.S. is funded by the Vanier Canada Graduate Scholarship and the University of British Columbia 4-Year Doctoral Fellowship (4YF). We would also like to acknowledge the Biomedical Research Centre Sequencing Core (BRC-Seq) for their assistance with RNA-Sequencing and processing.

REFERENCES

  1. 1.↵
    O’Shea JJ, Plenge R. JAK and STAT signaling molecules in immunoregulation and immune-mediated disease. Immunity 2012;36:542–50.
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.↵
    Luo Y, Alexander M, Gadina M, O’Shea JJ, Meylan F, Schwartz DM. JAK-STAT signaling in human disease: From genetic syndromes to clinical inhibition. J Allergy Clin Immunol 2021;148:911–25.
    OpenUrl
  3. 3.↵
    Jatiani SS, Baker SJ, Silverman LR, Reddy EP. Jak/STAT pathways in cytokine signaling and myeloproliferative disorders: approaches for targeted therapies. Genes Cancer 2010;1:979–93.
    OpenUrlCrossRefPubMed
  4. 4.↵
    Villarino AV, Kanno Y, Ferdinand JR, O’Shea JJ. Mechanisms of Jak/STAT signaling in immunity and disease. J Immunol 2015;194:21–7.
    OpenUrlAbstract/FREE Full Text
  5. 5.
    Ihle JN. Cytokine receptor signalling. Nature 1995;377:591–4.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.
    Mogensen TH. IRF and STAT Transcription Factors - From Basic Biology to Roles in Infection, Protective Immunity, and Primary Immunodeficiencies. Front Immunol 2018;9:3047.
    OpenUrlCrossRefPubMed
  7. 7.↵
    Briscoe J, Guschin D, Rogers NC, et al. JAKs, STATs and signal transduction in response to the interferons and other cytokines. Philos Trans R Soc Lond B Biol Sci 1996;351:167–71.
    OpenUrlCrossRefPubMed
  8. 8.↵
    Goenka S, Kaplan MH. Transcriptional regulation by STAT6. Immunologic research 2011;50:87–96.
    OpenUrlCrossRefPubMed
  9. 9.↵
    Takeda K, Tanaka T, Shi W, et al. Essential role of Stat6 in IL-4 signalling. Nature 1996;380:627–30.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.
    Villarino AV, Kanno Y, O’Shea JJ. Mechanisms and consequences of Jak-STAT signaling in the immune system. Nat Immunol 2017;18:374–84.
    OpenUrlCrossRefPubMed
  11. 11.↵
    Villarino AV, Gadina M, O’Shea JJ, Kanno Y. SnapShot: Jak-STAT Signaling II. Cell 2020;181:1696–.e1.
    OpenUrlCrossRef
  12. 12.↵
    Mjosberg J, Bernink J, Golebski K, et al. The transcription factor GATA3 is essential for the function of human type 2 innate lymphoid cells. Immunity 2012;37:649–59.
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    Junttila IS. Tuning the Cytokine Responses: An Update on Interleukin (IL)-4 and IL-13 Receptor Complexes. Front Immunol 2018;9:888.
    OpenUrlCrossRef
  14. 14.↵
    Spolski R, Gromer D, Leonard WJ. The gamma c family of cytokines: fine-tuning signals from IL-2 and IL-21 in the regulation of the immune response. F1000Res 2017;6:1872.
    OpenUrl
  15. 15.↵
    Kaplan MH, Schindler U, Smiley ST, Grusby MJ. Stat6 is required for mediating responses to IL-4 and for development of Th2 cells. Immunity 1996;4:313–9.
    OpenUrlCrossRefPubMedWeb of Science
  16. 16.↵
    Saatian B, Rezaee F, Desando S, et al. Interleukin-4 and interleukin-13 cause barrier dysfunction in human airway epithelial cells. Tissue Barriers 2013;1:e24333.
    OpenUrl
  17. 17.↵
    Mannon P, Reinisch W. Interleukin 13 and its role in gut defence and inflammation. Gut 2012;61:1765–73.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    Oetjen LK, Mack MR, Feng J, et al. Sensory Neurons Co-opt Classical Immune Signaling Pathways to Mediate Chronic Itch. Cell 2017;171:217–28 e13.
    OpenUrlCrossRef
  19. 19.↵
    Bacharier LB, Maspero JF, Katelaris CH, et al. Dupilumab in Children with Uncontrolled Moderate-to-Severe Asthma. N Engl J Med 2021;385:2230–40.
    OpenUrl
  20. 20.↵
    Beck LA, Thaci D, Hamilton JD, et al. Dupilumab treatment in adults with moderate-to-severe atopic dermatitis. N Engl J Med 2014;371:130–9.
    OpenUrlCrossRefPubMed
  21. 21.↵
    Tomkinson A, Kanehiro A, Rabinovitch N, Joetham A, Cieslewicz G, Gelfand EW. The failure of STAT6-deficient mice to develop airway eosinophilia and airway hyperresponsiveness is overcome by interleukin-5. Am J Respir Crit Care Med 1999;160:1283–91.
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.
    Akimoto T, Numata F, Tamura M, et al. Abrogation of bronchial eosinophilic inflammation and airway hyperreactivity in signal transducers and activators of transcription (STAT)6-deficient mice. J Exp Med 1998;187:1537–42.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    Kuperman D, Schofield B, Wills-Karp M, Grusby MJ. Signal transducer and activator of transcription factor 6 (Stat6)-deficient mice are protected from antigen-induced airway hyperresponsiveness and mucus production. J Exp Med 1998;187:939–48.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    Lyons JJ, Milner JD. Primary atopic disorders. J Exp Med 2018;215:1009–22.
    OpenUrlAbstract/FREE Full Text
  25. 25.
    Milner JD. Primary Atopic Disorders. Annu Rev Immunol 2020;38:785–808.
    OpenUrl
  26. 26.↵
    Vaseghi-Shanjani M, Smith KL, Sara RJ, et al. Inborn errors of immunity manifesting as atopic disorders. J Allergy Clin Immunol 2021;148:1130–9.
    OpenUrl
  27. 27.↵
    Asano T, Khourieh J, Zhang P, et al. Human STAT3 variants underlie autosomal dominant hyper-IgE syndrome by negative dominance. J Exp Med 2021;218.
  28. 28.↵
    Fabre A, Marchal S, Barlogis V, et al. Clinical Aspects of STAT3 Gain-of-Function Germline Mutations: A Systematic Review. J Allergy Clin Immunol Pract 2019;7:1958–69 e9.
    OpenUrl
  29. 29.↵
    Ma CA, Xi L, Cauff B, et al. Somatic STAT5b gain-of-function mutations in early onset nonclonal eosinophilia, urticaria, dermatitis, and diarrhea. Blood 2017;129:650–3.
    OpenUrlFREE Full Text
  30. 30.↵
    Kanai T, Jenks J, Nadeau KC. The STAT5b Pathway Defect and Autoimmunity. Front Immunol 2012;3:234.
    OpenUrlCrossRefPubMed
  31. 31.↵
    Del Bel KL, Ragotte RJ, Saferali A, et al. JAK1 gain-of-function causes an autosomal dominant immune dysregulatory and hypereosinophilic syndrome. J Allergy Clin Immunol 2017;139:2016–20 e5.
    OpenUrlCrossRefPubMed
  32. 32.↵
    Gruber CN, Calis JJA, Buta S, et al. Complex Autoinflammatory Syndrome Unveils Fundamental Principles of JAK1 Kinase Transcriptional and Biochemical Function. Immunity 2020;53:672–84 e11.
    OpenUrl
  33. 33.↵
    Uzel G, Sampaio EP, Lawrence MG, et al. Dominant gain-of-function STAT1 mutations in FOXP3 wild-type immune dysregulation-polyendocrinopathy-enteropathy-X-linked-like syndrome. J Allergy Clin Immunol 2013;131:1611–23.
    OpenUrlCrossRefPubMed
  34. 34.↵
    Sharma M, Fu MP, Lu HY, et al. Human germline biallelic complete NFAT1 deficiency causes the triad of progressive joint contractures, osteochondromas, and susceptibility to B cell malignancy. medRxiv 2022:2022.01.30.22269378.
  35. 35.↵
    Fung SY, Lu HY, Sharma M, et al. MALT1-Dependent Cleavage of HOIL1 Modulates Canonical NF-kappaB Signaling and Inflammatory Responsiveness. Front Immunol 2021;12:749794.
    OpenUrl
  36. 36.↵
    Kutner RH, Zhang XY, Reiser J. Production, concentration and titration of pseudotyped HIV-1-based lentiviral vectors. Nat Protoc 2009;4:495–505.
    OpenUrlCrossRefPubMedWeb of Science
  37. 37.↵
    Lu HY, Sharma M, Sharma AA, et al. Mechanistic understanding of the combined immunodeficiency in complete human CARD11 deficiency. J Allergy Clin Immunol 2021.
  38. 38.↵
    Li J, Rodriguez JP, Niu F, et al. Structural basis for DNA recognition by STAT6. Proc Natl Acad Sci U S A 2016;113:13015–20.
    OpenUrlAbstract/FREE Full Text
  39. 39.↵
    Chen Y, Lun AT, Smyth GK. From reads to genes to pathways: differential expression analysis of RNA-Seq experiments using Rsubread and the edgeR quasi-likelihood pipeline. F1000Res 2016;5:1438.
    OpenUrlPubMed
  40. 40.↵
    Ritchie ME, Phipson B, Wu D, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res 2015;43:e47.
    OpenUrlCrossRefPubMed
  41. 41.↵
    Kulpa DA, Talla A, Brehm JH, et al. Differentiation into an Effector Memory Phenotype Potentiates HIV-1 Latency Reversal in CD4(+) T Cells. J Virol 2019;93.
  42. 42.↵
    Elo LL, Jarvenpaa H, Tuomela S, et al. Genome-wide profiling of interleukin-4 and STAT6 transcription factor regulation of human Th2 cell programming. Immunity 2010;32:852–62.
    OpenUrlCrossRefPubMedWeb of Science
  43. 43.↵
    Waterhouse A, Bertoni M, Bienert S, et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res 2018;46:W296–W303.
    OpenUrlCrossRefPubMed
  44. 44.↵
    Pettersen EF, Goddard TD, Huang CC, et al. UCSF Chimera--a visualization system for exploratory research and analysis. J Comput Chem 2004;25:1605–12.
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    Yildiz M, Li H, Bernard D, et al. Activating STAT6 mutations in follicular lymphoma. Blood 2015;125:668–79.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    Zamo A, Pischimarov J, Horn H, Ott G, Rosenwald A, Leich E. The exomic landscape of t(14;18)-negative diffuse follicular lymphoma with 1p36 deletion. Br J Haematol 2018;180:391–4.
    OpenUrl
  47. 47.↵
    Rentzsch P, Schubach M, Shendure J, Kircher M. CADD-Splice-improving genome-wide variant effect prediction using deep learning-derived splice scores. Genome Med 2021;13:31.
    OpenUrlCrossRefPubMed
  48. 48.↵
    Sim NL, Kumar P, Hu J, Henikoff S, Schneider G, Ng PC. SIFT web server: predicting effects of amino acid substitutions on proteins. Nucleic Acids Res 2012;40:W452–7.
    OpenUrlCrossRefPubMedWeb of Science
  49. 49.↵
    Adzhubei IA, Schmidt S, Peshkin L, et al. A method and server for predicting damaging missense mutations. Nat Methods 2010;7:248–9.
    OpenUrlCrossRefPubMedWeb of Science
  50. 50.↵
    Mikita T, Campbell D, Wu P, Williamson K, Schindler U. Requirements for interleukin-4-induced gene expression and functional characterization of Stat6. Mol Cell Biol 1996;16:5811–20.
    OpenUrlAbstract/FREE Full Text
  51. 51.↵
    Kim Y, Kwon EK, Jeon JH, et al. Activation of the STAT6 transcription factor in Jurkat T-cells by the herpesvirus saimiri Tip protein. J Gen Virol 2012;93:330–40.
    OpenUrlCrossRefPubMed
  52. 52.
    Goenka S, Kaplan MH. Transcriptional regulation by STAT6. Immunol Res 2011;50:87–96.
    OpenUrlCrossRefPubMed
  53. 53.
    Bettigole SE, Lis R, Adoro S, et al. The transcription factor XBP1 is selectively required for eosinophil differentiation. Nat Immunol 2015;16:829–37.
    OpenUrlCrossRefPubMed
  54. 54.
    Frischmeyer-Guerrerio PA, Guerrerio AL, Oswald G, et al. TGFbeta receptor mutations impose a strong predisposition for human allergic disease. Sci Transl Med 2013;5:195ra94.
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    Tate JG, Bamford S, Jubb HC, et al. COSMIC: the Catalogue Of Somatic Mutations In Cancer. Nucleic Acids Res 2019;47:D941–D7.
    OpenUrlCrossRefPubMed
  56. 56.↵
    Ritz O, Guiter C, Castellano F, et al. Recurrent mutations of the STAT6 DNA binding domain in primary mediastinal B-cell lymphoma. Blood 2009;114:1236–42.
    OpenUrlAbstract/FREE Full Text
  57. 57.↵
    Bissonnette R, Papp KA, Poulin Y, et al. Topical tofacitinib for atopic dermatitis: a phase IIa randomized trial. Br J Dermatol 2016;175:902–11.
    OpenUrlPubMed
  58. 58.↵
    Kim BS, Howell MD, Sun K, et al. Treatment of atopic dermatitis with ruxolitinib cream (JAK1/JAK2 inhibitor) or triamcinolone cream. J Allergy Clin Immunol 2020;145:572–82.
    OpenUrl
  59. 59.↵
    Ginhoux F, Schultze JL, Murray PJ, Ochando J, Biswas SK. New insights into the multidimensional concept of macrophage ontogeny, activation and function. Nat Immunol 2016;17:34–40.
    OpenUrlCrossRefPubMed
  60. 60.↵
    Sloka S, Silva C, Wang J, Yong VW. Predominance of Th2 polarization by vitamin D through a STAT6-dependent mechanism. J Neuroinflammation 2011;8:56.
    OpenUrlCrossRefPubMed
  61. 61.↵
    Galli SJ, Tsai M. IgE and mast cells in allergic disease. Nat Med 2012;18:693–704.
    OpenUrlCrossRefPubMed
  62. 62.↵
    Matsukura S, Stellato C, Georas SN, et al. Interleukin-13 upregulates eotaxin expression in airway epithelial cells by a STAT6-dependent mechanism. Am J Respir Cell Mol Biol 2001;24:755–61.
    OpenUrlCrossRefPubMedWeb of Science
  63. 63.↵
    Medoff BD, Seung E, Hong S, et al. CD11b+ myeloid cells are the key mediators of Th2 cell homing into the airway in allergic inflammation. J Immunol 2009;182:623–35.
    OpenUrlAbstract/FREE Full Text
  64. 64.↵
    Chandesris MO, Azarine A, Ong KT, et al. Frequent and widespread vascular abnormalities in human signal transducer and activator of transcription 3 deficiency. Circ Cardiovasc Genet 2012;5:25–34.
    OpenUrlAbstract/FREE Full Text
  65. 65.
    Dadak M, Jacobs R, Skuljec J, et al. Gain-of-function STAT1 mutations are associated with intracranial aneurysms. Clin Immunol 2017;178:79–85.
    OpenUrl
  66. 66.↵
    Toubiana J, Okada S, Hiller J, et al. Heterozygous STAT1 gain-of-function mutations underlie an unexpectedly broad clinical phenotype. Blood 2016;127:3154–64.
    OpenUrlAbstract/FREE Full Text
  67. 67.↵
    Liao M, Xu J, Clair AJ, Ehrman B, Graham LM, Eagleton MJ. Local and systemic alterations in signal transducers and activators of transcription (STAT) associated with human abdominal aortic aneurysms. J Surg Res 2012;176:321–8.
    OpenUrlCrossRefPubMed
  68. 68.↵
    Daniel C, Salvekar A, Schindler U. A gain-of-function mutation in STAT6. J Biol Chem 2000;275:14255–9.
    OpenUrlAbstract/FREE Full Text
  69. 69.↵
    DaSilva-Arnold SC, Thyagarajan A, Seymour LJ, et al. Phenotyping acute and chronic atopic dermatitis-like lesions in Stat6VT mice identifies a role for IL-33 in disease pathogenesis. Arch Dermatol Res 2018;310:197–207.
    OpenUrl
  70. 70.↵
    Bruns HA, Schindler U, Kaplan MH. Expression of a constitutively active Stat6 in vivo alters lymphocyte homeostasis with distinct effects in T and B cells. J Immunol 2003;170:3478–87.
    OpenUrlAbstract/FREE Full Text
  71. 71.↵
    Freeman AF, Milner JD. The Child with Elevated IgE and Infection Susceptibility. Curr Allergy Asthma Rep 2020;20:65.
    OpenUrl
  72. 72.↵
    Minegishi Y. Hyper-IgE syndrome, 2021 update. Allergol Int 2021;70:407–14.
    OpenUrl
  73. 73.↵
    Bergerson JRE, Freeman AF. An Update on Syndromes with a Hyper-IgE Phenotype. Immunol Allergy Clin North Am 2019;39:49–61.
    OpenUrlPubMed
  74. 74.↵
    Biggs CM, Keles S, Chatila TA. DOCK8 deficiency: Insights into pathophysiology, clinical features and management. Clin Immunol 2017;181:75–82.
    OpenUrlCrossRef
  75. 75.↵
    Chu EY, Freeman AF, Jing H, et al. Cutaneous manifestations of DOCK8 deficiency syndrome. Arch Dermatol 2012;148:79–84.
    OpenUrlCrossRefPubMed
  76. 76.
    Sakaue S, Kanai M, Tanigawa Y, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nature Genetics 2021;53:1415–24.
    OpenUrl
  77. 77.
    Zhu Z, Guo Y, Shi H, et al. Shared genetic and experimental links between obesity-related traits and asthma subtypes in UK Biobank. Journal of Allergy and Clinical Immunology 2020;145:537–49.
    OpenUrl
  78. 78.
    Pividori M, Schoettler N, Nicolae DL, Ober C, Im HK. Shared and distinct genetic risk factors for childhood-onset and adult-onset asthma: genome-wide and transcriptome-wide studies. The Lancet Respiratory Medicine 2019;7:509–22.
    OpenUrl
  79. 79.
    Ferreira MAR, Mathur R, Vonk JM, et al. Genetic Architectures of Childhood-and Adult-Onset Asthma Are Partly Distinct. The American Journal of Human Genetics 2019;104:665–84.
    OpenUrlCrossRef
  80. 80.
    Zhu Z, Lee PH, Chaffin MD, et al. A genome-wide cross-trait analysis from UK Biobank highlights the shared genetic architecture of asthma and allergic diseases. Nature Genetics 2018;50:857–64.
    OpenUrl
  81. 81.
    Johansson Å, Rask-Andersen M, Karlsson T, Ek WE. Genome-wide association analysis of 350 000 Caucasians from the UK Biobank identifies novel loci for asthma, hay fever and eczema. Human Molecular Genetics 2019;28:4022–41.
    OpenUrlCrossRefPubMed
  82. 82.
    Zhu Z, Zhu X, Liu C-L, et al. Shared genetics of asthma and mental health disorders: a large-scale genome-wide cross-trait analysis. European Respiratory Journal 2019;54:1901507.
    OpenUrlAbstract/FREE Full Text
  83. 83.
    Demenais F, Margaritte-Jeannin P, Barnes KC, et al. Multiancestry association study identifies new asthma risk loci that colocalize with immune-cell enhancer marks. Nature Genetics 2018;50:42–53.
    OpenUrl
  84. 84.
    Pickrell JK, Berisa T, Liu JZ, Ségurel L, Tung JY, Hinds DA. Detection and interpretation of shared genetic influences on 42 human traits. Nature Genetics 2016;48:709–17.
    OpenUrlCrossRefPubMed
  85. 85.
    Shrine N, Portelli MA, John C, et al. Moderate-to-severe asthma in individuals of European ancestry: a genome-wide association study. The Lancet Respiratory Medicine 2019;7:20–34.
    OpenUrl
  86. 86.
    Valette K, Li Z, Bon-Baret V, et al. Prioritization of candidate causal genes for asthma in susceptibility loci derived from UK Biobank. Communications Biology 2021;4.
  87. 87.
    Olafsdottir TA, Theodors F, Bjarnadottir K, et al. Eighty-eight variants highlight the role of T cell regulation and airway remodeling in asthma pathogenesis. Nature Communications 2020;11.
  88. 88.
    Han Y, Jia Q, Jahani PS, et al. Genome-wide analysis highlights contribution of immune system pathways to the genetic architecture of asthma. Nature Communications 2020;11.
  89. 89.
    Astle WJ, Elding H, Jiang T, et al. The Allelic Landscape of Human Blood Cell Trait Variation and Links to Common Complex Disease. Cell 2016;167:1415–29.e19.
    OpenUrlCrossRefPubMed
  90. 90.
    Vuckovic D, Bao EL, Akbari P, et al. The Polygenic and Monogenic Basis of Blood Traits and Diseases. Cell 2020;182:1214–31.e11.
    OpenUrlCrossRefPubMed
  91. 91.
    Kachuri L, Jeon S, Dewan AT, et al. Genetic determinants of blood-cell traits influence susceptibility to childhood acute lymphoblastic leukemia. The American Journal of Human Genetics 2021;108:1823–35.
    OpenUrl
  92. 92.
    Chen M-H, Raffield LM, Mousas A, et al. Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations. Cell 2020;182:1198–213.e14.
    OpenUrlCrossRefPubMed
  93. 93.
    Kichaev G, Bhatia G, Loh P-R, et al. Leveraging Polygenic Functional Enrichment to Improve GWAS Power. The American Journal of Human Genetics 2019;104:65–75.
    OpenUrlCrossRefPubMed
  94. 94.
    Ferreira MAR, Vonk JM, Baurecht H, et al. Age-of-onset information helps identify 76 genetic variants associated with allergic disease. PLOS Genetics 2020;16:e1008725.
    OpenUrl
  95. 95.
    Ferreira MA, Vonk JM, Baurecht H, et al. Shared genetic origin of asthma, hay fever and eczema elucidates allergic disease biology. Nature Genetics 2017;49:1752–7.
    OpenUrlCrossRefPubMed
  96. 96.
    Tanaka N, Koido M, Suzuki A, et al. Eight novel susceptibility loci and putative causal variants in atopic dermatitis. Journal of Allergy and Clinical Immunology 2021;148:1293–306.
    OpenUrl
  97. 97.
    Granada M, Wilk JB, Tuzova M, et al. A genome-wide association study of plasma total IgE concentrations in the Framingham Heart Study. Journal of Allergy and Clinical Immunology 2012;129:840–5.e21.
    OpenUrlCrossRefPubMedWeb of Science
  98. 98.
    Daya M, Cox C, Acevedo N, et al. Multiethnic genome-wide and HLA association study of total serum IgE level. Journal of Allergy and Clinical Immunology 2021;148:1589–95.
    OpenUrl
  99. 99.
    Waage J, Standl M, Curtin JA, et al. Genome-wide association and HLA fine-mapping studies identify risk loci and genetic pathways underlying allergic rhinitis. Nature Genetics 2018;50:1072–80.
    OpenUrl
  100. 100.
    Bønnelykke K, Matheson MC, Pers TH, et al. Meta-analysis of genome-wide association studies identifies ten loci influencing allergic sensitization. Nature Genetics 2013;45:902–6.
    OpenUrlCrossRefPubMed
  101. 101.
    Sleiman PMA, Wang M-L, Cianferoni A, et al. GWAS identifies four novel eosinophilic esophagitis loci. Nature Communications 2014;5:5593.
    OpenUrl
Back to top
Posted April 28, 2022.
Download PDF
Data/Code
Email
Human Germline Heterozygous Gain-of-Function STAT6 Variants Cause Severe Allergic Disease
Mehul Sharma, Henry Y. Lu, Maryam Vaseghi-Shanjani, Kate L. Del Bel, Oriol Fornes, Robin van der Lee, Phillip A. Richmond, Susan Lin, Joshua Dalmann, Jessica J. Lee, Allison Matthews, Géraldine Blanchard-Rohner, Clara D M van Karnebeek, H. Melanie Bedford, Wyeth W. Wasserman, Michael Seear, Margaret L. McKinnon, Hanan Ahmed, Stuart E. Turvey
medRxiv 2022.04.25.22274265; doi: https://doi.org/10.1101/2022.04.25.22274265
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area