Large-Scale Exome Sequencing Study Implicates Both Developmental and Functional Changes in the Neurobiology of Autism.

We present the largest exome sequencing study of autism spectrum disorder (ASD) to date (n = 35,584 total samples, 11,986 with ASD). Using an enhanced analytical framework to integrate de novo and case-control rare variation, we identify 102 risk genes at a false discovery rate of 0.1 or less. Of these genes, 49 show higher frequencies of disruptive de novo variants in individuals ascertained to have severe neurodevelopmental delay, whereas 53 show higher frequencies in individuals ascertained to have ASD; comparing ASD cases with mutations in these groups reveals phenotypic differences. Expressed early in brain development, most risk genes have roles in regulation of gene expression or neuronal communication (i.e., mutations effect neurodevelopmental and neurophysiological changes), and 13 fall within loci recurrently hit by copy number variants. In cells from the human cortex, expression of risk genes is enriched in excitatory and inhibitory neuronal lineages, consistent with multiple paths to an excitatory-inhibitory imbalance underlying ASD.

Loss of function of OTUD7A in the schizophrenia- associated 15q13.3 deletion impairs synapse development and function in human neurons.

Identifying causative gene(s) within disease-associated large genomic regions of copy-number variants (CNVs) is challenging. Here, by targeted sequencing of genes within schizophrenia (SZ)-associated CNVs in 1,779 SZ cases and 1,418 controls, we identified three rare putative loss-of-function (LoF) mutations in OTU deubiquitinase 7A (OTUD7A) within the 15q13.3 deletion in cases but none in controls. To tie OTUD7A LoF with any SZ-relevant cellular phenotypes, we modeled the OTUD7A LoF mutation, rs757148409, in human induced pluripotent stem cell (hiPSC)-derived induced excitatory neurons (iNs) by CRISPR-Cas9 engineering. The mutant iNs showed a ∼50% decrease in OTUD7A expression without undergoing nonsense-mediated mRNA decay. The mutant iNs also exhibited marked reduction of dendritic complexity, density of synaptic proteins GluA1 and PSD-95, and neuronal network activity. Congruent with the neuronal phenotypes in mutant iNs, our transcriptomic analysis showed that the set of OTUD7A LoF-downregulated genes was enriched for those relating to synapse development and function and was associated with SZ and other neuropsychiatric disorders. These results suggest that OTUD7A LoF impairs synapse development and neuronal function in human neurons, providing mechanistic insight into the possible role of OTUD7A in driving neuropsychiatric phenotypes associated with the 15q13.3 deletion.

Identification of PSMB5 as a genetic modifier of fragile X-associated tremor/ataxia syndrome.

Fragile X­associated tremor/ataxia syndrome (FXTAS) is a debilitating late-onset neurodegenerative disease in premutation carriers of the expanded CGG repeat in FMR1 that presents with a spectrum of neurological manifestations, such as gait ataxia, intention tremor, and parkinsonism [P. J. Hagerman, R. J. Hagerman, Ann. N. Y. Acad. Sci. 1338, 58­70 (2015); S. Jacquemont et al., JAMA 291, 460­469 (2004)]. Here, we performed whole-genome sequencing (WGS) on male premutation carriers (CGG55­200) and prioritized candidate variants to screen for candidate genetic modifiers using a Drosophila model of FXTAS. We found 18 genes that genetically modulate CGG-associated neurotoxicity in Drosophila, such as Prosbeta5 (PSMB5), pAbp (PABPC1L), e(y)1 (TAF9), and CG14231 (OSGEPL1). Among them, knockdown of Prosbeta5 (PSMB5) suppressed CGG-associated neurodegeneration in the fly as well as in N2A cells. Interestingly, an expression quantitative trait locus variant in PSMB5, PSMB5rs11543947-A, was found to be associated with decreased expression of PSMB5 and delayed onset of FXTAS in human FMR1 premutation carriers. Finally, we demonstrate evidence that PSMB5 knockdown results in suppression of CGG neurotoxicity via both the RAN translation and RNA-mediated toxicity mechanisms, thereby presenting a therapeutic strategy for FXTAS.

