Identification of intergenic trans-regulatory RNAs containing a disease-linked SNP sequence and targeting cell cycle progression/differentiation pathways in multiple common human disorders

Meta-analysis of genomic coordinates of SNP variations identified in genome-wide association studies (GWAS) of up to 712,253 samples (comprising 221,158 disease cases, 322,862 controls, and 168,233 case/control subjects of obesity GWAS) reveals that 39% of SNPs associated with 22 common human disorders are located within intergenic regions. Chromatin-state maps based on H3K4me3-H3K36me3 signatures show that many intergenic disease-linked SNPs are located within the boundaries of the K4-K36 domains, suggesting that SNP-harboring genomic regions are transcribed. Here we report identification of 13 trans-regulatory RNAs (transRNAs) 100 to 200 nucleotides in length containing intergenic SNP sequences associated with Crohn's disease, rheumatoid arthritis, type 1 diabetes, vitiligo, and multiple types of epithelial malignancies (prostate, breast, ovarian, and colorectal cancers). We demonstrate that NALP1 loci intergenic SNP sequence, rs2670660, is expressed in human cells and may contribute to clinical manifestations of autoimmune and autoimflammatory phenotypes by generating distinct allelic variants of transRNAs. Stable expression of allele-specific sense and anti-sense variants of transRNAs markedly alters cellular behavior, affect cell cycle progression, and interfere with monocyte/macrophage transdifferentiation. On a molecular level, forced expression of allele-specific sense and anti-sense variants of transRNAs asserts allele-specific genome-wide effects on abundance of hundreds microRNAs and mRNAs. Using lentiviral gene transfer, microarray and Q-RT-PCR technologies, we identify rs2670660 allele-specific gene expression signatures (GES) which appear useful for detecting the activated states of innate immunity/inflammasome pathways in ~700 clinical samples from 185 control subjects and 350 patients diagnosed with 9 common human disorders, including Crohn's disease, ulcerative colitis, rheumatoid arthritis, Huntington's disease, autism, Alzheimer's disease, obesity, prostate and breast cancers. Microarray analysis of clinical samples demonstrates that rs2670660 allele-specific GES are engaged in patients' peripheral blood mononuclear cells (PBMC) which encounter pathological conditions in coherent tissues of a human body during immune surveillance and homeostasis monitoring. These data indicate that expression of transRNAs encoded by specific intergenic sequences can trigger activation of innate immunity/inflammasome pathways and contribute to clinical development of autoinflammatory and autoimmune syndromes. Documented in this work single-base substitution-driven molecular and biological antagonisms of intergenic SNP-containing transRNAs suggest a guiding mechanism of selection and retention of phenotype-compatible intergenic variations during evolution. According to this model, random genetic variations which generate transRNAs asserting antagonistic phenotype-altering effects compared to ancestral alleles will be selected and retained as SNP variants.     

WW-domain-containing oxidoreductase is associated with low plasma HDL-C levels

Low serum HDL-cholesterol (HDL-C) is a major risk factor for coronary artery disease. We performed targeted genotyping of a 12.4 Mb linked region on 16q to test for association with low HDL-C by using a regional-tag SNP strategy. We identified one SNP, rs2548861, in the WW-domain-containing oxidoreductase (WWOX) gene with region-wide significance for low HDL-C in dyslipidemic families of Mexican and European descent and in low-HDL-C cases and controls of European descent (p = 6.9 × 10⁻⁷). We extended our investigation to the population level by using two independent unascertained population-based Finnish cohorts, the cross-sectional METSIM cohort of 4,463 males and the prospective Young Finns cohort of 2,265 subjects. The combined analysis provided p = 4 × 10⁻⁴ to 2 × 10⁻⁵. Importantly, in the prospective cohort, we observed a significant longitudinal association of rs2548861 with HDL-C levels obtained at four different time points over 21 years (p = 0.003), and the T risk allele explained 1.5% of the variance in HDL-C levels. The rs2548861 resides in a highly conserved region in intron 8 of WWOX. Results from our in vitro reporter assay and electrophoretic mobility-shift assay demonstrate that this region functions as a cis-regulatory element whose associated rs2548861 SNP has a specific allelic effect and that the region forms an allele-specific DNA-nuclear-factor complex. In conclusion, analyses of 9,798 subjects show significant association between HDL-C and a WWOX variant with an allele-specific cis-regulatory function.     

