Revisiting the non-Gaucher-GBA-E326K carrier state: Is it sufficient to increase Parkinson's disease risk?

BACKGROUND: GBA variants are the most common genetic risk factors for Parkinson's disease (PD) world-wide, and can be found in up to 20% of Ashkenazi PD patients. The E326K variant, which is not considered a Gaucher's disease causing mutation, was recently shown to increase the risk for PD. Since E326K is a common variant among Europeans, Finnish and Ashkenazi (2.4, 8.6 and 1.2% carrier rate, respectively), we aimed to refine its involvement in PD. METHODS: 1200 consecutively recruited PD patients of a full Ashkenazi origin were genotyped for 10 GBA variants, the LRRK2-G2019S and the SMPD1-L302P. Alleles' frequencies were compared to controls, composed of 378 elderly healthy individuals and the non-neuro gnomAD Ashkenazi database. Odds-Ratio (OR) and age-at-motor-symptom-onset (AAO) were also calculated for all genotypes. RESULTS: All allelic variations tested had significant allelic ORs, demonstrating a wide range (1.86-12.84). The lowest allelic OR was observed for E326K (p = .013). Forty-five patients (of 1200, 3.75%) had at least two mutations (of the 12 tested), compared to 2 (0.53%) among 378 controls (p = .0013). Of the E326K carrier patients, 37% (10/27) carried additional mutations and the genotypic OR for individuals who carried only the E326K variant was 1.07. It did not reach statistical significance even when simulating the expected carrier frequency of E326K in 100,000 Ashkenazi controls (p = .39). In addition, an additive effect was demonstrated for risk in carriers of two mutations, the LRRK2-G2019S and a mild-GBA mutation (N370S or R496H), compared to carriers of only one mutation in one of these genes (simulated OR 11.79 compared to 7.58 and 2.49, respectively). An additive effect was also suggested for earlier AAO (5.0 years earlier than in non-carriers, compared to 3.1 and 2.2 years, respectively). CONCLUSIONS: Compared to previous studies, we demonstrate here a higher frequency of PD patients that carry two mutations. The GBA-E326K is more likely to affect PD risk when accompanied by another mutation, and an additive effect on risk and earlier AAO was proposed for carriers of LRRK2/mild-GBA double mutations. Altogether, these data support an oligogenic approach to PD genetics.

Functional confirmation that the R1488* variant in SCN9A results in complete loss-of-function of Nav1.7.

BACKGROUND: Individuals with an extremely rare inherited condition, termed Congenital Insensitivity to Pain (CIP), do not feel pain in response to noxious stimuli. Variants in SCN9A, encoding the transmembrane voltage-gated sodium channel Nav1.7, have previously been reported in subjects with CIP accompanied by anosmia, which are typically transmitted in a recessive pattern. Functional characterisations of some of these SCN9A mutations show that they result in complete loss-of-function of Nav1.7. METHODS: In a consanguineous family we performed whole exome sequencing of three members who have a diagnosis of CIP and one unaffected family member. The functional effects of the segregating variant in SCN9A were determined using patch clamp electrophysiology in human embryonic kidney (HEK) 293 cells transfected with the variant. RESULTS: We found that each CIP subject was homozygous for a putatively nonsense variant, R1488*, in SCN9A. This variant was reported elsewhere in a subject with CIP, though the functional effect was not determined. Using electrophysiology, we confirm that this variant results in a complete loss-of-function of Nav1.7. CONCLUSIONS: We confirm through electrophysiological analysis that this R1488* variant in SCN9A results in complete loss-of-function of Nav1.7, which is consistent with reports on other variants in this gene in subjects with CIP.

Whole exome re-sequencing implicates CCDC38 and cilia structure and function in resistance to smoking related airflow obstruction.

