Bipolar disorder susceptibility region on chromosome 3q29 not confirmed in a case-control association study

We identified a bipolar disorder (BPD) susceptibility region on chromosome 3q29 in a genome-wide linkage scan (Bailer et al. 2002 (Biol Psychiatry 52: 40), NPL-score 4.09) and follow-up linkage analysis (Schosser et al. 2004 (J Psychiatr Res 38(3): 357), NPL-scores >3 with five markers). These findings were supported by further fine-mapping of this region (Schosser et al. 2007 (Eur Neuropsychopharmacol 17(6-7): 501)), finding NPL-scores >3.9 with SNPs (single nucleotide polymorphisms) spanning a region of 3.46 Mbp in BPD families. Since genetic association studies are more powerful than linkage studies for detecting susceptibility genes of small effect size, we aimed to replicate these findings in an independent case-control sample collected in London (UK) and Vienna (Austria).

STATISTICAL METHODS FOR RARE AND COMMON VARIANT ASSOCIATION STUDIES

Complex human diseases are a major challenge for biological research. The goal of my research is to develop effective methods for biostatistics in order to create more opportunities for the prevention and cure of human diseases. This dissertation proposes statistical technologies that have the ability of being adapted to sequencing data in family-based designs, and that account for joint effects as well as gene-gene and gene-environment interactions in the GWA studies. The framework includes statistical methods for rare and common variant association studies. Although next-generation DNA sequencing technologies have made rare variant association studies feasible, the development of powerful statistical methods for rare variant association studies is still underway. Chapter 2 demonstrates two adaptive weighting methods for rare variant association studies based on family data for quantitative traits. The results show that both proposed methods are robust to population stratification, robust to the direction and magnitude of the effects of causal variants, and more powerful than the methods using weights suggested by Madsen and Browning [2009]. In Chapter 3, I extended the previously proposed test for Testing the effect of an Optimally Weighted combination of variants (TOW) [Sha et al., 2012] for unrelated individuals to TOW &ndash F, TOW for Family &ndash based design. Simulation results show that TOW &ndash F can control for population stratification in wide range of population structures including spatially structured populations, is robust to the directions of effect of causal variants, and is relatively robust to percentage of neutral variants. In GWA studies, this dissertation consists of a two &ndash locus joint effect analysis and a two-stage approach accounting for gene &ndash gene and gene &ndash environment interaction. Chapter 4 proposes a novel two &ndash stage approach, which is promising to identify joint effects, especially for monotonic models. The proposed approach outperforms a single &ndash marker method and a regular two &ndash stage analysis based on the two &ndash locus genotypic test. In Chapter 5, I proposed a gene &ndash based two &ndash stage approach to identify gene &ndash gene and gene &ndash environment interactions in GWA studies which can include rare variants. The two &ndash stage approach is applied to the GAW 17 dataset to identify the interaction between KDR gene and smoking status.

Association of retinoic acid receptor genes with meningomyelocele.

BACKGROUND: Neural tube defects (NTDs) occur in as many as 0.5-2 per 1000 live births in the United States. One of the most common and severe neural tube defects is meningomyelocele (MM) resulting from failed closure of the caudal end of the neural tube. MM has been induced by retinoic acid teratogenicity in rodent models. We hypothesized that genetic variants influencing retinoic acid (RA) induction via retinoic acid receptors (RARs) may be associated with risk for MM. METHODS: We analyzed 47 single nucleotide polymorphisms (SNPs) that span across the three retinoic acid receptor genes using the SNPlex genotyping platform. Our cohort consisted of 610 MM families. RESULTS: One variant in the RARA gene (rs12051734), three variants in the RARB gene (rs6799734, rs12630816, rs17016462), and a single variant in the RARG gene (rs3741434) were found to be statistically significant at p < 0.05. CONCLUSION: RAR genes were associated with risk for MM. For all associated SNPs, the rare allele conferred a protective effect for MM susceptibility.

Association of folate receptor (FOLR1, FOLR2, FOLR3) and reduced folate carrier (SLC19A1) genes with meningomyelocele.