Multiple Functional Risk Variants in a SMAD7 Enhancer Implicate a Colorectal Cancer Risk Haplotype

Genome-wide association studies (GWAS) of colorectal cancer (CRC) have led to the identification of a number of common variants associated with modest risk. Several risk variants map within the vicinity of TGFβ/BMP signaling pathway genes, including rs4939827 within an intron of SMAD7 at 18q21.1. A previous study implicated a novel SNP (novel 1 or rs58920878) as a functional variant within an enhancer element in SMAD7 intron 4. In this study, we show that four SNPs including novel 1 (rs6507874, rs6507875, rs8085824, and rs58920878) in linkage disequilibrium (LD) with the index SNP rs4939827 demonstrate allele-specific enhancer effects in a large, multi-component enhancer of SMAD7. All four SNPs demonstrate allele-specific protein binding to nuclear extracts of CRC cell lines. Furthermore, some of the risk-associated alleles correlate with increased expression of SMAD7 in normal colon tissues. Finally, we show that the enhancer is responsive to BMP4 stimulation. Taken together, we propose that the associated CRC risk at 18q21.1 is due to four functional variants that regulate SMAD7 expression and potentially perturb a BMP negative feedback loop in TGFβ/BMP signaling pathways.     

Genetic Predisposition to In Situ and Invasive Lobular Carcinoma of the Breast

Invasive lobular breast cancer (ILC) accounts for 10–15% of all invasive breast carcinomas. It is generally ER positive (ER+) and often associated with lobular carcinoma in situ (LCIS). Genome-wide association studies have identified more than 70 common polymorphisms that predispose to breast cancer, but these studies included predominantly ductal (IDC) carcinomas. To identify novel common polymorphisms that predispose to ILC and LCIS, we pooled data from 6,023 cases (5,622 ILC, 401 pure LCIS) and 34,271 controls from 36 studies genotyped using the iCOGS chip. Six novel SNPs most strongly associated with ILC/LCIS in the pooled analysis were genotyped in a further 516 lobular cases (482 ILC, 36 LCIS) and 1,467 controls. These analyses identified a lobular-specific SNP at 7q34 (rs11977670, OR (95%CI) for ILC = 1.13 (1.09–1.18), P = 6.0×10⁻¹⁰; P-het for ILC vs IDC ER+ tumors = 1.8×10⁻⁴). Of the 75 known breast cancer polymorphisms that were genotyped, 56 were associated with ILC and 15 with LCIS at P<0.05. Two SNPs showed significantly stronger associations for ILC than LCIS (rs2981579/10q26/FGFR2, P-het = 0.04 and rs889312/5q11/MAP3K1, P-het = 0.03); and two showed stronger associations for LCIS than ILC (rs6678914/1q32/LGR6, P-het = 0.001 and rs1752911/6q14, P-het = 0.04). In addition, seven of the 75 known loci showed significant differences between ER+ tumors with IDC and ILC histology, three of these showing stronger associations for ILC (rs11249433/1p11, rs2981579/10q26/FGFR2 and rs10995190/10q21/ZNF365) and four associated only with IDC (5p12/rs10941679; rs2588809/14q24/RAD51L1, rs6472903/8q21 and rs1550623/2q31/CDCA7). In conclusion, we have identified one novel lobular breast cancer specific predisposition polymorphism at 7q34, and shown for the first time that common breast cancer polymorphisms predispose to LCIS. We have shown that many of the ER+ breast cancer predisposition loci also predispose to ILC, although there is some heterogeneity between ER+ lobular and ER+ IDC tumors. These data provide evidence for overlapping, but distinct etiological pathways within ER+ breast cancer between morphological subtypes.     

A Retrospective Observational Study of the Relationship between Single Nucleotide Polymorphisms Associated with the Risk of Developing Colorectal Cancer and Survival

