Jason D. Hughes

Finding DNA regulatory motifs within unaligned noncoding sequences clustered by whole-genome mRNA quantitation

Whole-genome mRNA quantitation can be used to identify the genes that are most responsive to environmental or genotypic change. By searching for mutually similar DNA elements among the upstream non-coding DNA sequences of these genes, we can identify candidate regulatory motifs and corresponding candidate sets of coregulated genes. We have tested this strategy by applying it to three extensively studied regulatory systems in the yeast Saccharomyces cerevisiae: galactose response, heat shock, and mating type. Galactose-response data yielded the known binding site of Gal4, and six of nine genes known to be induced by galactose. Heat shock data yielded the cell-cycle activation motif, which is known to mediate cell-cycle dependent activation, and a set of genes coding for all four nucleosomal proteins. Mating type α and a data yielded all of the four relevant DNA motifs and most of the known a- and α-specific genes.

A computational analysis of whole-genome expression data reveals chromosomal domains of gene expression

Chromosome correlation maps display correlations between the expression patterns of genes on the same chromosome. Using these maps, we show here that adjacent pairs of genes, as well as nearby non-adjacent pairs of genes, show correlated expression independent of their orientation. We present specific examples of adjacent pairs with highly correlated expression patterns, in which the promoter of only one of the two genes contains an upstream activating sequence (UAS) known to be associated with that expression pattern. Finally, we show that genes with similar functions tend to occur in adjacent positions along the chromosomes. Our results suggest that, in certain chromosomal expression domains, an UAS can affect the transcription of genes that are not immediately downstream from it.

Early-onset lymphoproliferation and autoimmunity caused by germline STAT3 gain-of-function mutations.

Germline loss-of-function mutations in the transcription factor signal transducer and activator of transcription 3 (STAT3) cause immunodeficiency, whereas somatic gain-of-function mutations in STAT3 are associated with large granular lymphocytic leukemic, myelodysplastic syndrome, and aplastic anemia. Recently, germline mutations in STAT3 have also been associated with autoimmune disease. Here, we report on 13 individuals from 10 families with lymphoproliferation and early-onset solid-organ autoimmunity associated with 9 different germline heterozygous mutations in STAT3. Patients exhibited a variety of clinical features, with most having lymphadenopathy, autoimmune cytopenias, multiorgan autoimmunity (lung, gastrointestinal, hepatic, and/or endocrine dysfunction), infections, and short stature. Functional analyses demonstrate that these mutations confer a gain-of-function in STAT3 leading to secondary defects in STAT5 and STAT1 phosphorylation and the regulatory T-cell compartment. Treatment targeting a cytokine pathway that signals through STAT3 led to clinical improvement in 1 patient, suggesting a potential therapeutic option for such patients. These results suggest that there is a broad range of autoimmunity caused by germline STAT3 gain-of-function mutations, and that hematologic autoimmunity is a major component of this newly described disorder. Some patients for this study were enrolled in a trial registered at www.clinicaltrials.gov as #NCT00001350.

Patients with LRBA deficiency show CTLA4 loss and immune dysregulation responsive to abatacept therapy

