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.

Investigating the potential of single-cell DNA methylation data to detect allele-specific methylation and imprinting.

Allele-specific methylation (ASM) is an epigenetic modification whereby one parental allele becomes methylated and the other unmethylated at a specific locus. ASM is most often driven by the presence of nearby heterozygous variants that influence methylation, but also occurs somatically in the context of genomic imprinting. In this study, we investigate ASM using publicly available single-cell reduced representation bisulfite sequencing (scRRBS) data on 608 B cells sampled from six healthy B cell samples and 1,230 cells from 11 chronic lymphocytic leukemia (CLL) samples. We developed a likelihood-based criterion to test whether a CpG exhibited ASM, based on the distributions of methylated and unmethylated reads both within and across cells. Applying our likelihood ratio test, 65,998 CpG sites exhibited ASM in healthy B cell samples according to a Bonferroni criterion (p < 8.4 × 10-9), and 32,862 CpG sites exhibited ASM in CLL samples (p < 8.5 × 10-9). We also called ASM at the sample level. To evaluate the accuracy of our method, we called heterozygous variants from the scRRBS data, which enabled variant-based calls of ASM within each cell. Comparing sample-level ASM calls to the variant-based measures of ASM, we observed a positive predictive value of 76%-100% across samples. We observed high concordance of ASM across samples and an overrepresentation of ASM in previously reported imprinted genes and genes with imprinting binding motifs. Our study demonstrates that single-cell bisulfite sequencing is a potentially powerful tool to investigate ASM, especially as studies expand to increase the number of samples and cells sequenced.

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 control of the human brain proteome.

We generated an online brain pQTL resource for 7,376 proteins through the analysis of genetic and proteomic data derived from post-mortem samples of the dorsolateral prefrontal cortex of 330 older adults. The identified pQTLs tend to be non-synonymous variation, are over-represented among variants associated with brain diseases, and replicate well (77%) in an independent brain dataset. Comparison to a large study of brain eQTLs revealed that about 75% of pQTLs are also eQTLs. In contrast, about 40% of eQTLs were identified as pQTLs. These results are consistent with lower pQTL mapping power and greater evolutionary constraint on protein abundance. The latter is additionally supported by observations of pQTLs with large effects' tending to be rare, deleterious, and associated with proteins that have evidence for fewer protein-protein interactions. Mediation analyses using matched transcriptomic and proteomic data provided additional evidence that pQTL effects are often, but not always, mediated by mRNA. Specifically, we identified roughly 1.6 times more mRNA-mediated pQTLs than mRNA-independent pQTLs (550 versus 341). Our pQTL resource provides insight into the functional consequences of genetic variation in the human brain and a basis for novel investigations of genetics and disease.

Whole-genome sequencing of African Americans implicates differential genetic architecture in inflammatory bowel disease.

Whether or not populations diverge with respect to the genetic contribution to risk of specific complex diseases is relevant to understanding the evolution of susceptibility and origins of health disparities. Here, we describe a large-scale whole-genome sequencing study of inflammatory bowel disease encompassing 1,774 affected individuals and 1,644 healthy control Americans with African ancestry (African Americans). Although no new loci for inflammatory bowel disease are discovered at genome-wide significance levels, we identify numerous instances of differential effect sizes in combination with divergent allele frequencies. For example, the major effect at PTGER4 fine maps to a single credible interval of 22 SNPs corresponding to one of four independent associations at the locus in European ancestry individuals but with an elevated odds ratio for Crohn disease in African Americans. A rare variant aggregate analysis implicates Ca2+-binding neuro-immunomodulator CALB2 in ulcerative colitis. Highly significant overall overlap of common variant risk for inflammatory bowel disease susceptibility between individuals with African and European ancestries was observed, with 41 of 241 previously known lead variants replicated and overall correlations in effect sizes of 0.68 for combined inflammatory bowel disease. Nevertheless, subtle differences influence the performance of polygenic risk scores, and we show that ancestry-appropriate weights significantly improve polygenic prediction in the highest percentiles of risk. The median amount of variance explained per locus remains the same in African and European cohorts, providing evidence for compensation of effect sizes as allele frequencies diverge, as expected under a highly polygenic model of disease.

POIROT: a powerful test for parent-of-origin effects in unrelated samples leveraging multiple phenotypes.