BackgroundThere is variability in clinical outcome for patients with apparently the same stage colorectal cancer (CRC). Single nucleotide polymorphisms (SNPs) mapping to chromosomes 1q41, 3q26.2, 6p21, 8q23.3, 8q24.21, 10p14, 11q13, 11q23.1, 12q13.13, 14q22, 14q22.2, 15q13.3, 16q22.1, 18q21.1, 19q13.11, 20p12, 20p12.3, 20q13.33 and Xp22 have robustly been shown to be associated with the risk of developing CRC. Since germline variation can also influence patient outcome the relationship between these SNPs and patient survivorship from CRC was examined.MethodsAll enrolled into the National Study of Colorectal Cancer Genetics (NSCCG) were genotyped for 1q41, 3q26.2, 6p21, 8q23.3, 8q24.21, 10p14, 11q13, 11q23.1, 12q13.13, 14q22, 14q22.2, 15q13.3, 16q22.1, 18q21.1, 19q13.11, 20p12, 20p12.3, 20q13.33 and xp22 SNPs. Linking this information to the National Cancer Data Repository allowed patient genotype to be related to survival.ResultsThe linked dataset consisted of 4,327 individuals. 14q22.22 genotype defined by the SNP rs4444235 showed a significant association with overall survival. Specifically, the C allele was associated with poorer observed survival (per allele hazard ratio 1.13, 95% confidence interval 1.05–1.22, P = 0.0015).ConclusionThe CRC susceptibility SNP rs4444235 also appears to exert an influence in modulating patient survival and warrants further evaluation as a potential prognostic marker.     

Evaluating Genome-Wide Association Study-Identified Breast Cancer Risk Variants in African-American Women

Genome-wide association studies (GWAS), conducted mostly in European or Asian descendants, have identified approximately 67 genetic susceptibility loci for breast cancer. Given the large differences in genetic architecture between the African-ancestry genome and genomes of Asians and Europeans, it is important to investigate these loci in African-ancestry populations. We evaluated index SNPs in all 67 breast cancer susceptibility loci identified to date in our study including up to 3,300 African-American women (1,231 cases and 2,069 controls), recruited in the Southern Community Cohort Study (SCCS) and the Nashville Breast Health Study (NBHS). Seven SNPs were statistically significant (P≤0.05) with the risk of overall breast cancer in the same direction as previously reported: rs10069690 (5p15/TERT), rs999737 (14q24/RAD51L1), rs13387042 (2q35/TNP1), rs1219648 (10q26/FGFR2), rs8170 (19p13/BABAM1), rs17817449 (16q12/FTO), and rs13329835 (16q23/DYL2). A marginally significant association (P<0.10) was found for three additional SNPs: rs1045485 (2q33/CASP8), rs4849887 (2q14/INHBB), and rs4808801 (19p13/ELL). Three additional SNPs, including rs1011970 (9p21/CDKN2A/2B), rs941764 (14q32/CCDC88C), and rs17529111 (6q14/FAM46A), showed a significant association in analyses conducted by breast cancer subtype. The risk of breast cancer was elevated with an increasing number of risk variants, as measured by quintile of the genetic risk score, from 1.00 (reference), to 1.75 (1.30–2.37), 1.56 (1.15–2.11), 2.02 (1.50–2.74) and 2.63 (1.96–3.52), respectively, (P = 7.8×10^–10). Results from this study highlight the need for large genetic studies in AAs to identify risk variants impacting this population.     

The 12p13.33/RAD52 Locus and Genetic Susceptibility to Squamous Cell Cancers of Upper Aerodigestive Tract

Genetic variants located within the 12p13.33/RAD52 locus have been associated with lung squamous cell carcinoma (LUSC). Here, within 5,947 UADT cancers and 7,789 controls from 9 different studies, we found rs10849605, a common intronic variant in RAD52, to be also associated with upper aerodigestive tract (UADT) squamous cell carcinoma cases (OR = 1.09, 95% CI: 1.04–1.15, p = 6x10⁻⁴). We additionally identified rs10849605 as a RAD52 cis-eQTL inUADT(p = 1x10⁻³) and LUSC (p = 9x10⁻⁴) tumours, with the UADT/LUSC risk allele correlated with increased RAD52 expression levels. The 12p13.33 locus, encompassing rs10849605/RAD52, was identified as a significant somatic focal copy number amplification in UADT(n = 374, q-value = 0.075) and LUSC (n = 464, q-value = 0.007) tumors and correlated with higher RAD52 tumor expression levels (p = 6x10⁻⁴⁸ and p = 3x10⁻²⁹ in UADT and LUSC, respectively). In combination, these results implicate increased RAD52 expression in both genetic susceptibility and tumorigenesis of UADT and LUSC tumors.     

Polymorphisms on 8q24 are associated with lung cancer risk and survival in han chinese