Trafficking from bedside to bench Typically in translational research, a discovery in cell or molecular biology is later exploited to improve patient care. Occasionally, information flows in the opposite direction. Lo et al. found that patients with an autoimmune disorder caused by deficiency of a protein called LRBA responded dramatically to the drug abatacept (see the Perspective by Sansom). Abatacept contains a segment of a potent inhibitory immune receptor, CTLA4. Experiments prompted by this observation revealed the relationship between the two proteins: LRBA controls the intracellular trafficking and degradation of CTLA4. This information may further improve patient care, because other clinically approved drugs have the desired mechanism of action with potentially fewer side effects. Science, this issue p. 436; see also p. 377 A rare autoimmune disorder is caused by aberrant degradation of a potent inhibitory immune receptor. [Also see Perspective by Sansom] Mutations in the LRBA gene (encoding the lipopolysaccharide-responsive and beige-like anchor protein) cause a syndrome of autoimmunity, lymphoproliferation, and humoral immune deficiency. The biological role of LRBA in immunologic disease is unknown. We found that patients with LRBA deficiency manifested a dramatic and sustained improvement in response to abatacept, a CTLA4 (cytotoxic T lymphocyte antigen-4)–immunoglobulin fusion drug. Clinical responses and homology of LRBA to proteins controlling intracellular trafficking led us to hypothesize that it regulates CTLA4, a potent inhibitory immune receptor. We found that LRBA colocalized with CTLA4 in endosomal vesicles and that LRBA deficiency or knockdown increased CTLA4 turnover, which resulted in reduced levels of CTLA4 protein in FoxP3+ regulatory and activated conventional T cells. In LRBA-deficient cells, inhibition of lysosome degradation with chloroquine prevented CTLA4 loss. These findings elucidate a mechanism for CTLA4 trafficking and control of immune responses and suggest therapies for diseases involving the CTLA4 pathway.

Heterozygous splice mutation in PIK3R1 causes human immunodeficiency with lymphoproliferation due to dominant activation of PI3K

Lucas et al. identify humans with a gain-of-function mutation in PIK3R1, encoding the p85α subunit of PI3K. The splice site mutation causes in-frame skipping of exon 11, resulting in altered p85α association with p110δ that stabilizes the catalytic subunit but fails to properly inhibit catalytic activity. The patients have immunodeficiency and lymphoproliferation with skewing of CD8+ T cells toward terminally differentiated and senescent effector cells that have shortened telomeres.

Elevated basal serum tryptase identifies a multisystem disorder associated with increased TPSAB1 copy number

Elevated basal serum tryptase levels are present in 4–6% of the general population, but the cause and relevance of such increases are unknown. Previously, we described subjects with dominantly inherited elevated basal serum tryptase levels associated with multisystem complaints including cutaneous flushing and pruritus, dysautonomia, functional gastrointestinal symptoms, chronic pain, and connective tissue abnormalities, including joint hypermobility. Here we report the identification of germline duplications and triplications in the TPSAB1 gene encoding α-tryptase that segregate with inherited increases in basal serum tryptase levels in 35 families presenting with associated multisystem complaints. Individuals harboring alleles encoding three copies of α-tryptase had higher basal serum levels of tryptase and were more symptomatic than those with alleles encoding two copies, suggesting a gene-dose effect. Further, we found in two additional cohorts (172 individuals) that elevated basal serum tryptase levels were exclusively associated with duplication of α-tryptase–encoding sequence in TPSAB1, and affected individuals reported symptom complexes seen in our initial familial cohort. Thus, our findings link duplications in TPSAB1 with irritable bowel syndrome, cutaneous complaints, connective tissue abnormalities, and dysautonomia.

Recurrent rhinovirus infections in a child with inherited MDA5 deficiency

MDA5 is a cytosolic sensor of double-stranded RNA (ds)RNA including viral byproducts and intermediates. We studied a child with life-threatening, recurrent respiratory tract infections, caused by viruses including human rhinovirus (HRV), influenza virus, and respiratory syncytial virus (RSV). We identified in her a homozygous missense mutation in IFIH1 that encodes MDA5. Mutant MDA5 was expressed but did not recognize the synthetic MDA5 agonist/(ds)RNA mimic polyinosinic-polycytidylic acid. When overexpressed, mutant MDA5 failed to drive luciferase activity from the IFNB1 promoter or promoters containing ISRE or NF-&kgr;B sequence motifs. In respiratory epithelial cells or fibroblasts, wild-type but not knockdown of MDA5 restricted HRV infection while increasing IFN-stimulated gene expression and IFN-&bgr;/&lgr;. However, wild-type MDA5 did not restrict influenza virus or RSV replication. Moreover, nasal epithelial cells from the patient, or fibroblasts gene-edited to express mutant MDA5, showed increased replication of HRV but not influenza or RSV. Thus, human MDA5 deficiency is a novel inborn error of innate and/or intrinsic immunity that causes impaired (ds)RNA sensing, reduced IFN induction, and susceptibility to the common cold virus.

