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.

A novel unconventional T cell population enriched in Crohn's disease.

OBJECTIVE: One of the current hypotheses to explain the proinflammatory immune response in IBD is a dysregulated T cell reaction to yet unknown intestinal antigens. As such, it may be possible to identify disease-associated T cell clonotypes by analysing the peripheral and intestinal T-cell receptor (TCR) repertoire of patients with IBD and controls. DESIGN: We performed bulk TCR repertoire profiling of both the TCR alpha and beta chains using high-throughput sequencing in peripheral blood samples of a total of 244 patients with IBD and healthy controls as well as from matched blood and intestinal tissue of 59 patients with IBD and disease controls. We further characterised specific T cell clonotypes via single-cell RNAseq. RESULTS: We identified a group of clonotypes, characterised by semi-invariant TCR alpha chains, to be significantly enriched in the blood of patients with Crohn's disease (CD) and particularly expanded in the CD8+ T cell population. Single-cell RNAseq data showed an innate-like phenotype of these cells, with a comparable gene expression to unconventional T cells such as mucosal associated invariant T and natural killer T (NKT) cells, but with distinct TCRs. CONCLUSIONS: We identified and characterised a subpopulation of unconventional Crohn-associated invariant T (CAIT) cells. Multiple evidence suggests these cells to be part of the NKT type II population. The potential implications of this population for CD or a subset thereof remain to be elucidated, and the immunophenotype and antigen reactivity of CAIT cells need further investigations in future studies.

Genome-wide analysis of 944 133 individuals provides insights into the etiology of haemorrhoidal disease.

OBJECTIVE: Haemorrhoidal disease (HEM) affects a large and silently suffering fraction of the population but its aetiology, including suspected genetic predisposition, is poorly understood. We report the first genome-wide association study (GWAS) meta-analysis to identify genetic risk factors for HEM to date. DESIGN: We conducted a GWAS meta-analysis of 218 920 patients with HEM and 725 213 controls of European ancestry. Using GWAS summary statistics, we performed multiple genetic correlation analyses between HEM and other traits as well as calculated HEM polygenic risk scores (PRS) and evaluated their translational potential in independent datasets. Using functional annotation of GWAS results, we identified HEM candidate genes, which differential expression and coexpression in HEM tissues were evaluated employing RNA-seq analyses. The localisation of expressed proteins at selected loci was investigated by immunohistochemistry. RESULTS: We demonstrate modest heritability and genetic correlation of HEM with several other diseases from the GI, neuroaffective and cardiovascular domains. HEM PRS validated in 180 435 individuals from independent datasets allowed the identification of those at risk and correlated with younger age of onset and recurrent surgery. We identified 102 independent HEM risk loci harbouring genes whose expression is enriched in blood vessels and GI tissues, and in pathways associated with smooth muscles, epithelial and endothelial development and morphogenesis. Network transcriptomic analyses highlighted HEM gene coexpression modules that are relevant to the development and integrity of the musculoskeletal and epidermal systems, and the organisation of the extracellular matrix. CONCLUSION: HEM has a genetic component that predisposes to smooth muscle, epithelial and connective tissue dysfunction.

Protein-coding variants contribute to the risk of atopic dermatitis and skin-specific gene expression.

BACKGROUND: Fifteen percent of atopic dermatitis (AD) liability-scale heritability could be attributed to 31 susceptibility loci identified by using genome-wide association studies, with only 3 of them (IL13, IL-6 receptor [IL6R], and filaggrin [FLG]) resolved to protein-coding variants. OBJECTIVE: We examined whether a significant portion of unexplained AD heritability is further explained by low-frequency and rare variants in the gene-coding sequence. METHODS: We evaluated common, low-frequency, and rare protein-coding variants using exome chip and replication genotype data of 15,574 patients and 377,839 control subjects combined with whole-transcriptome data on lesional, nonlesional, and healthy skin samples of 27 patients and 38 control subjects. RESULTS: An additional 12.56% (SE, 0.74%) of AD heritability is explained by rare protein-coding variation. We identified docking protein 2 (DOK2) and CD200 receptor 1 (CD200R1) as novel genome-wide significant susceptibility genes. Rare coding variants associated with AD are further enriched in 5 genes (IL-4 receptor [IL4R], IL13, Janus kinase 1 [JAK1], JAK2, and tyrosine kinase 2 [TYK2]) of the IL13 pathway, all of which are targets for novel systemic AD therapeutics. Multiomics-based network and RNA sequencing analysis revealed DOK2 as a central hub interacting with, among others, CD200R1, IL6R, and signal transducer and activator of transcription 3 (STAT3). Multitissue gene expression profile analysis for 53 tissue types from the Genotype-Tissue Expression project showed that disease-associated protein-coding variants exert their greatest effect in skin tissues. CONCLUSION: Our discoveries highlight a major role of rare coding variants in AD acting independently of common variants. Further extensive functional studies are required to detect all potential causal variants and to specify the contribution of the novel susceptibility genes DOK2 and CD200R1 to overall disease susceptibility.

