Rare,Low-Frequency,and Common Variants in the Protein-Coding Sequence of Biological Candidate Genes from GWASs Contribute to Risk of Rheumatoid Arthritis

The extent to which variants in the protein-coding sequence of genes contribute to risk of rheumatoid arthritis (RA) is unknown. In this study, we addressed this issue by deep exon sequencing and large-scale genotyping of 25 biological candidate genes located within RA risk loci discovered by genome-wide association studies (GWASs). First, we assessed the contribution of rare coding variants in the 25 genes to the risk of RA in a pooled sequencing study of 500 RA cases and 650 controls of European ancestry. We observed an accumulation of rare nonsynonymous variants exclusive to RA cases in IL2RA and IL2RB (burden test: p = 0.007 and p = 0.018, respectively). Next, we assessed the aggregate contribution of low-frequency and common coding variants to the risk of RA by dense genotyping of the 25 gene loci in 10,609 RA cases and 35,605 controls. We observed a strong enrichment of coding variants with a nominal signal of association with RA (p < 0.05) after adjusting for the best signal of association at the loci (p_enrichment = 6.4 × 10⁻⁴). For one locus containing CD2, we found that a missense variant, rs699738 (c.798C>A [p.His266Gln]), and a noncoding variant, rs624988, reside on distinct haplotypes and independently contribute to the risk of RA (p = 4.6 × 10⁻⁶). Overall, our results indicate that variants (distributed across the allele-frequency spectrum) within the protein-coding portion of a subset of biological candidate genes identified by GWASs contribute to the risk of RA. Further, we have demonstrated that very large sample sizes will be required for comprehensively identifying the independent alleles contributing to the missing heritability of RA.     

Deep Resequencing of GWAS Loci Identifies Rare Variants in CARD9, IL23R and RNF186 That Are Associated with Ulcerative Colitis

Genome-wide association studies and follow-up meta-analyses in Crohn''s disease (CD) and ulcerative colitis (UC) have recently identified 163 disease-associated loci that meet genome-wide significance for these two inflammatory bowel diseases (IBD). These discoveries have already had a tremendous impact on our understanding of the genetic architecture of these diseases and have directed functional studies that have revealed some of the biological functions that are important to IBD (e.g. autophagy). Nonetheless, these loci can only explain a small proportion of disease variance (∼14% in CD and 7.5% in UC), suggesting that not only are additional loci to be found but that the known loci may contain high effect rare risk variants that have gone undetected by GWAS. To test this, we have used a targeted sequencing approach in 200 UC cases and 150 healthy controls (HC), all of French Canadian descent, to study 55 genes in regions associated with UC. We performed follow-up genotyping of 42 rare non-synonymous variants in independent case-control cohorts (totaling 14,435 UC cases and 20,204 HC). Our results confirmed significant association to rare non-synonymous coding variants in both IL23R and CARD9, previously identified from sequencing of CD loci, as well as identified a novel association in RNF186. With the exception of CARD9 (OR = 0.39), the rare non-synonymous variants identified were of moderate effect (OR = 1.49 for RNF186 and OR = 0.79 for IL23R). RNF186 encodes a protein with a RING domain having predicted E3 ubiquitin-protein ligase activity and two transmembrane domains. Importantly, the disease-coding variant is located in the ubiquitin ligase domain. Finally, our results suggest that rare variants in genes identified by genome-wide association in UC are unlikely to contribute significantly to the overall variance for the disease. Rather, these are expected to help focus functional studies of the corresponding disease loci.     

Investigation of Rare Single-Nucleotide PCDH15 Variants in Schizophrenia and Autism Spectrum Disorders

