Common variants in left/right asymmetry genes and pathways are associated with relative hand skill

Humans display structural and functional asymmetries in brain organization, strikingly with respect to language and handedness. The molecular basis of these asymmetries is unknown. We report a genome-wide association study meta-analysis for a quantitative measure of relative hand skill in individuals with dyslexia [reading disability (RD)] (n = 728). The most strongly associated variant, rs7182874 (P = 8.68×10-9), is located in PCSK6, further supporting an association we previously reported. We also confirmed the specificity of this association in individuals with RD; the same locus was not associated with relative hand skill in a general population cohort (n = 2,666). As PCSK6 is known to regulate NODAL in the development of left/right (LR) asymmetry in mice, we developed a novel approach to GWAS pathway analysis, using gene-set enrichment to test for an over-representation of highly associated variants within the orthologs of genes whose disruption in mice yields LR asymmetry phenotypes. Four out of 15 LR asymmetry phenotypes showed an over-representation (FDR≤5%). We replicated three of these phenotypes; situs inversus, heterotaxia, and double outlet right ventricle, in the general population cohort (FDR≤5%). Our findings lead us to propose that handedness is a polygenic trait controlled in part by the molecular mechanisms that establish LR body asymmetry early in development. © 2013 Brandler et al.

Entropy Based Approach to Finding Interacting Genes Responsible for Complex Human Disease

2000 Mathematics Subject Classification: 62P10, 92D10, 92D30, 94A17, 62L10.

The handedness-associated PCSK6 locus spans an intronic promoter regulating novel transcripts

We recently reported the association of the PCSK6 gene with handedness through a quantitative genome-wide association study (GWAS; P < 0.5 × 10(-8)) for a relative hand skill measure in individuals with dyslexia. PCSK6 activates Nodal, a morphogen involved in regulating left-right body axis determination. Therefore, the GWAS data suggest that the biology underlying the patterning of structural asymmetries may also contribute to behavioural laterality, e.g. handedness. The association is further supported by an independent study reporting a variable number tandem repeat (VNTR) within the same PCSK6 locus to be associated with degree of handedness in a general population cohort. Here, we have conducted a functional analysis of the PCSK6 locus combining further genetic analysis, in silico predictions and molecular assays. We have shown that the previous GWAS signal was not tagging a VNTR effect, suggesting that the two markers have independent effects. We demonstrated experimentally that one of the top GWAS-associated markers, rs11855145, directly alters the binding site for a nuclear factor. Furthermore, we have shown that the predicted regulatory region adjacent to rs11855415 acts as a bidirectional promoter controlling the expression of novel RNA transcripts. These include both an antisense long non-coding RNA (lncRNA) and a short PCSK6 isoform predicted to be coding. This is the first molecular characterization of a handedness-associated locus that supports the role of common variants in non-coding sequences in influencing complex phenotypes through gene expression regulation.

Further evidence for a parent-of-origin effect at the NOP9 locus on language-related phenotypes

Background - Specific language impairment (SLI) is a common neurodevelopmental disorder, observed in 5–10 % of children. Family and twin studies suggest a strong genetic component, but relatively few candidate genes have been reported to date. A recent genome-wide association study (GWAS) described the first statistically significant association specifically for a SLI cohort between a missense variant (rs4280164) in the NOP9 gene and language-related phenotypes under a parent-of-origin model. Replications of these findings are particularly challenging because the availability of parental DNA is required. Methods - We used two independent family-based cohorts characterised with reading- and language-related traits: a longitudinal cohort (n = 106 informative families) including children with language and reading difficulties and a nuclear family cohort (n = 264 families) selected for dyslexia. Results - We observed association with language-related measures when modelling for parent-of-origin effects at the NOP9 locus in both cohorts: minimum P = 0.001 for phonological awareness with a paternal effect in the first cohort and minimum P = 0.0004 for irregular word reading with a maternal effect in the second cohort. Allelic and parental trends were not consistent when compared to the original study. Conclusions - A parent-of-origin effect at this locus was detected in both cohorts, albeit with different trends. These findings contribute in interpreting the original GWAS report and support further investigations of the NOP9 locus and its role in language-related traits. A systematic evaluation of parent-of-origin effects in genetic association studies has the potential to reveal novel mechanisms underlying complex traits.