Chronic obstructive pulmonary disease (COPD) is a leading cause of global morbidity and mortality and, whilst smoking remains the single most important risk factor, COPD risk is heritable. Of 26 independent genomic regions showing association with lung function in genome-wide association studies, eleven have been reported to show association with airflow obstruction. Although the main risk factor for COPD is smoking, some individuals are observed to have a high forced expired volume in 1 second (FEV1) despite many years of heavy smoking. We hypothesised that these "resistant smokers" may harbour variants which protect against lung function decline caused by smoking and provide insight into the genetic determinants of lung health. We undertook whole exome re-sequencing of 100 heavy smokers who had healthy lung function given their age, sex, height and smoking history and applied three complementary approaches to explore the genetic architecture of smoking resistance. Firstly, we identified novel functional variants in the "resistant smokers" and looked for enrichment of these novel variants within biological pathways. Secondly, we undertook association testing of all exonic variants individually with two independent control sets. Thirdly, we undertook gene-based association testing of all exonic variants. Our strongest signal of association with smoking resistance for a non-synonymous SNP was for rs10859974 (P = 2.34 × 10(-4)) in CCDC38, a gene which has previously been reported to show association with FEV1/FVC, and we demonstrate moderate expression of CCDC38 in bronchial epithelial cells. We identified an enrichment of novel putatively functional variants in genes related to cilia structure and function in resistant smokers. Ciliary function abnormalities are known to be associated with both smoking and reduced mucociliary clearance in patients with COPD. We suggest that genetic influences on the development or function of cilia in the bronchial epithelium may affect growth of cilia or the extent of damage caused by tobacco smoke.

Genes contributing to pain sensitivity in the normal population: an exome sequencing study.

Sensitivity to pain varies considerably between individuals and is known to be heritable. Increased sensitivity to experimental pain is a risk factor for developing chronic pain, a common and debilitating but poorly understood symptom. To understand mechanisms underlying pain sensitivity and to search for rare gene variants (MAF<5%) influencing pain sensitivity, we explored the genetic variation in individuals' responses to experimental pain. Quantitative sensory testing to heat pain was performed in 2,500 volunteers from TwinsUK (TUK): exome sequencing to a depth of 70× was carried out on DNA from singletons at the high and low ends of the heat pain sensitivity distribution in two separate subsamples. Thus in TUK1, 101 pain-sensitive and 102 pain-insensitive were examined, while in TUK2 there were 114 and 96 individuals respectively. A combination of methods was used to test the association between rare variants and pain sensitivity, and the function of the genes identified was explored using network analysis. Using causal reasoning analysis on the genes with different patterns of SNVs by pain sensitivity status, we observed a significant enrichment of variants in genes of the angiotensin pathway (Bonferroni corrected p = 3.8×10(-4)). This pathway is already implicated in animal models and human studies of pain, supporting the notion that it may provide fruitful new targets in pain management. The approach of sequencing extreme exome variation in normal individuals has provided important insights into gene networks mediating pain sensitivity in humans and will be applicable to other common complex traits.

The VUKA family program: piloting a family-based psychosocial intervention to promote health and mental health among HIV infected early adolescents in South Africa.

An increasing number of adolescents born with HIV in South Africa are on antiretroviral treatment and have to confront complex issues related to coping with a chronic, stigmatizing and transmittable illness. Very few evidence-based mental health and health promotion programs for this population exist in South Africa. This study builds on a previous collaboratively designed and developmentally timed family-based intervention for early adolescents (CHAMP). The study uses community-based participatory approach as part of formative research to evaluate a pilot randomized control trial at two hospitals. The paper reports on the development, feasibility, and acceptability of the VUKA family-based program and its short-term impact on a range of psychosocial variables for HIV + preadolescents and their caregivers. A 10-session intervention of approximately 3-month duration was delivered to 65 preadolescents aged 10-13 years and their families. VUKA participants were noted to improve on all dimensions, including mental health, youth behavior, HIV treatment knowledge, stigma, communication, and adherence to medication. VUKA shows promise as a family-based mental and HIV prevention program for HIV + preadolescents and which could be delivered by trained lay staff.

Large candidate gene association study reveals genetic risk factors and therapeutic targets for fibromyalgia.