Both schizophrenia (SCZ) and autism spectrum disorders (ASD) are neuropsychiatric disorders with overlapping genetic etiology. Protocadherin 15 (PCDH15), which encodes a member of the cadherin super family that contributes to neural development and function, has been cited as a risk gene for neuropsychiatric disorders. Recently, rare variants of large effect have been paid attention to understand the etiopathology of these complex disorders. Thus, we evaluated the impacts of rare, single-nucleotide variants (SNVs) in PCDH15 on SCZ or ASD. First, we conducted coding exon-targeted resequencing of PCDH15 with next-generation sequencing technology in 562 Japanese patients (370 SCZ and 192 ASD) and detected 16 heterozygous SNVs. We then performed association analyses on 2,096 cases (1,714 SCZ and 382 ASD) and 1,917 controls with six novel variants of these 16 SNVs. Of these six variants, four (p.R219K, p.T281A, p.D642N, c.3010-1G>C) were ultra-rare variants (minor allele frequency < 0.0005) that may increase disease susceptibility. Finally, no statistically significant association between any of these rare, heterozygous PCDH15 point variants and SCZ or ASD was found. Our results suggest that a larger sample size of resequencing subjects is necessary to detect associations between rare PCDH15 variants and neuropsychiatric disorders.     

Exome Sequencing of 75 Individuals from Multiply Affected Coeliac Families and Large Scale Resequencing Follow Up

Coeliac disease (CeD) is a highly heritable common autoimmune disease involving chronic small intestinal inflammation in response to dietary wheat. The human leukocyte antigen (HLA) region, and 40 newer regions identified by genome wide association studies (GWAS) and dense fine mapping, account for ∼40% of the disease heritability. We hypothesized that in pedigrees with multiple individuals with CeD rare [minor allele frequency (MAF) <0.5%] mutations of larger effect size (odds ratios of ∼ 2–5) might exist. We sequenced the exomes of 75 coeliac individuals of European ancestry from 55 multiply affected families. We selected interesting variants and genes for further follow up using a combination of: an assessment of shared variants between related subjects, a model-free linkage test, and gene burden tests for multiple, potentially causal, variants. We next performed highly multiplexed amplicon resequencing of all RefSeq exons from 24 candidate genes selected on the basis of the exome sequencing data in 2,248 unrelated coeliac cases and 2,230 controls. 1,335 variants with a 99.9% genotyping call rate were observed in 4,478 samples, of which 939 were present in coding regions of 24 genes (Ti/Tv 2.99). 91.7% of coding variants were rare (MAF <0.5%) and 60% were novel. Gene burden tests performed on rare functional variants identified no significant associations (p<1×10⁻³) in the resequenced candidate genes. Our strategy of sequencing multiply affected families with deep follow up of candidate genes has not identified any new CeD risk mutations.     

Association and Mutation Analyses of 16p11.2 Autism Candidate Genes

Background

Autism is a complex childhood neurodevelopmental disorder with a strong genetic basis. Microdeletion or duplication of a ∼500–700-kb genomic rearrangement on 16p11.2 that contains 24 genes represents the second most frequent chromosomal disorder associated with autism. The role of common and rare 16p11.2 sequence variants in autism etiology is unknown.

Methodology/Principal Findings

To identify common 16p11.2 variants with a potential role in autism, we performed association studies using existing data generated from three microarray platforms: Affymetrix 5.0 (777 families), Illumina 550 K (943 families), and Affymetrix 500 K (60 families). No common variants were identified that were significantly associated with autism. To look for rare variants, we performed resequencing of coding and promoter regions for eight candidate genes selected based on their known expression patterns and functions. In total, we identified 26 novel variants in autism: 13 exonic (nine non-synonymous, three synonymous, and one untranslated region) and 13 promoter variants. We found a significant association between autism and a coding variant in the seizure-related gene SEZ6L2 (12/1106 autism vs. 3/1161 controls; p = 0.018). Sez6l2 expression in mouse embryos was restricted to the spinal cord and brain. SEZ6L2 expression in human fetal brain was highest in post-mitotic cortical layers, hippocampus, amygdala, and thalamus. Association analysis of SEZ6L2 in an independent sample set failed to replicate our initial findings.

Conclusions/Significance

We have identified sequence variation in at least one candidate gene in 16p11.2 that may represent a novel genetic risk factor for autism. However, further studies are required to substantiate these preliminary findings.     