BACKGROUND: Meningomyelocele (MM) results from lack of closure of the neural tube during embryologic development. Periconceptional folic acid supplementation is a modifier of MM risk in humans, leading toan interest in the folate transport genes as potential candidates for association to MM. METHODS: This study used the SNPlex Genotyping (ABI, Foster City, CA) platform to genotype 20 single polymorphic variants across the folate receptor genes (FOLR1, FOLR2, FOLR3) and the folate carrier gene (SLC19A1) to assess their association to MM. The study population included 329 trio and 281 duo families. Only cases with MM were included. Genetic association was assessed using the transmission disequilibrium test in PLINK. RESULTS: A variant in the FOLR2 gene (rs13908), three linked variants in the FOLR3 gene (rs7925545, rs7926875, rs7926987), and two variants in the SLC19A1 gene (rs1888530 and rs3788200) were statistically significant for association to MM in our population. CONCLUSION: This study involved the analyses of selected single nucleotide polymorphisms across the folate receptor genes and the folate carrier gene in a large population sample. It provided evidence that the rare alleles of specific single nucleotide polymorphisms within these genes appear to be statistically significant for association to MM in the patient population that was tested.

Genome-wide association studies based on sequence-derived genotypes reveal new QTL associated with conformation and performance traits in the Franches-Montagnes horse breed

To identify novel quantitative trait loci (QTL) within horses, we performed genome-wide association studies (GWAS) based on sequence-level genotypes for conformation and performance traits in the Franches-Montagnes (FM) horse breed. Sequence-level genotypes of FM horses were derived by re-sequencing 30 key founders and imputing 50K data of genotyped horses. In total, we included 1077 FM horses genotyped for ~4 million SNPs and their respective de-regressed breeding values of the traits in the analysis. Based on this dataset, we identified a total of 14 QTL associated with 18 conformation traits and one performance trait. Therefore, our results suggest that the application of sequence-derived genotypes increases the power to identify novel QTL which were not identified previously based on 50K SNP chip data.

The information bottleneck method for genome-wide association studies

In population studies, most current methods focus on identifying one outcome-related SNP at a time by testing for differences of genotype frequencies between disease and healthy groups or among different population groups. However, testing a great number of SNPs simultaneously has a problem of multiple testing and will give false-positive results. Although, this problem can be effectively dealt with through several approaches such as Bonferroni correction, permutation testing and false discovery rates, patterns of the joint effects by several genes, each with weak effect, might not be able to be determined. With the availability of high-throughput genotyping technology, searching for multiple scattered SNPs over the whole genome and modeling their joint effect on the target variable has become possible. Exhaustive search of all SNP subsets is computationally infeasible for millions of SNPs in a genome-wide study. Several effective feature selection methods combined with classification functions have been proposed to search for an optimal SNP subset among big data sets where the number of feature SNPs far exceeds the number of observations. ^ In this study, we take two steps to achieve the goal. First we selected 1000 SNPs through an effective filter method and then we performed a feature selection wrapped around a classifier to identify an optimal SNP subset for predicting disease. And also we developed a novel classification method-sequential information bottleneck method wrapped inside different search algorithms to identify an optimal subset of SNPs for classifying the outcome variable. This new method was compared with the classical linear discriminant analysis in terms of classification performance. Finally, we performed chi-square test to look at the relationship between each SNP and disease from another point of view. ^ In general, our results show that filtering features using harmononic mean of sensitivity and specificity(HMSS) through linear discriminant analysis (LDA) is better than using LDA training accuracy or mutual information in our study. Our results also demonstrate that exhaustive search of a small subset with one SNP, two SNPs or 3 SNP subset based on best 100 composite 2-SNPs can find an optimal subset and further inclusion of more SNPs through heuristic algorithm doesn't always increase the performance of SNP subsets. Although sequential forward floating selection can be applied to prevent from the nesting effect of forward selection, it does not always out-perform the latter due to overfitting from observing more complex subset states. ^ Our results also indicate that HMSS as a criterion to evaluate the classification ability of a function can be used in imbalanced data without modifying the original dataset as against classification accuracy. Our four studies suggest that Sequential Information Bottleneck(sIB), a new unsupervised technique, can be adopted to predict the outcome and its ability to detect the target status is superior to the traditional LDA in the study. ^ From our results we can see that the best test probability-HMSS for predicting CVD, stroke,CAD and psoriasis through sIB is 0.59406, 0.641815, 0.645315 and 0.678658, respectively. In terms of group prediction accuracy, the highest test accuracy of sIB for diagnosing a normal status among controls can reach 0.708999, 0.863216, 0.639918 and 0.850275 respectively in the four studies if the test accuracy among cases is required to be not less than 0.4. On the other hand, the highest test accuracy of sIB for diagnosing a disease among cases can reach 0.748644, 0.789916, 0.705701 and 0.749436 respectively in the four studies if the test accuracy among controls is required to be at least 0.4. ^ A further genome-wide association study through Chi square test shows that there are no significant SNPs detected at the cut-off level 9.09451E-08 in the Framingham heart study of CVD. Study results in WTCCC can only detect two significant SNPs that are associated with CAD. In the genome-wide study of psoriasis most of top 20 SNP markers with impressive classification accuracy are also significantly associated with the disease through chi-square test at the cut-off value 1.11E-07. ^ Although our classification methods can achieve high accuracy in the study, complete descriptions of those classification results(95% confidence interval or statistical test of differences) require more cost-effective methods or efficient computing system, both of which can't be accomplished currently in our genome-wide study. We should also note that the purpose of this study is to identify subsets of SNPs with high prediction ability and those SNPs with good discriminant power are not necessary to be causal markers for the disease.^

