Mutation and association analysis of the PVR and PVRL2 genes in patients with non-syndromic cleft lip and palate.

Orofacial clefts (OFC; MIM 119530) are among the most common major birth defects. Here, we carried out mutation screening of the PVR and PVRL2 genes, which are both located at an OFC linkage region at 19q13 (OFC3) and are closely related to PVRL1, which has been associated with both syndromic and non-syndromic cleft lip and palate (nsCLP). We screened a total of 73 nsCLP patients and 105 non-cleft controls from the USA for variants in PVR and PVRL2, including all exons and encompassing all isoforms. We identified four variants in PVR and five in PVRL2. One non-synonymous PVR variant, A67T, was more frequent among nsCLP patients than among normal controls, but this difference did not achieve statistical significance.

A Multi-Breed Genome-Wide Association Analysis for Canine Hypothyroidism Identifies a Shared Major Risk Locus on CFA12.

Hypothyroidism is a complex clinical condition found in both humans and dogs, thought to be caused by a combination of genetic and environmental factors. In this study we present a multi-breed analysis of predisposing genetic risk factors for hypothyroidism in dogs using three high-risk breeds-the Gordon Setter, Hovawart and the Rhodesian Ridgeback. Using a genome-wide association approach and meta-analysis, we identified a major hypothyroidism risk locus shared by these breeds on chromosome 12 (p = 2.1x10-11). Further characterisation of the candidate region revealed a shared ~167 kb risk haplotype (4,915,018-5,081,823 bp), tagged by two SNPs in almost complete linkage disequilibrium. This breed-shared risk haplotype includes three genes (LHFPL5, SRPK1 and SLC26A8) and does not extend to the dog leukocyte antigen (DLA) class II gene cluster located in the vicinity. These three genes have not been identified as candidate genes for hypothyroid disease previously, but have functions that could potentially contribute to the development of the disease. Our results implicate the potential involvement of novel genes and pathways for the development of canine hypothyroidism, raising new possibilities for screening, breeding programmes and treatments in dogs. This study may also contribute to our understanding of the genetic etiology of human hypothyroid disease, which is one of the most common endocrine disorders in humans.

Comparison of two statistical methods for rare variant association analysis

Studies have shown that rare genetic variants have stronger effects in predisposing common diseases, and several statistical methods have been developed for association studies involving rare variants. In order to better understand how these statistical methods perform, we seek to compare two recently developed rare variant statistical methods (VT and C-alpha) on 10,000 simulated re-sequencing data sets with disease status and the corresponding 10,000 simulated null data sets. The SLC1A1 gene has been suggested to be associated with diastolic blood pressure (DBP) in previous studies. In the current study, we applied VT and C-alpha methods to the empirical re-sequencing data for the SLC1A1 gene from 300 whites and 200 blacks. We found that VT method obtains higher power and performs better than C-alpha method with the simulated data we used. The type I errors were well-controlled for both methods. In addition, both VT and C-alpha methods suggested no statistical evidence for the association between the SLC1A1 gene and DBP. Overall, our findings provided an important comparison of the two statistical methods for future reference and provided preliminary and pioneer findings on the association between the SLC1A1 gene and blood pressure.^

Nonlinear dimension reduction in pathway-based genome-wide association analysis

Pathway based genome wide association study evolves from pathway analysis for microarray gene expression and is under rapid development as a complementary for single-SNP based genome wide association study. However, it faces new challenges, such as the summarization of SNP statistics to pathway statistics. The current study applies the ridge regularized Kernel Sliced Inverse Regression (KSIR) to achieve dimension reduction and compared this method to the other two widely used methods, the minimal-p-value (minP) approach of assigning the best test statistics of all SNPs in each pathway as the statistics of the pathway and the principal component analysis (PCA) method of utilizing PCA to calculate the principal components of each pathway. Comparison of the three methods using simulated datasets consisting of 500 cases, 500 controls and100 SNPs demonstrated that KSIR method outperformed the other two methods in terms of causal pathway ranking and the statistical power. PCA method showed similar performance as the minP method. KSIR method also showed a better performance over the other two methods in analyzing a real dataset, the WTCCC Ulcerative Colitis dataset consisting of 1762 cases, 3773 controls as the discovery cohort and 591 cases, 1639 controls as the replication cohort. Several immune and non-immune pathways relevant to ulcerative colitis were identified by these methods. Results from the current study provided a reference for further methodology development and identified novel pathways that may be of importance to the development of ulcerative colitis.^

Cladistic association analysis of Y chromosome effects on alcohol dependence and related personality traits