OBJECTIVE: Fibromyalgia (FM) represents a complex disorder that is characterized by widespread pain and tenderness and is frequently accompanied by additional somatic and cognitive/affective symptoms. Genetic risk factors are known to contribute to the etiology of the syndrome. The aim of this study was to examine >350 genes for association with FM, using a large-scale candidate gene approach. METHODS: The study group comprised 496 patients with FM (cases) and 348 individuals with no chronic pain (controls). Genotyping was performed using a dedicated gene array chip, the Pain Research Panel, which assays variants characterizing >350 genes known to be involved in the biologic pathways relevant to nociception, inflammation, and mood. Association testing was performed using logistic regression. RESULTS: Significant differences in allele frequencies between cases and controls were observed for 3 genes: GABRB3 (rs4906902; P = 3.65 × 10(-6)), TAAR1 (rs8192619; P = 1.11 × 10(-5)), and GBP1 (rs7911; P = 1.06 × 10(-4)). These 3 genes and 7 other genes with suggestive evidence for association were examined in a second, independent cohort of patients with FM and control subjects who were genotyped using the Perlegen 600K platform. Evidence of association in the replication cohort was observed for TAAR1, RGS4, CNR1, and GRIA4. CONCLUSION: Variation in these 4 replicated genes may serve as a basis for development of new diagnostic approaches, and the products of these genes may contribute to the pathophysiology of FM and represent potential targets for therapeutic action.

A combined functional annotation score for non-synonymous variants.

AIMS: Next-generation sequencing has opened the possibility of large-scale sequence-based disease association studies. A major challenge in interpreting whole-exome data is predicting which of the discovered variants are deleterious or neutral. To address this question in silico, we have developed a score called Combined Annotation scoRing toOL (CAROL), which combines information from 2 bioinformatics tools: PolyPhen-2 and SIFT, in order to improve the prediction of the effect of non-synonymous coding variants. METHODS: We used a weighted Z method that combines the probabilistic scores of PolyPhen-2 and SIFT. We defined 2 dataset pairs to train and test CAROL using information from the dbSNP: 'HGMD-PUBLIC' and 1000 Genomes Project databases. The training pair comprises a total of 980 positive control (disease-causing) and 4,845 negative control (non-disease-causing) variants. The test pair consists of 1,959 positive and 9,691 negative controls. RESULTS: CAROL has higher predictive power and accuracy for the effect of non-synonymous variants than each individual annotation tool (PolyPhen-2 and SIFT) and benefits from higher coverage. CONCLUSION: The combination of annotation tools can help improve automated prediction of whole-genome/exome non-synonymous variant functional consequences.

Genome-wide study of longitudinal brain imaging measures of multiple sclerosis progression across six clinical trials.

While the genetics of MS risk susceptibility are well-described, and recent progress has been made on the genetics of disease severity, the genetics of disease progression remain elusive. We therefore investigated the genetic determinants of MS progression on longitudinal brain MRI: change in brain volume (BV) and change in T2 lesion volume (T2LV), reflecting progressive tissue loss and increasing disease burden, respectively. We performed genome-wide association studies of change in BV (N = 3401) and change in T2LV (N = 3513) across six randomized clinical trials from Biogen and Roche/Genentech: ADVANCE, ASCEND, DECIDE, OPERA I & II, and ORATORIO. Analyses were adjusted for randomized treatment arm, age, sex, and ancestry. Results were pooled in a meta-analysis, and were evaluated for enrichment of MS risk variants. Variant colocalization and cell-specific expression analyses were performed using published cohorts. The strongest peaks were in PTPRD (rs77321193-C/A, p = 3.9 × 10-7) for BV change, and NEDD4L (rs11398377-GC/G, p = 9.3 × 10-8) for T2LV change. Evidence of colocalization was observed for NEDD4L, and both genes showed increased expression in neuronal and/or glial populations. No association between MS risk variants and MRI outcomes was observed. In this unique, precompetitive industry partnership, we report putative regions of interest in the neurodevelopmental gene PTPRD, and the ubiquitin ligase gene NEDD4L. These findings are distinct from known MS risk genetics, indicating an added role for genetic progression analyses and informing drug discovery.

