Stephen L. Guthery

Genome-wide meta-analysis increases to 71 the number of confirmed Crohn's disease susceptibility loci

We undertook a meta-analysis of six Crohns disease genome-wide association studies (GWAS) comprising 6,333 affected individuals (cases) and 15,056 controls and followed up the top association signals in 15,694 cases, 14,026 controls and 414 parent-offspring trios. We identified 30 new susceptibility loci meeting genome-wide significance (P < 5 × 10−8). A series of in silico analyses highlighted particular genes within these loci and, together with manual curation, implicated functionally interesting candidate genes including SMAD3, ERAP2, IL10, IL2RA, TYK2, FUT2, DNMT3A, DENND1B, BACH2 and TAGAP. Combined with previously confirmed loci, these results identify 71 distinct loci with genome-wide significant evidence for association with Crohns disease.

Common variants at five new loci associated with early-onset inflammatory bowel disease

The inflammatory bowel diseases (IBD) Crohns disease and ulcerative colitis are common causes of morbidity in children and young adults in the western world. Here we report the results of a genome-wide association study in early-onset IBD involving 3,426 affected individuals and 11,963 genetically matched controls recruited through international collaborations in Europe and North America, thereby extending the results from a previous study of 1,011 individuals with early-onset IBD. We have identified five new regions associated with early-onset IBD susceptibility, including 16p11 near the cytokine gene IL27 (rs8049439, P = 2.41 × 10−9), 22q12 (rs2412973, P = 1.55 × 10−9), 10q22 (rs1250550, P = 5.63 × 10−9), 2q37 (rs4676410, P = 3.64 × 10−8) and 19q13.11 (rs10500264, P = 4.26 × 10−10). Our scan also detected associations at 23 of 32 loci previously implicated in adult-onset Crohns disease and at 8 of 17 loci implicated in adult-onset ulcerative colitis, highlighting the close pathogenetic relationship between early- and adult-onset IBD.

Loci on 20q13 and 21q22 are associated with pediatric-onset inflammatory bowel disease

Inflammatory bowel disease (IBD) is a common inflammatory disorder with complex etiology that involves both genetic and environmental triggers, including but not limited to defects in bacterial clearance, defective mucosal barrier and persistent dysregulation of the immune response to commensal intestinal bacteria. IBD is characterized by two distinct phenotypes: Crohns disease (CD) and ulcerative colitis (UC). Previously reported GWA studies have identified genetic variation accounting for a small portion of the overall genetic susceptibility to CD and an even smaller contribution to UC pathogenesis. We hypothesized that stratification of IBD by age of onset might identify additional genes associated with IBD. To that end, we carried out a GWA analysis in a cohort of 1,011 individuals with pediatric-onset IBD and 4,250 matched controls. We identified and replicated significantly associated, previously unreported loci on chromosomes 20q13 (rs2315008[T] and rs4809330[A]; P = 6.30 × 10−8 and 6.95 × 10−8, respectively; odds ratio (OR) = 0.74 for both) and 21q22 (rs2836878[A]; P = 6.01 × 10−8; OR = 0.73), located close to the TNFRSF6B and PSMG1 genes, respectively.

Diverse Genome-wide Association Studies Associate the IL12/IL23 Pathway with Crohn Disease

Previous genome-wide association (GWA) studies typically focus on single-locus analysis, which may not have the power to detect the majority of genuinely associated loci. Here, we applied pathway analysis using Affymetrix SNP genotype data from the Wellcome Trust Case Control Consortium (WTCCC) and uncovered significant association between Crohn Disease (CD) and the IL12/IL23 pathway, harboring 20 genes (p = 8 x 10(-5)). Interestingly, the pathway contains multiple genes (IL12B and JAK2) or homologs of genes (STAT3 and CCR6) that were recently identified as genuine susceptibility genes only through meta-analysis of several GWA studies. In addition, the pathway contains other susceptibility genes for CD, including IL18R1, JUN, IL12RB1, and TYK2, which do not reach genome-wide significance by single-marker association tests. The observed pathway-specific association signal was subsequently replicated in three additional GWA studies of European and African American ancestry generated on the Illumina HumanHap550 platform. Our study suggests that examination beyond individual SNP hits, by focusing on genetic networks and pathways, is important to unleashing the true power of GWA studies.

Dense genotyping of immune-related disease regions identifies 14 new susceptibility loci for juvenile idiopathic arthritis