A genome-wide association scan of nonsynonymous SNPs identifies a susceptibility variant for Crohn disease in ATG16L1.

We performed a genome-wide association study of 19,779 nonsynonymous SNPs in 735 individuals with Crohn disease and 368 controls. A total of 7,159 of these SNPs were informative. We followed up on all 72 SNPs with P 0.4), these data suggest that the underlying biological process may be specific to Crohn disease.

XIAP variants in male Crohn's disease.

OBJECTIVE: The genetic basis of inflammatory bowel disease (IBD) is incompletely understood. The aim of this study was to identify rare genetic variants involved in the pathogenesis of IBD. DESIGN: Exome sequencing and immunological profiling were performed in a patient with early onset Crohn's disease (CD). The coding region of the gene encoding X-linked inhibitor of apoptosis protein (XIAP) was sequenced in samples of 275 paediatric IBD and 1047 adult-onset CD patients. XIAP genotyping was performed in samples of 2680 IBD patients and 2864 healthy controls. Functional effects of the variants identified were investigated in primary cells and cultured cell lines. RESULTS: Our results demonstrate the frequent occurrence of private variants in XIAP in about four percent of male patients with paediatric-onset CD. While XIAP mutations are known to be associated with the primary immunodeficiency (PID) X-linked lymphoproliferative disease type 2 (XLP2), CD patients described here exhibited intestinal inflammation in the absence of XLP2 and harboured a spectrum of mutations partially distinct from that observed in XLP2. The majority of XIAP variants identified was associated with a selective defect in NOD1/2 signalling, impaired NOD1/2-mediated activation of NF-κB, and altered NF-κB-dependent cytokine production. CONCLUSIONS: This study reveals the unanticipated, frequent occurrence of XIAP variants in male paediatric-onset CD. The link between XIAP and NOD1/2, and the association of XIAP variants with XLP2, support the concept of PID in a subset of IBD patients. Moreover, these studies provide a rationale for the implementation of XIAP sequencing in clinical diagnostics in male patients with severe CD.

Association between variants of PRDM1 and NDP52 and Crohn's disease, based on exome sequencing and functional studies.

BACKGROUND & AIMS: Genome-wide association studies (GWAS) have identified 140 Crohn's disease (CD) susceptibility loci. For most loci, the variants that cause disease are not known and the genes affected by these variants have not been identified. We aimed to identify variants that cause CD through detailed sequencing, genetic association, expression, and functional studies. METHODS: We sequenced whole exomes of 42 unrelated subjects with CD and 5 healthy subjects (controls) and then filtered single nucleotide variants by incorporating association results from meta-analyses of CD GWAS and in silico mutation effect prediction algorithms. We then genotyped 9348 subjects with CD, 2868 subjects with ulcerative colitis, and 14,567 control subjects and associated variants analyzed in functional studies using materials from subjects and controls and in vitro model systems. RESULTS: We identified rare missense mutations in PR domain-containing 1 (PRDM1) and associated these with CD. These mutations increased proliferation of T cells and secretion of cytokines on activation and increased expression of the adhesion molecule L-selectin. A common CD risk allele, identified in GWAS, correlated with reduced expression of PRDM1 in ileal biopsy specimens and peripheral blood mononuclear cells (combined P = 1.6 × 10(-8)). We identified an association between CD and a common missense variant, Val248Ala, in nuclear domain 10 protein 52 (NDP52) (P = 4.83 × 10(-9)). We found that this variant impairs the regulatory functions of NDP52 to inhibit nuclear factor κB activation of genes that regulate inflammation and affect the stability of proteins in Toll-like receptor pathways. CONCLUSIONS: We have extended the results of GWAS and provide evidence that variants in PRDM1 and NDP52 determine susceptibility to CD. PRDM1 maps adjacent to a CD interval identified in GWAS and encodes a transcription factor expressed by T and B cells. NDP52 is an adaptor protein that functions in selective autophagy of intracellular bacteria and signaling molecules, supporting the role of autophagy in the pathogenesis of CD.

Intragenic allele pyramiding combines different specificities of wheat Pm3 resistance alleles.