Genetic Predictors of Clinical Outcomes in Non-Small Cell Lung Cancer Patients

Lung cancer is the leading cause of cancer-related mortality in the US. Emerging evidence has shown that host genetic factors can interact with environmental exposures to influence patient susceptibility to the diseases as well as clinical outcomes, such as survival and recurrence. We aimed to identify genetic prognostic markers for non-small cell lung cancer (NSCLC), a major (85%) subtype of lung cancer, and also in other subgroups. With the fast evolution of genotyping technology, genetic association studies have went through candidate gene approach, to pathway-based approach, to the genome wide association study (GWAS). Even in the era of GWAS, pathway-based approach has its own advantages on studying cancer clinical outcomes: it is cost-effective, requiring a smaller sample size than GWAS easier to identify a validation population and explore gene-gene interactions. In the current study, we adopted pathway-based approach focusing on two critical pathways - miRNA and inflammation pathways. MicroRNAs (miRNA) post-transcriptionally regulate around 30% of human genes. Polymorphisms within miRNA processing pathways and binding sites may influence patients’ prognosis through altered gene regulation. Inflammation plays an important role in cancer initiation and progression, and also has shown to impact patients’ clinical outcomes. We first evaluated 240 single nucleotide polymorphisms (SNPs) in miRNA biogenesis genes and predicted binding sites in NSCLC patients to determine associations with clinical outcomes in early-stage (stage I and II) and late-stage (stage III and IV) lung cancer patients, respectively. First, in 535 early-stage patients, after correcting multiple comparisons, FZD4:rs713065 (hazard ratio [HR]:0.46, 95% confidence interval [CI]:0.32-0.65) showed a significant inverse association with survival in early stage surgery-only patients. SP1:rs17695156 (HR:2.22, 95% CI:1.44-3.41) and DROSHA:rs6886834 (HR:6.38, 95% CI:2.49-16.31) conferred increased risk of progression in the all patients and surgery-only populations, respectively. FAS:rs2234978 was significantly associated with improved survival in all patients (HR:0.59, 95% CI:0.44-0.77) and in the surgery plus chemotherapy populations (HR:0.19, 95% CI:0.07-0.46).. Functional genomics analysis demonstrated that this variant creates a miR-651 binding site resulting in altered miRNA regulation of FAS, providing biological plausibility for the observed association. We then analyzed these associations in 598 late-stage patients. After multiple comparison corrections, no SNPs remained significant in the late stage group, while the top SNP NAT1:rs15561 (HR=1.98, 96%CI=1.32-2.94) conferred a significantly increased risk of death in the chemotherapy subgroup. To test the hypothesis that genetic variants in the inflammation-related pathways may be associated with survival in NSCLC patients, we first conducted a three-stage study. In the discovery phase, we investigated a comprehensive panel of 11,930 inflammation-related SNPs in three independent lung cancer populations. A missense SNP (rs2071554) in HLA-DOB was significantly associated with poor survival in the discovery population (HR: 1.46, 95% CI: 1.02-2.09), internal validation population (HR: 1.51, 95% CI: 1.02-2.25), and external validation (HR: 1.52, 95% CI: 1.01-2.29) population. Rs2900420 in KLRK1 was significantly associated with a reduced risk for death in the discovery (HR: 0.76, 95% CI: 0.60-0.96) and internal validation (HR: 0.77, 95% CI: 0.61-0.99) populations, and the association reached borderline significance in the external validation population (HR: 0.80, 95% CI: 0.63-1.02). We also evaluated these inflammation-related SNPs in NSCLC patients in never smokers. Lung cancer in never smokers has been increasingly recognized as distinct disease from that in ever-smokers. A two-stage study was performed using a discovery population from MD Anderson (411 patients) and a validation population from Mayo Clinic (311 patients). Three SNPs (IL17RA:rs879576, BMP8A:rs698141, and STK:rs290229) that were significantly associated with survival were validated (pCD74:rs1056400 and CD38:rs10805347) were borderline significant (p=0.08) in the Mayo Clinic population. In the combined analysis, IL17RA:rs879576 resulted in a 40% reduction in the risk for death (p=4.1 × 10-5 [p=0.61, heterogeneity test]). We also validated a survival tree created in MD Anderson population in the Mayo Clinic population. In conclusion, our results provided strong evidence that genetic variations in specific pathways that examined (miRNA and inflammation pathways) influenced clinical outcomes in NSCLC patients, and with further functional studies, the novel loci have potential to be translated into clinical use.

