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.

Transethnic analysis of the human leukocyte antigen region for ulcerative colitis reveals not only shared but also ethnicity-specific disease associations.

Inflammatory bowel disease (IBD) is a chronic inflammatory disease of the gut. Genetic association studies have identified the highly variable human leukocyte antigen (HLA) region as the strongest susceptibility locus for IBD and specifically DRB1*01:03 as a determining factor for ulcerative colitis (UC). However, for most of the association signal such as delineation could not be made because of tight structures of linkage disequilibrium within the HLA. The aim of this study was therefore to further characterize the HLA signal using a transethnic approach. We performed a comprehensive fine mapping of single HLA alleles in UC in a cohort of 9272 individuals with African American, East Asian, Puerto Rican, Indian and Iranian descent and 40 691 previously analyzed Caucasians, additionally analyzing whole HLA haplotypes. We computationally characterized the binding of associated HLA alleles to human self-peptides and analyzed the physicochemical properties of the HLA proteins and predicted self-peptidomes. Highlighting alleles of the HLA-DRB1*15 group and their correlated HLA-DQ-DR haplotypes, we not only identified consistent associations (regarding effects directions/magnitudes) across different ethnicities but also identified population-specific signals (regarding differences in allele frequencies). We observed that DRB1*01:03 is mostly present in individuals of Western European descent and hardly present in non-Caucasian individuals. We found peptides predicted to bind to risk HLA alleles to be rich in positively charged amino acids. We conclude that the HLA plays an important role for UC susceptibility across different ethnicities. This research further implicates specific features of peptides that are predicted to bind risk and protective HLA proteins.

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.

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.

Early Onset Granulomatous Colitis Associated with a Mutation in NCF4 Resolved with Hematopoietic Stem Cell Transplantation.

A 4-year-old boy presented with perianal abscess and granulomatous colitis, which led the diagnosis of Crohn's disease. He became refractory to all available therapies and required colectomy. Targeted sequencing revealed a deleterious variant in NCF4, causing severe neutrophil dysfunction. He underwent hematopoietic stem cell transplantation (HSCT) with an excellent outcome.

Genome-Wide Association Study Identifies African-Specific Susceptibility Loci in African Americans With Inflammatory Bowel Disease.

BACKGROUND & AIMS: The inflammatory bowel diseases (IBD) ulcerative colitis (UC) and Crohn's disease (CD) cause significant morbidity and are increasing in prevalence among all populations, including African Americans. More than 200 susceptibility loci have been identified in populations of predominantly European ancestry, but few loci have been associated with IBD in other ethnicities. METHODS: We performed 2 high-density, genome-wide scans comprising 2345 cases of African Americans with IBD (1646 with CD, 583 with UC, and 116 inflammatory bowel disease unclassified) and 5002 individuals without IBD (controls, identified from the Health Retirement Study and Kaiser Permanente database). Single-nucleotide polymorphisms (SNPs) associated at P < 5.0 × 10-8 in meta-analysis with a nominal evidence (P < .05) in each scan were considered to have genome-wide significance. RESULTS: We detected SNPs at HLA-DRB1, and African-specific SNPs at ZNF649 and LSAMP, with associations of genome-wide significance for UC. We detected SNPs at USP25 with associations of genome-wide significance for IBD. No associations of genome-wide significance were detected for CD. In addition, 9 genes previously associated with IBD contained SNPs with significant evidence for replication (P < 1.6 × 10-6): ADCY3, CXCR6, HLA-DRB1 to HLA-DQA1 (genome-wide significance on conditioning), IL12B,PTGER4, and TNC for IBD; IL23R, PTGER4, and SNX20 (in strong linkage disequilibrium with NOD2) for CD; and KCNQ2 (near TNFRSF6B) for UC. Several of these genes, such as TNC (near TNFSF15), CXCR6, and genes associated with IBD at the HLA locus, contained SNPs with unique association patterns with African-specific alleles. CONCLUSIONS: We performed a genome-wide association study of African Americans with IBD and identified loci associated with UC in only this population; we also replicated IBD, CD, and UC loci identified in European populations. The detection of variants associated with IBD risk in only people of African descent demonstrates the importance of studying the genetics of IBD and other complex diseases in populations beyond those of European ancestry.

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.

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.