Rare and common variants in ROM1 and PRPH2 genes trans-modify Stargardt/ABCA4 disease

Over 1,500 variants in the ABCA4 locus cause phenotypes ranging from severe, early-onset retinal degeneration to very late-onset maculopathies. The resulting ABCA4/Stargardt disease is the most prevalent Mendelian eye disorder, although its underlying clinical heterogeneity, including penetrance of many alleles, are not well-understood. We hypothesized that a share of this complexity is explained by trans-modifiers, i.e., variants in unlinked loci, which are currently unknown. We sought to identify these by performing exome sequencing in a large cohort for a rare disease of 622 cases and compared variation in seven genes known to clinically phenocopy ABCA4 disease to cohorts of ethnically matched controls. We identified a significant enrichment of variants in 2 out of the 7 genes. Moderately rare, likely functional, variants, at the minor allele frequency (MAF) <0.005 and CADD>25, were enriched in ROM1, where 1.3% of 622 patients harbored a ROM1 variant compared to 0.3% of 10,865 controls (p = 2.41E04; OR 3.81 95% CI [1.77; 8.22]). More importantly, analysis of common variants (MAF>0.1) identified a frequent haplotype in PRPH2, tagged by the p.Asp338 variant with MAF = 0.21 in the matched general population that was significantly increased in the patient cohort, MAF 0.25, p = 0.0014. Significant differences were also observed between ABCA4 disease subgroups. In the late-onset subgroup, defined by the hypomorphic p.Asn1868Ile variant and including c.4253+43G>A, the allele frequency for the PRPH2 p.Asp338 variant was 0.15 vs 0.27 in the remaining cohort, p = 0.00057. Known functional data allowed suggesting a mechanism by which the PRPH2 haplotype influences the ABCA4 disease penetrance. These associations were replicated in an independent cohort of 408 patients. The association was highly statistically significant in the combined cohorts of 1,030 cases, p = 4.00E-05 for all patients and p = 0.00014 for the hypomorph subgroup, suggesting a substantial trans-modifying role in ABCA4 disease for both rare and common variants in two unlinked loci.     

An Excess of Deleterious Variants in VEGF-A Pathway Genes in Down-Syndrome-Associated Atrioventricular Septal Defects

About half of people with trisomy 21 have a congenital heart defect (CHD), whereas the remainder have a structurally normal heart, demonstrating that trisomy 21 is a significant risk factor but is not causal for abnormal heart development. Atrioventricular septal defects (AVSD) are the most commonly occurring heart defects in Down syndrome (DS), and ∼65% of all AVSD is associated with DS. We used a candidate-gene approach among individuals with DS and complete AVSD (cases = 141) and DS with no CHD (controls = 141) to determine whether rare genetic variants in genes involved in atrioventricular valvuloseptal morphogenesis contribute to AVSD in this sensitized population. We found a significant excess (p < 0.0001) of variants predicted to be deleterious in cases compared to controls. At the most stringent level of filtering, we found potentially damaging variants in nearly 20% of cases but fewer than 3% of controls. The variants with the highest probability of being damaging in cases only were found in six genes: COL6A1, COL6A2, CRELD1, FBLN2, FRZB, and GATA5. Several of the case-specific variants were recurrent in unrelated individuals, occurring in 10% of cases studied. No variants with an equal probability of being damaging were found in controls, demonstrating a highly specific association with AVSD. Of note, all of these genes are in the VEGF-A pathway, even though the candidate genes analyzed in this study represented numerous biochemical and developmental pathways, suggesting that rare variants in the VEGF-A pathway might contribute to the genetic underpinnings of AVSD in humans.     

Low-Frequency Coding Variants at 6p21.33 and 20q11.21 Are Associated with Lung Cancer Risk in Chinese Populations