Association between Y chromosome haplotype variation and alcohol dependence and related personality traits was investigated in a large sample of psychiatrically diagnosed Finnish males. Haplotypes were constructed for 359 individuals using alleles at eight loci (seven microsatellite loci and a nucleotide substitution in the DYZ3 alphoid satellite locus). A cladogram linking the 102 observed haplotype configurations was constructed by using parsimony with a single-step mutation model. Then, a series of contingency tables nested according to the cladogram hierarchy were used to test for association between Y haplotype and alcohol dependence. Finally, using only alcohol-dependent subjects, we tested for association between Y haplotype and personality variables postulated to define subtypes of alcoholism—antisocial personality disorder, novelty seeking, harm avoidance, and reward dependence. Significant association with alcohol dependence was observed at three Y haplotype clades, with significance levels of P = 0.002, P = 0.020, and P = 0.010. Within alcohol-dependent subjects, no relationship was revealed between Y haplotype and antisocial personality disorder, novelty seeking, harm avoidance, or reward dependence. These results demonstrate, by using a fully objective association design, that differences among Y chromosomes contribute to variation in vulnerability to alcohol dependence. However, they do not demonstrate an association between Y haplotype and the personality variables thought to underlie the subtypes of alcoholism.

Effects of dopamine receptor D4 variation on alcohol and tobacco use and on novelty seeking: Multivariate linkage and association analysis

The dopamine D4 receptor gene contains a polymorphic sequence consisting of a variable number of 48-base-pair (bp) repeats, and there have been a number of reports that this polymorphism is associated with variation in novelty seeking or in substance abuse and addictive behaviors. In this study we have assessed the linkage and association of DRD4 genotype with novelty seeking, alcohol use, and smoking in a sample of 377 dizygotic twin pairs and 15 single twins recruited from the Australian Twin Registry (ATR). We found no evidence of linkage or association of the DRD4 locus with any of the phenotypes. We made use of repeated measures for some phenotypes to increase power by multivariate genetic analysis, but allelic effects were still non-significant. Specifically, it has been suggested that the DRD4 7-repeat allele is associated with increased novelty seeking in males but we found no evidence for this, despite considerable power to do so. We conclude that DRD4 variation does not have an effect on use of alcohol and the problems that arise from it, on smoking, or on novelty seeking behavior. (C) 2003 Wiley-Liss, Inc.

Genome-wide association analysis identifies novel loci for chronotype in 100,420 individuals from the UK Biobank

Our sleep timing preference, or chronotype, is a manifestation of our internal biological clock. Variation in chronotype has been linked to sleep disorders, cognitive and physical performance, and chronic disease. Here we perform a genome-wide association study of self-reported chronotype within the UK Biobank cohort (n=100,420). We identify 12 new genetic loci that implicate known components of the circadian clock machinery and point to previously unstudied genetic variants and candidate genes that might modulate core circadian rhythms or light-sensing pathways. Pathway analyses highlight central nervous and ocular systems and fear-response-related processes. Genetic correlation analysis suggests chronotype shares underlying genetic pathways with schizophrenia, educational attainment and possibly BMI. Further, Mendelian randomization suggests that evening chronotype relates to higher educational attainment. These results not only expand our knowledge of the circadian system in humans but also expose the influence of circadian characteristics over human health and life-history variables such as educational attainment.

Flexible strategies for association analysis with genomic phenotypes

Thesis (Ph.D.)--University of Washington, 2016-08

Genome-wide analysis in German shepherd dogs reveals association of a locus on CFA 27 with atopic dermatitis

Humans and dogs are both affected by the allergic skin disease atopic dermatitis (AD), caused by an interaction between genetic and environmental factors. The German shepherd dog (GSD) is a high-risk breed for canine AD (CAD). In this study, we used a Swedish cohort of GSDs as a model for human AD. Serum IgA levels are known to be lower in GSDs compared to other breeds. We detected significantly lower IgA levels in the CAD cases compared to controls (p = 1.1 × 10(-5)) in our study population. We also detected a separation within the GSD cohort, where dogs could be grouped into two different subpopulations. Disease prevalence differed significantly between the subpopulations contributing to population stratification (λ = 1.3), which was successfully corrected for using a mixed model approach. A genome-wide association analysis of CAD was performed (n cases = 91, n controls = 88). IgA levels were included in the model, due to the high correlation between CAD and low IgA levels. In addition, we detected a correlation between IgA levels and the age at the time of sampling (corr = 0.42, p = 3.0 × 10(-9)), thus age was included in the model. A genome-wide significant association was detected on chromosome 27 (praw = 3.1 × 10(-7), pgenome = 0.03). The total associated region was defined as a ~1.5-Mb-long haplotype including eight genes. Through targeted re-sequencing and additional genotyping of a subset of identified SNPs, we defined 11 smaller haplotype blocks within the associated region. Two blocks showed the strongest association to CAD. The ~209-kb region, defined by the two blocks, harbors only the PKP2 gene, encoding Plakophilin 2 expressed in the desmosomes and important for skin structure. Our results may yield further insight into the genetics behind both canine and human AD.

Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

