Distinguishing Pleiotropy from Linked QTL Between Milk Production Traits

Total Page:16

File Type:pdf, Size:1020Kb

Distinguishing Pleiotropy from Linked QTL Between Milk Production Traits Cai et al. Genet Sel Evol (2020) 52:19 https://doi.org/10.1186/s12711-020-00538-6 Genetics Selection Evolution RESEARCH ARTICLE Open Access Distinguishing pleiotropy from linked QTL between milk production traits and mastitis resistance in Nordic Holstein cattle Zexi Cai1*†, Magdalena Dusza2†, Bernt Guldbrandtsen1, Mogens Sandø Lund1 and Goutam Sahana1 Abstract Background: Production and health traits are central in cattle breeding. Advances in next-generation sequencing technologies and genotype imputation have increased the resolution of gene mapping based on genome-wide association studies (GWAS). Thus, numerous candidate genes that afect milk yield, milk composition, and mastitis resistance in dairy cattle are reported in the literature. Efect-bearing variants often afect multiple traits. Because the detection of overlapping quantitative trait loci (QTL) regions from single-trait GWAS is too inaccurate and subjective, multi-trait analysis is a better approach to detect pleiotropic efects of variants in candidate genes. However, large sample sizes are required to achieve sufcient power. Multi-trait meta-analysis is one approach to deal with this prob- lem. Thus, we performed two multi-trait meta-analyses, one for three milk production traits (milk yield, protein yield and fat yield), and one for milk yield and mastitis resistance. Results: For highly correlated traits, the power to detect pleiotropy was increased by multi-trait meta-analysis com- pared with the subjective assessment of overlapping of single-trait QTL confdence intervals. Pleiotropic efects of lead single nucleotide polymorphisms (SNPs) that were detected from the multi-trait meta-analysis were confrmed by bivariate association analysis. The previously reported pleiotropic efects of variants within the DGAT1 and MGST1 genes on three milk production traits, and pleiotropic efects of variants in GHR on milk yield and fat yield were con- frmed. Furthermore, our results suggested that variants in KCTD16, KCNK18 and ENSBTAG00000023629 had pleiotropic efects on milk production traits. For milk yield and mastitis resistance, we identifed possible pleiotropic efects of variants in two genes, GC and DGAT1. Conclusions: Multi-trait meta-analysis improves our ability to detect pleiotropic interactions between milk produc- tion traits and identifes variants with pleiotropic efects on milk production traits and mastitis resistance. In particular, this should contribute to better understand the biological mechanisms that underlie the unfavorable genetic correla- tion between milk yield and mastitis. Background Intense selection for increased milk yield has negative Holstein is an important cattle breed in the Danish dairy consequences on the udder health of cows [1]. Unfavora- production and much efort has gone in the genetic ble genetic correlations between milk production and improvement of its milk production and functional traits. clinical mastitis (from 0.21 to 0.55) have been reported [2]. A genetic correlation between two traits could be due to the pleiotropic action of genetic variants or the corre- *Correspondence: [email protected] †Zexi Cai and Magdalena Dusza authors contributed equally to this study lation (i.e., linkage disequilibrium (LD)) between causal 1 Faculty of Technical Sciences, Center for Quantitative Genetics variants. Te identifcation of a quantitative trait locus and Genomics, Aarhus University, 8830 Tjele, Denmark (QTL) that afects simultaneously milk yield and udder Full list of author information is available at the end of the article health can help reveal some of the genetic basis of the © The Author(s) 2020. This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creat iveco mmons .org/publi cdoma in/ zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data. Cai et al. Genet Sel Evol (2020) 52:19 Page 2 of 15 genetic connection between milk production and masti- the animals were imputed to the high-density (HD) level. tis resistance. In combination with specifc genetic tests, Next, the imputed HD genotypes were imputed to the this information can contribute to reduce the unfavora- WGS level using a multi-breed reference of 1228 animals: ble correlated response on mastitis due to selection that 1148 WGS from Run4 of the 1000 Bull Genomes Pro- focused on improving milk production traits by diferen- ject (288 Holstein, 56 Red, and 61 Jersey, as well as 743 tially weighting variants based on their favorable or unfa- individuals from various breeds) [8] and 80 animals from vorable efects on the two traits. Aarhus University (23 HOL, 30 RDC, and 27 JER). Impu- One application of genome-wide association studies tation to HD genotypes was done with the IMPUTE2 (GWAS) is to detect pleiotropic efects for the QTL iden- v2.3.1 software [9], and imputation to the whole-genome tifed from single-trait analysis. If a genomic region is sig- level with the Minimac2 software [10]. Te average accu- nifcant for two or more traits, it may be due to causal racy (r2-values from Minimac2) was 0.85 for across- variants that are in LD and afect individual traits (link- breed imputation. Imputation accuracy was previously age), or that these traits are afected by the same vari- reported by Wu et al. [11]. ant (pleiotropy). Te number of segregating variants in Before phasing and imputation, we fltered the 54 k and a population is large, but fnite. Te proportion of the HD genotypes based on an SNP call rate higher than 85% segregating variants that are associated with the genetic and an animal call rate higher than 90%. Te imputed variation of complex traits is unknown. However, traits sequence data included 22,751,039 bi-allelic variants. For often appear to be associated with the same or closely- each breed, SNPs with a minor allele frequency (MAF) linked variants in the genome [3, 4], which strongly sug- lower than 1% or with a highly signifcant deviation from gests that, at least some of the underlying causal variants, Hardy–Weinberg proportions (p < 1.0−6) were excluded. afect several traits. Terefore, the primary aim of this After quality fltering, 16,503,508 SNPs remained for analysis was to determine whether QTL associated with analysis. more than one trait were indeed pleiotropic. We used our previous GWAS results (summary statistics) of milk pro- Single‑SNP association analysis for a trait duction traits [3] and mastitis resistance [4] to perform a Te GWAS summary statistics were from two previous multi-trait meta-analysis for scanning lead SNPs associ- association analyses [3, 4] and, here, we provide a brief ated with three milk production traits or with milk yield description of the GWAS method used. Te genetic rela- and mastitis resistance for pleiotropy. In combination tionship matrix (GRM) was estimated using imputed HD with a bivariate analysis, we examined the possible pleio- genotypes by the GCTA software [12]. We followed the tropic nature of the QTL identifed. leave-one-chromosome-out approach [13] to build a kin- ship matrix that was specifc to the analysis of each chro- Methods mosome by leaving out markers on that chromosome to Animals and phenotypes avoid loss of power due to double-counting (ftting the We used de-regressed proofs (DRP) based on estimated SNP as a fxed efect for testing associations and as a ran- breeding values (EBV) [5, 6] for milk, fat, protein yields dom efect as part of the GRM) [14]. and mastitis resistance (udder health index, which is an First, we performed a single-SNP GWAS analysis using index for clinical mastitis from frst to third lactation) for GCTA [12] for each chromosome. A Bonferroni mul- about 5000 Nordic Holstein (HOL) bulls. Nordic Cattle tiple-testing correction was applied to control for false- Genetic Evaluation (https ://www.nordi cebv.info/) pro- positive associations: a SNP was signifcant if its test vided the EBV. probability p-value, pM , was less than 0.05/ M , where M is the number of SNPs. Tis corresponds to a trait-wise Genotype and sequence data nominal type 1 error-rate of 5%. Here, the signifcance pM = 8.5 M ≈ 15.36 We performed an association analysis using imputed threshold value was − log10( ) with whole-genome sequence (WGS) data. All bulls (~ 5000) million SNPs. We identifed the lead SNPs for each inde- were genotyped with the BovineSNP50 BeadChip SNP pendent QTL signal on a chromosome by iteratively ft- array (54 k) versions 1 or 2 (Illumina Inc., San Diego, ting the allele dosages of the lead SNPs identifed in the CA). Imputation to WGS variants was described ear- previous runs as covariate (for details see [3, 4]). lier by Iso-Touru et al. [7]. Briefy, 54k genotypes were imputed to WGS variants by a 2-step approach. First, Genetic variance explained by the identifed QTL using a multi-breed reference of 3383 animals (1222 We used GCTA [12] to estimate the genetic variance Holstein (HOL), 1326 Nordic Red cattle (RDC) and 835 explained by all the identifed QTL together for each trait.
Recommended publications
  • PARSANA-DISSERTATION-2020.Pdf
    DECIPHERING TRANSCRIPTIONAL PATTERNS OF GENE REGULATION: A COMPUTATIONAL APPROACH by Princy Parsana A dissertation submitted to The Johns Hopkins University in conformity with the requirements for the degree of Doctor of Philosophy Baltimore, Maryland July, 2020 © 2020 Princy Parsana All rights reserved Abstract With rapid advancements in sequencing technology, we now have the ability to sequence the entire human genome, and to quantify expression of tens of thousands of genes from hundreds of individuals. This provides an extraordinary opportunity to learn phenotype relevant genomic patterns that can improve our understanding of molecular and cellular processes underlying a trait. The high dimensional nature of genomic data presents a range of computational and statistical challenges. This dissertation presents a compilation of projects that were driven by the motivation to efficiently capture gene regulatory patterns in the human transcriptome, while addressing statistical and computational challenges that accompany this data. We attempt to address two major difficulties in this domain: a) artifacts and noise in transcriptomic data, andb) limited statistical power. First, we present our work on investigating the effect of artifactual variation in gene expression data and its impact on trans-eQTL discovery. Here we performed an in-depth analysis of diverse pre-recorded covariates and latent confounders to understand their contribution to heterogeneity in gene expression measurements. Next, we discovered 673 trans-eQTLs across 16 human tissues using v6 data from the Genotype Tissue Expression (GTEx) project. Finally, we characterized two trait-associated trans-eQTLs; one in Skeletal Muscle and another in Thyroid. Second, we present a principal component based residualization method to correct gene expression measurements prior to reconstruction of co-expression networks.
    [Show full text]
  • University of Copenhagen, Copenhagen, Denmark, Citation: Pereira J, Johnson WE, O’Brien SJ, 8
    Evolutionary genomics and adaptive evolution of the hedgehog gene family (Shh, Ihh and Dhh) in vertebrates Pereira, Joana; Johnson, Warren E.; O'Brien, Stephen J.; Jarvis, Erich D; Zhang, Guojie; Gilbert, M. Thomas P.; Vasconcelos, Vitor; Antunes, Agostinho Published in: PloS one DOI: 10.1371/journal.pone.0074132 Publication date: 2014 Document version Publisher's PDF, also known as Version of record Document license: CC BY Citation for published version (APA): Pereira, J., Johnson, W. E., O'Brien, S. J., Jarvis, E. D., Zhang, G., Gilbert, M. T. P., Vasconcelos, V., & Antunes, A. (2014). Evolutionary genomics and adaptive evolution of the hedgehog gene family (Shh, Ihh and Dhh) in vertebrates. PloS one, 9(12), [e74132]. https://doi.org/10.1371/journal.pone.0074132 Download date: 24. sep.. 2021 RESEARCH ARTICLE Evolutionary Genomics and Adaptive Evolution of the Hedgehog Gene Family (Shh, Ihh and Dhh) in Vertebrates Joana Pereira1¤, Warren E. Johnson2, Stephen J. O’Brien3,4, Erich D. Jarvis5, Guojie Zhang6, M. Thomas P. Gilbert7, Vitor Vasconcelos1,8, Agostinho Antunes1,8* 1. CIIMAR/CIMAR, Interdisciplinary Centre of Marine and Environmental Research, University of Porto, Porto, Portugal, 2. Smithsonian Conservation Biology Institute, National Zoological Park, Front Royal, Virginia, United States of America, 3. Theodosius Dobzhansky Center for Genome Bioinformatics, St. Petersburg State University, St. Petersburg, Russia, 4. Oceanographic Center, N. Ocean Drive, Nova Southeastern University, Ft. Lauderdale, Florida, United States of America, 5. Howard Hughes Medical Institute, Department of Neurobiology, Duke University Medical Center, Durham, North Carolina, United States of America, 6. BGI-Shenzhen, Beishan Industrial Zoon, Yantian District, Shenzhen, China, 7.
    [Show full text]
  • Analysis of Gene Expression Data for Gene Ontology
    ANALYSIS OF GENE EXPRESSION DATA FOR GENE ONTOLOGY BASED PROTEIN FUNCTION PREDICTION A Thesis Presented to The Graduate Faculty of The University of Akron In Partial Fulfillment of the Requirements for the Degree Master of Science Robert Daniel Macholan May 2011 ANALYSIS OF GENE EXPRESSION DATA FOR GENE ONTOLOGY BASED PROTEIN FUNCTION PREDICTION Robert Daniel Macholan Thesis Approved: Accepted: _______________________________ _______________________________ Advisor Department Chair Dr. Zhong-Hui Duan Dr. Chien-Chung Chan _______________________________ _______________________________ Committee Member Dean of the College Dr. Chien-Chung Chan Dr. Chand K. Midha _______________________________ _______________________________ Committee Member Dean of the Graduate School Dr. Yingcai Xiao Dr. George R. Newkome _______________________________ Date ii ABSTRACT A tremendous increase in genomic data has encouraged biologists to turn to bioinformatics in order to assist in its interpretation and processing. One of the present challenges that need to be overcome in order to understand this data more completely is the development of a reliable method to accurately predict the function of a protein from its genomic information. This study focuses on developing an effective algorithm for protein function prediction. The algorithm is based on proteins that have similar expression patterns. The similarity of the expression data is determined using a novel measure, the slope matrix. The slope matrix introduces a normalized method for the comparison of expression levels throughout a proteome. The algorithm is tested using real microarray gene expression data. Their functions are characterized using gene ontology annotations. The results of the case study indicate the protein function prediction algorithm developed is comparable to the prediction algorithms that are based on the annotations of homologous proteins.
    [Show full text]
  • Universidade Estadual De Campinas Instituto De Biologia
    UNIVERSIDADE ESTADUAL DE CAMPINAS INSTITUTO DE BIOLOGIA VERÔNICA APARECIDA MONTEIRO SAIA CEREDA O PROTEOMA DO CORPO CALOSO DA ESQUIZOFRENIA THE PROTEOME OF THE CORPUS CALLOSUM IN SCHIZOPHRENIA CAMPINAS 2016 1 VERÔNICA APARECIDA MONTEIRO SAIA CEREDA O PROTEOMA DO CORPO CALOSO DA ESQUIZOFRENIA THE PROTEOME OF THE CORPUS CALLOSUM IN SCHIZOPHRENIA Dissertação apresentada ao Instituto de Biologia da Universidade Estadual de Campinas como parte dos requisitos exigidos para a obtenção do Título de Mestra em Biologia Funcional e Molecular na área de concentração de Bioquímica. Dissertation presented to the Institute of Biology of the University of Campinas in partial fulfillment of the requirements for the degree of Master in Functional and Molecular Biology, in the area of Biochemistry. ESTE ARQUIVO DIGITAL CORRESPONDE À VERSÃO FINAL DA DISSERTAÇÃO DEFENDIDA PELA ALUNA VERÔNICA APARECIDA MONTEIRO SAIA CEREDA E ORIENTADA PELO DANIEL MARTINS-DE-SOUZA. Orientador: Daniel Martins-de-Souza CAMPINAS 2016 2 Agência(s) de fomento e nº(s) de processo(s): CNPq, 151787/2F2014-0 Ficha catalográfica Universidade Estadual de Campinas Biblioteca do Instituto de Biologia Mara Janaina de Oliveira - CRB 8/6972 Saia-Cereda, Verônica Aparecida Monteiro, 1988- Sa21p O proteoma do corpo caloso da esquizofrenia / Verônica Aparecida Monteiro Saia Cereda. – Campinas, SP : [s.n.], 2016. Orientador: Daniel Martins de Souza. Dissertação (mestrado) – Universidade Estadual de Campinas, Instituto de Biologia. 1. Esquizofrenia. 2. Espectrometria de massas. 3. Corpo caloso.
    [Show full text]
  • Evolution, Dynamic Expression Changes and Regulatory Characteristics of Gene Families Involved in the Glycerophosphate Pathway O
    www.nature.com/scientificreports OPEN Evolution, dynamic expression changes and regulatory characteristics of gene families Received: 2 October 2018 Accepted: 14 August 2019 involved in the glycerophosphate Published: xx xx xxxx pathway of triglyceride synthesis in chicken (Gallus gallus) Liyu Yang1, Ziming Liu1, Kepeng Ou2, Taian Wang1, Zhuanjian Li1,3,4, Yadong Tian1,3,4, Yanbin Wang1,3,4, Xiangtao Kang1,3,4, Hong Li1,3,4 & Xiaojun Liu1,3,4 It is well documented that four gene families, including the glycerophosphate acyltransferases (GPATs), acylglycerophosphate acyltransferases (AGPATs), lipid phosphate phosphohydrolases (LPINs) and diacylglycerol acyltransferases (DGATs), are involved in the glycerophosphate pathway of de novo triglyceride (TG) biosynthesis in mammals. However, no systematic analysis has been conducted to characterize the gene families in poultry. In this study, the sequences of gene family members in the glycerophosphate pathway were obtained by screening the public databases. The phylogenetic tree, gene structures and conserved motifs of the corresponding proteins were evaluated. Dynamic expression changes of the genes at diferent developmental stages were analyzed by qRT-PCR. The regulatory characteristics of the genes were analyzed by in vivo experiments. The results showed that the GPAT, AGPAT and LPIN gene families have 2, 7 and 2 members, respectively, and they were classifed into 2, 4 and 2 cluster respectively based on phylogenetic analysis. All of the genes except AGPAT1 were extensively expressed in various tissues. Estrogen induction upregulated the expression of GPAM and AGPAT2, downregulated the expression of AGPAT3, AGPAT9, LPIN1 and LPIN2, and had no efect on the expression of the other genes. These fndings provide a valuable resource for further investigation of lipid metabolism in liver of chicken.
    [Show full text]
  • Chemoprevention with Cyclooxygenase and Epidermal Growth Factor Receptor Inhibitors in Familial Adenomatous Polyposis Patients
    Published OnlineFirst November 6, 2017; DOI: 10.1158/1940-6207.CAPR-17-0130 Research Article Cancer Prevention Research Chemoprevention with Cyclooxygenase and Epidermal Growth Factor Receptor Inhibitors in Familial Adenomatous Polyposis Patients: mRNA Signatures of Duodenal Neoplasia Don A. Delker1, Austin C. Wood2, Angela K. Snow2, N. Jewel Samadder1,2, Wade S. Samowitz2,3, Kajsa E. Affolter2,3, Kenneth M. Boucher1,2, Lisa M. Pappas2, Inge J. Stijleman2, Priyanka Kanth1, Kathryn R. Byrne1, Randall W. Burt1,2, Philip S. Bernard2,3, and Deborah W. Neklason1,2 Abstract To identify gene expression biomarkers and pathways tar- pared with paired baseline normal for patients on placebo geted by sulindac and erlotinib given in a chemoprevention and drug show that pathways activated in polyp growth and trial with a significant decrease in duodenal polyp burden at proliferation are blocked by this drug combination. Directly 6months(P < 0.001) in familial adenomatous polyposis comparing polyp gene expression between patients on drug (FAP) patients, we biopsied normal and polyp duodenal and placebo also identified innate immune response genes tissuesfrompatientsondrugversusplaceboandanalyzed (IL12 and IFNg) preferentially expressed in patients on drug. the RNA expression. RNA sequencing was performed on Gene expression analyses from tissue obtained at endpoint of biopsies from the duodenum of FAP patients obtained at the trial demonstrated inhibition of the cancer pathways baseline and 6-month endpoint endoscopy. Ten FAP patients COX2/PGE2, EGFR, and WNT. These findings provide molec- on placebo and 10 on sulindac and erlotinib were selected for ular evidence that the drug combination of sulindac and analysis. Purity of biopsied polyp tissue was calculated from erlotinib reached the intended tissue and was on target for RNA expression data.
    [Show full text]
  • Molecular and Physiological Basis for Hair Loss in Near Naked Hairless and Oak Ridge Rhino-Like Mouse Models: Tracking the Role of the Hairless Gene
    University of Tennessee, Knoxville TRACE: Tennessee Research and Creative Exchange Doctoral Dissertations Graduate School 5-2006 Molecular and Physiological Basis for Hair Loss in Near Naked Hairless and Oak Ridge Rhino-like Mouse Models: Tracking the Role of the Hairless Gene Yutao Liu University of Tennessee - Knoxville Follow this and additional works at: https://trace.tennessee.edu/utk_graddiss Part of the Life Sciences Commons Recommended Citation Liu, Yutao, "Molecular and Physiological Basis for Hair Loss in Near Naked Hairless and Oak Ridge Rhino- like Mouse Models: Tracking the Role of the Hairless Gene. " PhD diss., University of Tennessee, 2006. https://trace.tennessee.edu/utk_graddiss/1824 This Dissertation is brought to you for free and open access by the Graduate School at TRACE: Tennessee Research and Creative Exchange. It has been accepted for inclusion in Doctoral Dissertations by an authorized administrator of TRACE: Tennessee Research and Creative Exchange. For more information, please contact [email protected]. To the Graduate Council: I am submitting herewith a dissertation written by Yutao Liu entitled "Molecular and Physiological Basis for Hair Loss in Near Naked Hairless and Oak Ridge Rhino-like Mouse Models: Tracking the Role of the Hairless Gene." I have examined the final electronic copy of this dissertation for form and content and recommend that it be accepted in partial fulfillment of the requirements for the degree of Doctor of Philosophy, with a major in Life Sciences. Brynn H. Voy, Major Professor We have read this dissertation and recommend its acceptance: Naima Moustaid-Moussa, Yisong Wang, Rogert Hettich Accepted for the Council: Carolyn R.
    [Show full text]
  • Potassium Channels in Epilepsy
    Downloaded from http://perspectivesinmedicine.cshlp.org/ on September 28, 2021 - Published by Cold Spring Harbor Laboratory Press Potassium Channels in Epilepsy Ru¨diger Ko¨hling and Jakob Wolfart Oscar Langendorff Institute of Physiology, University of Rostock, Rostock 18057, Germany Correspondence: [email protected] This review attempts to give a concise and up-to-date overview on the role of potassium channels in epilepsies. Their role can be defined from a genetic perspective, focusing on variants and de novo mutations identified in genetic studies or animal models with targeted, specific mutations in genes coding for a member of the large potassium channel family. In these genetic studies, a demonstrated functional link to hyperexcitability often remains elusive. However, their role can also be defined from a functional perspective, based on dy- namic, aggravating, or adaptive transcriptional and posttranslational alterations. In these cases, it often remains elusive whether the alteration is causal or merely incidental. With 80 potassium channel types, of which 10% are known to be associated with epilepsies (in humans) or a seizure phenotype (in animals), if genetically mutated, a comprehensive review is a challenging endeavor. This goal may seem all the more ambitious once the data on posttranslational alterations, found both in human tissue from epilepsy patients and in chronic or acute animal models, are included. We therefore summarize the literature, and expand only on key findings, particularly regarding functional alterations found in patient brain tissue and chronic animal models. INTRODUCTION TO POTASSIUM evolutionary appearance of voltage-gated so- CHANNELS dium (Nav)andcalcium (Cav)channels, Kchan- nels are further diversified in relation to their otassium (K) channels are related to epilepsy newer function, namely, keeping neuronal exci- Psyndromes on many different levels, ranging tation within limits (Anderson and Greenberg from direct control of neuronal excitability and 2001; Hille 2001).
    [Show full text]
  • Supplementary Information Integrative Analyses of Splicing in the Aging Brain: Role in Susceptibility to Alzheimer’S Disease
    Supplementary Information Integrative analyses of splicing in the aging brain: role in susceptibility to Alzheimer’s Disease Contents 1. Supplementary Notes 1.1. Religious Orders Study and Memory and Aging Project 1.2. Mount Sinai Brain Bank Alzheimer’s Disease 1.3. CommonMind Consortium 1.4. Data Availability 2. Supplementary Tables 3. Supplementary Figures Note: Supplementary Tables are provided as separate Excel files. 1. Supplementary Notes 1.1. Religious Orders Study and Memory and Aging Project Gene expression data1. Gene expression data were generated using RNA- sequencing from Dorsolateral Prefrontal Cortex (DLPFC) of 540 individuals, at an average sequence depth of 90M reads. Detailed description of data generation and processing was previously described2 (Mostafavi, Gaiteri et al., under review). Samples were submitted to the Broad Institute’s Genomics Platform for transcriptome analysis following the dUTP protocol with Poly(A) selection developed by Levin and colleagues3. All samples were chosen to pass two initial quality filters: RNA integrity (RIN) score >5 and quantity threshold of 5 ug (and were selected from a larger set of 724 samples). Sequencing was performed on the Illumina HiSeq with 101bp paired-end reads and achieved coverage of 150M reads of the first 12 samples. These 12 samples will serve as a deep coverage reference and included 2 males and 2 females of nonimpaired, mild cognitive impaired, and Alzheimer's cases. The remaining samples were sequenced with target coverage of 50M reads; the mean coverage for the samples passing QC is 95 million reads (median 90 million reads). The libraries were constructed and pooled according to the RIN scores such that similar RIN scores would be pooled together.
    [Show full text]
  • A Dominant-Negative Mutation in the TRESK Potassium Channel Is Linked to Familial Migraine with Aura
    LETTERS A dominant-negative mutation in the TRESK potassium channel is linked to familial migraine with aura Ronald G Lafrenière1,2,13, M Zameel Cader3,4,13, Jean-François Poulin2, Isabelle Andres-Enguix5, Maryse Simoneau2, Namrata Gupta2, Karine Boisvert2, François Lafrenière2, Shannon McLaughlan2, Marie-Pierre Dubé6, Martin M Marcinkiewicz7, Sreeram Ramagopalan8, Olaf Ansorge9, Bernard Brais1, Jorge Sequeiros10, Jose Maria Pereira-Monteiro11, Lyn R Griffiths12, Stephen J Tucker5, George Ebers8 & Guy A Rouleau1,2 Migraine with aura is a common, debilitating, recurrent of the trigeminal ganglion afferents5 and may lead to central sensiti- headache disorder associated with transient and reversible focal zation. Considerable insights into the pathogenesis of migraine have neurological symptoms1. A role has been suggested for the come from the investigation of the rare autosomal dominant subtype two-pore domain (K2P) potassium channel, TWIK-related spinal of migraine with aura, familial hemiplegic migraine. Three suscepti- cord potassium channel (TRESK, encoded by KCNK18), in pain bility genes (CACNA1A, ATP1A2 and SCNA1), which encode either pathways and general anaesthesia2. We therefore examined ion channels or ion transport proteins, have so far been identified6–8, whether TRESK is involved in migraine by screening the and it is likely that mutations in these genes reduce the threshold for KCNK18 gene in subjects diagnosed with migraine. Here we CSD9,10. However, such mutations are not found in typical migraine report a frameshift mutation, F139WfsX24, which segregates with aura, suggesting that other ion channels may be involved. perfectly with typical migraine with aura in a large pedigree. K2P channels are expressed throughout the central nervous system We also identified prominent TRESK expression in migraine- and have a key role in controlling neuronal resting membrane potential salient areas such as the trigeminal ganglion.
    [Show full text]
  • Dual Proteome-Scale Networks Reveal Cell-Specific Remodeling of the Human Interactome
    bioRxiv preprint doi: https://doi.org/10.1101/2020.01.19.905109; this version posted January 19, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. Dual Proteome-scale Networks Reveal Cell-specific Remodeling of the Human Interactome Edward L. Huttlin1*, Raphael J. Bruckner1,3, Jose Navarrete-Perea1, Joe R. Cannon1,4, Kurt Baltier1,5, Fana Gebreab1, Melanie P. Gygi1, Alexandra Thornock1, Gabriela Zarraga1,6, Stanley Tam1,7, John Szpyt1, Alexandra Panov1, Hannah Parzen1,8, Sipei Fu1, Arvene Golbazi1, Eila Maenpaa1, Keegan Stricker1, Sanjukta Guha Thakurta1, Ramin Rad1, Joshua Pan2, David P. Nusinow1, Joao A. Paulo1, Devin K. Schweppe1, Laura Pontano Vaites1, J. Wade Harper1*, Steven P. Gygi1*# 1Department of Cell Biology, Harvard Medical School, Boston, MA, 02115, USA. 2Broad Institute, Cambridge, MA, 02142, USA. 3Present address: ICCB-Longwood Screening Facility, Harvard Medical School, Boston, MA, 02115, USA. 4Present address: Merck, West Point, PA, 19486, USA. 5Present address: IQ Proteomics, Cambridge, MA, 02139, USA. 6Present address: Vor Biopharma, Cambridge, MA, 02142, USA. 7Present address: Rubius Therapeutics, Cambridge, MA, 02139, USA. 8Present address: RPS North America, South Kingstown, RI, 02879, USA. *Correspondence: [email protected] (E.L.H.), [email protected] (J.W.H.), [email protected] (S.P.G.) #Lead Contact: [email protected] bioRxiv preprint doi: https://doi.org/10.1101/2020.01.19.905109; this version posted January 19, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder.
    [Show full text]
  • The Endocytic Membrane Trafficking Pathway Plays a Major Role
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by University of Liverpool Repository RESEARCH ARTICLE The Endocytic Membrane Trafficking Pathway Plays a Major Role in the Risk of Parkinson’s Disease Sara Bandres-Ciga, PhD,1,2 Sara Saez-Atienzar, PhD,3 Luis Bonet-Ponce, PhD,4 Kimberley Billingsley, MSc,1,5,6 Dan Vitale, MSc,7 Cornelis Blauwendraat, PhD,1 Jesse Raphael Gibbs, PhD,7 Lasse Pihlstrøm, MD, PhD,8 Ziv Gan-Or, MD, PhD,9,10 The International Parkinson’s Disease Genomics Consortium (IPDGC), Mark R. Cookson, PhD,4 Mike A. Nalls, PhD,1,11 and Andrew B. Singleton, PhD1* 1Molecular Genetics Section, Laboratory of Neurogenetics, National Institute on Aging, National Institutes of Health, Bethesda, Maryland, USA 2Instituto de Investigación Biosanitaria de Granada (ibs.GRANADA), Granada, Spain 3Transgenics Section, Laboratory of Neurogenetics, National Institute on Aging, National Institutes of Health, Bethesda, Maryland, USA 4Cell Biology and Gene Expression Section, Laboratory of Neurogenetics, National Institute on Aging, National Institutes of Health, Bethesda, Maryland, USA 5Department of Molecular and Clinical Pharmacology, Institute of Translational Medicine, University of Liverpool, Liverpool, United Kingdom 6Department of Pathophysiology, University of Tartu, Tartu, Estonia 7Computational Biology Group, Laboratory of Neurogenetics, National Institute on Aging, National Institutes of Health, Bethesda, Maryland, USA 8Department of Neurology, Oslo University Hospital, Oslo, Norway 9Department of Neurology and Neurosurgery, Department of Human Genetics, McGill University, Montréal, Quebec, Canada 10Department of Neurology and Neurosurgery, Montreal Neurological Institute, McGill University, Montréal, Quebec, Canada 11Data Tecnica International, Glen Echo, Maryland, USA ABSTRACT studies, summary-data based Mendelian randomization Background: PD is a complex polygenic disorder.
    [Show full text]