Genome-wide association studies have successfully identified a subset of common variants associated with lung cancer risk. However, these variants explain only a fraction of lung cancer heritability. It has been proposed that low-frequency or rare variants might have strong effects and contribute to the missing heritability. To assess the role of low-frequency or rare variants in lung cancer development, we analyzed exome chips representing 1,348 lung cancer subjects and 1,998 control subjects during the discovery stage and subsequently evaluated promising associations in an additional 4,699 affected subjects and 4,915 control subjects during the replication stages. Single-variant and gene-based analyses were carried out for coding variants with a minor allele frequency less than 0.05. We identified three low-frequency missense variants in BAT2 (rs9469031, c.1544C>T [p.Pro515Leu]; odds ratio [OR] = 0.55, p = 1.28 × 10⁻¹⁰), FKBPL (rs200847762, c.410C>T [p.Pro137Leu]; OR = 0.25, p = 9.79 × 10⁻¹²), and BPIFB1 (rs6141383, c.850G>A [p.Val284Met]; OR = 1.72, p = 1.79 × 10⁻⁷); these variants were associated with lung cancer risk. rs9469031 in BAT2 and rs6141383 in BPIFB1 were also associated with the age of onset of lung cancer (p = 0.001 and 0.006, respectively). BAT2 and FKBPL at 6p21.33 and BPIFB1 at 20q11.21 were differentially expressed in lung tumors and paired normal tissues. Gene-based analysis revealed that FKBPL, in which two independent variants were identified, might account for the association with lung cancer risk at 6p21.33. Our results highlight the important role low-frequency variants play in lung cancer susceptibility and indicate that candidate genes at 6p21.33 and 20q11.21 are potentially biologically relevant to lung carcinogenesis.     

Novel Crohn disease locus identified by genome-wide association maps to a gene desert on 5p13.1 and modulates expression of PTGER4

To identify novel susceptibility loci for Crohn disease (CD), we undertook a genome-wide association study with more than 300,000 SNPs characterized in 547 patients and 928 controls. We found three chromosome regions that provided evidence of disease association with p-values between 10⁻⁶ and 10⁻⁹. Two of these (IL23R on Chromosome 1 and CARD15 on Chromosome 16) correspond to genes previously reported to be associated with CD. In addition, a 250-kb region of Chromosome 5p13.1 was found to contain multiple markers with strongly suggestive evidence of disease association (including four markers with p < 10⁻⁷). We replicated the results for 5p13.1 by studying 1,266 additional CD patients, 559 additional controls, and 428 trios. Significant evidence of association (p < 4 × 10⁻⁴) was found in case/control comparisons with the replication data, while associated alleles were over-transmitted to affected offspring (p < 0.05), thus confirming that the 5p13.1 locus contributes to CD susceptibility. The CD-associated 250-kb region was saturated with 111 SNP markers. Haplotype analysis supports a complex locus architecture with multiple variants contributing to disease susceptibility. The novel 5p13.1 CD locus is contained within a 1.25-Mb gene desert. We present evidence that disease-associated alleles correlate with quantitative expression levels of the prostaglandin receptor EP4, PTGER4, the gene that resides closest to the associated region. Our results identify a major new susceptibility locus for CD, and suggest that genetic variants associated with disease risk at this locus could modulate cis-acting regulatory elements of PTGER4.     

Genome-wide association scan shows genetic variants in the FTO gene are associated with obesity-related traits