Some plant resistance genes occur as allelic series, with each member conferring specific resistance against a subset of pathogen races. In wheat, there are 17 alleles of the Pm3 gene. They encode nucleotide-binding (NB-ARC) and leucine-rich-repeat (LRR) domain proteins, which mediate resistance to distinct race spectra of powdery mildew. It is not known if specificities from different alleles can be combined to create resistance genes with broader specificity. Here, we used an approach based on avirulence analysis of pathogen populations to characterize the molecular basis of Pm3 recognition spectra. A large survey of mildew races for avirulence on the Pm3 alleles revealed that Pm3a has a resistance spectrum that completely contains that of Pm3f, but also extends towards additional races. The same is true for the Pm3b and Pm3c gene pair. The molecular analysis of these allelic pairs revealed a role of the NB-ARC protein domain in the efficiency of effector-dependent resistance. Analysis of the wild-type and chimeric Pm3 alleles identified single residues in the C-terminal LRR motifs as the main determinant of allele specificity. Variable residues of the N-terminal LRRs are necessary, but not sufficient, to confer resistance specificity. Based on these data, we constructed a chimeric Pm3 gene by intragenic allele pyramiding of Pm3d and Pm3e that showed the combined resistance specificity and, thus, a broader recognition spectrum compared with the parental alleles. Our findings support a model of stepwise evolution of Pm3 recognition specificities.

Variants in a novel epidermal collagen gene (COL29A1) are associated with atopic dermatitis.

Atopic dermatitis (AD) is a common chronic inflammatory skin disorder and a major manifestation of allergic disease. AD typically presents in early childhood often preceding the onset of an allergic airway disease, such as asthma or hay fever. We previously mapped a susceptibility locus for AD on Chromosome 3q21. To identify the underlying disease gene, we used a dense map of microsatellite markers and single nucleotide polymorphisms, and we detected association with AD. In concordance with the linkage results, we found a maternal transmission pattern. Furthermore, we demonstrated that the same families contribute to linkage and association. We replicated the association and the maternal effect in a large independent family cohort. A common haplotype showed strong association with AD (p = 0.000059). The associated region contained a single gene, COL29A1, which encodes a novel epidermal collagen. COL29A1 shows a specific gene expression pattern with the highest transcript levels in skin, lung, and the gastrointestinal tract, which are the major sites of allergic disease manifestation. Lack of COL29A1 expression in the outer epidermis of AD patients points to a role of collagen XXIX in epidermal integrity and function, the breakdown of which is a clinical hallmark of AD.

Missense variants in NOX1 and p22phox in a case of very-early-onset inflammatory bowel disease are functionally linked to NOD2.

Whole-genome and whole-exome sequencing of individual patients allow the study of rare and potentially causative genetic variation. In this study, we sequenced DNA of a trio comprising a boy with very-early-onset inflammatory bowel disease (veoIBD) and his unaffected parents. We identified a rare, X-linked missense variant in the NAPDH oxidase NOX1 gene (c.C721T, p.R241C) in heterozygous state in the mother and in hemizygous state in the patient. We discovered that, in addition, the patient was homozygous for a common missense variant in the CYBA gene (c.T214C, p.Y72H). CYBA encodes the p22phox protein, a cofactor for NOX1. Functional assays revealed reduced cellular ROS generation and antibacterial capacity of NOX1 and p22phox variants in intestinal epithelial cells. Moreover, the identified NADPH oxidase complex variants affected NOD2-mediated immune responses, and p22phox was identified as a novel NOD2 interactor. In conclusion, we detected missense variants in a veoIBD patient that disrupt the host response to bacterial challenges and reduce protective innate immune signaling via NOD2. We assume that the patient's individual genetic makeup favored disturbed intestinal mucosal barrier function.

Structure-function analysis of the NB-ARC domain of plant disease resistance proteins.

Resistance (R) proteins in plants are involved in pathogen recognition and subsequent activation of innate immune responses. Most resistance proteins contain a central nucleotide-binding domain. This so-called NB-ARC domain consists of three subdomains: NB, ARC1, and ARC2. The NB-ARC domain is a functional ATPase domain, and its nucleotide-binding state is proposed to regulate activity of the R protein. A highly conserved methionine-histidine-aspartate (MHD) motif is present at the carboxy-terminus of ARC2. An extensive mutational analysis of the MHD motif in the R proteins I-2 and Mi-1 is reported. Several novel autoactivating mutations of the MHD invariant histidine and conserved aspartate were identified. The combination of MHD mutants with autoactivating hydrolysis mutants in the NB subdomain showed that the autoactivation phenotypes are not additive. This finding indicates an important regulatory role for the MHD motif in the control of R protein activity. To explain these observations, a three-dimensional model of the NB-ARC domain of I-2 was built, based on the APAF-1 template structure. The model was used to identify residues important for I-2 function. Substitution of the selected residues resulted in the expected distinct phenotypes. Based on the model, it is proposed that the MHD motif fulfils the same function as the sensor II motif found in AAA+ proteins (ATPases associated with diverse cellular activities)-co-ordination of the nucleotide and control of subdomain interactions. The presented 3D model provides a framework for the formulation of hypotheses on how mutations in the NB-ARC exert their effects.