The burden of rare protein-truncating genetic variants on human lifespan.

Genetic predisposition has been shown to contribute substantially to the age at which we die. Genome-wide association studies (GWASs) have linked more than 20 loci to phenotypes related to human lifespan1. However, little is known about how lifespan is impacted by gene loss of function. Through whole-exome sequencing of 352,338 UK Biobank participants of European ancestry, we assessed the relevance of protein-truncating variant (PTV) gene burden on individual and parental survival. We identified four exome-wide significant (P < 4.2 × 10-7) human lifespan genes, BRCA1, BRCA2, ATM and TET2. Gene and gene-set, PTV-burden, phenome-wide association studies support known roles of these genes in cancer to impact lifespan at the population level. The TET2 PTV burden was associated with a lifespan through somatic mutation events presumably due to clonal hematopoiesis. The overlap between PTV burden and common variant-based lifespan GWASs was modest, underscoring the value of exome sequencing in well-powered biobank cohorts to complement GWASs for identifying genes underlying complex traits.

OmicsView: Omics data analysis through interactive visual analytics.

With advances in NGS technologies, transcriptional profiling of human tissue across many diseases is becoming more routine, leading to the generation of petabytes of data deposited in public repositories. There is a need for bench scientists with little computational expertise to be able to access and mine this data to understand disease pathology, identify robust biomarkers of disease and the effect of interventions (in vivo or in vitro). To this end we release an open source analytics and visualization platform for expression data called OmicsView, http://omicsview.org. This platform comes preloaded with 1000 s of samples across many disease areas and normal tissue, including the GTEx database, all processed with a harmonized pipeline. We demonstrate the power and ease-of-use of the platform by means of a Crohn's disease data mining exercise where we can quickly uncover disease pathology and identify strong biomarkers of disease and response to treatment.

Pairwise effects between lipid GWAS genes modulate lipid plasma levels and cellular uptake.

Complex traits are characterized by multiple genes and variants acting simultaneously on a phenotype. However, studying the contribution of individual pairs of genes to complex traits has been challenging since human genetics necessitates very large population sizes, while findings from model systems do not always translate to humans. Here, we combine genetics with combinatorial RNAi (coRNAi) to systematically test for pairwise additive effects (AEs) and genetic interactions (GIs) between 30 lipid genome-wide association studies (GWAS) genes. Gene-based burden tests from 240,970 exomes show that in carriers with truncating mutations in both, APOB and either PCSK9 or LPL ("human double knock-outs") plasma lipid levels change additively. Genetics and coRNAi identify overlapping AEs for 12 additional gene pairs. Overlapping GIs are observed for TOMM40/APOE with SORT1 and NCAN. Our study identifies distinct gene pairs that modulate plasma and cellular lipid levels primarily via AEs and nominates putative drug target pairs for improved lipid-lowering combination therapies.

Neurology-related protein biomarkers are associated with cognitive ability and brain volume in older age.

Identifying biological correlates of late life cognitive function is important if we are to ascertain biomarkers for, and develop treatments to help reduce, age-related cognitive decline. Here, we investigated the associations between plasma levels of 90 neurology-related proteins (Olink® Proteomics) and general fluid cognitive ability in the Lothian Birth Cohort 1936 (LBC1936, N = 798), Lothian Birth Cohort 1921 (LBC1921, N = 165), and the INTERVAL BioResource (N = 4451). In the LBC1936, 22 of the proteins were significantly associated with general fluid cognitive ability (ß between -0.11 and -0.17). MRI-assessed total brain volume partially mediated the association between 10 of these proteins and general fluid cognitive ability. In an age-matched subsample of INTERVAL, effect sizes for the 22 proteins, although smaller, were all in the same direction as in LBC1936. Plasma levels of a number of neurology-related proteins are associated with general fluid cognitive ability in later life, mediated by brain volume in some cases.

Increased estrogen rather than decreased androgen action is associated with longer androgen receptor CAG repeats.