The obesity epidemic is responsible for a substantial economic burden in developed countries and is a major risk factor for type 2 diabetes and cardiovascular disease. The disease is the result not only of several environmental risk factors, but also of genetic predisposition. To take advantage of recent advances in gene-mapping technology, we executed a genome-wide association scan to identify genetic variants associated with obesity-related quantitative traits in the genetically isolated population of Sardinia. Initial analysis suggested that several SNPs in the FTO and PFKP genes were associated with increased BMI, hip circumference, and weight. Within the FTO gene, rs9930506 showed the strongest association with BMI (p = 8.6 ×10⁻⁷), hip circumference (p = 3.4 × 10⁻⁸), and weight (p = 9.1 × 10⁻⁷). In Sardinia, homozygotes for the rare “G” allele of this SNP (minor allele frequency = 0.46) were 1.3 BMI units heavier than homozygotes for the common “A” allele. Within the PFKP gene, rs6602024 showed very strong association with BMI (p = 4.9 × 10⁻⁶). Homozygotes for the rare “A” allele of this SNP (minor allele frequency = 0.12) were 1.8 BMI units heavier than homozygotes for the common “G” allele. To replicate our findings, we genotyped these two SNPs in the GenNet study. In European Americans (N = 1,496) and in Hispanic Americans (N = 839), we replicated significant association between rs9930506 in the FTO gene and BMI (p-value for meta-analysis of European American and Hispanic American follow-up samples, p = 0.001), weight (p = 0.001), and hip circumference (p = 0.0005). We did not replicate association between rs6602024 and obesity-related traits in the GenNet sample, although we found that in European Americans, Hispanic Americans, and African Americans, homozygotes for the rare “A” allele were, on average, 1.0–3.0 BMI units heavier than homozygotes for the more common “G” allele. In summary, we have completed a whole genome–association scan for three obesity-related quantitative traits and report that common genetic variants in the FTO gene are associated with substantial changes in BMI, hip circumference, and body weight. These changes could have a significant impact on the risk of obesity-related morbidity in the general population.     

Genetic Variation in the Familial Mediterranean Fever Gene (MEFV) and Risk for Crohn's Disease and Ulcerative Colitis

Background and Aims

The familial Mediterranean fever (FMF) gene (MEFV) encodes pyrin, a major regulator of the inflammasome platform controlling caspase-1 activation and IL-1β processing. Pyrin has been shown to interact with the gene product of NLRP3, NALP3/cryopyrin, also an important active member of the inflammasome. The NLRP3 region was recently reported to be associated with Crohn''s disease (CD) susceptibility. We therefore sought to evaluate MEFV as an inflammatory bowel disease (IBD) susceptibility gene.

Methodology and Results

MEFV colonic mucosal gene expression was significantly increased in experimental colitis mice models (TNBS p<0.0003; DSS p<0.006), in biopsies from CD (p<0.02) and severe ulcerative colitis (UC) patients (p<0.008). Comprehensive genetic screening of the MEFV region in the Belgian exploratory sample set (440 CD trios, 137 UC trios, 239 CD cases, 96 UC cases, and 107 healthy controls) identified SNPs located in the MEFV 5′ haplotype block that were significantly associated with UC (rs224217; p = 0.003; A allele frequency: 56% cases, 45% controls), while no CD associations were observed. Sequencing and subsequent genotyping of variants located in this associated haplotype block identified three synonymous variants (D102D/rs224225, G138G/rs224224, A165A/rs224223) and one non-synonymous variant (R202Q/rs224222) located in MEFV exon 2 that were significantly associated with UC (rs224222: p = 0.0005; A allele frequency: 32% in cases, 23% in controls). No consistent associations were observed in additional Canadian (256 CD trios, 91 UC trios) and Scottish (495 UC, 370 controls) sample sets. We note that rs224222 showed marginal association (p = 0.012; G allele frequency: 82% in cases, 70% in controls) in the Canadian sample, but with a different risk allele. None of the NLRP3 common variants were associated with UC in the Belgian-Canadian UC samples and no significant interactions were observed between NLRP3 and MEFV that could explain the observed flip-flop of the rs224222 risk allele.

Conclusion

The differences in association levels observed between the sample sets may be a consequence of distinct founder effects or of the relative small sample size of the cohorts evaluated in this study. However, the results suggest that common variants in the MEFV region do not contribute to CD and UC susceptibility.     

In-Silico Analysis of Inflammatory Bowel Disease (IBD) GWAS Loci to Novel Connections

