Impact of common regulatory single-nucleotide variants on gene expression profiles in whole blood

Genome-wide association studies (GWASs) have uncovered susceptibility loci for a large number of complex traits. Functional interpretation of candidate genes identified by GWAS and confident assignment of the causal variant still remains a major challenge. Expression quantitative trait (eQTL) mapping has facilitated identification of risk loci for quantitative traits and might allow prioritization of GWAS candidate genes. One major challenge of eQTL studies is the need for larger sample numbers and replication. The aim of this study was to evaluate the robustness and reproducibility of whole-blood eQTLs in humans and test their value in the identification of putative functional variants involved in the etiology of complex traits. In the current study, we performed comphrehensive eQTL mapping from whole blood. The discovery sample included 322 Caucasians from a general population sample (KORA F3). We identified 363 cis and 8 trans eQTLs after stringent Bonferroni correction for multiple testing. Of these, 98.6% and 50% of cis and trans eQTLs, respectively, could be replicated in two independent populations (KORA F4 (n=740) and SHIP-TREND (n=653)). Furthermore, we identified evidence of regulatory variation for SNPs previously reported to be associated with disease loci (n=59) or quantitative trait loci (n=20), indicating a possible functional mechanism for these eSNPs. Our data demonstrate that eQTLs in whole blood are highly robust and reproducible across studies and highlight the relevance of whole-blood eQTL mapping in prioritization of GWAS candidate genes in humans.     

Association of GWAS loci with PD in China

Genome‐wide association studies (GWAS) have identified numerous single‐nucleotide polymorphisms (SNPs) at four loci (SNCA, PARK16, LRRK2, BST1) that can modulate the risk of Parkinson's disease (PD). The strength of these associations has yet to be clarified in Mainland China. Ethnic specific effect is an important consideration in GWAS analysis. Using a case–control methodology, we genotyped multiple SNPs at these four loci to investigate their association with risk of PD in Mainland China. A total of 1,146 study subjects comprising 636 patients with PD and 510 unrelated healthy controls were recruited. The minor alleles at SNPs rs894278, rs1994090, rs2046932, rs4698412, and rs7304279 were found to be significantly higher in cases than in controls, while the minor alleles were found to significantly reduce the risk of developing PD at SNPs rs823128, rs823156, rs6532194, rs1191532, and rs16856139. These associations remained after taking into considerations the effects of age and gender. We showed that multiple SNPs at LRRK2 and SNCA increase risk of PD, while PARK16 SNPs are associated with a lower risk of PD in China. Our study findings will contribute to further research using GWAS‐linked data and research on ethnic specific effect of common variants. © 2011 Wiley‐Liss, Inc.     

The CELSR1 polymorphisms rs6007897 and rs4044210 are associated with ischaemic stroke in Chinese Han population

Background: Recently, CELSR1 was identified by genome-wide association studies (GWAS) as a susceptibility gene for ischaemic stroke (IS) in Japanese individuals.

Aim: The goal was to examine whether CELSR1 variants are associated with IS in the Chinese Han population.

Subjects and methods: This study genotyped two single nucleotide polymorphisms (SNPs) of CELSR1, rs6007897 and rs4044210, in a Chinese sample of 569 IS cases and 581 controls and assessed their genotype and allele associations with IS.