MOTIVATION: There is widespread interest in identifying genetic variants that exhibit parent-of-origin effects (POEs) wherein the effect of an allele on phenotype expression depends on its parental origin. POEs can arise from different phenomena including genomic imprinting and have been documented for many complex traits. Traditional tests for POEs require family data to determine parental origins of transmitted alleles. As most genome-wide association studies (GWAS) sample unrelated individuals (where allelic parental origin is unknown), the study of POEs in such datasets requires sophisticated statistical methods that exploit genetic patterns we anticipate observing when POEs exist. We propose a method to improve discovery of POE variants in large-scale GWAS samples that leverages potential pleiotropy among multiple correlated traits often collected in such studies. Our method compares the phenotypic covariance matrix of heterozygotes to homozygotes based on a Robust Omnibus Test. We refer to our method as the Parent of Origin Inference using Robust Omnibus Test (POIROT) of multiple quantitative traits. RESULTS: Through simulation studies, we compared POIROT to a competing univariate variance-based method which considers separate analysis of each phenotype. We observed POIROT to be well-calibrated with improved power to detect POEs compared to univariate methods. POIROT is robust to non-normality of phenotypes and can adjust for population stratification and other confounders. Finally, we applied POIROT to GWAS data from the UK Biobank using BMI and two cholesterol phenotypes. We identified 338 genome-wide significant loci for follow-up investigation. AVAILABILITY AND IMPLEMENTATION: The code for this method is available at https://github.com/staylorhead/POIROT-POE.