Chromosome 8q24 is commonly amplified in many types of cancer, particularly lung cancer. Polymorphisms in this region are associated with risk of different cancers. To investigate the relationship between three single nucleotide polymorphisms (SNPs) (rs1447295, rs16901979 and rs6983267) on 8q24 and lung cancer risk, we conducted an association study in two Han Chinese populations: one population was from Zhejiang Province (576 case patients and 576 control subjects), whereas the other was from Fujian Province (576 case patients and 576 control subjects). We found that rs6983267 was significantly associated with an increased risk of lung cancer in both populations. Compared with the TT genotype, the GG genotype was associated with a significant 1.555-fold increased risk of lung cancer [95% confidence interval (CI) 1.218–1.986, P = 4.0×10⁻⁴]. This effect was more pronounced in never-smokers [odds ratio (OR) = 2.366, 95% CI 1.605–3.488, P = 1.4×10⁻⁵]. Analyses stratified by histology revealed that rs6983267 GG genotype was most associated with patients with other histological types (OR = 3.012, 95% CI 1.675–5.417, P = 2.3×10⁻⁴). The AA genotype of rs1447295 was associated with increased risk for adenocarcinoma compared with the CC genotype (OR = 2.260, 95% CI 1.174–4.353, P = 0.015). Furthermore, the GG genotype of rs6983267 was associated with worse survival in the Zhejiang population (hazard ratio (HR) = 1.646, 95% CI 1.099–2.464, P = 0.016). No association was observed for rs16901979. These results suggest that genetic variations on 8q24 may play significant roles in the development and progression of lung cancer.     

Allele-specific hybridization using streptavidin-coated magnetic beads for species identification, S genotyping, and SNP analysis in plants

Dot-blot hybridization has been successfully used for the construction of single nucleotide polymorphism (SNP)-based linkage maps, quantitative trait locus analysis, marker-assisted selection, and the identification of species and cultivars. This method is, however, time-consuming, even for a small number of plant samples. We propose a method in which streptavidin-coated magnetic beads replace the nylon membrane for immobilization of the PCR products and are hybridized with allele-specific oligonucleotide probes and a digoxigenin-labeled oligonucleotide hybridized with the allele-specific oligonucleotide probe. After amplification of plant DNA by PCR with the biotinylated primers, those oligonucleotide probes having species-specific or allele-specific sequences were mixed together with the digoxigenin-labeled oligonucleotide and the streptavidin-coated magnetic beads at a temperature suitable for each probe. Species-specific internal transcribed spacer 1 (ITS1) sequences and allele-specific sequences of the hypervariable region I of S-locus receptor kinase (SRK) specifically detected ITS1 sequences and SRK alleles in Brassica species, respectively. SNPs were also successfully analyzed by using allele-specific oligonucleotide probes and competitive oligonucleotides. In the SNP analysis, PCR products were indirectly captured by magnetic beads. SNP alleles of eight cultivars each of Brassica rapa and Raphanus sativus were analyzed using streptavidin-coated magnetic beads. The genotyping results corresponded well with those of dot-blot-SNP analysis. Although allele-specific hybridization using streptavidin-coated magnetic beads is somewhat costly, it is easier and more rapid than dot-blot hybridization. This method is suitable for the analysis of a small number of plant samples with a large number of DNA markers.     

MiRNA-Related Genetic Variations Associated with Radiotherapy-Induced Toxicities in Patients with Locally Advanced Non–Small Cell Lung Cancer

Severe radiation-induced toxicities limit treatment efficacy and compromise outcomes of lung cancer. We aimed to identify microRNA-related genetic variations as biomarkers for the prediction of radiotherapy-induced acute toxicities. We genotyped 233 SNPs (161 in microRNA binding site and 72 in processing gene) and analyzed their associations with pneumonitis and esophagitis in 167 stage III NSCLC patients received definitive radiation therapy. Sixteen and 11 SNPs were associated with esophagitis and pneumonitis, respectively. After multiple comparison correction, RPS6KB2:rs10274, SMO:rs1061280, SMO:rs1061285 remained significantly associated with esophagitis, while processing gene DGCR8:rs720014, DGCR8:rs3757, DGCR8:rs1633445 remained significantly associated with pneumonitis. Patients with the AA genotype of RPS6KB2:rs10274 had an 81% reduced risk of developing esophagitis (OR: 0.19, 95% CI: 0.07–0.51, p = 0.001, q = 0.06). Patients with the AG+GG genotype of SMO:rs1061280 had an 81% reduced risk of developing esophagitis (OR: 0.19, 95% CI: 0.07–0.53, p = 0.001, q = 0.06). Patients with the GG+GA genotype of DGCR8:rs720014 had a 3.54-fold increased risk of pneumonitis (OR: 3.54, 95% CI: 1.65–7.61, p <0.05, q <0.1). Significantly cumulative effects of the top SNPs were observed for both toxicities (P-trend <0.001). Using bioinformatics tools, we found that the genotype of rs10274 was associated with altered expression of the RPS6KB2 gene. Gene-based analysis showed DGCR8 (p = 0.010) and GEMIN4 (p = 0.039) were the top genes associated with the risk of developing pneumonitis. Our results provide strong evidence that microRNA-related genetic variations contribute to the development of radiotherapy-induced acute esophagitis and pneumonitis and could thus serve as biomarkers to help accurately predict radiotherapy-induced toxicity in NSCLC patients.     