Genetic contributors to risk of schizophrenia in the presence of a 22q11.2 deletion.

Schizophrenia occurs in about one in four individuals with 22q11.2 deletion syndrome (22q11.2DS). The aim of this International Brain and Behavior 22q11.2DS Consortium (IBBC) study was to identify genetic factors that contribute to schizophrenia, in addition to the ~20-fold increased risk conveyed by the 22q11.2 deletion. Using whole-genome sequencing data from 519 unrelated individuals with 22q11.2DS, we conducted genome-wide comparisons of common and rare variants between those with schizophrenia and those with no psychotic disorder at age ≥25 years. Available microarray data enabled direct comparison of polygenic risk for schizophrenia between 22q11.2DS and independent population samples with no 22q11.2 deletion, with and without schizophrenia (total n = 35,182). Polygenic risk for schizophrenia within 22q11.2DS was significantly greater for those with schizophrenia (padj = 6.73 × 10-6). Novel reciprocal case-control comparisons between the 22q11.2DS and population-based cohorts showed that polygenic risk score was significantly greater in individuals with psychotic illness, regardless of the presence of the 22q11.2 deletion. Within the 22q11.2DS cohort, results of gene-set analyses showed some support for rare variants affecting synaptic genes. No common or rare variants within the 22q11.2 deletion region were significantly associated with schizophrenia. These findings suggest that in addition to the deletion conferring a greatly increased risk to schizophrenia, the risk is higher when the 22q11.2 deletion and common polygenic risk factors that contribute to schizophrenia in the general population are both present.

Inherited genetic susceptibility to acute lymphoblastic leukemia in Down syndrome.

Children with Down syndrome (DS) have a 20-fold increased risk of acute lymphoblastic leukemia (ALL) and distinct somatic features, including CRLF2 rearrangement in ∼50% of cases; however, the role of inherited genetic variation in DS-ALL susceptibility is unknown. We report the first genome-wide association study of DS-ALL, comprising a meta-analysis of 4 independent studies, with 542 DS-ALL cases and 1192 DS controls. We identified 4 susceptibility loci at genome-wide significance: rs58923657 near IKZF1 (odds ratio [OR], 2.02; Pmeta = 5.32 × 10-15), rs3731249 in CDKN2A (OR, 3.63; Pmeta = 3.91 × 10-10), rs7090445 in ARID5B (OR, 1.60; Pmeta = 8.44 × 10-9), and rs3781093 in GATA3 (OR, 1.73; Pmeta = 2.89 × 10-8). We performed DS-ALL vs non-DS ALL case-case analyses, comparing risk allele frequencies at these and other established susceptibility loci (BMI1, PIP4K2A, and CEBPE) and found significant association with DS status for CDKN2A (OR, 1.58; Pmeta = 4.1 × 10-4). This association was maintained in separate regression models, both adjusting for and stratifying on CRLF2 overexpression and other molecular subgroups, indicating an increased penetrance of CDKN2A risk alleles in children with DS. Finally, we investigated functional significance of the IKZF1 risk locus, and demonstrated mapping to a B-cell super-enhancer, and risk allele association with decreased enhancer activity and differential protein binding. IKZF1 knockdown resulted in significantly higher proliferation in DS than non-DS lymphoblastoid cell lines. Our findings demonstrate a higher penetrance of the CDKN2A risk locus in DS and serve as a basis for further biological insights into DS-ALL etiology.

Bayesian Pathway Analysis for Complex Interactions.