Rare variant modifier analysis identifies variants in SEC24D associated with orofacial cleft subtypes.

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p = 6.86 [Formula: see text] 10-7). Of the 15 variants within SEC24D, 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p = 0.0009), and created binding sites for 23 TFs, including Pax6 (p = 6.12 [Formula: see text] 10-5) and Pax9 (p = 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study indicates that rare genetic variation may contribute to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC.

Genome-wide Enrichment of De Novo Coding Mutations in Orofacial Cleft Trios.

Although de novo mutations (DNMs) are known to increase an individual's risk of congenital defects, DNMs have not been fully explored regarding orofacial clefts (OFCs), one of the most common human birth defects. Therefore, whole-genome sequencing of 756 child-parent trios of European, Colombian, and Taiwanese ancestry was performed to determine the contributions of coding DNMs to an individual's OFC risk. Overall, we identified a significant excess of loss-of-function DNMs in genes highly expressed in craniofacial tissues, as well as genes associated with known autosomal dominant OFC syndromes. This analysis also revealed roles for zinc-finger homeobox domain and SOX2-interacting genes in OFC etiology.

Rare variants found in clinical gene panels illuminate the genetic and allelic architecture of orofacial clefting.

PURPOSE: Orofacial clefts (OFCs) are common birth defects including cleft lip, cleft lip and palate, and cleft palate. OFCs have heterogeneous etiologies, complicating clinical diagnostics because it is not always apparent if the cause is Mendelian, environmental, or multifactorial. Sequencing is not currently performed for isolated or sporadic OFCs; therefore, we estimated the diagnostic yield for 418 genes in 841 cases and 294 controls. METHODS: We evaluated 418 genes using genome sequencing and curated variants to assess their pathogenicity using American College of Medical Genetics criteria. RESULTS: 9.04% of cases and 1.02% of controls had "likely pathogenic" variants (P < .0001), which was almost exclusively driven by heterozygous variants in autosomal genes. Cleft palate (17.6%) and cleft lip and palate (9.09%) cases had the highest yield, whereas cleft lip cases had a 2.80% yield. Out of 39 genes with likely pathogenic variants, 9 genes, including CTNND1 and IRF6, accounted for more than half of the yield (4.64% of cases). Most variants (61.8%) were "variants of uncertain significance", occurring more frequently in cases (P = .004), but no individual gene showed a significant excess of variants of uncertain significance. CONCLUSION: These results underscore the etiological heterogeneity of OFCs and suggest sequencing could reduce the diagnostic gap in OFCs.

Racial and ethnic diversity of classic and clinical variant galactosemia in the United States.

Classic and clinical variant galactosemia (CG/CVG) are allelic, autosomal recessive disorders that result from deficiency of galactose-1-P uridylyltransferase (GALT). CG/CVG has been reported globally among patients of diverse ancestries, but most large studies of outcomes have included, almost exclusively, patients categorized as White or Caucasian. As a first step to explore whether the cohorts studied are representative of the CG/CVG population at large, we sought to define the racial and ethnic makeup of CG/CVG newborns in a diverse population with essentially universal newborn screening (NBS) for galactosemia: the United States (US). First, we estimated the predicted racial and ethnic distribution of CG/CVG by combining the reported demographics of US newborns from 2016 to 2018 with predicted homozygosity or compound heterozygosity of pathogenic, or likely pathogenic, GALT alleles from the relevant ancestral groups. Incorporating some simplifying assumptions, we predicted that of US newborns diagnosed with CG/CVG, 65% should be White (non-Hispanic), 23% should be Black (non-Hispanic), 10% should be Hispanic, and 2% should be Asian (non-Hispanic). Next, we calculated the observed racial and ethnic distribution of US newborns diagnosed with CG/CVG using available de-identified data from state NBS programs from 2016 to 2018. Of the 235 newborns in this cohort, 41 were categorized as other or unknown. Of the remaining 194, 66% were White (non-Hispanic or ethnicity unknown), 16% were Black (non-Hispanic or ethnicity unknown),15% were Hispanic, and 2% were Asian (non-Hispanic or ethnicity unknown). This observed distribution was statistically indistinguishable from the predicted distribution. To the limits of our study, these data confirm the racial and ethnic diversity of newborns with CG/CVG in the US, demonstrate an approach for estimating CG/CVG racial and ethnic diversity in other populations, and raise the troubling possibility that current understanding of long-term outcomes in CG/CVG may be skewed by ascertainment bias of the cohorts studied.

LDL cholesterol is associated with higher AD neuropathology burden independent of APOE.

OBJECTIVE: APOE is a strong risk factor for Alzheimer's disease (AD) and associated with higher low-density lipoprotein cholesterol (LDL-C) levels. Moreover, LDL-C is associated with the development of clinically ascertained AD; however, whether this association is present with the underlying neuropathological manifestations of AD or whether it is independent of the effect of APOE is unknown and is the focus of this paper. METHODS: Individuals in the Religious Orders Study/Memory and Ageing Project cohorts with longitudinal measures of blood lipids and detailed autopsies were studied. We modelled the relationship between blood lipids and 12 age-related brain pathologies using a linear mixed model adjusted for potential confounding factors and stratified by APOE genotype with overall significance determined by meta-analysis. Blood lipids considered were LDL-C, high-density lipoprotein cholesterol and triglycerides. Brain pathologies included AD pathology measured by silver staining (Braak stage, a modified Consortium to Establish a Registry for Alzheimer's Disease [CERAD] score and global AD pathology) and immunohistochemistry (beta-amyloid and neurofibrillary tangles) as well as cerebral microinfarct, cerebral macroinfarct, cerebral amyloid angiopathy, cerebral atherosclerosis, hippocampal sclerosis, TDP-43 cytoplasmic inclusions and Lewy bodies. RESULTS: 559 participants (69.1% female) had complete data for analysis. They were followed for a median of 7 years and a median of 3 years prior to dementia onset. LDL-C was associated with all measures of AD neuropathology (neurofibrillary tangles, beta-amyloid, Braak stage, modified CERAD score and global AD pathology) and cerebral amyloid angiopathy independent of APOE after adjusting for age, sex, cholesterol-lowering medication use, body mass index, smoking and education at false discovery rate (FDR) p-value <0.05. CONCLUSIONS: These findings implicate LDL-C in the pathophysiology of AD independent of APOE and suggest LDL-C is a modifiable risk factor for AD.

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.

Blood-Derived DNA Methylation Signatures of Crohn's Disease and Severity of Intestinal Inflammation.