Population-based meta-analysis in Caucasians confirms association with COL5A1 and ZNF469 but not COL8A2 with central corneal thickness

Central corneal thickness (CCT) has become an endophenotype of major interest for the genetically complex disorder glaucoma. CCT has a high heritability, and thin CCT is an independent risk factor for the diagnosis and progression of open-angle glaucoma. Genome-wide association studies thus provide genetic loci associated with CCT and potentially related to open-angle glaucoma. The distribution of CCT and prevalence of glaucoma in population-based studies have demonstrated ethnic differences suggesting ethnic-dependent variations in the genetic determinants of CCT. We conducted a genome-wide association study in Caucasians (n?=?3,931) from the Gutenberg Health Study (Germany) followed by replication of 30 genome-wide significant SNPs or SNPs of interest (P?<?10^?5) in the Rotterdam Study (The Netherlands, n?=?1,418). In a combined analysis, we confirmed quantitative trait loci on chromosomes 9q34 and 16q24 for association with CCT. On chromosome 16q24, the locus is located in an intergenic region near the ZNF469 gene (top SNP: rs9938149, P?=?1.45?×?10^?12). ZNF469 missense mutation is involved in a syndrome with very thin cornea (brittle cornea syndrome). The second locus on chromosome 9q34 represents the intergenic region between the RXRA and COL5A1 gene (top SNP: rs3132306, P?=?2.71?×?10^?10). Collagen type 5 determines the diameter of the corneal collagen fibrils. In our Caucasian population-based GWA study, we reinforce the involvement of collagen-related genes influencing CCT in Caucasians. We could not confirm the collagen type 8 locus on chromosome 1 as reported in Asian studies.     

Identification of potential functional variants and genes at 18q21.1 associated with the carcinogenesis of colorectal cancer

Genome-wide association studies (GWAS) have identified more than 160 susceptibility loci for colorectal cancer (CRC). The effects of these variants, particularly their mechanisms, however, remain unclear. In this study, a comprehensive functional annotation of CRC-related GWAS signals was firstly conducted to identify the potential causal variants. We found that the SNP rs7229639 in intron 3 of SMAD7 at 18q21.1 might serve as a putative functional variant in CRC. The SNP rs7229639 is located in a region with evidence of regulatory potential. Dual-luciferase reporter assays revealed that three other SNPs (rs77544449, rs60385309 and rs72917785), in strong linkage disequilibrium (LD) with rs7229639, exhibited allele-specific enhancer activity, of which one of the target genes may conceivably be LIPG, as suggested by eQTL association data and Hi-C data. We also verified that LIPG promoted malignancy of CRC cells in vitro, with supporting clinical data indicating that LIPG is upregulated and correlated with a poor prognosis in CRC. Finally, pitavastatin was observed to exhibit an anti-CRC activity and modest inhibition of LIPG mRNA levels. Collectively, our data suggest that these functional variants at 18q21.1 are involved in the pathogenesis of CRC by modulating enhancer activity, and possibly LIPG expression, thus indicating a promising therapeutic target for CRC. The results of functional annotation in our investigation could also serve as an inventory for CRC susceptibility SNPs and offer guides for post-GWAS downstream functional studies.     

Analysis of TNFAIP3, a feedback inhibitor of nuclear factor-κB and the neighbor intergenic 6q23 region in rheumatoid arthritis susceptibility

Introduction

Genome-wide association studies of rheumatoid arthritis (RA) have identified an association of the disease with a 6q23 region devoid of genes. TNFAIP3, an RA candidate gene, flanks this region, and polymorphisms in both the TNFAIP3 gene and the intergenic region are associated with systemic lupus erythematosus. We hypothesized that there is a similar association with RA, including polymorphisms in TNFAIP3 and the intergenic region.