Genome-wide association studies (GWASs) for many complex diseases, including inflammatory bowel disease (IBD), produced hundreds of disease-associated loci—the majority of which are noncoding. The number of GWAS loci is increasing very rapidly, but the process of translating single nucleotide polymorphisms (SNPs) from these loci to genomic medicine is lagging. In this study, we investigated 4,734 variants from 152 IBD associated GWAS loci (IBD associated 152 lead noncoding SNPs identified from pooled GWAS results + 4,582 variants in strong linkage-disequilibrium (LD) (r² ≥0.8) for EUR population of 1K Genomes Project) using four publicly available bioinformatics tools, e.g. dbPSHP, CADD, GWAVA, and RegulomeDB, to annotate and prioritize putative regulatory variants. Of the 152 lead noncoding SNPs, around 11% are under strong negative selection (GERP++ RS ≥2); and ~30% are under balancing selection (Tajima’s D score >2) in CEU population (1K Genomes Project)—though these regions are positively selected (GERP++ RS <0) in mammalian evolution. The analysis of 4,734 variants using three integrative annotation tools produced 929 putative functional SNPs, of which 18 SNPs (from 15 GWAS loci) are in concordance with all three classifiers. These prioritized noncoding SNPs may contribute to IBD pathogenesis by dysregulating the expression of nearby genes. This study showed the usefulness of integrative annotation for prioritizing fewer functional variants from a large number of GWAS markers.     

Identification of Environmental Factors Associated with Inflammatory Bowel Disease in a Southwestern Highland Region of China: A Nested Case-Control Study

BackgroundThe aim of this study was to examine environmental factors associated with inflammatory bowel disease (IBD) in Yunnan Province, a southwestern highland region of China.MethodsIn this nested case-control study, newly diagnosed ulcerative colitis (UC) cases in 2 cities in Yunnan Province and Crohn’s disease (CD) cases in 16 cities in Yunnan Province were recruited between 2008 and 2013. Controls were matched by geography, sex and age at a ratio of 1:4. Data were collected using the designed questionnaire. Conditional logistic regression models were used to estimate adjusted odds ratios (ORs).ResultsA total of 678 UC and 102 CD cases were recruited. For UC, various factors were associated with an increased risk of developing UC: dietary habits, including frequent irregular meal times; consumption of fried foods, salty foods and frozen dinners; childhood factors, including intestinal infectious diseases and frequent use of antibiotics; and other factors, such as mental labor, high work stress, use of non-aspirin non-steroidal anti-inflammatory drugs and allergies (OR > 1, p < 0.05). Other factors showed a protective effect: such as consumption of fruits, current smoking, physical activity, and drinking tea (OR < 1, p < 0.05). For CD, appendectomy and irregular meal times increased the disease risk (OR >1, p < 0.05), whereas physical activity may have reduced this risk (OR < 1, p < 0.05).ConclusionsThis study is the first nested case-control study to analyze the association between environmental factors and IBD onset in a southwestern highland region of China. Certain dietary habits, lifestyles, allergies and childhood factors may play important roles in IBD, particularly UC.     

Genetic Variants in Genes of the Inflammatory Response in Association with Infective Endocarditis

Aims

Inflammation in infective endocarditis (IE) is a complex network including interactions of inflammatory cytokines and other components of host response. Certainly, any variation in this network could influence susceptibility or disease progression of IE. In this study, 14 single nucleotide variants (SNVs) in genes coding for interleukin-1β, interleukin-6, interleukin-10, toll–like receptor-4, tumor necrosis factor-α, selectin E and intercellular adhesion molecule-1 were analyzed for an association with susceptibility to IE and correlated with disease-related laboratory parameters. Furthermore, the occurrence of SNVs was examined to elucidate pathogen-dependent associations.

Methods and Results

The distribution of SNVs was determined in IE-patients and healthy blood donors by RFLP analysis. White blood cells (WBC) were counted using flow cytometry, concentration of C-reactive protein and procalcitonin was measured immunologically. Interleukin-6 c.471+870G>A genotypes differed significantly between IE patients and controls. The frequency of the heterozygote genotype GA was considerably higher in the patient group (68.9% vs. 43.8%, P_c<0.0003). Interleukin-6 c.-237 minor allele frequency was increased in patients, although not statistically significant. Additionally, we detected a potential relation between interleukin-1β c.315C>T and IE. Pathogen-dependent analysis showed no significantly associated subgroup in relation to IE susceptibility, but gave hints towards alterations regarding Enterococcus-caused IE cases. Patients with genotype selectin-E c.-19 GT tend to have higher preoperative WBC counts than patients with genotype GG. We further showed an association between two interleukin-1β SNVs and laboratory biomarkers.