Compound heterozygous mutations in IL10RA combined with a complement factor properdin mutation in infantile-onset inflammatory bowel disease.

OBJECTIVES: Inflammatory bowel diseases (IBDs) are chronic and multifactorial diseases resulting from a complex interaction of host genetic factors and environmental stimuli. Although many genome-wide association studies have identified host genetic factors associated with IBD, rare Mendelian forms of IBD have been reported in patients with very early onset forms. Therefore, this study aimed to identify genetic variants associated with infantile-onset IBD. PARTICIPANTS AND METHODS: We obtained genomic DNA from whole blood samples of a male patient with infantile-onset IBD and nonconsanguineous Korean parents. Whole-exome sequencing was performed using trio samples. Then, we analyzed the data using susceptibility genes for monogenic forms of IBD and various immunodeficiencies and protein structural analysis. RESULTS: The patient who presented with oral aphthous ulcers at the age of 14 days suffered from severe colitis and was refractory to medical treatment. Compound heterozygous mutations in IL10RA (p.R101W; p.T179T) were found in the patient. In addition, a hemizygous mutation in complement factor properdin (CFP) (p.L456V) located on the X-chromosome was detected, inherited from the patient's mother. Protein structural modeling suggested impaired properdin subunit interactions by p.L456V that may hamper protein oligomerization required for complement activation. CONCLUSION: This study identified compound heterozygous mutations in IL10RA combined with a hemizygous CFP mutation in infantile-onset IBD by using whole-exome sequencing. CFP p.L456V may exacerbate symptoms of infantile-onset IBD by disturbing oligomerization of properdin.

Comparative analysis of protein structure alignments.

BACKGROUND: Several methods are currently available for the comparison of protein structures. These methods have been analysed regarding the performance in the identification of structurally/evolutionary related proteins, but so far there has been less focus on the objective comparison between the alignments produced by different methods. RESULTS: We analysed and compared the structural alignments obtained by different methods using three sets of pairs of structurally related proteins. The first set corresponds to 355 pairs of remote homologous proteins according to the SCOP database (ASTRAL40 set). The second set was derived from the SISYPHUS database and includes 69 protein pairs (SISY set). The third set consists of 40 pairs that are challenging to align (RIPC set). The alignment of pairs of this set requires indels of considerable number and size and some of the proteins are related by circular permutations, show extensive conformational variability or include repetitions. Two standard methods (CE and DALI) were applied to align the proteins in the ASTRAL40 set. The extent of structural similarity identified by both methods is highly correlated and the alignments from the two methods agree on average in more than half of the aligned positions. CE, DALI, as well as four additional methods (FATCAT, MATRAS, Calpha-match and SHEBA) were then compared using the SISY and RIPC sets. The accuracy of the alignments was assessed by comparison to reference alignments. The alignments generated by the different methods on average match more than half of the reference alignments in the SISY set. The alignments obtained in the more challenging RIPC set tend to differ considerably and match reference alignments less successfully than the SISY set alignments. CONCLUSION: The alignments produced by different methods tend to agree to a considerable extent, but the agreement is lower for the more challenging pairs. The results for the comparison to reference alignments are encouraging, but also indicate that there is still room for improvement.

Targeted Resequencing and Functional Testing Identifies Low-Frequency Missense Variants in the Gene Encoding GARP as Significant Contributors to Atopic Dermatitis Risk.

Gene-mapping studies have consistently identified a susceptibility locus for atopic dermatitis and other inflammatory diseases on chromosome band 11q13.5, with the strongest association observed for a common variant located in an intergenic region between the two annotated genes C11orf30 and LRRC32. Using a targeted resequencing approach we identified low-frequency and rare missense mutations within the LRRC32 gene encoding the protein GARP, a receptor on activated regulatory T cells that binds latent transforming growth factor-ß. Subsequent association testing in more than 2,000 atopic dermatitis patients and 2,000 control subjects showed a significant excess of these LRRC32 variants in individuals with atopic dermatitis. Structural protein modeling and bioinformatic analysis predicted a disruption of protein transport upon these variants, and overexpression assays in CD4+CD25- T cells showed a significant reduction in surface expression of the mutated protein. Consistently, flow cytometric (FACS) analyses of different T-cell subtypes obtained from atopic dermatitis patients showed a significantly reduced surface expression of GARP and a reduced conversion of CD4+CD25- T cells into regulatory T cells, along with lower expression of latency-associated protein upon stimulation in carriers of the LRRC32 A407T variant. These results link inherited disturbances of transforming growth factor-ß signaling with atopic dermatitis risk.