Methods

To test this hypothesis, we selected tag-single nucleotide polymorphisms (SNPs) in both loci. They were analyzed in 1,651 patients with RA and 1,619 control individuals of Spanish ancestry.

Results

Weak evidence of association was found both in the 6q23 intergenic region and in the TNFAIP3 locus. The rs582757 SNP and a common haplotype in the TNFAIP3 locus exhibited association with RA. In the intergenic region, two SNPs were associated, namely rs609438 and rs13207033. The latter was only associated in patients with anti-citrullinated peptide antibodies. Overall, statistical association was best explained by the interdependent contribution of SNPs from the two loci TNFAIP3 and the 6q23 intergenic region.

Conclusions

Our data are consistent with the hypothesis that several RA genetic factors exist in the 6q23 region, including polymorphisms in the TNFAIP3 gene, like that previously described for systemic lupus erythematosus.     

Identification of a Sex-Linked SNP Marker in the Salmon Louse (Lepeophtheirus salmonis) Using RAD Sequencing

The salmon louse (Lepeophtheirus salmonis (Krøyer, 1837)) is a parasitic copepod that can, if untreated, cause considerable damage to Atlantic salmon (Salmo salar Linnaeus, 1758) and incurs significant costs to the Atlantic salmon mariculture industry. Salmon lice are gonochoristic and normally show sex ratios close to 1:1. While this observation suggests that sex determination in salmon lice is genetic, with only minor environmental influences, the mechanism of sex determination in the salmon louse is unknown. This paper describes the identification of a sex-linked Single Nucleotide Polymorphism (SNP) marker, providing the first evidence for a genetic mechanism of sex determination in the salmon louse. Restriction site-associated DNA sequencing (RAD-seq) was used to isolate SNP markers in a laboratory-maintained salmon louse strain. A total of 85 million raw Illumina 100 base paired-end reads produced 281,838 unique RAD-tags across 24 unrelated individuals. RAD marker Lsa101901 showed complete association with phenotypic sex for all individuals analysed, being heterozygous in females and homozygous in males. Using an allele-specific PCR assay for genotyping, this SNP association pattern was further confirmed for three unrelated salmon louse strains, displaying complete association with phenotypic sex in a total of 96 genotyped individuals. The marker Lsa101901 was located in the coding region of the prohibitin-2 gene, which showed a sex-dependent differential expression, with mRNA levels determined by RT-qPCR about 1.8-fold higher in adult female than adult male salmon lice. This study’s observations of a novel sex-linked SNP marker are consistent with sex determination in the salmon louse being genetic and following a female heterozygous system. Marker Lsa101901 provides a tool to determine the genetic sex of salmon lice, and could be useful in the development of control strategies.     

Comprehensive resequence analysis of a 97 kb region of chromosome 10q11.2 containing the MSMB gene associated with prostate cancer

Genome-wide association studies of prostate cancer have identified single nucleotide polymorphism (SNP) markers in a region of chromosome 10q11.2, harboring the microseminoprotein-β (MSMB) gene. Both the gene product of MSMB, the prostate secretory protein 94 (PSP94) and its binding protein (PSPBP), have been previously investigated as serum biomarkers for prostate cancer progression. Recent functional work has shown that different alleles of the significantly associated SNP in the promoter of MSMB found to be associated with prostate cancer risk, rs10993994, can influence its expression in tumors and in vitro studies. Since it is plausible that additional variants in this region contribute to the risk of prostate cancer, we have used next-generation sequencing technology to resequence a ~97-kb region that includes the area surrounding MSMB (chr10: 51,168,025–51,265,101) in 36 prostate cancer cases, 26 controls of European origin, and 8 unrelated CEPH individuals in order to identify additional variants to investigate in functional studies. We identified 241 novel polymorphisms within this region, including 142 in the 51-kb block of linkage disequilibrium (LD) that contains rs10993994 and the proximal promoter of MSMB. No sites were observed to be polymorphic within the exons of MSMB.     

The NOD2 p.Leu1007fsX1008 Mutation (rs2066847) Is a Stronger Predictor of the Clinical Course of Crohn's Disease than the FOXO3A Intron Variant rs12212067

Background