Conclusion

This study shows genetic predispositions for the establishment of IE. Furthermore, correlation of SNVs with disease-related biomarkers suggests a role of genetic variants regarding the inflammatory response in IE.     

Genetic Association Analysis of the Functional c.714T>G Polymorphism and Mucosal Expression of Dectin-1 in Inflammatory Bowel Disease

Background

Dectin-1 is a pattern recognition receptor (PRR) expressed by myeloid cells that specifically recognizes β-1,3 glucan, a polysaccharide and major component of the fungal cell wall. Upon activation, dectin-1 signaling converges, similar to NOD2, on the adaptor molecule CARD9 which is associated with inflammatory bowel disease (IBD). An early stop codon polymorphism (c.714T>G) in DECTIN-1 results in a loss-of-function (p.Y238X) and impaired cytokine responses, including TNF-α, interleukin (IL)-1β and IL-17 upon in vitro stimulation with Candida albicans or β-glucan. The aim of the present study was to test the hypothesis that the DECTIN-1 c.714T>G (p.Y238X) polymorphism is associated with lower disease susceptibility or severity in IBD and to investigate the level of dectin-1 expression in inflamed and non-inflamed colon tissue of IBD patients.

Methodology

Paraffin embedded tissue samples from non-inflamed and inflamed colon of IBD patients and from diverticulitis patients were immunohistochemically stained for dectin-1 and related to CD68 macrophage staining. Genomic DNA of IBD patients (778 patients with Crohn''s disease and 759 patients with ulcerative colitis) and healthy controls (n = 772) was genotyped for the c.714T>G polymorphism and genotype-phenotype interactions were investigated.

Principal Findings

Increased expression of dectin-1 was observed in actively inflamed colon tissue, as compared to non-inflamed tissue of the same patients. Also an increase in dectin-1 expression was apparent in diverticulitis tissue. No statistically significant difference in DECTIN-1 c.714T>G allele frequencies was observed between IBD patients and healthy controls. Furthermore, no differences in clinical characteristics could be observed related to DECTIN-1 genotype, neither alone, nor stratified for NOD2 genotype.

Conclusions

Our data demonstrate that dectin-1 expression is elevated on macrophages, neutrophils, and other immune cells involved in the inflammatory reaction in IBD. The DECTIN-1 c.714T>G polymorphism however, is not a major susceptibility factor for developing IBD.     

Genome-Wide Analyses of Exonic Copy Number Variants in a Family-Based Study Point to Novel Autism Susceptibility Genes

The genetics underlying the autism spectrum disorders (ASDs) is complex and remains poorly understood. Previous work has demonstrated an important role for structural variation in a subset of cases, but has lacked the resolution necessary to move beyond detection of large regions of potential interest to identification of individual genes. To pinpoint genes likely to contribute to ASD etiology, we performed high density genotyping in 912 multiplex families from the Autism Genetics Resource Exchange (AGRE) collection and contrasted results to those obtained for 1,488 healthy controls. Through prioritization of exonic deletions (eDels), exonic duplications (eDups), and whole gene duplication events (gDups), we identified more than 150 loci harboring rare variants in multiple unrelated probands, but no controls. Importantly, 27 of these were confirmed on examination of an independent replication cohort comprised of 859 cases and an additional 1,051 controls. Rare variants at known loci, including exonic deletions at NRXN1 and whole gene duplications encompassing UBE3A and several other genes in the 15q11–q13 region, were observed in the course of these analyses. Strong support was likewise observed for previously unreported genes such as BZRAP1, an adaptor molecule known to regulate synaptic transmission, with eDels or eDups observed in twelve unrelated cases but no controls (p=2.3×10⁻⁵). Less is known about MDGA2, likewise observed to be case-specific (p=1.3×10⁻⁴). But, it is notable that the encoded protein shows an unexpectedly high similarity to Contactin 4 (BLAST E-value=3×10⁻³⁹), which has also been linked to disease. That hundreds of distinct rare variants were each seen only once further highlights complexity in the ASDs and points to the continued need for larger cohorts.     