CONTEXT: The individual variability in the waning androgenic-anabolic functions of aging men may be influenced by the CAG repeat polymorphism in exon 1 of the androgen receptor (AR), affecting androgen sensitivity. However, findings on its phenotypic effects are inconclusive. OBJECTIVE: The aim was to investigate the relationships between health status, various reproductive hormones, and the AR CAG repeat length. DESIGN: We conducted a multinational prospective cohort observational study with cross-sectional baseline data. SETTING: This was a population survey of community-dwelling men. PARTICIPANTS: Men (40-79 yr old; n = 3369) were randomly recruited from centers in eight European countries; CAG repeat analysis was performed in 2878 men. MAIN OUTCOME MEASURES: We measured the correlations of the CAG repeat length with selected endocrine, metabolic, and phenotypic parameters related to aging and sex hormone action. RESULTS: Only minor differences were found in CAG repeat lengths between the eight European countries. They showed significant positive association with total, free, and bioavailable levels of testosterone (T) and estradiol. FSH but not LH correlated inversely with CAG repeat length. Significant associations were found with bone ultrasound parameters at the calcaneus. Negative correlation was found with triglycerides, but not with other blood lipids or with anthropometry, blood pressure, hemoglobin, insulin sensitivity, or sexual and prostatic functions. CONCLUSIONS: The AR CAG repeat length correlates significantly with serum T and estradiol of aging men. Weaker transcriptional activity of the AR with longer CAG-encoded polyglutamine repeats appears to be totally or nearly totally compensated for by higher T levels. The residual phenotypic correlations may reflect differences in estrogen levels/actions after aromatization of the higher T levels.

Discovering in vivo cytokine-eQTL interactions from a lupus clinical trial.

BACKGROUND: Cytokines are critical to human disease and are attractive therapeutic targets given their widespread influence on gene regulation and transcription. Defining the downstream regulatory mechanisms influenced by cytokines is central to defining drug and disease mechanisms. One promising strategy is to use interactions between expression quantitative trait loci (eQTLs) and cytokine levels to define target genes and mechanisms. RESULTS: In a clinical trial for anti-IL-6 in patients with systemic lupus erythematosus, we measure interferon (IFN) status, anti-IL-6 drug exposure, and whole blood genome-wide gene expression at three time points. We show that repeat transcriptomic measurements increases the number of cis eQTLs identified compared to using a single time point. We observe a statistically significant enrichment of in vivo eQTL interactions with IFN status and anti-IL-6 drug exposure and find many novel interactions that have not been previously described. Finally, we find transcription factor binding motifs interrupted by eQTL interaction SNPs, which point to key regulatory mediators of these environmental stimuli and therefore potential therapeutic targets for autoimmune diseases. In particular, genes with IFN interactions are enriched for ISRE binding site motifs, while those with anti-IL-6 interactions are enriched for IRF4 motifs. CONCLUSIONS: This study highlights the potential to exploit clinical trial data to discover in vivo eQTL interactions with therapeutically relevant environmental variables.

Phenome-wide association studies across large population cohorts support drug target validation.

Phenome-wide association studies (PheWAS) have been proposed as a possible aid in drug development through elucidating mechanisms of action, identifying alternative indications, or predicting adverse drug events (ADEs). Here, we select 25 single nucleotide polymorphisms (SNPs) linked through genome-wide association studies (GWAS) to 19 candidate drug targets for common disease indications. We interrogate these SNPs by PheWAS in four large cohorts with extensive health information (23andMe, UK Biobank, FINRISK, CHOP) for association with 1683 binary endpoints in up to 697,815 individuals and conduct meta-analyses for 145 mapped disease endpoints. Our analyses replicate 75% of known GWAS associations (P < 0.05) and identify nine study-wide significant novel associations (of 71 with FDR < 0.1). We describe associations that may predict ADEs, e.g., acne, high cholesterol, gout, and gallstones with rs738409 (p.I148M) in PNPLA3 and asthma with rs1990760 (p.T946A) in IFIH1. Our results demonstrate PheWAS as a powerful addition to the toolkit for drug discovery.

Religious Beliefs and Depression: Psychosocial Factors Affecting HIV Treatment Outcomes in South Africa.