Functional Evaluation of Causal Mutations Identified in Human Genetic Studies

Human genetics has been experiencing a wave of genetic discoveries thanks to the development of several technologies, such as genome-wide association studies (GWAS), whole-exome sequencing, and whole genome sequencing. Despite the massive genetic discoveries of new variants associated with human diseases, several key challenges emerge following the genetic discovery. GWAS is known to be good at identifying the locus associated with the patient phenotype. However, the actually causal variants responsible for the phenotype are often elusive. Another challenge in human genetics is that even the causal mutations are already known, the underlying biological effect might remain largely ambiguous. Functional evaluation plays a key role to solve these key challenges in human genetics both to identify causal variants responsible for the phenotype, and to further develop the biological insights from the disease-causing mutations.

We adopted various methods to characterize the effects of variants identified in human genetic studies, including patient genetic and phenotypic data, RNA chemistry, molecular biology, virology, and multi-electrode array and primary neuronal culture systems. Chapter 1 is a broader introduction for the motivation and challenges for functional evaluation in human genetic studies, and the background of several genetics discoveries, such as hepatitis C treatment response, in which we performed functional characterization.

Chapter 2 focuses on the characterization of causal variants following the GWAS study for hepatitis C treatment response. We characterized a non-coding SNP (rs4803217) of IL28B (IFNL3) in high linkage disequilibrium (LD) with the discovery SNP identified in the GWAS. In this chapter, we used inter-disciplinary approaches to characterize rs4803217 on RNA structure, disease association, and protein translation.

Chapter 3 describes another avenue of functional characterization following GWAS focusing on the novel transcripts and proteins identified near the IL28B (IFNL3) locus. It has been recently speculated that this novel protein, which was named IFNL4, may affect the HCV treatment response and clearance. In this chapter, we used molecular biology, virology, and patient genetic and phenotypic data to further characterize and understand the biology of IFNL4. The efforts in chapter 2 and 3 provided new insights to the candidate causal variant(s) responsible for the GWAS for HCV treatment response, however, more evidence is still required to make claims for the exact causal roles of these variants for the GWAS association.

Chapter 4 aims to characterize a mutation already known to cause a disease (seizure) in a mouse model. We demonstrate the potential use of multi-electrode array (MEA) system for the functional characterization and drug testing on mutations found in neurological diseases, such as seizure. Functional characterization in neurological diseases is relatively challenging and available systematic tools are relatively limited. This chapter shows an exploratory research and example to establish a system for the broader use for functional characterization and translational opportunities for mutations found in neurological diseases.

Overall, this dissertation spans a range of challenges of functional evaluations in human genetics. It is expected that the functional characterization to understand human mutations will become more central in human genetics, because there are still many biological questions remaining to be answered after the explosion of human genetic discoveries. The recent advance in several technologies, including genome editing and pluripotent stem cells, is also expected to make new tools available for functional studies in human diseases.

The Genetics of Success: How Single-Nucleotide Polymorphisms Associated With Educational Attainment Relate to Life-Course Development.

A previous genome-wide association study (GWAS) of more than 100,000 individuals identified molecular-genetic predictors of educational attainment. We undertook in-depth life-course investigation of the polygenic score derived from this GWAS using the four-decade Dunedin Study (N = 918). There were five main findings. First, polygenic scores predicted adult economic outcomes even after accounting for educational attainments. Second, genes and environments were correlated: Children with higher polygenic scores were born into better-off homes. Third, children's polygenic scores predicted their adult outcomes even when analyses accounted for their social-class origins; social-mobility analysis showed that children with higher polygenic scores were more upwardly mobile than children with lower scores. Fourth, polygenic scores predicted behavior across the life course, from early acquisition of speech and reading skills through geographic mobility and mate choice and on to financial planning for retirement. Fifth, polygenic-score associations were mediated by psychological characteristics, including intelligence, self-control, and interpersonal skill. Effect sizes were small. Factors connecting DNA sequence with life outcomes may provide targets for interventions to promote population-wide positive development.