We used the Immunochip array to analyze 2,816 individuals with juvenile idiopathic arthritis (JIA), comprising the most common subtypes (oligoarticular and rheumatoid factor–negative polyarticular JIA), and 13,056 controls. We confirmed association of 3 known JIA risk loci (the human leukocyte antigen (HLA) region, PTPN22 and PTPN2) and identified 14 loci reaching genome-wide significance (P < 5 × 10−8) for the first time. Eleven additional new regions showed suggestive evidence of association with JIA (P < 1 × 10−6). Dense mapping of loci along with bioinformatics analysis refined the associations to one gene in each of eight regions, highlighting crucial pathways, including the interleukin (IL)-2 pathway, in JIA disease pathogenesis. The entire Immunochip content, the HLA region and the top 27 loci (P < 1 × 10−6) explain an estimated 18, 13 and 6% of the risk of JIA, respectively. In summary, this is the largest collection of JIA cases investigated so far and provides new insight into the genetic basis of this childhood autoimmune disease.

Comparative genetic analysis of inflammatory bowel disease and type 1 diabetes implicates multiple loci with opposite effects

Inflammatory bowel disease, including Crohns disease (CD) and ulcerative colitis (UC), and type 1 diabetes (T1D) are autoimmune diseases that may share common susceptibility pathways. We examined known susceptibility loci for these diseases in a cohort of 1689 CD cases, 777 UC cases, 989 T1D cases and 6197 shared control subjects of European ancestry, who were genotyped by the Illumina HumanHap550 SNP arrays. We identified multiple previously unreported or unconfirmed disease associations, including known CD loci (ICOSLG and TNFSF15) and T1D loci (TNFAIP3) that confer UC risk, known UC loci (HERC2 and IL26) that confer T1D risk and known UC loci (IL10 and CCNY) that confer CD risk. Additionally, we show that T1D risk alleles residing at the PTPN22, IL27, IL18RAP and IL10 loci protect against CD. Furthermore, the strongest risk alleles for T1D within the major histocompatibility complex (MHC) confer strong protection against CD and UC; however, given the multi-allelic nature of the MHC haplotypes, sequencing of the MHC locus will be required to interpret this observation. These results extend our current knowledge on genetic variants that predispose to autoimmunity, and suggest that many loci involved in autoimmunity may be under a balancing selection due to antagonistic pleiotropic effect. Our analysis implies that variants with opposite effects on different diseases may facilitate the maintenance of common susceptibility alleles in human populations, making autoimmune diseases especially amenable to genetic dissection by genome-wide association studies.

Increased Effectiveness of Early Therapy With Anti-Tumor Necrosis Factor-α vs an Immunomodulator in Children With Crohn's Disease

BACKGROUND & AIMS Standard therapy for children newly diagnosed with Crohns disease (CD) includes early administration of immunomodulators after initial treatment with corticosteroids. We compared the effectiveness of early (≤3 mo after diagnosis) treatment with an anti-tumor necrosis factor (TNF)α with that of an immunomodulator in attaining clinical remission and facilitating growth of pediatric patients. METHODS We analyzed data from the RISK study, an observational research program that enrolled patients younger than age 17 diagnosed with inflammatory (nonpenetrating, nonstricturing) CD from 2008 through 2012 at 28 pediatric gastroenterology centers in North America. Patients were managed by physician dictate. From 552 children (median age, 11.8 y; 61% male; 63% with pediatric CD activity index scores >30; and median C-reactive protein level 5.6-fold the upper limit of normal), we used propensity score methodology to identify 68 triads of patients matched for baseline characteristics who were treated with early anti-TNFα therapy, early immunomodulator, or no early immunotherapy. We evaluated relationships among therapies, corticosteroid and surgery-free remission (pediatric CD activity index scores, ≤10), and growth at 1 year for 204 children. Treatment after 3 months was a covariate. RESULTS Early treatment with anti-TNFα was superior to early treatment with an immunomodulator (85.3% vs 60.3% in remission; relative risk, 1.41; 95% confidence interval [CI], 1.14-1.75; P = .0017), whereas early immunomodulator therapy was no different than no early immunotherapy (60.3% vs 54.4% in remission; relative risk, 1.11; 95% CI, 0.83-1.48; P = .49) in achieving remission at 1 year. Accounting for therapy after 3 months, early treatment with anti-TNFα remained superior to early treatment with an immunomodulator (relative risk, 1.51; 95% CI, 1.20-1.89; P = .0004), whereas early immunomodulator therapy was no different than no early immunotherapy (relative risk, 1.00; 95% CI, 0.75-1.34; P = .99). The mean height z-score increased compared with baseline only in the early anti-TNFα group. CONCLUSIONS In children newly diagnosed with comparably severe CD, early monotherapy with anti-TNFα produced better overall clinical and growth outcomes at 1 year than early monotherapy with an immunomodulator. Further data will be required to best identify children most likely to benefit from early treatment with anti-TNFα therapy.