Vitamin D status and bone mineral density in African American children with Crohn disease.

BACKGROUND: Vitamin D deficiency and low bone mineral density (BMD) are complications of inflammatory bowel disease. Vitamin D deficiency is more prevalent among African Americans compared with whites. There are little data comparing differences in serum 25-hydroxyvitamin D (25OHD) concentrations and BMD between African American and white children with Crohn disease (CD). METHODS: We compared serum 25OHD concentrations of African American children with CD (n = 52) to white children with CD (n = 64) and healthy African American controls (n = 40). We also analyzed BMD using dual-energy x-ray absorptiometry results from our pediatric CD population. RESULTS: African American children with CD had lower serum 25OHD concentrations (16.1 [95% confidence interval, CI 14.5-17.9] ng/mL) than whites with CD (22.3 [95% CI 20.2-24.6] ng/mL; P < 0.001). African Americans with CD and controls exhibited similar serum 25OHD concentration (16.1 [95% CI 14.5-17.9] vs 16.3 [95% CI 14.4-18.4] ng/mL; NS). African Americans with CD exhibited no difference in serum 25OHD concentration when controlling for seasonality, disease severity, and surgical history, although serum 25OHD concentration was significantly decreased in overweight children (body mass index ≥85%, P = 0.003). Multiple regression analysis demonstrated that obese African American girls with CD had the lowest serum 25OHD concentrations (9.6 [95% CI 6.8-13.5] ng/mL). BMD was comparable between African American and white children with CD (z score -0.4 ± 0.9 vs -0.7 ± 1.2; NS). CONCLUSIONS: African American children with CD are more likely to have vitamin D deficiency compared with white children with CD, but have similar BMD. CD disease severity and history of surgery do not affect serum 25OHD concentrations among African American children with CD. African American children have low serum 25OHD concentrations, independent of CD, compared with white children. Future research should focus on how race affects vitamin D status and BMD in children with CD.

Methylation quantitative trait loci are largely consistent across disease states in Crohn's disease.

Recently, we identified 1,189 CpG sites whose DNA methylation level in blood associated with Crohn's disease. Here, we examined associations between DNA methylation and genetic variants to identify methylation quantitative trait loci across disease states in (1) 402 blood samples from 164 newly diagnosed pediatric Crohn's disease patients taken at 2 time points (diagnosis and follow-up), and 74 non-inflammatory bowel disease controls, (2) 780 blood samples from a non-Crohn's disease adult population, and (3) 40 ileal biopsies (17 Crohn's disease cases and 23 non-inflammatory bowel disease controls) from group (1). Genome-wide DNAm profiling and genotyping were performed using the Illumina MethylationEPIC and Illumina Multi-Ethnic arrays. SNP-CpG associations were identified via linear models adjusted for age, sex, disease status, disease subtype, estimated cell proportions, and genotype-based principal components. In total, we observed 535,448 SNP-CpG associations between 287,881 SNPs and 12,843 CpG sites (P < 8.21 × 10-14). Associations were highly consistent across different ages, races, disease states, and tissue types, suggesting that the majority of these methylation quantitative trait loci participate in common gene regulation. However, genes near CpGs associated with inflammatory bowel disease SNPs were enriched for 18 KEGG pathways relevant to inflammatory bowel disease-linked immune function and inflammatory responses. We observed suggestive evidence for a small number of tissue-specific associations and disease-specific associations in ileum, though larger studies will be needed to confirm these results. Our study concludes that the vast majority of blood-derived methylation quantitative trait loci are common across individuals, though a subset may be involved in processes related to Crohn's disease. Independent cohort studies will be required to validate these findings.

Eicosatetraynoic Acid and Butyrate Regulate Human Intestinal Organoid Mitochondrial and Extracellular Matrix Pathways Implicated in Crohn's Disease Strictures.