Autophagy receptor CALCOCO2/NDP52 takes center stage in Crohn disease.

To advance understanding of the complex genetics of Crohn disease (CD) we sequenced 42 whole exomes of patients with CD and five healthy control individuals, resulting in identification of a missense mutation in the autophagy receptor calcium binding and coiled-coil domain 2 (CALCOCO2/NDP52) gene. Protein domain modeling and functional studies highlight the potential role of this mutation in controlling NFKB signaling downstream of toll-like receptor (TLR) pathways. We summarize our recent findings and discuss the role of autophagy as a major modulator of proinflammatory signaling in the context of chronic inflammation.

Dense genotyping of immune-related disease regions identifies nine new risk loci for primary sclerosing cholangitis.

Primary sclerosing cholangitis (PSC) is a severe liver disease of unknown etiology leading to fibrotic destruction of the bile ducts and ultimately to the need for liver transplantation. We compared 3,789 PSC cases of European ancestry to 25,079 population controls across 130,422 SNPs genotyped using the Immunochip. We identified 12 genome-wide significant associations outside the human leukocyte antigen (HLA) complex, 9 of which were new, increasing the number of known PSC risk loci to 16. Despite comorbidity with inflammatory bowel disease (IBD) in 72% of the cases, 6 of the 12 loci showed significantly stronger association with PSC than with IBD, suggesting overlapping yet distinct genetic architectures for these two diseases. We incorporated association statistics from 7 diseases clinically occurring with PSC in the analysis and found suggestive evidence for 33 additional pleiotropic PSC risk loci. Together with network analyses, these findings add to the genetic risk map of PSC and expand on the relationship between PSC and other immune-mediated diseases.

Structure collisions between interacting proteins.

Protein-protein interactions take place at defined binding interfaces. One protein may bind two or more proteins at different interfaces at the same time. So far it has been commonly accepted that non-overlapping interfaces allow a given protein to bind other proteins simultaneously while no collisions occur between the binding protein structures. To test this assumption, we performed a comprehensive analysis of structural protein interactions to detect potential collisions. Our results did not indicate cases of biologically relevant collisions in the Protein Data Bank of protein structures. However, we discovered a number of collisions that originate from alternative protein conformations or quaternary structures due to different experimental conditions.

Genome-wide association study identifies a psoriasis susceptibility locus at TRAF3IP2.

Psoriasis is a multifactorial skin disease characterized by epidermal hyperproliferation and chronic inflammation, the most common form of which is psoriasis vulgaris (PsV). We present a genome-wide association analysis of 2,339,118 SNPs in 472 PsV cases and 1,146 controls from Germany, with follow-up of the 147 most significant SNPs in 2,746 PsV cases and 4,140 controls from three independent replication panels. We identified an association at TRAF3IP2 on 6q21 and genotyped two SNPs at this locus in two additional replication panels (the combined discovery and replication panels consisted of 6,487 cases and 8,037 controls; combined P = 2.36 × 10⁻¹° for rs13210247 and combined P = 1.24 × 10⁻¹6 for rs33980500). About 15% of psoriasis cases develop psoriatic arthritis (PsA). A stratified analysis of our datasets including only PsA cases (1,922 cases compared to 8,037 controls, P = 4.57 × 10⁻¹² for rs33980500) suggested that TRAF3IP2 represents a shared susceptibility for PsV and PsA. TRAF3IP2 encodes a protein involved in IL-17 signaling and which interacts with members of the Rel/NF-κB transcription factor family.

Genome-wide association study for ulcerative colitis identifies risk loci at 7q22 and 22q13 (IL17REL).

We performed a genome-wide association analysis of 1,897,764 SNPs in 1,043 German ulcerative colitis (UC) cases and 1,703 controls. We discovered new associations at chromosome 7q22 (rs7809799) and at chromosome 22q13 in IL17REL (rs5771069) and confirmed these associations in six replication panels (2,539 UC cases and 5,428 controls) from different regions of Europe (overall study sample P(rs7809799) = 8.81 x 10(-11) and P(rs5771069) = 4.21 x 10(-8), respectively).