Evaluation of Five Interleukin Genes for Association with End-Stage Renal Disease in White Europeans

Background: Genetic variation within interleukin genes has been reported to be associated with end-stage renal disease (ESRD). These findings have not been consistently replicated. No study has yet reported the comprehensive investigation of IL1A, IL1B, IL1RN, IL6 and IL10 genes. Methods: 664 kidney transplant recipients (cases) and 577 kidney donors (controls) were genotyped to establish if common variants in interleukin genes are associated with ESRD. Single nucleotide polymorphism (SNP) genotype data for each gene were downloaded for a northern and western European population from the International HapMap Project. Haploview was used to visualize linkage disequilibrium and select tag SNPs. Thirty SNPs were genotyped using MassARRAY (R) iPLEX Gold technology and data were analyzed using the chi(2) test for trend. Independent replication was conducted in 1,269 individuals with similar phenotypic characteristics. Results: Investigating all common variants in IL1A, IL1B, IL1RN, IL6 and IL10 genes revealed a statistically significant association (rs452204 p(empirical) = 0.02) with one IL1RN variant and ESRD. This IL1RN SNP tags three other variants, none of which have previously been reported to be associated with renal disease. Independent replication in a separate transplant population of comparable size did not confirm the original observation. Conclusions: Common variants in these five candidate interleukin genes are not major risk factors for ESRD in white Europeans. Copyright (C) 2010 S. Karger AG, Basel

The role of human demographic history in determining the distribution and frequency of transferase-deficient galactosaemia mutations

Classical or transferase-deficient galactosaemia is an inherited metabolic disorder caused by mutation in the human Galactose-1-phosphate uridyl transferase (GALT) gene. Of some 170 causative mutations reported, fewer than 10% are observed in more than one geographic region or ethnic group. To better understand the population history of the common GALT mutations, we have established a haplotyping system for the GALT locus incorporating eight single nucleotide polymorphisms and three short tandem repeat markers. We analysed haplotypes associated with the three most frequent GALT gene mutations, Q188R, K285N and Duarte-2 (D2), and estimated their age. Haplotype diversity, in conjunction with measures of genetic diversity and of linkage disequilibrium, indicated that Q188R and K285N are European mutations. The Q188R mutation arose in central Europe within the last 20 000 years, with its observed east-west cline of increasing relative allele frequency possibly being due to population expansion during the re-colonization of Europe by Homo sapiens in the Mesolithic age. K285N was found to be a younger mutation that originated in Eastern Europe and is probably more geographically restricted as it arose after all major European population expansions. The D2 variant was found to be an ancient mutation that originated before the expansion of Homo sapiens out of Africa. Heredity (2010) 104, 148-154; doi:10.1038/hdy.2009.84; published online 29 July 2009

Replication of association between schizophrenia and ZNF804A in the Irish Case-Control Study of Schizophrenia sample

A recent genome-wide association study reported association between schizophrenia and the ZNF804A gene on chromosome 2q32.1. We attempted to replicate these findings in our Irish Case-Control Study of Schizophrenia (ICCSS) sample (N=1021 cases, 626 controls). Following consultation with the original investigators, we genotyped three of the most promising single-nucleotide polymorphisms (SNPs) from the Cardiff study. We replicate association with rs1344706 (trend test one-tailed P=0.0113 with the previously associated A allele) in ZNF804A. We detect no evidence of association with rs6490121 in NOS1 (one-tailed P=0.21), and only a trend with rs9922369 in RGRIP1L (one-tailed P=0.0515). On the basis of these results, we completed genotyping of 11 additional linkage disequilibrium-tagging SNPs in ZNF804A. Of 12 SNPs genotyped, 11 pass quality control criteria and 4 are nominally associated, with our most significant evidence of association at rs7597593 (P=0.0013) followed by rs1344706. We observe no evidence of differential association in ZNF804A on the basis of family history or sex of case. The associated SNP rs1344706 lies in approximately 30 bp of conserved mammalian sequence, and the associated A allele is predicted to maintain binding sites for the brain-expressed transcription factors MYT1l and POU3F1/OCT-6. In controls, expression is significantly increased from the A allele of rs1344706 compared with the C allele. Expression is increased in schizophrenic cases compared with controls, but this difference does not achieve statistical significance. This study replicates the original reported association of ZNF804A with schizophrenia and suggests that there is a consistent link between the A allele of rs1344706, increased expression of ZNF804A and risk for schizophrenia.