BACKGROUND: Perturbagen analysis of Crohn's disease (CD) ileal gene expression data identified small molecules including eicosatetraynoic acid (ETYA), which may exert an antifibrotic effect. We developed a patient-specific human intestinal organoid (HIO) model system to test small molecule regulation of mitochondrial and wound-healing functions implicated in stricturing behavior. METHODS: HIOs were made from CD induced pluripotent stem cells with and without a loss-of-function haplotype in the DUOX2 gene implicated in ileal homeostasis and characterized under basal conditions and following exposure to butyrate and ETYA using RNA sequencing, flow cytometry, and immunofluorescent and polarized light microscopy. Mitochondrial activity was measured using high-resolution respirometry and tissue stiffness using atomic force microscopy. RESULTS: HIOs expressed core mitochondrial and extracellular matrix (ECM) genes and enriched biologic functions implicated in CD ileal strictures; ECM gene expression was suppressed by both butyrate and ETYA, with butyrate also suppressing genes regulating epithelial proliferation. Consistent with this, butyrate, but not ETYA, exerted a profound effect on HIO epithelial mitochondrial function, reactive oxygen species production, and cellular abundance. Butyrate and ETYA suppressed HIO expression of alpha smooth muscle actin expressed by myofibroblasts, type I collagen, and collagen protein abundance. HIOs exhibited tissue stiffness comparable to normal human ileum; this was reduced by chronic ETYA exposure in HIOs carrying the DUOX2 loss-of-function haplotype. CONCLUSIONS: ETYA regulates ECM genes implicated in strictures and suppresses collagen content and tissue stiffness in an HIO model. HIOs provide a platform to test personalized therapeutics, including small molecules prioritized by perturbagen analysis.

A subset of pediatric Crohn's disease patients develop intestinal strictures requiring surgery. The microbial metabolite butyrate and eicosatetraynoic acid regulate pathways implicated in stricture formation in a human intestinal organoid model system, which may be used to test new therapies.

Predicting adverse drug reactions of two-drug combinations using structural and transcriptomic drug representations to train an artificial neural network.

Adverse drug reactions (ADRs) are pharmacological events triggered by drug interactions with various sources of origin including drug-drug interactions. While there are many computational studies that explore models to predict ADRs originating from single drugs, only a few of them explore models that predict ADRs from drug combinations. Further, as far as we know, none of them have developed models using transcriptomic data, specifically the LINCS L1000 drug-induced gene expression data to predict ADRs for drug combinations. In this study, we use the TWOSIDES database as a source of ADRs originating from two-drug combinations. 34,549 common drug pairs between these two databases were used to train an artificial neural network (ANN), to predict 243 ADRs that were induced by at least 10% of the drug pairs. Our model predicts the occurrence of these ADRs with an average accuracy of 82% across a multifold cross-validation.

Site- and Taxa-Specific Disease-Associated Oral Microbial Structures Distinguish Inflammatory Bowel Diseases.

BACKGROUND: The gut and oral microbiome have independently been shown to be associated with inflammatory bowel disease (IBD). However, it is not known to what extent gut and oral microbial disease markers converge in terms of their composition in IBD. Further, the spatial and temporal variation within the oral microenvironments of IBD remain to be elucidated. PATIENTS AND METHODS: We used a prospectively recruited cohort of patients with IBD (n = 47) and unrelated healthy control patients (n = 18) to examine the spatial and temporal distribution of microbiota within the various oral microenvironments, represented by saliva, tongue, buccal mucosa, and plaque, and compared them with stool. Microbiome characterization was performed using 16S rRNA gene sequencing. RESULTS: The oral microbiome displayed IBD-associated dysbiosis, in a site- and taxa-specific manner. Plaque samples depicted a relatively severe degree of dysbiosis, and the disease-associated dysbiotic bacterial groups were predominantly the members of the phylum Firmicutes. Our 16S rRNA gene analyses show that oral microbiota can distinguish patients with IBD from healthy control patients, with salivary microbiota performing the best, closely matched by stool and other oral sites. Longitudinal profiles of microbial composition suggest that some taxa are more consistently perturbed than others, preferentially in a site-dependent fashion. CONCLUSIONS: Collectively, these data indicate the potential of using oral microbial profiles in screening and monitoring patients with IBD. Furthermore, these results support the importance of spatial and longitudinal microbiome sampling to interpret disease-associated dysbiotic states and eventually to gain insights into disease pathogenesis.