Very recently, a sub-analysis of genome-wide association scans revealed that the non-coding single nucleotide polymorphism (SNP) rs12212067 in the FOXO3A gene is associated with a milder course of Crohn''s disease (CD) (Cell 2013;155:57–69). The aim of our study was to evaluate the clinical value of the SNP rs12212067 in predicting the severity of CD by correlating CD patient genotype status with the most relevant complications of CD such as stenoses, fistulas, and CD-related surgery.

Methodology/Principal Findings

We genotyped 550 CD patients for rs12212067 (FOXO3A) and the three common CD-associated NOD2 mutations rs2066844, rs2066847, and rs2066847 and performed genotype-phenotype analyses.

Results

No significant phenotypic differences were found between the wild-type genotype TT of the FOXO3A SNP rs12212067 and the minor genotypes TG and GG independently from NOD2 variants. The allele frequency of the minor G allele was 12.7%. Age at diagnosis, disease duration, body mass index, surgery rate, stenoses, fistula, need for immunosuppressive therapy, and disease course were not significantly different. In contrast, the NOD2 mutant p.Leu1007fsX1008 (rs2066847) was highly associated with penetrating CD (p = 0.01), the development of fistulas (p = 0.01) and stenoses (p = 0.01), and ileal disease localization (p = 0.03). Importantly, the NOD2 SNP rs2066847 was a strong separator between an aggressive and a mild course of CD (p = 2.99×10⁻⁵), while the FOXO3A SNP rs12212067 did not separate between mild and aggressive CD behavior in our cohort (p = 0.35). 96.2% of the homozygous NOD2 p.Leu1007fsX1008 carriers had an aggressive disease behavior compared to 69.3% of the patients with the NOD2 wild-type genotype (p = 0.007).

Conclusion/Significance

In clinical practice, the NOD2 variant p.Leu1007fsX1008 (rs2066847), in particular in homozygous form, is a much stronger marker for a severe clinical phenotype than the FOXO3A rs12212067 SNP for a mild disease course on an individual patient level despite its important impact on the inflammatory response of monocytes.     

A GWAS Study on Liver Function Test Using eMERGE Network Participants

Introduction

Liver enzyme levels and total serum bilirubin are under genetic control and in recent years genome-wide population-based association studies have identified different susceptibility loci for these traits. We conducted a genome-wide association study in European ancestry participants from the Electronic Medical Records and Genomics (eMERGE) Network dataset of patient medical records with available genotyping data in order to identify genetic contributors to variability in serum bilirubin levels and other liver function tests and to compare the effects between adult and pediatric populations.

Methods

The process of whole genome imputation of eMERGE samples with standard quality control measures have been described previously. After removing missing data and outliers based on principal components (PC) analyses, 3294 samples from European ancestry were used for the GWAS study. The association between each single nucleotide polymorphism (SNP) and total serum bilirubin and other liver function tests was tested using linear regression, adjusting for age, gender, site, platform and ancestry principal components (PC).

Results

Consistent with previous results, a strong association signal has been detected for UGT1A gene cluster (best SNP rs887829, beta = 0.15, p = 1.30x10^-118) for total serum bilirubin level. Indeed, in this region more than 176 SNPs (or indels) had p<10⁻⁸ spanning 150Kb on the long arm of chromosome 2q37.1. In addition, we found a similar level of magnitude in a pediatric group (p = 8.26x10^-47, beta = 0.17). Further imputation using sequencing data as a reference panel revealed association of other markers including known TA7 repeat indels (rs8175347) (p = 9.78x10^-117) and rs111741722 (p = 5.41x10^-119) which were in proxy (r2 = 0.99) with rs887829. Among rare variants, two Asian subjects homozygous for coding SNP rs4148323 (G71R) were identified. Additional known effects for total serum bilirubin were also confirmed including organic anion transporters SLCO1B1-SLCO1B3, TDRP and ZMYND8 at FDR<0.05 with no gene-gene interaction effects. Phenome-wide association studies (PheWAS) suggest a protective effect of TA7 repeat against cerebrovascular disease in an adult cohort (OR = 0.75, p = 0.0008). Among other liver function tests, we also confirmed the previous effect of the ABO blood group locus for variation in serum alkaline phosphatase (rs579459, p = 9.44x10^-15).

Conclusions

Taken together, our data present interesting findings with strong confirmation of previous effects by simply using the eMERGE electronic health record phenotyping. In addition, our findings indicate that similar to the adult population, the UGT1A1 is the main locus responsible for normal variation of serum bilirubin in pediatric populations.