A Validation Study of Vascular Cognitive Impairment Genetics Meta-Analysis Findings in an Independent Collaborative Cohort

Vascular cognitive impairment (VCI), including its severe form, vascular dementia (VaD), is the second most common form of dementia. The genetic etiology of sporadic VCI remains largely unknown. We previously conducted a systematic review and meta-analysis of all published genetic association studies of sporadic VCI prior to 6 July 2012, which demonstrated that APOE (ɛ4, ɛ2) and MTHFR (rs1801133) variants were associated with susceptibility for VCI. De novo genotyping was conducted in a new independent relatively large collaborative European cohort of VaD (nmax = 549) and elderly non-demented samples (nmax = 552). Where available, genotype data derived from Illumina's 610-quad array for 1210 GERAD1 control samples were also included in analyses of genes examined. Associations were tested using the Cochran-Armitage trend test: MTHFR rs1801133 (OR = 1.36, 95% CI 1.16-1.58, p = <0.0001), APOE rs7412 (OR = 0.62, 95% CI 0.42-0.90, p = 0.01), and APOE rs429358 (OR = 1.59, 95% CI 1.17-2.16, p = 0.003). Association was also observed with APOE epsilon alleles; ɛ4 (OR = 1.85, 95% CI 1.35-2.52, p = <0.0001) and ɛ2 (OR = 0.67, 95% CI 0.46-0.98, p = 0.03). Logistic Regression and Bonferroni correction in a subgroup of the cohort adjusted for gender, age, and population maintained the association of APOE rs429358 and ɛ4 allele.

Large national series of patients with Xq28 duplication involving MECP2: Delineation of brain MRI abnormalities in 30 affected patients.