Structure and complexity of a bacterial transcriptome.

Although gene expression has been studied in bacteria for decades, many aspects of the bacterial transcriptome remain poorly understood. Transcript structure, operon linkages, and information on absolute abundance all provide valuable insights into gene function and regulation, but none has ever been determined on a genome-wide scale for any bacterium. Indeed, these aspects of the prokaryotic transcriptome have been explored on a large scale in only a few instances, and consequently little is known about the absolute composition of the mRNA population within a bacterial cell. Here we report the use of a high-throughput sequencing-based approach in assembling the first comprehensive, single-nucleotide resolution view of a bacterial transcriptome. We sampled the Bacillus anthracis transcriptome under a variety of growth conditions and showed that the data provide an accurate and high-resolution map of transcript start sites and operon structure throughout the genome. Further, the sequence data identified previously nonannotated regions with significant transcriptional activity and enhanced the accuracy of existing genome annotations. Finally, our data provide estimates of absolute transcript abundance and suggest that there is significant transcriptional heterogeneity within a clonal, synchronized bacterial population. Overall, our results offer an unprecedented view of gene expression and regulation in a bacterial cell.

Combining microarray-based genomic selection (MGS) with the Illumina Genome Analyzer platform to sequence diploid target regions.

Novel methods of targeted sequencing of unique regions from complex eukaryotic genomes have generated a great deal of excitement, but critical demonstrations of these methods efficacy with respect to diploid genotype calling and experimental variation are lacking. To address this issue, we optimized microarray-based genomic selection (MGS) for use with the Illumina Genome Analyzer (IGA). A set of 202 fragments (304 kb total) contained within a 1.7 Mb genomic region on human chromosome X were MGS/IGA sequenced in ten female HapMap samples generating a total of 2.4 GB of DNA sequence. At a minimum coverage threshold of 5X, 93.9% of all bases and 94.9% of segregating sites were called, while 57.7% of bases (57.4% of segregating sites) were called at a 50X threshold. Data accuracy at known segregating sites was 98.9% at 5X coverage, rising to 99.6% at 50X coverage. Accuracy at homozygous sites was 98.7% at 5X sequence coverage and 99.5% at 50X coverage. Although accuracy at heterozygous sites was modestly lower, it was still over 92% at 5X coverage and increased to nearly 97% at 50X coverage. These data provide the first demonstration that MGS/IGA sequencing can generate the very high quality sequence data necessary for human genetics research. All sequences generated in this study have been deposited in NCBI Short Read Archive (http://www.ncbi.nlm.nih.gov/Traces/sra, Accession # SRA007913).

Neutrophil GM-CSF signaling in inflammatory bowel disease patients is influenced by non-coding genetic variants.

Neutrophil dysfunction and GM-CSF auto-antibodies are observed in pediatric and adult patients with Crohn's disease (CD). We associated damaging coding variants with low GM-CSF induced STAT5 stimulation index (GMSI) in pediatric CD patients and implicated variation of neutrophil GM-CSF signaling in cell function and disease complications. Because many CD patients with low GMSI do not carry damaging coding mutations, we sought to test the hypothesis that non-coding variants contribute to this phenotype. We enrolled, performed whole genome sequencing, and measured the GMSI in 77 CD and ulcerative colitis (UC) patients (24 low and 53 normal GMSI). We identified 4 non-coding variants (rs3808851, rs10974787, rs10974788 and rs10974789) in RCL1 significantly associated with variation of GMSI level (p < 0.011). They were validated in two independent cohorts with: RNAseq data (n = 50) and blood eQTL dataset (n = 31,684). These variants are in LD and affect expression of JAK2 (p 0.005 to 0.013), RCL1 (p 8.17E-13 to 2.98E-11) and AK3 (p 2.00E-68 to 3.03E-55) genes. Additionally, they influence proteins involved in differentiation of gut epithelium, inflammation, and immune system regulation. In summary, our study outlines the contribution of non-coding variants in neutrophil GM-CSF signaling and the potential importance of RCL1 and AK3 in CD pathogenesis.

Microarray-based mutation detection in the dystrophin gene.

Duchenne and Becker muscular dystrophies (DMD and BMD) are X-linked recessive neuromuscular disorders caused by mutations in the dystrophin gene affecting approximately 1 in 3,500 males. The human dystrophin gene spans>2,200 kb, or roughly 0.1% of the genome, and is composed of 79 exons. The mutational spectrum of disease-causing alleles, including exonic copy number variations (CNVs), is complex. Deletions account for approximately 65% of DMD mutations and 85% of BMD mutations. Duplications occur in approximately 6 to 10% of males with either DMD or BMD. The remaining 30 to 35% of mutations consist of small deletions, insertions, point mutations, or splicing mutations, most of which introduce a premature stop codon. Laboratory analysis of dystrophin can be used to confirm a clinical diagnosis of DMD, characterize the type of dystrophin mutation, and perform prenatal testing and carrier testing for females. Current dystrophin diagnostic assays involve a variety of methodologies, including multiplex PCR, Southern blot analysis, multiplex ligation-dependent probe amplification (MLPA), detection of virtually all mutations-SSCP (DOVAM-S), and single condition amplification/internal primer sequencing (SCAIP); however, these methods are time-consuming, laborious, and do not accurately detect duplication mutations in the dystrophin gene. Furthermore, carrier testing in females is often difficult when a related affected male is unavailable. Here we describe the development, design, validation, and implementation of a high-resolution comparative genomic hybridization (CGH) microarray-based approach capable of accurately detecting both deletions and duplications in the dystrophin gene. This assay can be readily adopted by clinical molecular testing laboratories and represents a rapid, cost-effective approach for screening a large gene, such as dystrophin.

Enhanced Contribution of HLA in Pediatric Onset Ulcerative Colitis.

Background: The genetic contributions to pediatric onset ulcerative colitis (UC), characterized by severe disease and extensive colonic involvement, are largely unknown. In adult onset UC, Genome Wide Association Study (GWAS) has identified numerous loci, most of which have a modest susceptibility risk (OR 0.84-1.14), with the exception of the human leukocyte antigen (HLA) region on Chromosome 6 (OR 3.59). Method: To study the genetic contribution to exclusive pediatric onset UC, a GWAS was performed on 466 cases with 2099 healthy controls using UK Biobank array. SNP2HLA was used to impute classical HLA alleles and their corresponding amino acids, and the results are compared with adult onset UC. Results: HLA explained the almost entire association signal, dominated with 191 single nucleotide polymorphisms (SNPs) (p = 5 x 10-8 to 5 x 10-10). Although very small effects, established SNPs in adult onset UC loci had similar direction and magnitude in pediatric onset UC. SNP2HLA imputation identified HLA-DRB1*0103 (odds ratio [OR] = 6.941, p = 1.92*10-13) as the most significant association for pediatric UC compared with adult onset UC (OR = 3.59). Further conditioning showed independent effects for HLA-DRB1*1301 (OR = 2.25, p = 7.92*10-9) and another SNP rs17188113 (OR = 0.48, p = 7.56*10-9). Two HLA-DRB1 causal alleles are shared with adult onset UC, while at least 2 signals are unique to pediatric UC. Subsequent stratified analyses indicated that HLA-DRB1*0103 has stronger association for extensive disease (E4: OR = 8.28, p = 4.66x10-10) and female gender (OR = 8.85, p = 4.82x10-13). Conclusion: In pediatric onset UC, the HLA explains almost the entire genetic associations. In addition, the HLA association is approximately twice as strong in pediatric UC compared with adults, due to a combination of novel and shared effects. We speculate the paramount importance of antigenic stimulation either by infectious or noninfectious stimuli as a causal event in pediatric UC onset.