BACKGROUND & AIMS: Crohn's disease is a relapsing and remitting inflammatory disorder with a variable clinical course. Although most patients present with an inflammatory phenotype (B1), approximately 20% of patients rapidly progress to complicated disease, which includes stricturing (B2), within 5 years. We analyzed DNA methylation patterns in blood samples of pediatric patients with Crohn's disease at diagnosis and later time points to identify changes that associate with and might contribute to disease development and progression. METHODS: We obtained blood samples from 164 pediatric patients (1-17 years old) with Crohn's disease (B1 or B2) who participated in a North American study and were followed for 5 years. Participants without intestinal inflammation or symptoms served as controls (n = 74). DNA methylation patterns were analyzed in samples collected at time of diagnosis and 1-3 years later at approximately 850,000 sites. We used genetic association and the concept of Mendelian randomization to identify changes in DNA methylation patterns that might contribute to the development of or result from Crohn's disease. RESULTS: We identified 1189 5'-cytosine-phosphate-guanosine-3' (CpG) sites that were differentially methylated between patients with Crohn's disease (at diagnosis) and controls. Methylation changes at these sites correlated with plasma levels of C-reactive protein. A comparison of methylation profiles of DNA collected at diagnosis of Crohn's disease vs during the follow-up period showed that, during treatment, alterations identified in methylation profiles at the time of diagnosis of Crohn's disease more closely resembled patterns observed in controls, irrespective of disease progression to B2. We identified methylation changes at 3 CpG sites that might contribute to the development of Crohn's disease. Most CpG methylation changes associated with Crohn's disease disappeared with treatment of inflammation and might be a result of Crohn's disease. CONCLUSIONS: Methylation patterns observed in blood samples from patients with Crohn's disease accompany acute inflammation; with treatment, these change to resemble methylation patterns observed in patients without intestinal inflammation. These findings indicate that Crohn's disease-associated patterns of DNA methylation observed in blood samples are a result of the inflammatory features of the disease and are less likely to contribute to disease development or progression.

Identification of Polycystic Kidney Disease 1 Like 1 Gene Variants in Children With Biliary Atresia Splenic Malformation Syndrome.

Biliary atresia (BA) is the most common cause of end-stage liver disease in children and the primary indication for pediatric liver transplantation, yet underlying etiologies remain unknown. Approximately 10% of infants affected by BA exhibit various laterality defects (heterotaxy) including splenic abnormalities and complex cardiac malformations-a distinctive subgroup commonly referred to as the biliary atresia splenic malformation (BASM) syndrome. We hypothesized that genetic factors linking laterality features with the etiopathogenesis of BA in BASM patients could be identified through whole-exome sequencing (WES) of an affected cohort. DNA specimens from 67 BASM subjects, including 58 patient-parent trios, from the National Institute of Diabetes and Digestive and Kidney Diseases-supported Childhood Liver Disease Research Network (ChiLDReN) underwent WES. Candidate gene variants derived from a prespecified set of 2,016 genes associated with ciliary dysgenesis and/or dysfunction or cholestasis were prioritized according to pathogenicity, population frequency, and mode of inheritance. Five BASM subjects harbored rare and potentially deleterious biallelic variants in polycystic kidney disease 1 like 1 (PKD1L1), a gene associated with ciliary calcium signaling and embryonic laterality determination in fish, mice, and humans. Heterozygous PKD1L1 variants were found in 3 additional subjects. Immunohistochemical analysis of liver from the one BASM subject available revealed decreased PKD1L1 expression in bile duct epithelium when compared to normal livers and livers affected by other noncholestatic diseases. Conclusion: WES identified biallelic and heterozygous PKD1L1 variants of interest in 8 BASM subjects from the ChiLDReN data set; the dual roles for PKD1L1 in laterality determination and ciliary function suggest that PKD1L1 is a biologically plausible, cholangiocyte-expressed candidate gene for the BASM syndrome.

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.

A Statistical Approach for Testing Cross-Phenotype Effects of Rare Variants.

Increasing empirical evidence suggests that many genetic variants influence multiple distinct phenotypes. When cross-phenotype effects exist, multivariate association methods that consider pleiotropy are often more powerful than univariate methods that model each phenotype separately. Although several statistical approaches exist for testing cross-phenotype effects for common variants, there is a lack of similar tests for gene-based analysis of rare variants. In order to fill this important gap, we introduce a statistical method for cross-phenotype analysis of rare variants using a nonparametric distance-covariance approach that compares similarity in multivariate phenotypes to similarity in rare-variant genotypes across a gene. The approach can accommodate both binary and continuous phenotypes and further can adjust for covariates. Our approach yields a closed-form test whose significance can be evaluated analytically, thereby improving computational efficiency and permitting application on a genome-wide scale. We use simulated data to demonstrate that our method, which we refer to as the Gene Association with Multiple Traits (GAMuT) test, provides increased power over competing approaches. We also illustrate our approach using exome-chip data from the Genetic Epidemiology Network of Arteriopathy.

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.