Two non-synonymous markers in PTPN21, identified by genome-wide association study data-mining and replication, are associated with schizophrenia

We conducted data-mining analyses of genome wide association (GWA) studies of the CATIE and MGS-GAIN datasets, and found 13 markers in the two physically linked genes, PTPN21 and EML5, showing nominally significant association with schizophrenia. Linkage disequilibrium (LD) analysis indicated that all 7 markers from PTPN21 shared high LD (r(2)>0.8), including rs2274736 and rs2401751, the two non-synonymous markers with the most significant association signals (rs2401751, P=1.10 × 10(-3) and rs2274736, P=1.21 × 10(-3)). In a meta-analysis of all 13 replication datasets with a total of 13,940 subjects, we found that the two non-synonymous markers are significantly associated with schizophrenia (rs2274736, OR=0.92, 95% CI: 0.86-0.97, P=5.45 × 10(-3) and rs2401751, OR=0.92, 95% CI: 0.86-0.97, P=5.29 × 10(-3)). One SNP (rs7147796) in EML5 is also significantly associated with the disease (OR=1.08, 95% CI: 1.02-1.14, P=6.43 × 10(-3)). These 3 markers remain significant after Bonferroni correction. Furthermore, haplotype conditioned analyses indicated that the association signals observed between rs2274736/rs2401751 and rs7147796 are statistically independent. Given the results that 2 non-synonymous markers in PTPN21 are associated with schizophrenia, further investigation of this locus is warranted.

A population-based association study of SNPs of GSTP1, MnSOD, GPX2 and Barrett's esophagus and esophageal adenocarcinoma