BACH2 immunodeficiency illustrates an association between super-enhancers and haploinsufficiency.

The transcriptional programs that guide lymphocyte differentiation depend on the precise expression and timing of transcription factors (TFs). The TF BACH2 is essential for T and B lymphocytes and is associated with an archetypal super-enhancer (SE). Single-nucleotide variants in the BACH2 locus are associated with several autoimmune diseases, but BACH2 mutations that cause Mendelian monogenic primary immunodeficiency have not previously been identified. Here we describe a syndrome of BACH2-related immunodeficiency and autoimmunity (BRIDA) that results from BACH2 haploinsufficiency. Affected subjects had lymphocyte-maturation defects that caused immunoglobulin deficiency and intestinal inflammation. The mutations disrupted protein stability by interfering with homodimerization or by causing aggregation. We observed analogous lymphocyte defects in Bach2-heterozygous mice. More generally, we observed that genes that cause monogenic haploinsufficient diseases were substantially enriched for TFs and SE architecture. These findings reveal a previously unrecognized feature of SE architecture in Mendelian diseases of immunity: heterozygous mutations in SE-regulated genes identified by whole-exome/genome sequencing may have greater significance than previously recognized.

ERBIN deficiency links STAT3 and TGF-β pathway defects with atopy in humans

Nonimmunological connective tissue phenotypes in humans are common among some congenital and acquired allergic diseases. Several of these congenital disorders have been associated with either increased TGF-&bgr; activity or impaired STAT3 activation, suggesting that these pathways might intersect and that their disruption may contribute to atopy. In this study, we show that STAT3 negatively regulates TGF-&bgr; signaling via ERBB2-interacting protein (ERBIN), a SMAD anchor for receptor activation and SMAD2/3 binding protein. Individuals with dominant-negative STAT3 mutations (STAT3mut) or a loss-of-function mutation in ERBB2IP (ERBB2IPmut) have evidence of deregulated TGF-&bgr; signaling with increased regulatory T cells and total FOXP3 expression. These naturally occurring mutations, recapitulated in vitro, impair STAT3–ERBIN–SMAD2/3 complex formation and fail to constrain nuclear pSMAD2/3 in response to TGF-&bgr;. In turn, cell-intrinsic deregulation of TGF-&bgr; signaling is associated with increased functional IL-4R&agr; expression on naive lymphocytes and can induce expression and activation of the IL-4/IL-4R&agr;/GATA3 axis in vitro. These findings link increased TGF-&bgr; pathway activation in ERBB2IPmut and STAT3mut patient lymphocytes with increased T helper type 2 cytokine expression and elevated IgE.

Estimation of kinship coefficient in structured and admixed populations using sparse sequencing data

Knowledge of biological relatedness between samples is important for many genetic studies. In large-scale human genetic association studies, the estimated kinship is used to remove cryptic relatedness, control for family structure, and estimate trait heritability. However, estimation of kinship is challenging for sparse sequencing data, such as those from off-target regions in target sequencing studies, where genotypes are largely uncertain or missing. Existing methods often assume accurate genotypes at a large number of markers across the genome. We show that these methods, without accounting for the genotype uncertainty in sparse sequencing data, can yield a strong downward bias in kinship estimation. We develop a computationally efficient method called SEEKIN to estimate kinship for both homogeneous samples and heterogeneous samples with population structure and admixture. Our method models genotype uncertainty and leverages linkage disequilibrium through imputation. We test SEEKIN on a whole exome sequencing dataset (WES) of Singapore Chinese and Malays, which involves substantial population structure and admixture. We show that SEEKIN can accurately estimate kinship coefficient and classify genetic relatedness using off-target sequencing data down sampled to ~0.15X depth. In application to the full WES dataset without down sampling, SEEKIN also outperforms existing methods by properly analyzing shallow off-target data (~0.75X). Using both simulated and real phenotypes, we further illustrate how our method improves estimation of trait heritability for WES studies.