Analyzing factors associated with virological failure (VF) may improve antiretroviral therapy (ART) outcomes for individuals living with HIV. The Risk Factors for Virological Failure (RFVF) study compared 158 cases with VF (viral load, VL, >1,000 copies/mL) and 300 controls with virological suppression (VL ≤1,000 copies/mL) after ≥5 months on their first ART regimen at McCord Hospital in Durban, South Africa between October 2010 and June 2012. RFVF participants completed a battery of various psychosocial measures. Using multivariate logistic regression stratified for gender, the association of various psychosocial factors with VF was assessed. It was found that not all factors were equally significant for both genders. The factors that were significantly associated with VF for both genders were younger age, shorter treatment duration and reporting depressive symptoms. The factors associated with VF that differed by gender were religious inactivity, having HIV+ family members, and status disclosure to friends.

TNF-857T, a genetic risk marker for acute anterior uveitis.

PURPOSE: To determine the association between 17 single nucleotide polymorphisms (SNPs) of tumor necrosis factor-alpha, lymphotoxin-alpha, and the TNF-receptors genes (TNF, LTA, and TNFRSF1A and -B) and idiopathic acute anterior uveitis (IAU) and to investigate their association with HLA-B27 and/or the development of visually significant complications. METHODS: Ninety-eight white patients in the United Kingdom were identified (by SL) from the uveitis clinics of Moorfields Eye Hospital (London, UK). Sequence-specific primers with 3' end mismatches were used to identify the presence of specific allelic variants by PCR amplification. RESULTS: There was a significant increase in the frequency of the TNF-857T allele in patients with IAU when compared with control subjects (15.3% vs. 7.3%, P = 0.0006). The frequency of haplotype 4, containing the T allele at nucleotide position -857, was also significantly increased in patients with IAU compared with control subjects (15.4% vs. 7.1%, P = 0.0003, OR 2.4, 95% confidence interval 1.4-4.0). In subgroup analysis, there were significant differences in the frequencies of the uncommon TNFRSF1A-201T and TNFRSF1A-1135T alleles between HLA-B27(+) patients with inflammation-related complications and those without complications (80.0% vs. 33.6%, P = 0.006; 80.0% vs. 36.6%, P = 0.01, respectively). CONCLUSIONS: A significant difference in the frequency of TNF-857T allele was found in patients with IAU. There was a trend toward the development of inflammation-related complications in HLA-B27(+) patients with IAU who were carriers of TNFRSF1A-201T or TNFRSF1A-1135T alleles. Genetic variations in these proinflammatory mediators and their receptors appear to influence the susceptibility and severity of the inflammatory response within the eyes of patients during the development of IAU.

Association of PADI4 and rheumatoid arthritis: a successful multidisciplinary approach.

The identification of functionally relevant polymorphisms of peptidylarginine deiminase 4, an enzyme that catalyzes the post-translational citrullination of proteins, as a rheumatoid arthritis gene is one of the most convincing success stories of complex disease gene mapping to date. In addition to an extensive single nucleotide polymorphism-based association study in a Japanese cohort, a range of techniques have been used to validate this finding.

Will the real disease gene please stand up?

A common dilemma arising in linkage studies of complex genetic diseases is the selection of positive signals, their follow-up with association studies and discrimination between true and false positive results. Several strategies for overcoming these issues have been devised. Using the Genetic Analysis Workshop 14 simulated dataset, we aimed to apply different analytical approaches and evaluate their performance in discerning real associations. We considered a) haplotype analyses, b) different methods adjusting for multiple testing, c) replication in a second dataset, and d) exhaustive genotyping of all markers in a sufficiently powered, large sample group. We found that haplotype-based analyses did not substantially improve over single-point analysis, although this may reflect the low levels of linkage disequilibrium simulated in the datasets provided. Multiple testing correction methods were in general found to be over-conservative. Replication of nominally positive results in a second dataset appears to be less stringent, resulting in the follow-up of false positives. Performing a comprehensive assay of all markers in a large, well-powered dataset appears to be the most effective strategy for complex disease gene identification.