Oxidative stress appears to be important in the pathogenesis of Barrett's esophagus (BE) and esophageal adenocarcinoma (EAC). Single-nucleotide polymorphisms (SNPs) of antioxidant enzyme genes may play a part in determining individual susceptibility to these diseases. The Factors Influencing the Barrett's Adenocarcinoma Relationship (FINBAR) study is a population-based, case-control study of BE and EAC in Ireland. DNA from EAC (n = 207), BE (> or =3 cm BE at endoscopy with specialized intestinal metaplasia on biopsy, n = 189) and normal population controls (n = 223) were analyzed. Several SNPs spanning the genes for glutathione S-transferase P1 (GSTP1), manganese superoxide dismutase (MnSOD) and glutathione peroxidase 2 (GPX2) were genotyped using multiplex polymerase chain reaction and SNaPshottrade mark. The chi(2) test was used to compare genotype and allele frequencies between case and control subjects. Linkage disequilibrium between SNPs was quantified using Lewontin's D' value and haplotype frequency estimates obtained using Haploview. Eleven SNPs were genotyped (six for GSTP1, three for MnSOD and two for GPX2); all were in Hardy-Weinberg equilibrium. None was significantly associated with EAC or BE even before Bonferroni correction. Odds ratios for EAC for individual SNPs ranged from 0.68 [95% confidence interval (CI) 0.43-1.08] to 1.25 (95% CI 0.73-2.16), and for BE from 0.84 (95% CI 0.52-1.30) to 1.30 (95% CI 0.85-1.97). SNPs in all three genes were in strong linkage disequilibrium (D' > 0.887) but haplotype analysis did not show any significant association with EAC or BE. SNPs involving the GSTP1, MnSOD and GPX2 genes were not associated with BE or EAC. Further studies aimed at identifying susceptibility genes should focus on different antioxidant genes or different pathways.

Dense fine-mapping study identifies new susceptibility loci for primary biliary cirrhosis

We genotyped 2,861 cases of primary biliary cirrhosis (PBC) from the UK PBC Consortium and 8,514 UK population controls across 196,524 variants within 186 known autoimmune risk loci. We identified 3 loci newly associated with PBC (at P <5 × 10(-8)), increasing the number of known susceptibility loci to 25. The most associated variant at 19p12 is a low-frequency nonsynonymous SNP in TYK2, further implicating JAK-STAT and cytokine signaling in disease pathogenesis. An additional five loci contained nonsynonymous variants in high linkage disequilibrium (LD; r(2) > 0.8) with the most associated variant at the locus. We found multiple independent common, low-frequency and rare variant association signals at five loci. Of the 26 independent non-human leukocyte antigen (HLA) signals tagged on the Immunochip, 15 have SNPs in B-lymphoblastoid open chromatin regions in high LD (r(2) > 0.8) with the most associated variant. This study shows how data from dense fine-mapping arrays coupled with functional genomic data can be used to identify candidate causal variants for functional follow-up.

Investigation of Genetic Variation in Scavenger Receptor Class B, Member 1 (SCARB1) and Association with Serum Carotenoids

Objective: To investigate association of scavenger receptor class B, member 1 (SCARB1) genetic variants with serum carotenoid levels of lutein (L) and zeaxanthin (Z) and macular pigment optical density (MPOD).
Design: A cross-sectional study of healthy adults aged 20 to 70.
Participants: We recruited 302 participants after local advertisement.
Methods: We measured MPOD by customized heterochromatic flicker photometry. Fasting blood samples were taken for serum L and Z measurement by high-performance liquid chromatography and lipoprotein analysis by spectrophotometric assay. Forty-seven single nucleotide polymorphisms (SNPs) across SCARB1 were genotyped using Sequenom technology. Association analyses were performed using PLINK to compare allele and haplotype means, with adjustment for potential confounding and correction for multiple comparisons by permutation testing. Replication analysis was performed in the TwinsUK and Carotenoids in Age-Related Eye Disease Study (CAREDS) cohorts.
Main Outcome Measures: Odds ratios for MPOD area, serum L and Z concentrations associated with genetic variations in SCARB1 and interactions between SCARB1 and gender.
Results: After multiple regression analysis with adjustment for age, body mass index, gender, high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, triglycerides, smoking, and dietary L and Z levels, 5 SNPs were significantly associated with serum L concentration and 1 SNP with MPOD (P<0.01). Only the association between rs11057841 and serum L withstood correction for multiple comparisons by permutation testing (P<0.01) and replicated in the TwinsUK cohort (P = 0.014). Independent replication was also observed in the CAREDS cohort with rs10846744 (P = 2×10-4), an SNP in high linkage disequilibrium with rs11057841 (r2 = 0.93). No interactions by gender were found. Haplotype analysis revealed no stronger association than obtained with single SNP analyses.
Conclusions: Our study has identified association between rs11057841 and serum L concentration (24% increase per T allele) in healthy subjects, independent of potential confounding factors. Our data supports further evaluation of the role for SCARB1 in the transport of macular pigment and the possible modulation of age-related macular degeneration risk through combating the effects of oxidative stress within the retina.
Financial Disclosure(s): Proprietary or commercial disclosures may be found after the references. Ophthalmology 2013;120:1632–1640 © 2013 by the American Academy of Ophthalmology.

Five authors reply

We read with interest the comments offered by Drs. Hughes and Bradley (1) on our systematic review (2). Four single nucleotide polymorphisms (SNPs), rs9332739 and rs547154 in the complement component 2 gene (C2) and rs4151667 and rs641153 in the complement factor B gene (CFB), were pooled. Hughes and Bradley point out that we omitted the most common variant, rs12614. In fact, rs12614 is in high linkage disequilibrium (LD) with rs641153, which was included, and the major allele of both of these SNPs is in the range of 90% (population code, CEU, in the International HapMap Project (http://hapmap.ncbi.nlm.nih.gov/)). Moreover, our review was initiated in September 2010, at which point only 4 studies had published associations with rs12614, whereas 14 studies (n = 11,378) were available for rs641153. While it is true that both SNPs are better analyzed as a haplotype, these data were simply not available for pooling.
Hughes and Bradley also point out that we obtained and pooled new data that were not previously published. While it is recommended that contact with authors be completed as part of a comprehensive meta-analysis, we acknowledge that these additional data were not previously published and peer reviewed and, hence, do not have the same level of transparency. However, given that sample collections often increase over time and that the instrumentation for genotyping is continually improving, we thought that it would be advantageous to use the most recent information; this is a subjective decision.
We also agree that the allele frequencies given by Kaur et al. (3) were exactly opposite to those expected and were suggestive of strand flipping. However, we specifically queried this with the lead author on 2 separate occasions and were assured it was not.
Hughes and Bradley do make an interesting suggestion that SNPs in high LD should be used as a gauge of genotyping quality in HuGE reviews. This is an interesting idea but difficult to put into practice as the r2 parameter they propose as a measure of LD has some unusual properties. Although r2 is a measure of LD, it is also linked to the allele frequency; even small differences in allele frequencies between 2 linked SNPs can reduce the r2 dramatically. Wray (4) explored these effects and found that, at a baseline allele frequency of 10%, even a difference in allele frequency between 2 SNPs as small as 2% can drop the r2 value below 0.8. This degree of allele frequency difference is consistent with what could be expected for sampling error. Furthermore, when we look at 2 linked dialleleic SNPs, giving 4 possible haplotypes, the absence of 1 haplotype dramatically reduces r2, despite the 2 loci being in high LD as measured by D'. In fact, this is the situation for rs12614 and rs641153, where the low frequency of 1 haplotype means that the r2 is 0.01 but the D' is 1.
Hughes and Bradley also suggest consideration of genotype call rate restrictions as an inclusion criterion for metaanalysis. This would be more appropriate when focusing on genetic variants per se, as considered within the context of a genome-wide association study or other specific genetic analysis where large numbers of SNPs are evaluated (5).
The concerns raised by Hughes and Bradley reflect the limited ability of a meta-analysis based on summary data to tease out inconsistencies best identified at the individual level. We agree that SNPs in LD should be evaluated, but this will not necessarily be straightforward. A move to make genetic data sets publicly available, as in the Database of Genotypes and Phenotypes (http://www.ncbi.nlm.nih.gov/ gap), is a step in the right direction for greater transparency.

Identification of a high-risk haplotype for the dystrobrevin binding protein 1 (DTNBP1) gene in the Irish study of high-density schizophrenia families

A recent report showed significant associations between several SNPs in a previously unknown EST cluster with schizophrenia. (1). The cluster was identified as the human dystrobrevin binding protein 1 gene (DTNBP1) by sequence database comparisons and homology with mouse DTNBP1. (2). However, the linkage disequilibrium (LD) among the SNPs in DTNBP1 as well as the pattern of significant SNP-schizophrenia association was complex. This raised several questions such as the number of susceptibility alleles that may be involved and the size of the region where the actual disease mutation(s) could be located. To address these questions, we performed different single-marker tests on the 12 previously studied and 2 new SNPs in DTNBP1 that were re-scored using an improved procedure, and performed a variety of haplotype analyses. The sample consisted of 268 Irish multiplex families selected for high density of schizophrenia. Results suggested a simple structure where the LD in the target region could be explained by 6 haplotypes that together accounted for 96% of haplotype diversity in the whole sample. From these six, a single high-risk haplotype was identified that showed a significant association with schizophrenia and explained the pattern of significant findings in the analyses with individual markers. This haplotype was 30 kb long, had a large effect, could be measured with two tag SNPs only, had a frequency of 6% in our sample, seemed to be of relatively recent origin in evolutionary terms, and was equally distributed over Ireland. Implications of these findings for follow-up and replication studies are discussed.

Variants in the catechol-o-methyltransferase (COMT) gene are associated with schizophrenia in Irish high-density families

The enzyme catechol-o-methyltransferase (COMT) transfers a methyl group from adenosylmethionine to catecholamines including the neurotransmitters dopamine, epinephrine and norepinephrine. This methylation results in the degradation of catecholamines. The involvement of the COMT gene in the metabolic pathway of these neurotransmitters has made it an attractive candidate gene for many psychiatric disorders. In this article, we reported our study of association of COMT with schizophrenia in Irish families with a high density of schizophrenia. Three single nucleotide polymorphisms (SNPs) were genotyped for the 274 such families and within-family transmission disequilibrium tests were performed. SNP rs4680, which is the functional Val/Met polymorphism, showed modest association with the disease by the TRANSMIT, FBAT and PDT programs, while the other two SNPs were negative. These SNPs showed lower level of LDs with each other in the Irish subjects than in Ashkenazi Jews. Haplotype analysis indicated that a haplotype, haplotype A-G-A for SNPs rs737865-rs4680-rs165599, was preferentially transmitted to the affected subjects. This was different from the reported G-G-G haplotype found in Ashkenazi Jews, but both haplotypes shared the Val allele. We concluded that COMT gene is associated with schizophrenia and carries a small but significant risk to the susceptibility in the Irish subjects.

Regulator of G-protein signaling 4 (RGS4) gene is associated with schizophrenia in Irish high density families

The regulator of the G-protein signaling 4 (RGS4) gene was shown to have a different expression pattern in schizophrenia patients in a microarray study. A family-based study subsequently implicated the association of this gene with schizophrenia. We replicated the study with our sample from the Irish Study of High Density Schizophrenia Families (ISHDSF). Single marker transmission disequilibrium tests (TDT) for the four core SNPs showed modest association for SNP 18 (using a narrow diagnostic approach with FBAT P = 0.044; with PDT P = 0.0073) and a trend for SNP 4 (with FBAT P = 0.1098; with PDT P = 0.0249). For SNP 1 and 7, alleles overtransmitted to affected subjects were the same as previously reported. Haplotype analyses suggested that haplotype G-G-G for SNP1-4-18, which is the most abundant haplotype (42.3%) in the Irish families, was associated with the disease (narrow diagnosis, FBAT P = 0.0061, PDT P = 0.0498). This was the same haplotype implicated in the original study. While P values were not corrected for multiple testing because of the clear prior hypothesis, these results could be interpreted as supporting evidence for the association between RGS4 and schizophrenia.

Association between the 5q31.1 gene neurogenin1 and schizophrenia

Multiple lines of evidence suggest that schizophrenia results from aberrant neurodevelopment. The neurogenin1 gene (neurog1) consists of a single 1,666 bp exon that encodes a basic helix-loop-helix (bHLH) transcription factor that causes neuronal differentiation and induces cortical and glutamatergic differentiation programs. Because of its function and its location in 5q31.1, which has been linked to schizophrenia in multiple samples, we tested it for association with the disorder. We sequenced neurog1 in 25 affected subjects from the Irish Study of High-Density Schizophrenia Families. We observed a 5'-UTR SNP at position -60, already present in databases as rs8192558, and tested it along with rs2344485, rs8192559, and rs2344484. Narrow, intermediate, and broad diagnostic definitions were used. The major alleles of rs8192558 and rs2344484 were over-transmitted to affected subjects using both Pedigree Disequilibrium Test (PDT) (0.01 <or = P <or = 0.06) and FBAT (0.02 <or = P <or = 0.07). A haplotype consisting of the major alleles of all four SNPs was significantly over-transmitted in FBAT to the broad definition (P = 0.049), with trend significance to the narrow and intermediate definitions, and with trend significance in PDT. In confirmatory tests using 657 cases and 411 controls, this haplotype was slightly but not significantly over-represented in cases (81% vs. 77%, P = 0.21). These results, along with a priori evidence for the involvement of neurog1 in neurodevelopment, suggest that variants in neurog1 might have a small effect on susceptibility to schizophrenia. This gene should be tested in additional and larger samples.

An inherited duplication at the gene p21 Protein-Activated Kinase 7 (PAK7) is a risk factor for psychosis

Identifying rare, highly penetrant risk mutations may be an important step in dissecting the molecular etiology of schizophrenia. We conducted a gene-based analysis of large (>100kb), rare copy number variants (CNVs) in the Wellcome Trust Case Control Consortium 2 (WTCCC2) schizophrenia sample of 1,564 cases and 1,748 controls all from Ireland, and further extended the analysis to include an additional 5,196 UK controls. We found association with duplications at chr20p12.2 (P=0.007) and evidence of replication in large independent European schizophrenia (P=0.052) and UK bipolar disorder case-control cohorts (P=0.047). A combined analysis of Irish/UK subjects including additional psychosis cases (schizophrenia and bipolar disorder) identified 22 carriers in 11,707 cases and 10 carriers in 21,204 controls (meta-analysis CMH P value=2x10(-4) (odds ratio (OR)=11.3, 95% CI=3.7, ∞)). Nineteen of the 22 cases and 8 of the 10 controls carried duplications starting at 9.68Mb with similar breakpoints across samples. By haplotype analysis and sequencing we identified a tandem ∼149kb duplication overlapping the gene p21 Protein-Activated Kinase 7 (PAK7, also called PAK5) which was in linkage disequilibrium with local haplotypes (P=2.5x10(-21)), indicative of a single ancestral duplication event. We confirmed the breakpoints in 8/8 carriers tested and found co-segregation of the duplication with illness in two additional family members of one of the affected probands. We demonstrate that PAK7 is developmentally co-expressed with another known psychosis risk gene (DISC1) suggesting a potential molecular mechanism involving aberrant synapse development and plasticity.

FLT1 Genetic Variation Predisposes to Neovascular AMD in Ethnically Diverse Populations and Alters Systemic FLT1 Expression

Purpose: Current understanding of the genetic risk factors for age-related macular degeneration (AMD) is not sufficiently predictive of the clinical course. The VEGF pathway is a key therapeutic target for treatment of neovascular AMD; however, risk attributable to genetic variation within pathway genes is unclear. We sought to identify single nucleotide polymorphisms (SNPs) associated with AMD within the VEGF pathway.
Methods: Using a tagSNP, direct sequencing and meta-analysis approach within four ethnically diverse cohorts, we identified genetic risk present in FLT1, though not within other VEGF pathway genes KDR, VEGFA, or VASH1. We used ChIP and ELISA in functional analysis.
Results: The FLT1 SNPs rs9943922, rs9508034, rs2281827, rs7324510, and rs9513115 were significantly associated with increased risk of neovascular AMD. Each association was more significant after meta-analysis than in any one of the four cohorts. All associations were novel, within noncoding regions of FLT1 that do not tag for coding variants in linkage disequilibrium. Analysis of soluble FLT1 demonstrated higher expression in unaffected individuals homozygous for the FLT1 risk alleles rs9943922 (P = 0.0086) and rs7324510 (P = 0.0057). In silico analysis suggests that these variants change predicted splice sites and RNA secondary structure, and have been identified in other neovascular pathologies. These data were supported further by murine chromatin immunoprecipitation demonstrating that FLT1 is a target of Nr2e3, a nuclear receptor gene implicated in regulating an AMD pathway.
Conclusions: Although exact variant functions are not known, these data demonstrate relevancy across ethnically diverse genetic backgrounds within our study and, therefore, hold potential for global efficacy.

Isolation and characterization of eight polymorphic microsatellite loci for Natterer's bat, Myotis nattereri (Vespertilionidae, Chiroptera)

Eight new microsatellite loci were isolated and characterized for the Natterer's bat Myotis nattereri from a microsatellite-enriched genomic library. The usefulness of these markers was assessed by screening a sample comprising 100 specimens collected from throughout the species range in Europe. Both moderately and highly polymorphic loci were identified with 3-17 alleles segregating per locus (mean 8.1 SE +/- A 0.048). No evidence for departure from HWE or linkage disequilibrium among loci was observed. These markers will provide a valuable addition to the molecular toolbox currently available for studies of population genetic structure, parentage and social organisation of M. nattereri and related species.