Xq28 duplications encompassing MECP2 have been described in male patients with a severe neurodevelopmental disorder associated with hypotonia and spasticity, severe learning disability, stereotyped movements, and recurrent pulmonary infections. We report on standardized brain magnetic resonance imaging (MRI) data of 30 affected patients carrying an Xq28 duplication involving MECP2 of various sizes (228 kb to 11.7 Mb). The aim of this study was to seek recurrent malformations and attempt to determine whether variations in imaging features could be explained by differences in the size of the duplications. We showed that 93% of patients had brain MRI abnormalities such as corpus callosum abnormalities (n = 20), reduced volume of the white matter (WM) (n = 12), ventricular dilatation (n = 9), abnormal increased hyperintensities on T2-weighted images involving posterior periventricular WM (n = 6), and vermis hypoplasia (n = 5). The occipitofrontal circumference varied considerably between >+2SD in five patients and <-2SD in four patients. Among the nine patients with dilatation of the lateral ventricles, six had a duplication involving L1CAM. The only patient harboring bilateral posterior subependymal nodular heterotopia also carried an FLNA gene duplication. We could not demonstrate a correlation between periventricular WM hyperintensities/delayed myelination and duplication of the IKBKG gene. We thus conclude that patients with an Xq28 duplication involving MECP2 share some similar but non-specific brain abnormalities. These imaging features, therefore, could not constitute a diagnostic clue. The genotype-phenotype correlation failed to demonstrate a relationship between the presence of nodular heterotopia, ventricular dilatation, WM abnormalities, and the presence of FLNA, L1CAM, or IKBKG, respectively, in the duplicated segment.

Partial agonists of the [alpha]3[beta]4[ast] neuronal nicotinic acetylcholine receptor reduce ethanol consumption and seeking in rats

Alcohol use disorders (AUDs) impact millions of individuals and there remain few effective treatment strategies. Despite evidence that neuronal nicotinic acetylcholine receptors (nAChRs) have a role in AUDs, it has not been established which subtypes of the nAChR are involved. Recent human genetic association studies have implicated the gene cluster CHRNA3-CHRNA5-CHRNB4 encoding the α3, α5, and β4 subunits of the nAChR in susceptibility to develop nicotine and alcohol dependence; however, their role in ethanol-mediated behaviors is unknown due to the lack of suitable and selective research tools. To determine the role of the α3, and β4 subunits of the nAChR in ethanol self-administration, we developed and characterized high-affinity partial agonists at α3β4 nAChRs, CP-601932, and PF-4575180. Both CP-601932 and PF-4575180 selectively decrease ethanol but not sucrose consumption and operant self-administration following long-term exposure. We show that the functional potencies of CP-601932 and PF-4575180 at α3β4 nAChRs correlate with their unbound rat brain concentrations, suggesting that the effects on ethanol self-administration are mediated via interaction with α3β4 nAChRs. Also varenicline, an approved smoking cessation aid previously shown to decrease ethanol consumption and seeking in rats and mice, reduces ethanol intake at unbound brain concentrations that allow functional interactions with α3β4 nAChRs. Furthermore, the selective α4β2(*) nAChR antagonist, DHβE, did not reduce ethanol intake. Together, these data provide further support for the human genetic association studies, implicating CHRNA3 and CHRNB4 genes in ethanol-mediated behaviors. CP-601932 has been shown to be safe in humans and may represent a potential novel treatment for AUDs.

A pharmacogenomic evaluation of migraine therapy