Fine-scaled human genetic structure revealed by SNP microarrays

We report an analysis of more than 240,000 loci genotyped using the Affymetrix SNP microarray in 554 individuals from 27 worldwide populations in Africa, Asia, and Europe. To provide a more extensive and complete sampling of human genetic variation, we have included caste and tribal samples from two states in South India, Daghestanis from eastern Europe, and the Iban from Malaysia. Consistent with observations made by Charles Darwin, our results highlight shared variation among human populations and demonstrate that much genetic variation is geographically continuous. At the same time, principal components analyses reveal discernible genetic differentiation among almost all identified populations in our sample, and in most cases, individuals can be clearly assigned to defined populations on the basis of SNP genotypes. All individuals are accurately classified into continental groups using a model-based clustering algorithm, but between closely related populations, genetic and self-classifications conflict for some individuals. The 250K data permitted high-level resolution of genetic variation among Indian caste and tribal populations and between highland and lowland Daghestani populations. In particular, upper-caste individuals from Tamil Nadu and Andhra Pradesh form one defined group, lower-caste individuals from these two states form another, and the tribal Irula samples form a third. Our results emphasize the correlation of genetic and geographic distances and highlight other elements, including social factors that have contributed to population structure.

Maximum-likelihood estimation of recent shared ancestry (ERSA)

Accurate estimation of recent shared ancestry is important for genetics, evolution, medicine, conservation biology, and forensics. Established methods estimate kinship accurately for first-degree through third-degree relatives. We demonstrate that chromosomal segments shared by two individuals due to identity by descent (IBD) provide much additional information about shared ancestry. We developed a maximum-likelihood method for the estimation of recent shared ancestry (ERSA) from the number and lengths of IBD segments derived from high-density SNP or whole-genome sequence data. We used ERSA to estimate relationships from SNP genotypes in 169 individuals from three large, well-defined human pedigrees. ERSA is accurate to within one degree of relationship for 97% of first-degree through fifth-degree relatives and 80% of sixth-degree and seventh-degree relatives. We demonstrate that ERSAs statistical power approaches the maximum theoretical limit imposed by the fact that distant relatives frequently share no DNA through a common ancestor. ERSA greatly expands the range of relationships that can be estimated from genetic data and is implemented in a freely available software package.

Variants in TNFAIP3, STAT4, and C12orf30 loci associated with multiple autoimmune diseases are also associated with juvenile idiopathic arthritis.

OBJECTIVE Subtypes of juvenile idiopathic arthritis (JIA) share phenotypic features with other autoimmune disorders. We investigated several genetic variants associated with rheumatoid arthritis (RA) and other autoimmune disorders for association with JIA to test the hypothesis that clinically distinct phenotypes share common genetic susceptibility factors. METHODS Cases were 445 children with JIA, and controls were 643 healthy adults. Using the TaqMan assay, subjects were genotyped for 8 single-nucleotide polymorphisms in 7 loci including rs10499194 and rs6920220 in the TNFAIP3 locus, rs6679677 in the RSBN1 locus, rs17696736 in the C12orf30 locus, rs3761847 in the TRAF1/C5 locus, rs2104286 in the IL2RA locus, rs7574865 in the STAT4 locus, and rs2542151 in the PTPN2 locus. Alleles and genotypes were analyzed for association with JIA and JIA subtypes. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated. RESULTS The strongest associations with JIA risk or protection were observed for TNFAIP3 variants rs10499194 (OR 0.74 [95% CI 0.61-0.91], P < 0.004) and rs6920220 (OR 1.30 [95% CI 1.05-1.61], P = 0.015). We also observed associations between JIA and both STAT4 (OR 1.24 [95% CI 1.02-1.51], P = 0.029) and C12orf30 (OR 1.20 [95% CI 1.01-1.43], P = 0.041) variants. The PTPN2 variant rs2542151 deviated from Hardy-Weinberg equilibrium and was excluded from analyses. Variants in IL2RA, TRAF1/C5, and RSBN1 were not associated with JIA. After stratification by JIA subtype, the TNFAIP3 and C12orf30 variants were associated with oligoarticular JIA, while the STAT4 variant was associated primarily with polyarticular JIA. CONCLUSION We have demonstrated associations between JIA and variants in the TNFAIP3, STAT4, and C12orf30 regions that have previously shown associations with other autoimmune diseases, including RA and systemic lupus erythematosus. Our results suggest that clinically distinct autoimmune phenotypes share common genetic susceptibility factors.