Whole Exome Sequencing Reveals Homozygous Mutations in RAI1, OTOF,and SLC26A4 Genes Associated with Nonsyndromic Hearing Loss in Altaian Families (South Siberia)

Hearing loss (HL) is one of the most common sensorineural disorders and several dozen genes contribute to its pathogenesis. Establishing a genetic diagnosis of HL is of great importance for clinical evaluation of deaf patients and for estimating recurrence risks for their families. Efforts to identify genes responsible for HL have been challenged by high genetic heterogeneity and different ethnic-specific prevalence of inherited deafness. Here we present the utility of whole exome sequencing (WES) for identifying candidate causal variants for previously unexplained nonsyndromic HL of seven patients from four unrelated Altaian families (the Altai Republic, South Siberia). The WES analysis revealed homozygous missense mutations in three genes associated with HL. Mutation c.2168A>G (SLC26A4) was found in one family, a novel mutation c.1111G>C (OTOF) was revealed in another family, and mutation c.5254G>A (RAI1) was found in two families. Sanger sequencing was applied for screening of identified variants in an ethnically diverse cohort of other patients with HL (n = 116) and in Altaian controls (n = 120). Identified variants were found only in patients of Altaian ethnicity (n = 93). Several lines of evidences support the association of homozygosity for discovered variants c.5254G>A (RAI1), c.1111C>G (OTOF), and c.2168A>G (SLC26A4) with HL in Altaian patients. Local prevalence of identified variants implies possible founder effect in significant number of HL cases in indigenous population of the Altai region. Notably, this is the first reported instance of patients with RAI1 missense mutation whose HL is not accompanied by specific traits typical for Smith-Magenis syndrome. Presumed association of RAI1 gene variant c.5254G>A with isolated HL needs to be proved by further experimental studies.     

Low Frequency Variants in the Exons Only Encoding Isoform A of HNF1A Do Not Contribute to Susceptibility to Type 2 Diabetes

Background

There is considerable interest in the hypothesis that low frequency, intermediate penetrance variants contribute to the proportion of Type 2 Diabetes (T2D) susceptibility not attributable to the common variants uncovered through genome-wide association approaches. Genes previously implicated in monogenic and multifactorial forms of diabetes are obvious candidates in this respect. In this study, we focussed on exons 8–10 of the HNF1A gene since rare, penetrant mutations in these exons (which are only transcribed in selected HNF1A isoforms) are associated with a later age of diagnosis of Maturity onset diabetes of the young (MODY) than mutations in exons 1–7. The age of diagnosis in the subgroup of HNF1A-MODY individuals with exon 8–10 mutations overlaps with that of early multifactorial T2D, and we set out to test the hypothesis that these exons might also harbour low-frequency coding variants of intermediate penetrance that contribute to risk of multifactorial T2D.

Methodology and Principal Findings

We performed targeted capillary resequencing of HNF1A exons 8–10 in 591 European T2D subjects enriched for genetic aetiology on the basis of an early age of diagnosis (≤45 years) and/or family history of T2D (≥1 affected sibling). PCR products were sequenced and compared to the published HNF1A sequence. We identified several variants (rs735396 [IVS9−24T>C], rs1169304 [IVS8+29T>C], c.1768+44C>T [IVS9+44C>T] and rs61953349 [c.1545G>A, p.T515T] but no novel non-synonymous coding variants were detected.

Conclusions and Significance

We conclude that low frequency, nonsynonymous coding variants in the terminal exons of HNF1A are unlikely to contribute to T2D-susceptibility in European samples. Nevertheless, the rationale for seeking low-frequency causal variants in genes known to contain rare, penetrant mutations remains strong and should motivate efforts to screen other genes in a similar fashion.