Genome-wide association studies in oesophageal adenocarcinoma and Barrett's oesophagus: a large-scale meta-analysis

Background: Esophageal adenocarcinoma (EA) is one of the fastest rising cancers in western countries. Barrett’s Esophagus (BE) is the premalignant precursor of EA. However, only a subset of BE patients develop EA, which complicates the clinical management in the absence of valid predictors. Genetic risk factors for BE and EA are incompletely understood. This study aimed to identify novel genetic risk factors for BE and EA.Methods: Within an international consortium of groups involved in the genetics of BE/EA, we performed the first meta-analysis of all genome-wide association studies (GWAS) available, involving 6,167 BE patients, 4,112 EA patients, and 17,159 representative controls, all of European ancestry, genotyped on Illumina high-density SNP-arrays, collected from four separate studies within North America, Europe, and Australia. Meta-analysis was conducted using the fixed-effects inverse variance-weighting approach. We used the standard genome-wide significant threshold of 5×10-8 for this study. We also conducted an association analysis following reweighting of loci using an approach that investigates annotation enrichment among the genome-wide significant loci. The entire GWAS-data set was also analyzed using bioinformatics approaches including functional annotation databases as well as gene-based and pathway-based methods in order to identify pathophysiologically relevant cellular pathways.Findings: We identified eight new associated risk loci for BE and EA, within or near the CFTR (rs17451754, P=4·8×10-10), MSRA (rs17749155, P=5·2×10-10), BLK (rs10108511, P=2·1×10-9), KHDRBS2 (rs62423175, P=3·0×10-9), TPPP/CEP72 (rs9918259, P=3·2×10-9), TMOD1 (rs7852462, P=1·5×10-8), SATB2 (rs139606545, P=2·0×10-8), and HTR3C/ABCC5 genes (rs9823696, P=1·6×10-8). A further novel risk locus at LPA (rs12207195, posteriori probability=0·925) was identified after re-weighting using significantly enriched annotations. This study thereby doubled the number of known risk loci. The strongest disease pathways identified (P<10-6) belong to muscle cell differentiation and to mesenchyme development/differentiation, which fit with current pathophysiological BE/EA concepts. To our knowledge, this study identified for the first time an EA-specific association (rs9823696, P=1·6×10-8) near HTR3C/ABCC5 which is independent of BE development (P=0·45).Interpretation: The identified disease loci and pathways reveal new insights into the etiology of BE and EA. Furthermore, the EA-specific association at HTR3C/ABCC5 may constitute a novel genetic marker for the prediction of transition from BE to EA. Mutations in CFTR, one of the new risk loci identified in this study, cause cystic fibrosis (CF), the most common recessive disorder in Europeans. Gastroesophageal reflux (GER) belongs to the phenotypic CF-spectrum and represents the main risk factor for BE/EA. Thus, the CFTR locus may trigger a common GER-mediated pathophysiology.

Heritability and genome-wide linkage analysis of migraine in the genetic isolate of Norfolk Island

Migraine is a common neurovascular disorder with a complex envirogenomic aetiology. In an effort to identify migraine susceptibility genes, we conducted a study of the isolated population of Norfolk Island, Australia. A large portion of the permanent inhabitants of Norfolk Island are descended from 18th Century English sailors involved in the infamous mutiny on the Bounty and their Polynesian consorts. In total, 600 subjects were recruited including a large pedigree of 377 individuals with lineage to the founders. All individuals were phenotyped for migraine using International Classification of Headache Disorders-II criterion. All subjects were genotyped for a genome-wide panel of microsatellite markers. Genotype and phenotype data for the pedigree were analysed using heritability and linkage methods implemented in the programme SOLAR. Follow-up association analysis was performed using the CLUMP programme. A total of 154 migraine cases (25%) were identified indicating the Norfolk Island population is high-risk for migraine. Heritability estimation of the 377-member pedigree indicated a significant genetic component for migraine (h2 = 0.53, P = 0.016). Linkage analysis showed peaks on chromosome 13q33.1 (P = 0.003) and chromosome 9q22.32 (P = 0.008). Association analysis of the key microsatellites in the remaining 223 unrelated Norfolk Island individuals showed evidence of association, which strengthen support for the linkage findings (P ≤ 0.05). In conclusion, a genome-wide linkage analysis and follow-up association analysis of migraine in the genetic isolate of Norfolk Island provided evidence for migraine susceptibility loci on chromosomes 9q22.22 and 13q33.1.

Association between a 19 bp deletion polymorphism at the dopamine beta-hydroxylase (DBH) locus and migraine with aura