Modern epidemiologic studies permit investigation of the complex pathways that mediate effects of social, behavioral, and molecular factors on health outcomes. Conventional analytical approaches struggle with high-dimensional data, leading to high likelihoods of both false-positive and false-negative inferences. Herein, we describe a novel Bayesian pathway analysis approach, the algorithm for learning pathway structure (ALPS), which addresses key limitations in existing approaches to complex data analysis. ALPS uses prior information about pathways in concert with empirical data to identify and quantify complex interactions within networks of factors that mediate an association between an exposure and an outcome. We illustrate ALPS through application to a complex gene-drug interaction analysis in the Predictors of Breast Cancer Recurrence (ProBe CaRe) Study, a Danish cohort study of premenopausal breast cancer patients (2002-2011), for which conventional analyses severely limit the quality of inference.

New phenotypes associated with 3q29 duplication syndrome: Results from the 3q29 registry.

3q29 duplication syndrome (3q29dup) is a rare genomic disorder caused by a 1.6 Mb duplication (GRCh38 chr3:195,998,000-197,623,000). Case reports indicate the 3q29dup is likely to be pathogenic, but the full range of manifestations is not well understood. We used the 3q29 registry (https://3q29.com) to ascertain 31 individuals with 3q29dup, the largest cohort ever surveyed in a systematic way. For comparison, we ascertained 117 individuals with the reciprocal 3q29 deletion and 64 typically developing controls. We used a custom medical and demographic questionnaire to assess physical and developmental phenotypes, and two standardized instruments, the Social Responsiveness Scale and Child Behavior Checklist/Adult Behavior Checklist, to assess social disability. Participants with 3q29dup report a high rate of problems in the first year of life (80.6%), including feeding problems (55%), failure to gain weight (42%), hypotonia (39%), and respiratory distress (29%). In early childhood, learning problems (71.0%) and seizures (25.8%) are common. Additionally, the rate of self-reported autism spectrum disorder diagnoses (39%) is substantially elevated compared to the general population, suggesting that the 3q29 duplication may be an autism susceptibility locus. This is the most comprehensive description of 3q29dup to date. Our findings can be used to develop evidence-based strategies for early intervention and management of 3q29dup.

PEMapper and PECaller provide a simplified approach to whole-genome sequencing.

The analysis of human whole-genome sequencing data presents significant computational challenges. The sheer size of datasets places an enormous burden on computational, disk array, and network resources. Here, we present an integrated computational package, PEMapper/PECaller, that was designed specifically to minimize the burden on networks and disk arrays, create output files that are minimal in size, and run in a highly computationally efficient way, with the single goal of enabling whole-genome sequencing at scale. In addition to improved computational efficiency, we implement a statistical framework that allows for a base by base error model, allowing this package to perform as well or better than the widely used Genome Analysis Toolkit (GATK) in all key measures of performance on human whole-genome sequences.

Genetic variants and pathways implicated in a pediatric inflammatory bowel disease cohort.

In the United States, approximately 5% of individuals with inflammatory bowel disease (IBD) are younger than 20 years old. Studies of pediatric cohorts can provide unique insights into genetic architecture of IBD, which includes Crohn's disease (CD) and ulcerative colitis (UC). Large genome-wide association studies have found more than 200 IBD-associated loci but explain a minority of disease variance for CD and UC. We sought to characterize the contribution of rare variants to disease development, comparing exome sequencing of 368 pediatric IBD patients to publicly available exome sequencing (dbGaP) and aggregate frequency data (ExAC). Using dbGaP data, we performed logistic regression for common variants and optimal unified association tests (SKAT-O) for rare, likely-deleterious variants. We further compared rare variants to ExAC counts with Fisher's exact tests. We did pathway enrichment analysis on the most significant genes from each comparison. Many variants overlapped with known IBD-associated genes (e.g. NOD2). Rare variants were enriched in CD-associated loci (p = 0.009) and showed suggestive enrichment in neutrophil function genes (p = 0.05). Pathway enrichment implicated immune-related pathways, especially cell killing and apoptosis. Variants in extracellular matrix genes also emerged as an important theme in our analysis.

Clinical and Genomic Correlates of Neutrophil Reactive Oxygen Species Production in Pediatric Patients With Crohn's Disease.

BACKGROUND & AIMS: Individuals with monogenic disorders of phagocyte function develop chronic colitis that resembles Crohn's disease (CD). We tested for associations between mutations in genes encoding reduced nicotinamide adenine dinucleotide phosphate (NADPH) oxidases, neutrophil function, and phenotypes of CD in pediatric patients. METHODS: We performed whole-exome sequence analysis to identify mutations in genes encoding NADPH oxidases (such as CYBA, CYBB, NCF1, NCF2, NCF4, RAC1, and RAC2) using DNA from 543 pediatric patients with inflammatory bowel diseases. Blood samples were collected from an additional 129 pediatric patients with CD and 26 children without IBD (controls); we performed assays for neutrophil activation, reactive oxygen species (ROS) production, and bacteria uptake and killing. Whole-exome sequence analysis was performed using DNA from 46 of the children with CD to examine associations with NADPH gene mutations; RNA sequence analyses were performed using blood cells from 46 children with CD to test for variations in neutrophil gene expression associated with ROS production. RESULTS: We identified 26 missense mutations in CYBA, CYBB, NCF1, NCF2, and NCF4. Patients with CD who carried mutations in these genes were 3-fold more likely to have perianal disease (P = .0008) and stricturing complications (P = .002) than children with CD without these mutations. Among patients with CD with none of these mutations, 9% had undergone abdominal surgery; among patients with mutations in these NADPH oxidase genes, 31% had undergone abdominal surgery (P = .0004). A higher proportion of neutrophils from children with CD had low ROS production (47%) than from controls (15%) among the 129 patients tested for ROS (P = .002). Minor alleles of the NADPH genes were detected in 7% of children with CD whose neutrophils produced normal levels of ROS vs 38% of children whose neutrophils produced low levels of ROS (P = .009). Neutrophils that produced low levels of ROS had specific alterations in genes that regulate glucose metabolism and antimicrobial responses. CONCLUSIONS: We identified missense mutations in genes that encode NADPH oxidases in children with CD; these were associated with a more aggressive disease course and reduced ROS production by neutrophils from the patients.

Omicseq: a web-based search engine for exploring omics datasets.

The development and application of high-throughput genomics technologies has resulted in massive quantities of diverse omics data that continue to accumulate rapidly. These rich datasets offer unprecedented and exciting opportunities to address long standing questions in biomedical research. However, our ability to explore and query the content of diverse omics data is very limited. Existing dataset search tools rely almost exclusively on the metadata. A text-based query for gene name(s) does not work well on datasets wherein the vast majority of their content is numeric. To overcome this barrier, we have developed Omicseq, a novel web-based platform that facilitates the easy interrogation of omics datasets holistically to improve 'findability' of relevant data. The core component of Omicseq is trackRank, a novel algorithm for ranking omics datasets that fully uses the numerical content of the dataset to determine relevance to the query entity. The Omicseq system is supported by a scalable and elastic, NoSQL database that hosts a large collection of processed omics datasets. In the front end, a simple, web-based interface allows users to enter queries and instantly receive search results as a list of ranked datasets deemed to be the most relevant. Omicseq is freely available at http://www.omicseq.org.

Characterization of genetic loci that affect susceptibility to inflammatory bowel diseases in African Americans.

BACKGROUND & AIMS: Inflammatory bowel disease (IBD) has familial aggregation in African Americans (AAs), but little is known about the molecular genetic susceptibility. Mapping studies using the Immunochip genotyping array expand the number of susceptibility loci for IBD in Caucasians to 163, but the contribution of the 163 loci and European admixture to IBD risk in AAs is unclear. We performed a genetic mapping study using the Immunochip to determine whether IBD susceptibility loci in Caucasians also affect risk in AAs and identify new associated loci. METHODS: We recruited AAs with IBD and without IBD (controls) from 34 IBD centers in the United States; additional controls were collected from 4 other Immunochip studies. Association and admixture loci were mapped for 1088 patients with Crohn's disease, 361 with ulcerative colitis, 62 with IBD type unknown, and 1797 controls; 130,241 autosomal single-nucleotide polymorphisms (SNPs) were analyzed. RESULTS: The strongest associations were observed between ulcerative colitis and HLA rs9271366 (P = 7.5 × 10(-6)), Crohn's disease and 5p13.1 rs4286721 (P = 3.5 × 10(-6)), and IBD and KAT2A rs730086 (P = 2.3 × 10(-6)). Additional suggestive associations (P < 4.2 × 10(-5)) were observed between Crohn's disease and IBD and African-specific SNPs in STAT5A and STAT3; between IBD and SNPs in IL23R, IL12B, and C2orf43; and between ulcerative colitis and SNPs near HDAC11 and near LINC00994. The latter 3 loci have not been previously associated with IBD, but require replication. Established Caucasian associations were replicated in AAs (P < 3.1 × 10(-4)) at NOD2, IL23R, 5p15.3, and IKZF3. Significant admixture (P < 3.9 × 10(-4)) was observed for 17q12-17q21.31 (IZKF3 through STAT3), 10q11.23-10q21.2, 15q22.2-15q23, and 16p12.2-16p12.1. Network analyses showed significant enrichment (false discovery rate <1 × 10(-5)) in genes that encode members of the JAK-STAT, cytokine, and chemokine signaling pathways, as well those involved in pathogenesis of measles. CONCLUSIONS: In a genetic analysis of 3308 AA IBD cases and controls, we found that many variants associated with IBD in Caucasians also showed association evidence with these diseases in AAs; we also found evidence for variants and loci not previously associated with IBD. The complex genetic factors that determine risk for or protection against IBD in different populations require further study.

Novel features of 3q29 deletion syndrome: Results from the 3q29 registry.

3q29 deletion syndrome is caused by a recurrent, typically de novo heterozygous 1.6 Mb deletion, but because incidence of the deletion is rare (1 in 30,000 births) the phenotype is not well described. To characterize the range of phenotypic manifestations associated with 3q29 deletion syndrome, we have developed an online registry (3q29deletion.org) for ascertainment of study subjects and phenotypic data collection via Internet-based survey instruments. We report here on data collected during the first 18 months of registry operation, from 44 patients. This is the largest cohort of 3q29 deletion carriers ever assembled and surveyed in a systematic way. Our data reveal that 28% of registry participants report neuropsychiatric phenotypes, including anxiety disorder, panic attacks, depression, bipolar disorder, and schizophrenia. Other novel findings include a high prevalence (64%) of feeding problems in infancy and reduced weight at birth for 3q29 deletion carriers (average reduction 13.9 oz (394 g), adjusted for gestational age and sex, P = 6.5e-07). We further report on the frequency of heart defects, autism, recurrent ear infections, gastrointestinal phenotypes, and dental phenotypes, among others. We also report on the expected timing of delayed developmental milestones. This is the most comprehensive description of the 3q29 deletion phenotype to date. These results are clinically actionable toward improving patient care for 3q29 deletion carriers, and can guide the expectations of physicians and parents. These data also demonstrate the value of patient-reported outcomes to reveal the full phenotypic spectrum of rare genomic disorders.

Genomic characterization of the Bacillus cereus sensu lato species: backdrop to the evolution of Bacillus anthracis.

The key genes required for Bacillus anthracis to cause anthrax have been acquired recently by horizontal gene transfer. To understand the genetic background for the evolution of B. anthracis virulence, we obtained high-redundancy genome sequences of 45 strains of the Bacillus cereus sensu lato (s.l.) species that were chosen for their genetic diversity within the species based on the existing multilocus sequence typing scheme. From the resulting data, we called more than 324,000 new genes representing more than 12,333 new gene families for this group. The core genome size for the B. cereus s.l. group was ∼1750 genes, with another 2150 genes found in almost every genome constituting the extended core. There was a paucity of genes specific and conserved in any clade. We found no evidence of recent large-scale gene loss in B. anthracis or for unusual accumulation of nonsynonymous DNA substitutions in the chromosome; however, several B. cereus genomes isolated from soil and not previously associated with human disease were degraded to various degrees. Although B. anthracis has undergone an ecological shift within the species, its chromosome does not appear to be exceptional on a macroscopic scale compared with close relatives.

Contribution of copy-number variation to Down syndrome-associated atrioventricular septal defects.

PURPOSE: The goal of this study was to identify the contribution of large copy-number variants to Down syndrome-associated atrioventricular septal defects, the risk for which in the trisomic population is 2,000-fold more as compared with that of the general disomic population. METHODS: Genome-wide copy-number variant analysis was performed on 452 individuals with Down syndrome (210 cases with complete atrioventricular septal defects; 242 controls with structurally normal hearts) using Affymetrix SNP 6.0 arrays, making this the largest heart study conducted to date on a trisomic background. RESULTS: Large, common copy-number variants with substantial effect sizes (OR > 2.0) do not account for the increased risk observed in Down syndrome-associated atrioventricular septal defects. By contrast, cases had a greater burden of large, rare deletions (P < 0.01) and intersected more genes (P < 0.007) as compared with controls. We also observed a suggestive enrichment of deletions intersecting ciliome genes in cases as compared with controls. CONCLUSION: Our data provide strong evidence that large, rare deletions increase the risk of Down syndrome-associated atrioventricular septal defects, whereas large, common copy-number variants do not appear to increase the risk of Down syndrome-associated atrioventricular septal defects. The genetic architecture of atrioventricular septal defects is complex and multifactorial in nature.

Variant ATRX syndrome with dysfunction of ATRX and MAGT1 genes.

A 0.8 kb intronic duplication in MAGT1 and a single base pair deletion in the last exon of ATRX were identified using a chromosome X-specific microarray and exome sequencing in a family with five males demonstrating intellectual disability (ID) and unusual skin findings (e.g., generalized pruritus). MAGT1 is an Mg²âº transporter previously associated with primary immunodeficiency and ID, whereas mutations in ATRX cause ATRX-ID syndrome. In patient cells, the function of ATRX was demonstrated to be abnormal based on altered RNA/protein expression, hypomethylation of rDNA, and abnormal cytokinesis. Dysfunction of MAGT1 was reflected in reduced RNA/protein expression and Mg²âº influx. The mutation in ATRX most likely explains the ID, whereas MAGT1 disruption could be linked to abnormal skin findings, as normal magnesium homeostasis is necessary for skin health. This work supports observations that multiple mutations collectively contribute to the phenotypic variability of syndromic ID, and emphasizes the importance of correlating clinical phenotype with genomic and cell function analyses.

Excess variants in AFF2 detected by massively parallel sequencing of males with autism spectrum disorder.

Autism spectrum disorder (ASD) is a heterogeneous disorder with substantial heritability, most of which is unexplained. ASD has a population prevalence of one percent and affects four times as many males as females. Patients with fragile X E (FRAXE) intellectual disability, which is caused by a silencing of the X-linked gene AFF2, display a number of ASD-like phenotypes. Duplications and deletions at the AFF2 locus have also been reported in cases with moderate intellectual disability and ASD. We hypothesized that other rare X-linked sequence variants at the AFF2 locus might contribute to ASD. We sequenced the AFF2 genomic region in 202 male ASD probands and found that 2.5% of males sequenced had missense mutations at highly conserved evolutionary sites. When compared with the frequency of missense mutations in 5545 X chromosomes from unaffected controls, we saw a statistically significant enrichment in patients with ASD (OR: 4.9; P < 0.014). In addition, we identified rare AFF2 3' UTR variants at conserved sites which alter gene expression in a luciferase assay. These data suggest that rare variation in AFF2 may be a previously unrecognized ASD susceptibility locus and may help explain some of the male excess of ASD.

Susceptibility to childhood-onset rheumatoid arthritis: investigation of a weighted genetic risk score that integrates cumulative effects of variants at five genetic loci.

OBJECTIVE: Children with childhood-onset rheumatoid arthritis (RA) include those with rheumatoid factor or anti-citrullinated protein antibody-positive juvenile idiopathic arthritis. To test the hypothesis that adult-onset RA-associated variants are also associated with childhood-onset RA, we investigated RA-associated variants at 5 loci in a cohort of patients with childhood-onset RA. We also assessed the cumulative association of these variants in susceptibility to childhood-onset RA using a weighted genetic risk score (wGRS). METHODS: A total of 155 children with childhood-onset RA and 684 healthy controls were genotyped for 5 variants in the PTPN22, TRAF1/C5, STAT4, and TNFAIP3 loci. High-resolution HLA-DRB1 genotypes were available for 149 cases and 373 controls. We tested each locus for association with childhood-onset RA via logistic regression. We also computed a wGRS for each subject, with weights based on the natural log of the published odds ratios (ORs) for the alleles investigated, and used logistic regression to test the wGRS for association with childhood-onset RA. RESULTS: Childhood-onset RA was associated with TNFAIP3 rs10499194 (OR 0.60 [95% confidence interval 0.44-0.83]), PTPN22 rs2476601 (OR 1.61 [95% confidence interval 1.11-2.31]), and STAT4 rs7574865 (OR 1.41 [95% confidence interval 1.06-1.87]) variants. The wGRS was significantly different between cases and controls (P < 2 × 10(-16) ). Individuals in the third to fifth quintiles of wGRS had a significantly increased disease risk compared to baseline (individuals in the first quintile). Higher wGRS was associated with increased risk of childhood-onset RA, especially among males. CONCLUSION: The magnitude and direction of the association between TNFAIP3, STAT4, and PTPN22 variants and childhood-onset RA are similar to those observed in RA, suggesting that adult-onset RA and childhood-onset RA share common genetic risk factors. Using a wGRS, we have demonstrated the cumulative association of RA-associated variants with susceptibility to childhood-onset RA.

Exome sequencing identifies a novel FOXP3 mutation in a 2-generation family with inflammatory bowel disease.

OBJECTIVES: Inflammatory bowel disease (IBD) is heritable, but a total of 163 variants commonly implicated in IBD pathogenesis account for only 25% of the heritability. Rare, highly penetrant genetic variants may also explain mendelian forms of IBD and some of the missing heritability. To test the hypothesis that rare loss-of-function mutations can be causative, we performed whole exome sequencing (WES) on 5 members of a 2-generation family of European ancestry presenting with an early-onset and atypical form of IBD. METHODS: WES was performed for all of the 5 family members; the mother and 3 male offspring were affected, whereas the father was unaffected. Mapping, annotation, and filtering criteria were used to reduce candidate variants. For functional testing we performed forkhead box P3 (FOXP3) staining and a T-cell suppression assay. RESULTS: We identified a novel missense variant in exon 6 of the X-linked FOXP3 gene. The c.694A>C substitution in FOXP3 results in a cysteine-to-glycine change at the protein position 232 that is completely conserved among all vertebrates. This variant (heterozygous in the mother and hemizygous in all 3 affected sons) did not impair FOXP3 protein expression, but significantly reduced the ability of the host's T regulatory cells to suppress an inappropriate autoimmune response. The variant results in a milder immune dysregulation, polyendocrinopathy, enteropathy, and X-linked phenotype with early-onset IBD. CONCLUSIONS: Our study illustrates the successful application of WES for making a definitive molecular diagnosis in a case of multiply affected families, with atypical IBD-like phenotype. Our results also have important implications for disease biology and disease-directed therapeutic development.