IBD risk loci are enriched in multigenic regulatory modules encompassing putative causative genes.

GWAS have identified >200 risk loci for Inflammatory Bowel Disease (IBD). The majority of disease associations are known to be driven by regulatory variants. To identify the putative causative genes that are perturbed by these variants, we generate a large transcriptome data set (nine disease-relevant cell types) and identify 23,650 cis-eQTL. We show that these are determined by ∼9720 regulatory modules, of which ∼3000 operate in multiple tissues and ∼970 on multiple genes. We identify regulatory modules that drive the disease association for 63 of the 200 risk loci, and show that these are enriched in multigenic modules. Based on these analyses, we resequence 45 of the corresponding 100 candidate genes in 6600 Crohn disease (CD) cases and 5500 controls, and show with burden tests that they include likely causative genes. Our analyses indicate that ≥10-fold larger sample sizes will be required to demonstrate the causality of individual genes using this approach.

Analysis of Genes Associated With Monogenic Primary Immunodeficiency Identifies Rare Variants in XIAP in Patients With Crohn's Disease.

BACKGROUND & AIMS: A few rare monogenic primary immunodeficiencies (PIDs) are characterized by chronic intestinal inflammation that resembles Crohn's disease (CD). We investigated whether 23 genes associated with 10 of these monogenic disorders contain common, low-frequency, or rare variants that increase risk for CD. METHODS: Common and low frequency variants in 1 Mb loci centered on the candidate genes were analyzed using meta-data corresponding to genotypes of approximately 17,000 patients with CD or without CD (controls) in Europe. The contribution of rare variants was assessed by high-throughput sequencing of 4750 individuals, including 660 early-onset and/or familial cases among the 2390 patients with CD. Variants were expressed from vectors in SW480 or HeLa cells and functions of their products were analyzed in immunofluorescence, luciferase, immunoprecipitation, and immunoblot assays. RESULTS: We reproduced the association of the interleukin 10 locus with CD (P = .007), although none of the significantly associated variants modified the coding sequence of interleukin 10. We found XIAP to be significantly enriched for rare coding mutations in patients with CD vs controls (P = .02). We identified 4 previously unreported missense variants associated with CD. Variants in XIAP cause the PID X-linked lymphoproliferative disease type 2, yet none of the carriers of these variants had all the clinical features of X-linked lymphoproliferative disease type 2. Identified XIAP variants S123N, R233Q, and P257A were associated with an impaired activation of NOD2 signaling after muramyl dipeptide stimulation. CONCLUSIONS: In a systematic analysis of variants in 23 PID-associated genes, we confirmed the association of variants in XIAP with CD. Further screenings for CD-associated variants and analyses of their functions could increase our understanding of the relationship between PID-associated genes and CD pathogenesis.

Fine-mapping inflammatory bowel disease loci to single-variant resolution.

Inflammatory bowel diseases are chronic gastrointestinal inflammatory disorders that affect millions of people worldwide. Genome-wide association studies have identified 200 inflammatory bowel disease-associated loci, but few have been conclusively resolved to specific functional variants. Here we report fine-mapping of 94 inflammatory bowel disease loci using high-density genotyping in 67,852 individuals. We pinpoint 18 associations to a single causal variant with greater than 95% certainty, and an additional 27 associations to a single variant with greater than 50% certainty. These 45 variants are significantly enriched for protein-coding changes (n = 13), direct disruption of transcription-factor binding sites (n = 3), and tissue-specific epigenetic marks (n = 10), with the last category showing enrichment in specific immune cells among associations stronger in Crohn's disease and in gut mucosa among associations stronger in ulcerative colitis. The results of this study suggest that high-resolution fine-mapping in large samples can convert many discoveries from genome-wide association studies into statistically convincing causal variants, providing a powerful substrate for experimental elucidation of disease mechanisms.

Characterization of bacteria in biopsies of colon and stools by high throughput sequencing of the V2 region of bacterial 16S rRNA gene in human.