Migraine is a debilitating neurological disorder, affecting 12% of Caucasian populations. It is well known that migraine has a strong genetic component, although the type and number of genes involved is unclear. Our previous work has investigated dopamine related migraine candidate genes and has reported a significant allelic association with migraine of a microsatellite localised to the promoter region of the dopamine beta-hydroxylase (DBH) gene. The present study performed an association analysis in a larger population of case-controls (275 unrelated Caucasian migraineurs versus 275 controls) examining two different genetic DBH polymorphisms (a functional insertion/deletion promoter and a coding SNP A444G polymorphism). Although no significant association was found for the SNP polymorphism, the results showed a significant association between the insertion/deletion variant and disease (chi(2)=8.92, P=0.011), in particular in migraine with aura (chi(2)=11.53, P=0.003) compared to the control group. Furthermore, the analysis of this polymorphism stratified by gender, revealed that male individuals with the homozygote deletion genotype had three times the risk of developing migraine, compared to females. The DBH insertion/deletion polymorphism is in linkage disequilibrium with the previously reported migraine associated DBH microsatellite and this insertion/deletion polymorphism is functional, which may explain a potential role in susceptibility to migraine.

Association and linkage analyses of restriction fragment length polymorphisms for the human renin and antithrombin III genes in essential hypertension

Essential hypertension is a highly hereditable disorder in which genetic influences predominate over environmental factors. The molecular genetic profiles which predispose to essential hypertension are not known. In rats with genetic hypertension, there is some recent evidence pointing to linkage of renin gene alleles with blood pressure. The genes for renin and antithrombin III belong to a conserved synteny group which, in humans, spans the q21.3-32.3 region of chromosome I and, in rats, is linkage group X on chromosome 13. The present study examined the association of particular human renin gene (REN) and antithrombin III gene (AT3) polymorphisms with essential hypertension by comparing the frequency of specific alleles for each of these genes in 50 hypertensive offspring of hypertensive parents and 91 normotensive offspring of normotensive parents. In addition, linkage relationships were examined in hypertensive pedigrees with multiple affected individuals. Alleles of a REN HindIII restriction fragment length polymorphism (RFLP) were detected using a genomic clone, λHR5, to probe Southern blots of HindIII-cut leucocyte DNA, and those for an AT3 Pstl RFLP were detected by phATIII 113 complementary DNA probe. The frequencies of each REN allele in the hypertensive group were 0.76 and 0.24 compared with 0.74 and 0.26 in the normotensive group. For AT3, hypertensive allele frequencies were 0.49 and 0.51 compared with normotensive values of 0.54 and 0.46. These differences were not significant by χ2 analysis (P > 0.2). Linkage analysis of a family (data from 16 family members, 10 of whom were hypertensive), informative for both markers, without an age-of-onset correction, and assuming dominant inheritance of hypertension, complete penetrance and a disease frequency of 20%, did not indicate linkage of REN with hypertension, but gave a positive, although not significant, logarithm of the odds for linkage score of 0.784 at a recombination fraction of 0 for AT3 linkage to hypertension. In conclusion, the present study could find no evidence for an association of a REN HindIII RFLP with essential hypertension or for a linkage of the locus defined by this RFLP in a family segregating for hypertension. In the case of an AT3 Pstl RFLP, although association analysis was negative, linkage analysis suggested possible involvement (odds of 6:1 in favour) of a gene located near the 1q23 locus with hypertension in one informative family.

Evidence for allelic association of the dopamine β-hydroxylase gene (DBH) with susceptibility to typical migraine

Migraine is a debilitating neurological disorder characterized by recurrent attacks of severe headache. The disorder is highly prevalent, affecting approximately 12% of Caucasian populations. It is well known that migraine has a strong genetic component, although the type and number of genes involved is not yet clear. However, the calcium channel gene, CACNA1A, on chromosome 19 contains mutations responsible for familial hemiplegic migraine, a rare and severe subtype of migraine. There is also evidence to suggest that serotonin- and dopamine-related genes may be involved in the pathogenesis of migraine. This study employed a linkage and association approach to investigate neurotransmitter-related migraine candidate genes. Polymorphisms within the dopamine beta-hydroxylase (DBH) gene, serotonin transporter gene (SERT), and dopamine receptor gene (DRD2) were tested in 177 unrelated Caucasian migraineurs and 182 control individuals. In addition, an independent sample of 82 families affected with migraine was examined. Unrelated case-control association analysis of a DBH intragenic dinucleotide polymorphism indicated altered allelic distribution between migraine and control groups (L2=16.53, P=0.019). Furthermore, the transmission/disequilibrium test, which was implemented on the family data, also indicated distortion of allele transmission for the same DBH marker (L2=4.44, P=0.035). Together, these results provide evidence for allelic association of the DBH gene with typical migraine susceptibility (Fisher's combined P value=0.006) and indicate that further research into the role of the DBH gene in the etiology of migraine is warranted.