Associations between smoking behavior-related alleles and the risk of melanoma.

Several studies have reported that cigarette smoking is inversely associated with the risk of melanoma. This study further tested whether incorporating genetic factors will provide another level of evaluation of mechanisms underlying the association between smoking and risk of melanoma. We investigated the association between SNPs selected from genome-wide association studies (GWAS) on smoking behaviors and risk of melanoma using 2,298 melanoma cases and 6,654 controls. Among 16 SNPs, three (rs16969968 [A], rs1051730 [A] and rs2036534 [C] in the 15q25.1 region) reached significance for association with melanoma risk in men (0.01 < = P values < = 0.02; 0.85 < = Odds Ratios (ORs) <= 1.20). There was association between the genetic scores based on the number of smoking behavior-risk alleles and melanoma risk with P-trend = 0.005 among HPFS. Further association with smoking behaviors indicating those three SNPs (rs16969968 [A], rs1051730 [A] and rs2036534 [C]) significantly associated with number of cigarettes smoked per day, CPD, with P = 0.009, 0.011 and 0.001 respectively. The SNPs rs215605 in the PDE1C gene and rs6265 in the BDNF gene significantly interacted with smoking status on melanoma risk (interaction P = 0.005 and P = 0.003 respectively). Our study suggests that smoking behavior-related SNPs are likely to play a role in melanoma development and the potential public health importance of polymorphisms in the CHRNA5-A3-B4 gene cluster. Further larger studies are warranted to validate the findings.

The Genetic Landscape of Renal Complications in Type 1 Diabetes

Diabetes is the leading cause of end stage renal disease. Despite evidence for a substantial heritability of diabetic kidney disease, efforts to identify genetic susceptibility variants have had limited success. We extended previous efforts in three dimensions, examining a more comprehensive set of genetic variants in larger numbers of subjects with type 1 diabetes characterized for a wider range of cross-sectional diabetic kidney disease phenotypes. In 2,843 subjects, we estimated that the heritability of diabetic kidney disease was 35% ( p=6x10-3 ). Genome-wide association analysis and replication in 12,540 individuals identified no single variants reaching stringent levels of significance and, despite excellent power, provided little independent confirmation of previously published associated variants. Whole exome sequencing in 997 subjects failed to identify any large-effect coding alleles of lower frequency influencing the risk of diabetic kidney disease. However, sets of alleles increasing body mass index ( p=2.2×10-5) and the risk of type 2 diabetes (p=6.1x10-4 ) were associated with the risk of diabetic kidney disease. We also found genome-wide genetic correlation between diabetic kidney disease and failure at smoking cessation ( p=1.1×10-4 ). Pathway analysis implicated ascorbate and aldarate metabolism ( p=9×10-6), and pentose and glucuronate interconversions ( p=3×10-6) in pathogenesis of diabetic kidney disease. These data provide further evidence for the role of genetic factors influencing diabetic kidney disease in those with type 1 diabetes and highlight some key pathways that may be responsible. Altogether these results reveal important biology behind the major cause of kidney disease.  

Novel Genetic Determinants of Diabetic Kidney Disease

Diabetic kidney disease (DKD) is a devastating diabetes complication, with known heritability not fully revealed by previous genetics studies. We performed the largest genome-wide association study of type 1 DKD to date, in a 13-cohort consortium of 15,590 individuals of European ancestry genotyped on the Illumina HumanCoreExome Beadchip, which allows exploration of coding variation in addition to genomic markers.

As prior work has shown that different characterizations of the DKD phenotype highlight distinct genetic associations, we investigated a spectrum of DKD definitions based on proteinuria and renal function criteria. Controls were DKD-free after a minimum of 15 years diabetes duration; cases had diabetes for at least 10 years prior to DKD diagnosis. We also performed a quantitative trait analysis of estimated glomerular filtration rate in all participants.