Migraine is a common idiopathic primary headache disorder with significant mental, physical and social health implications. Accompanying an intense unilateral pulsating head pain other characteristic migraine symptoms include nausea, emesis, phonophobia, photophobia and in approximately 20-30% of migraine cases, neurologic disturbances associated with the aura phase. Although selective serotonin (5-HT) receptor agonists (i.e., 5-HT(1B/1D)) are successful in alleviating migrainous symptoms in < or = 70% of known sufferers, for the remaining 30%, additional migraine abortive medications remain unsuccessful, not tested or yet to be identified. Genetic characterization of the migrainous disorder is making steady progress with an increasing number of genomic susceptibility loci now identified on chromosomes 1q, 4q, 5q, 6p, 11q, 14q, 15q, 17p, 18q, 19p and Xq. The 4q, 5q, 17p and 18q loci involve endophenotypic susceptibility regions for various migrainous symptoms. In an effort to develop individualized pharmacotherapeutics, the identification of these migraine endophenotypic loci may well be the catalyst needed to aid in this goal. In this review the authors discuss the present treatment of migraine, known genomic susceptibility regions and results from migraine (genetic) association studies. The authors also discuss pharmacogenomic considerations for more individualized migraine prophylactic treatments.

The α5 subunit regulates the expression and function of α4*-containing neuronal nicotinic acetylcholine receptors in the ventral-tegmental area

Human genetic association studies have shown gene variants in the α5 subunit of the neuronal nicotinic receptor (nAChR) influence both ethanol and nicotine dependence. The α5 subunit is an accessory subunit that facilitates α4* nAChRs assembly in vitro. However, it is unknown whether this occurs in the brain, as there are few research tools to adequately address this question. As the α4*-containing nAChRs are highly expressed in the ventral tegmental area (VTA) we assessed the molecular, functional and pharmacological roles of α5 in α4*-containing nAChRs in the VTA. We utilized transgenic mice α5+/+(α4YFP) and α5-/-(α4YFP) that allow the direct visualization and measurement of α4-YFP expression and the effect of the presence (α5+/+) and absence of α5 (-/-) subunit, as the antibodies for detecting the α4* subunits of the nAChR are not specific. We performed voltage clamp electrophysiological experiments to study baseline nicotinic currents in VTA dopaminergic neurons. We show that in the presence of the α5 subunit, the overall expression of α4 subunit is increased significantly by 60% in the VTA. Furthermore, the α5 subunit strengthens baseline nAChR currents, suggesting the increased expression of α4* nAChRs to be likely on the cell surface. While the presence of the α5 subunit blunts the desensitization of nAChRs following nicotine exposure, it does not alter the amount of ethanol potentiation of VTA dopaminergic neurons. Our data demonstrates a major regulatory role for the α5 subunit in both the maintenance of α4*-containing nAChRs expression and in modulating nicotinic currents in VTA dopaminergic neurons. Additionally, the α5α4* nAChR in VTA dopaminergic neurons regulates the effect of nicotine but not ethanol on currents. Together, the data suggest that the α5 subunit is critical for controlling the expression and functional role of a population of α4*-containing nAChRs in the VTA.

Genetics of ankylosing spondylitis

Purpose of Review Over the past 3 years, several new genes and gene deserts have been identified that are associated with ankylosing spondylitis (AS). The purpose of this review is to discuss the major findings of these studies, and the answers they provide and questions they raise about the pathogenesis of this common condition. Recent Findings: Five genes/genetic regions have now definitively been associated with AS [the major histocompatibility complex (MHC), IL23R, ERAP1, 2p15 and 21q22]. Strong evidence to support association with the disease has been demonstrated for the genes IL1R2, ANTXR2, TNFSF15, TNFR1 and a region on chromosome 16q including the gene TRADD. There is an overrepresentation of genes involved in Th17 lymphocyte differentiation/activation among genes associated with AS and the related diseases inflammatory bowel disease and psoriasis, pointing strongly to this pathway as playing a major causative role in the disease. Increasing information about differential association of HLA-B27 subtypes with disease suggests a hierarchy of strength of association of those alleles with AS, providing a useful test as to the validity of different potential mechanisms of association of HLA-B27 with AS. The mechanism underlying the association of the gene deserts, 2p15 and 21q22, suggests the involvement of noncoding RNA in AS etiopathogenesis. Summary: The increasing list of genes identified as being definitely involved in AS provides a useful platform for hypothesis-driven research in the field, providing a potential alternative route to determining the underlying mechanisms involved in the disease to research focusing on HLA-B27 alone.