Results: The results showed that rs6007897 and rs4044210 variants of CELSR1 were significantly (p?p?Conclusion: Taken together, the present study has proven for the first time that CELSR1 is a susceptibility gene for IS in the Chinese Han population, especially for LAA.     

Meta‐analysis of association between obsessive‐compulsive disorder and the 3′ region of neuronal glutamate transporter gene SLC1A1

The neuronal glutamate transporter gene SLC1A1 is a candidate gene for obsessive‐compulsive disorder (OCD) based on linkage studies and convergent evidence implicating glutamate in OCD etiology. The 3′ end of SLC1A1 is the only genomic region with consistently demonstrated OCD association, especially when analyzing male‐only probands. However, specific allele associations have not been consistently replicated, and recent OCD genome‐wide association and meta‐analysis studies have not incorporated all previously associated SLC1A1 SNPs. To clarify the nature of association between SLC1A1 and OCD, pooled analysis was performed on all available relevant raw study data, comprising a final sample of 815 trios, 306 cases and 634 controls. This revealed weak association between OCD and one of nine tested SLC1A1 polymorphisms (rs301443; uncorrected P = 0.046; non‐significant corrected P). Secondary analyses of male‐affecteds only (N = 358 trios and 133 cases) demonstrated modest association between OCD and a different SNP (rs12682807; uncorrected P = 0.012; non‐significant corrected P). Findings of this meta‐analysis are consistent with the trend of previous candidate gene studies in psychiatry and do not clarify the putative role of SLC1A1 in OCD pathophysiology. Nonetheless, it may be important to further examine the potential associations demonstrated in this amalgamated sample, especially since the SNPs with modest associations were not included in the more highly powered recent GWAS or in a past meta‐analysis including five SLC1A1 polymorphisms. This study underscores the need for much larger sample sizes in future genetic association studies and suggests that next‐generation sequencing may be beneficial in examining the potential role of rare variants in OCD. © 2013 Wiley Periodicals, Inc.     

A case–control genome wide association study of substance use disorder (SUD) identifies novel variants on chromosome 7p14.1 in patients from the United Arab Emirates (UAE)

Genome wide association studies (GWASs) have provided insights into the molecular basis of the disorder in different population. This study presents the first GWAS of substance use disorder (SUD) in patients from the United Arab Emirates (UAE). The aim was to identify genetic association(s) that may provide insights into the molecular basis of the disorder. The GWAS discovery cohort consisted of 512 (250 cases and 262 controls) male participants from the UAE. Controls with no prior history of SUD were available from the Emirates family registry. The replication cohort consisted of 520 (415 cases and 105 controls) Australian male Caucasian participants. The GWAS discovery samples were genotyped for 4.6 million single nucleotide polymorphism (SNP). The replication cohort was genotyped using TaqMan assay. The GWAS association analysis identified three potential SNPs rs118129027 (p‐value = 6.24 × 10^?8), rs74477937 (p‐value = 8.56 × 10^?8) and rs78707086 (p‐value = 8.55 × 10^?8) on ch7p14.1, that did not meet the GWAS significance threshold but were highly suggestive. In the replication cohort, the association of the three top SNPs did not reach statistical significance. In a meta‐analysis of the discovery and the replication cohorts, there were no strengthen evidence for association of the three SNPs. The top identified rs118129027 overlaps with a regulatory factor (enhancer) region that targets three neighboring genes LOC105375237, LOC105375240, and YAE1D1. The YAE1D1, which represents a potential locus that is involved in regulating translation initiation pathway. Novel associations that require further confirmation were identified, suggesting a new insight to the genetic basis of SUD.     

A network of dopaminergic gene variations implicated as risk factors for schizophrenia

We evaluated the hypothesis that dopaminergic polymorphisms are risk factors for schizophrenia (SZ). In stage I, we screened 18 dopamine-related genes in two independent US Caucasian samples: 150 trios and 328 cases/501 controls. The most promising associations were detected with SLC6A3 (alias DAT), DRD3, COMT and SLC18A2 (alias VMAT2). In stage II, we comprehensively evaluated these four genes by genotyping 68 SNPs in all 478 cases and 501 controls from stage I. Fifteen (23.1%) significant associations were found (p < or = 0.05). We sought epistasis between pairs of SNPs providing evidence of a main effect and observed 17 significant interactions (169 tests); 41.2% of significant interactions involved rs3756450 (5' near promoter) or rs464049 (intron 4) at SLC6A3. In stage III, we confirmed our findings by genotyping 65 SNPs among 659 Bulgarian trios. Both SLC6A3 variants implicated in the US interactions were overtransmitted in this cohort (rs3756450, p = 0.035; rs464049, p = 0.011). Joint analyses from stages II and III identified associations at all four genes (p(joint) < 0.05). We tested 29 putative interactions from stage II and detected replication between seven locus pairs (p < or = 0.05). Simulations suggested our stage II and stage III interaction results were unlikely to have occurred by chance (p = 0.008 and 0.001, respectively). In stage IV we evaluated rs464049 and rs3756450 for functional effects and found significant allele-specific differences at rs3756450 using electrophoretic mobility shift assays and dual-luciferase promoter assays. Our data suggest that a network of dopaminergic polymorphisms increase risk for SZ.     

A genome-wide association study identifies FSHR rs2300441 associated with follicle-stimulating hormone levels

Follicle-stimulating hormone (FSH) and luteinizing hormone (LH) play critical roles in female reproduction, while the underlying genetic basis is poorly understood. Genome-wide association studies (GWASs) of FSH and LH levels were conducted in 2590 Chinese females including 1882 polycystic ovary syndrome (PCOS) cases and 708 controls. GWAS for FSH level identified multiple variants at FSHR showing genome-wide significance with the top variant (rs2300441) located in the intron of FSHR. The A allele of rs2300441 led to a reduced level of FSH in the PCOS group (β = −.43, P = 6.70 × 10⁻¹⁴) as well as in the control group (β = −.35, P = 6.52 × 10⁻⁴). In the combined sample, this association was enhanced after adjusting for the PCOS status (before: β = −.38, P = 1.77 × 10⁻¹³; after: β = −.42, P = 3.33 × 10⁻¹⁶), suggesting the genetic effect is independent of the PCOS status. The rs2300441 explained sevenfold higher proportion of the FSH variance than the total variance explained by the two previously reported FSHR missense variants (rs2300441 R² = 1.40% vs rs6166 R² = 0.17%, rs6165 R² = 0.03%). GWAS for LH did not identify any genome-wide significant associations. In conclusion, we identified genome-wide significant association between variants in FSHR and circulating FSH first, with the top associated variant rs2300441 might be a primary contributor at the population level.     

Population‐based identity‐by‐descent mapping combined with exome sequencing to detect rare risk variants for schizophrenia

Genome‐wide association studies (GWASs) are highly effective at identifying common risk variants for schizophrenia. Rare risk variants are also important contributors to schizophrenia etiology but, with the exception of large copy number variants, are difficult to detect with GWAS. Exome and genome sequencing, which have accelerated the study of rare variants, are expensive so alternative methods are needed to aid detection of rare variants. Here we re‐analyze an Irish schizophrenia GWAS dataset (n = 3,473) by performing identity‐by‐descent (IBD) mapping followed by exome sequencing of individuals identified as sharing risk haplotypes to search for rare risk variants in coding regions. We identified 45 rare haplotypes (>1 cM) that were significantly more common in cases than controls. By exome sequencing 105 haplotype carriers, we investigated these haplotypes for functional coding variants that could be tested for association in independent GWAS samples. We identified one rare missense variant in PCNT but did not find statistical support for an association with schizophrenia in a replication analysis. However, IBD mapping can prioritize both individual samples and genomic regions for follow‐up analysis but genome rather than exome sequencing may be more effective at detecting risk variants on rare haplotypes.     

Replication study of genome‐wide associated SNPs with late‐onset Alzheimer's disease

Late‐onset Alzheimer's disease (LOAD) is a multifactorial disease with the potential involvement of multiple genes. Four recent genome‐wide association studies (GWAS) have found variants showing significant association with LOAD on chromosomes 6, 10, 11, 12, 14, 18, 19, and on the X chromosome. We examined a total of 12 significant SNPs from these studies to determine if the results could be replicated in an independent large case–control sample. We genotyped these 12 SNPs as well the E2/E3/E4 APOE polymorphisms in up to 993 Caucasian Americans with LOAD and up to 976 age‐matched healthy Caucasian Americans. We found no statistically significant associations between the 12 SNPs and the risk of AD. Stratification by APOE*4 carrier status also failed to reveal statistically significant associations. Additional analyses were performed to examine potential associations between the 12 SNPs and age‐at‐onset (AAO) and disease duration among AD cases. Significant associations were observed between AAO and ZNF224/rs3746319 (P = 0.002) and KCNMA1/rs16934131 (P = 0.0066). KCNMA1/rs16934131 also demonstrated statistically significant association with disease duration (P = 0.0002). Although we have been unable to replicate the reported GWAS association with AD risk in our sample, we have identified two new associations with AAO and disease duration that need to be confirmed in additional studies. © 2011 Wiley‐Liss, Inc.     

Evaluation of the relationships of the WBP1L gene with schizophrenia and the general psychopathology scale based on a case–control study

WBP1L is a target of microRNA 137 (miR‐137) and has been considered a candidate gene for schizophrenia (SCZ). To investigate the relationships between WBP1L and SCZ and its related symptom scales, a total of 5,993 Chinese Han subjects, including 2,128 SCZ patients and 3,865 controls, were enrolled. In addition, an independent sample set for replication study including 1,052 SCZ patients and 2,124 controls were also recruited. Thirty‐two tag single nucleotide polymorphisms (SNPs) located within gene region of WBP1L were selected for genotyping and analyzing. The expression quantitative trait loci (eQTL) effects for the targeted SNPs were investigated with gene expression data from multiple human tissues. Rs4147157 (OR = 0.84, p = 1.51 × 10^?5) and rs284854 (OR = 1.14, p = 7.00 × 10^?4) were significantly associated with SCZ disease status and these association signals were replicated in our replication sample. A significant association was identified between rs4147157 and the general (β = ?.66, p = .001) and total (β = ?.8, p = .0042) scores of positive and negative syndrome scale scores in SCZ patients. Both SNPs were significant eQTL for genes around WBP1L in human brain tissues including ARL3 and AS3MT. To conclude, SNPs rs4147157 and rs284854 were associated with SCZ in the Chinese Han population. Additionally, rs4147157 was significantly associated with specific symptom features of SCZ.     

Impact of rare variants in ARHGAP29 to the etiology of oral clefts: role of loss‐of‐function vs missense variants

Non‐syndromic cleft lip with or without cleft palate (NSCL/P) is a prevalent, complex congenital malformation. Genome‐wide association studies (GWAS) on NSCL/P have consistently identified association for the 1p22 region, in which ARHGAP29 has emerged as the main candidate gene. ARHGAP29 re‐sequencing studies in NSCL/P patients have identified rare variants; however, their clinical impact is still unclear. In this study we identified 10 rare variants in ARHGAP29, including five missense, one in‐frame deletion, and four loss‐of‐function (LoF) variants, in a cohort of 188 familial NSCL/P cases. A significant mutational burden was found for LoF (Sequence Kernel Association Test, p = 0.0005) but not for missense variants in ARHGAP29, suggesting that only LoF variants contribute to the etiology of NSCL/P. Penetrance was estimated as 59%, indicating that heterozygous LoF variants in ARHGAP29 confer a moderate risk to NSCL/P. The GWAS hits in IRF6 (rs642961) and 1p22 (rs560426 and rs4147811) do not seem to contribute to the penetrance of the phenotype, based on co‐segregation analysis. Our data show that rare variants leading to haploinsufficiency of ARHGAP29 represent an important etiological clefting mechanism, and genetic testing for this gene might be taken into consideration in genetic counseling of familial cases.     

Evaluation of risk loci for schizophrenia derived from genome‐wide association studies in a German population

In the genome‐wide association study (GWAS) on schizophrenia [O'Donovan et al. (2008); Nat Genet 40:1053–1055] a UK‐sample of 479 cases with DSM‐IV schizophrenia was genotyped in comparison to control subjects with follow up of 12 putative loci in international replication sets of approximately 15,000 cases and controls. In these cohorts and a combined bipolar and schizophrenia UK‐sample, six single nucleotide polymorphisms (SNPs) supported association, with the strongest evidence for SNP‐marker rs1344706 at the zinc finger ZNF804A locus on chromosome 2q32.1 (P = 1.61 × 10^?7). We attempted replication of these findings in a German population of 2,154 individuals (632 with affective disorders, 937 with schizophrenia, and 585 controls), but found none of the GWAS risk alleles significantly associated with psychosis. Particularly rs1344706, initially surpassing the genome‐wide significance level in an extended phenotype of schizophrenia and affective disorder, produced consistently negative results. At the ZNF804A locus estimated Odds ratios reached 1.08 (0.93–1.26 95% CI) for the schizophrenia sample and 1.04 (0.90–1.20 95% CI) for the combined set of cases with schizophrenia and affective disorder. The main limitation of our study may be the reduced power of the sample size, but our data may be useful for future meta‐analysis of GWA data sets. Although GWAS have proven extraordinary successful in identifying susceptibility genes for complex genetic disorders, the hypothesis of common genetic variants in the complex group of the schizophrenic psychoses with small effect size but relatively high frequency is still put to further scrutiny. © 2011 Wiley‐Liss, Inc.     

A genome-wide association study identifies new genes associated with developmental dysplasia of the hip

Developmental dysplasia of the hip (DDH) is one of the most common congenital malformations and covers a spectrum of hip disorders from mild dysplasia to irreducible dislocation. The pathological mechanisms of DDH are poorly understood, which hampers the development of diagnostic tools and treatments. To gain insight into its disease mechanism, we explored the potential biological processes that underlie DDH by integrating pathway analysis tools and performing a genome-wide association study (GWAS). A total of 406 DDH-associated genes (P < 0.001) were identified by our GWAS using a Chinese Han cohort consisting of 386 DDH cases and 500 healthy controls (Set A). We verified the significant loci (P < 10⁻⁵) in another Chinese Han cohort consisting of 574 DDH patients and 569 healthy controls (Set B). An intronic Single Nucleotide Polymorphism (SNP) (rs61930502) showed significant association in Set A and Set B (P = 2.65 × 10⁻⁷ and 2.0 × 10⁻⁴, respectively). The minor allele, rs61930502-A, which tended to prevent DDH showed a dominant effect. Heat shock 70 kDa protein 8 (HSPA8) showed the most direct interactions with other proteins which were coded by DDH-associated genes in the protein-protein interaction analysis. Interestingly, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis suggested a relation between DDH and the genes involved in type II diabetes mellitus pathway (P = 0.0067). Our genetic and protein interaction evidence could open avenues for future studies of DDH.     

The CHRNA5/CHRNA3/CHRNB4 Nicotinic Receptor Regulome: Genomic Architecture,Regulatory Variants,and Clinical Associations

Functionally related genes often cluster into a genome region under coordinated regulation, forming a local regulome. To understand regulation of the CHRNA5/CHRNA3/CHRNB4 nicotinic receptor gene cluster, we integrate large‐scale RNA expression data (brain and peripheral) from GTEx (Genotype Tissue Expression), clinical associations (GRASP), and linkage disequilibrium data (1000 Genomes) to find candidate SNPs representing independent regulatory variants. CHRNA3, CHRNA5, CHRNB4 mRNAs, and a well‐expressed CHRNA5 antisense RNA (RP11‐650L12.2) are co‐expressed in many human tissues, suggesting common regulatory elements. The CHRNA5 enhancer haplotype tagged by rs880395 not only increases CHRNA5 mRNA expression in all tissues, but also enhances RP11‐650L12.2 and CHRNA3 expression, suggesting DNA looping to multiple promoters. However, in nucleus accumbens and putamen, but not other brain regions, CHRNA3 expression associates uniquely with a haplotype tagged by rs1948 (located in the CHRNB4 3′UTR). Haplotype/diplotype analysis of rs880395 and rs1948 plus rs16969968 (a nonsynonymous CHRNA5 risk variant) in GWAS (COGEND, UW‐TTURC, SAGE) yields a nicotine dependence risk profile only partially captured by rs16969968 alone. An example of local gene clusters, this nicotinic regulome is controlled by complex genetic variation, with broad implications for interpreting GWAS.     

A non-synonymous variant in SLC30A8 is not associated with type 1 diabetes in the Danish population

Genome-wide association scans in type 2 diabetes (T2D) have identified a risk variant, rs13266634 (Arg325Trp), in SLC30A8 on chromosome 8. SLC30A8 encodes a β-cell specific zinc-ion transporter and rs13266634 has been shown to affect insulin secretion. Recently, autoantibodies for Slc30A8 with high predictive value were demonstrated in individuals with type 1 diabetes (T1D), making this gene an interesting T1D candidate gene. We genotyped rs13266634 in 3008 cases and controls and 246 families from Denmark. Association to T1D could not be demonstrated.     

Common variants in DLG1 locus are associated with non‐syndromic cleft lip with or without cleft palate

Non‐syndromic cleft lip with or without cleft palate (nsCL/P) is a common craniofacial anomaly with a complex and heterogeneous aetiology. Knowledge regarding specific genetic factors underlying this birth defect is still not well understood. Therefore, we conducted an independent replication analysis for the top‐associated variants located within the DLG1 locus at chromosome 3q29, which was identified as a novel cleft‐susceptibility locus in our genome‐wide association study (GWAS). Mega‐analysis of the pooled individual data from the GWAS and replication study confirmed that common DLG1 variants are associated with the risk of nsCL/P. Two single nucleotide polymorphisms (SNPs), rs338217 and rs7649443, were statistically significant even at the genome‐wide level (P_trend = 9.70E?10 and P_trend = 8.96E?09, respectively). Three other SNPs, rs9826379, rs6805920 and rs6583202, reached a suggestive genome‐wide significance threshold (P_trend < 1.00E?05). The location of the strongest individual SNP in the intronic sequence of the gene encoding DLG1 antisense RNA suggests that the true causal variant implicated in the risk of nsCL/P may affect the DLG1 gene expression level rather than structure of the encoded protein. In conclusion, we identified a novel cleft‐susceptibility locus at chromosome 3q29 with a DLG1 as a novel candidate gene for this common craniofacial anomaly.