Our top finding was a missense mutation in COL4A3, rs55703767 (Asp326Tyr); the minor allele is common in Europeans (20%) and East Asians (13%) but not Africans (2%). This SNP had a genome-wide significant association with traditionally defined DKD (macroalbuminuria or end-stage renal disease [ESRD], (OR= 0.79, P=1.9×10-9), and a suggestive association with macroalbuminuria (OR= 0.79, P=1.6×10-6) and ESRD (OR= 0.79, P=4.5×10-5) individually. Though its PolyPhen score is 0.3 (benign), this SNP has been implicated as a splice site disruptor.

The COL4A3 gene encodes the alpha 3 subunit of Type IV collagen, the major structural component of basement membranes. Pathogenic mutations in COL4A3 have been identified in thin basement membrane nephropathy, familial focal segmental glomerulosclerosis, and Alport syndrome. A proxy (r2=0.6) for rs55703767 had no significant associations in the CKDGen consortium, suggesting its pathogenicity occurs solely in the setting of hyperglycemia.

By significantly increasing sample size we have discovered a novel locus underlying DKD risk, paving the way for better understanding of pathology, prevention, and treatment.

Pleiotropic Analysis of Lung Cancer and Blood Triglycerides

Epidemiologically related traits may share genetic risk factors, and pleiotropic analysis could identify individual loci associated with these traits. Because of their shared epidemiological associations, we conducted pleiotropic analysis of genome-wide association studies of lung cancer (12 160 lung cancer case patients and 16 838 control subjects) and cardiovascular disease risk factors (blood lipids from 188 577 subjects, type 2 diabetes from 148 821 subjects, body mass index from 123 865 subjects, and smoking phenotypes from 74 053 subjects). We found that 6p22.1 (rs6904596, ZNF184) was associated with both lung cancer (P = 5.50x10(-6)) and blood triglycerides (P = 1.39x10(-5)). We replicated the association in 6097 lung cancer case patients and 204 657 control subjects (P = 2.40 × 10(-4)) and in 71 113 subjects with triglycerides data (P = .01). rs6904596 reached genome-wide significance in lung cancer meta-analysis (odds ratio = 1.15, 95% confidence interval = 1.10 to 1.21 ,: Pcombined = 5.20x10(-9)). The large sample size provided by the lipid GWAS data and the shared genetic risk factors between the two traits contributed to the uncovering of a hitherto unidentified genetic locus for lung cancer.

Evidence for three genetic loci involved in both anorexia nervosa risk and variation of body mass index

The maintenance of normal body weight is disrupted in patients with anorexia nervosa (AN) for prolonged periods of time. Prior to the onset of AN, premorbid body mass index (BMI) spans the entire range from underweight to obese. After recovery, patients have reduced rates of overweight and obesity. As such, loci involved in body weight regulation may also be relevant for AN and vice versa. Our primary analysis comprised a cross-trait analysis of the 1000 single-nucleotide polymorphisms (SNPs) with the lowest P-values in a genome-wide association meta-analysis (GWAMA) of AN (GCAN) for evidence of association in the largest published GWAMA for BMI (GIANT). Subsequently we performed sex-stratified analyses for these 1000 SNPs. Functional ex vivo studies on four genes ensued. Lastly, a look-up of GWAMA-derived BMI-related loci was performed in the AN GWAMA. We detected significant associations (P-values <5 × 10(-5), Bonferroni-corrected P<0.05) for nine SNP alleles at three independent loci. Interestingly, all AN susceptibility alleles were consistently associated with increased BMI. None of the genes (chr. 10: CTBP2, chr. 19: CCNE1, chr. 2: CARF and NBEAL1; the latter is a region with high linkage disequilibrium) nearest to these SNPs has previously been associated with AN or obesity. Sex-stratified analyses revealed that the strongest BMI signal originated predominantly from females (chr. 10 rs1561589; Poverall: 2.47 × 10(-06)/Pfemales: 3.45 × 10(-07)/Pmales: 0.043). Functional ex vivo studies in mice revealed reduced hypothalamic expression of Ctbp2 and Nbeal1 after fasting. Hypothalamic expression of Ctbp2 was increased in diet-induced obese (DIO) mice as compared with age-matched lean controls. We observed no evidence for associations for the look-up of BMI-related loci in the AN GWAMA. A cross-trait analysis of AN and BMI loci revealed variants at three chromosomal loci with potential joint impact. The chromosome 10 locus is particularly promising given that the association with obesity was primarily driven by females. In addition, the detected altered hypothalamic expression patterns of Ctbp2 and Nbeal1 as a result of fasting and DIO implicate these genes in weight regulation.

Rare coding variants in the phospholipase D3 gene confer risk for Alzheimer's disease.

Genome-wide association studies (GWAS) have identified several risk variants for late-onset Alzheimer's disease (LOAD)1, 2. These common variants have replicable but small effects on LOAD risk and generally do not have obvious functional effects. Low-frequency coding variants, not detected by GWAS, are predicted to include functional variants with larger effects on risk. To identify low-frequency coding variants with large effects on LOAD risk, we carried out whole-exome sequencing (WES) in 14 large LOAD families and follow-up analyses of the candidate variants in several large LOAD case–control data sets. A rare variant in PLD3 (phospholipase D3; Val232Met) segregated with disease status in two independent families and doubled risk for Alzheimer’s disease in seven independent case–control series with a total of more than 11,000 cases and controls of European descent. Gene-based burden analyses in 4,387 cases and controls of European descent and 302 African American cases and controls, with complete sequence data for PLD3, reveal that several variants in this gene increase risk for Alzheimer’s disease in both populations. PLD3 is highly expressed in brain regions that are vulnerable to Alzheimer’s disease pathology, including hippocampus and cortex, and is expressed at significantly lower levels in neurons from Alzheimer’s disease brains compared to control brains. Overexpression of PLD3 leads to a significant decrease in intracellular amyloid-β precursor protein (APP) and extracellular Aβ42 and Aβ40 (the 42- and 40-residue isoforms of the amyloid-β peptide), and knockdown of PLD3 leads to a significant increase in extracellular Aβ42 and Aβ40. Together, our genetic and functional data indicate that carriers of PLD3 coding variants have a twofold increased risk for LOAD and that PLD3 influences APP processing. This study provides an example of how densely affected families may help to identify rare variants with large effects on risk for disease or other complex traits.

HNF1B variants associate with promoter methylation and regulate gene networks activated in prostate and ovarian cancer

Two independent regions within HNF1B are consistently identified in prostate and ovarian cancer genome-wide association studies (GWAS); their functional roles are unclear. We link prostate cancer (PC) risk SNPs rs11649743 and rs3760511 with elevated HNF1B gene expression and allele-specific epigenetic silencing, and outline a mechanism by which common risk variants could effect functional changes that increase disease risk: functional assays suggest that HNF1B is a pro-differentiation factor that suppresses epithelial-to-mesenchymal transition (EMT) in unmethylated, healthy tissues. This tumor-suppressor activity is lost when HNF1B is silenced by promoter methylation in the progression to PC. Epigenetic inactivation of HNF1B in ovarian cancer also associates with known risk SNPs, with a similar impact on EMT. This represents one of the first comprehensive studies into the pleiotropic role of a GWAS-associated transcription factor across distinct cancer types, and is the first to describe a conserved role for a multi-cancer genetic risk factor.

Eph kinases and their ligands ephrins act in concert with sex hormones in regulating blood pressure

Les Erythropoietin-producing hepatocyte (EPH) sont la plus grande famille de récepteurs tyrosine kinase. Leurs ligands, les éphrines (EFNs), sont aussi des molécules exprimées à la surface cellulaire. Les EPH/EFNs sont impliqués dans de nombreux processus biologiques. L'hypertension artérielle (PA) est une maladie chronique qui, aujourd'hui, est devenue un problème médical critique dans le monde entier et un enjeu de santé publique. La découverte de nouvelles thérapeutiques de l'hypertension sont d'une grande importance pour la santé publique. Jusqu’à tout récemment, il existe seulement quelques études concernant le rôle de l’axe EPH/EFNs sur la fonction des cellules musculaires lisses vasculaires (CMLV). Dans nos études précédentes, nous avons montré qu'EPHB6 et EFNB1, de concert avec les hormones sexuelles, régulent la PA. Dans la présente étude, nous avons constaté que les différents membres de la famille EPH/EFN peuvent réguler soit positivement, soit négativement, la contractilité des CMLV et la PA: tandis que EPHB4 et EFNB2 appartiennent à la première catégorie, EFNB1, EFNB3 et EPHB6 appartiennent à la deuxième. In vivo, des souris males, mais non pas des femelles, porteuses d’une mutation EPHB4 (KO) spécifique du muscle lisse présentent une PA diminuée, comparée aux souris témoins (WT). Les CMLV de souris EPHB4 KO, en présence de testostérone, ont montré une contractilité réduite lors de la stimulation par la phényléphrine (PE). Au niveau moléculaire, la phosphorylation de la protéine kinase II dépendante de Ca2+/calmoduline et de la kinase de la chaine légère de la myosine (CLM) est augmentée, tandis que la phosphorylation de la kinase de la CLM est réduite dans les CMLV KO lors de la stimulation par PE, par rapport au WT CMLV. Cela fournit une base moléculaire à la réduction de la PA et de la contractilité des CMLV chez les souris EPHB4 KO. EFNB2 est le ligand majeur de l’EPHB4. Comme attendu, les souris EFNB2 KO spécifique du muscle lisse avaient un phénotype de PA semblable, quoique non identique, aux souris EPHB4 KO. Les souris mâles EFNB2 KO, mais pas femelles, sous régime régulier ou riche en sel, présentent une PA réduite, par rapport à leurs homologues WT. Au niveau cellulaire, les CMLV des souris KO ont montré une contractilité réduite lors de la stimulation par PE par rapport aux témoins WT. Une région de l’acide aminé (aa) 313 à l’aa 331 dans la partie intracellulaire d’EFNB2 est essentielle pour la signalisation inverse qui régule la contractilité des CMLV, selon des études de mutation-délétion. Dans une étude de génétique humaine, nous avons identifié, dans le gène EFNB2, six SNP qui étaient associées significativement au risque d'hypertension artérielle, de façon dépendante du sexe, ce qui corrobore nos résultats chez les souris. En revanche, la délétion du gène EFNB3 (KO) chez les souris femelles aboutit à une PA élevée et à une augmentation des résistances des petites artères in vivo, améliore la contractilité des petites artères ex-vivo et augmente la contractilité des CMLV in vitro. Les souris mâles KO ont une PA normale, mais la castration conduit à une augmentation significative de la PA dans les souris KO, mais pas dans les souris WT. Les CMLV des souris KO femelles ont montré une phosphorylation accrue de la CLM et une phosphorylation réduite de la kinase de la CLM, ce qui fournit à nouveau une base moléculaire aux phénotypes de PA et de contractilité des CMLV observés. Ce changement de signalisation est attribuable à une protéine adaptatrice Grip1. En effet, dans une étude d'association pan génomique par le Consortium International pour la Pression Sanguine, un SNP dans le gène GRIP1 a approché le seuil de significativité de la valeur p pour son association avec la pression diastolique. Nos recherches, pour la première fois, ont révélé que EPH/EFNs sont de nouveaux composants dans le système de régulation de la PA. Les membres de la famille EPH/EFN peuvent agir comme des forces Yin et Yang pour régler finement le tonus des vaisseaux pour assurer l'homéostasie de la PA et de sa régulation. Ces effets de EPH/EFNs dépendent du sexe et des niveaux d’hormones sexuelles. À partir de ces nouvelles connaissances, nous pourrions développer une nouvelle thérapie personnalisée pour l’hypertension artérielle, utilisant des antagonistes d'hormones sexuelles ou des thérapies de remplacement d'hormones sexuelles, selon les niveaux d'hormones sexuelles des patients et les mutations dans les gènes de l'EPH/EFN.

IL10 Low-Frequency Variants in Behçet's Disease Patients

To explain the missing heritability after the genome-wide association studies era, sequencing studies allow the identification of low-frequency variants with a stronger effect on disease risk. Common variants in the interleukin 10 gene (IL10) have been consistently associated with Behçet's disease (BD) and the goal of this study is to investigate the role of low-frequency IL10 variants in BD susceptibility.