BACKGROUND: The characterization of the human intestinal microflora and their interactions with the host have been identified as key components in the study of intestinal disorders such as inflammatory bowel diseases. High-throughput sequencing has enabled culture-independent studies to deeply analyze bacteria in the gut. It is possible with this technology to systematically analyze links between microbes and the genetic constitution of the host, such as DNA polymorphisms and methylation, and gene expression. METHODS AND FINDINGS: In this study the V2 region of the bacterial 16S ribosomal RNA (rRNA) gene using 454 pyrosequencing from seven anatomic regions of human colon and two types of stool specimens were analyzed. The study examined the number of reads needed to ascertain differences between samples, the effect of DNA extraction procedures and PCR reproducibility, and differences between biopsies and stools in order to design a large scale systematic analysis of gut microbes. It was shown (1) that sequence coverage lower than 1,000 reads influenced quantitative and qualitative differences between samples measured by UniFrac distances. Distances between samples became stable after 1,000 reads. (2) Difference of extracted bacteria was observed between the two DNA extraction methods. In particular, Firmicutes Bacilli were not extracted well by one method. (3) Quantitative and qualitative difference in bacteria from ileum to rectum colon were not observed, but there was a significant positive trend between distances within colon and quantitative differences. Between sample type, biopsies or stools, quantitative and qualitative differences were observed. CONCLUSIONS: Results of human colonic bacteria analyzed using high-throughput sequencing were highly dependent on the experimental design, especially the number of sequence reads, DNA extraction method, and sample type.

A relict bank vole lineage highlights the biogeographic history of the Pyrenean region in Europe.

The Pyrenean region exhibits high levels of endemism suggesting a major contribution to the phylogeography of European species. But, to date, the role of the Pyrenees and surrounding areas as a glacial refugium for temperate species remains poorly explored. In the current study, we investigated the biogeographic role of the Pyrenean region through the analyses of genetic polymorphism and morphology of a typical forest-dwelling small mammal, the bank vole (Myodes glareolus). Analyses of the mitochondrial cytochrome b gene and the third upper molar (M(3)) show a complex phylogeographic structure in the Pyrenean region with at least three distinct lineages: the Western European, Spanish and Basque lineages. The Basque lineage in the northwestern (NW) Pyrenees was identified as a new clearly differentiated and geographically localized bank vole lineage in Europe. The average M(3) shape of Basque bank voles suggests morphological differentiation but also restricted genetic exchanges with other populations. Our genetic and morphological results as well as palaeo-environmental and fossils records support the hypothesis of a new glacial refugium in Europe situated in the NW Pyrenees. The permissive microclimatic conditions that prevailed for a long time in this region may have allowed the survival of temperate species, including humans. Moreover, local differentiation around the Pyrenees is favoured by the opportunity for populations to track the shift of the vegetation belt in altitude rather than in latitude. The finding of the Basque lineage is in agreement with the high level of endemic taxa reported in the NW Pyrenees.

A northern glacial refugium for bank voles (Clethrionomys glareolus).

There is controversy and uncertainty on how far north there were glacial refugia for temperate species during the Pleistocene glaciations and in the extent of the contribution of such refugia to present-day populations. We examined these issues using phylogeographic analysis of a European woodland mammal, the bank vole (Clethrionomys glareolus). A Bayesian coalescence analysis indicates that a bank vole population survived the height of the last glaciation (approximately 25,000-10,000 years B.P.) in the vicinity of the Carpathians, a major central European mountain chain well north of the Mediterranean areas typically regarded as glacial refugia for temperate species. Parameter estimates from the fitted isolation with migration model show that the divergence of the Carpathian population started at least 22,000 years ago, and it was likely followed by only negligible immigration from adjacent regions, suggesting the persistence of bank voles in the Carpathians through the height of the last glaciation. On the contrary, there is clear evidence for gene flow out of the Carpathians, demonstrating the contribution of the Carpathian population to the colonization of Europe after the Pleistocene. These findings are consistent with data from animal and plant fossils recovered in the Carpathians and provide the clearest phylogeographic evidence to date of a northern glacial refugium for temperate species in Europe.