Genome Wide Association of Chronic Kidney Disease Progression: the CRIC Study (Author List and Affiliations Listed at End of Document)

Total Page:16

File Type:pdf, Size:1020Kb

Genome Wide Association of Chronic Kidney Disease Progression: the CRIC Study (Author List and Affiliations Listed at End of Document) SUPPLEMENTARY MATERIALS Genome Wide Association of Chronic Kidney Disease Progression: The CRIC Study (Author list and affiliations listed at end of document) Genotyping information page 2 Molecular pathway analysis information page 3 Replication cohort acknowledgments page 4 Supplementary Table 1. AA top hit region gene function page 5-6 Supplementary Table 2. EA top hit region gene function page 7 Supplementary Table 3. GSA pathway results page 8 Supplementary Table 4. Number of molecular interaction based on top candidate gene molecular networks page 9 Supplementary Table 5. Results of top gene marker association in AA, based on EA derived candidate gene regions page 10 Supplementary Table 6. Results of top gene marker association in EA, based on AA derived candidate gene regions page 11 Supplementary Table 7. EA Candidate SNP look up page 12 Supplementary Table 8. AA Candidate SNP look up page 13 Supplementary Table 9. Replication cohorts page 14 Supplementary Table 10. Replication cohort study characteristics page 15 Supplementary Figure 1a-b. Boxplot of eGFR decline in AA and EA page 16 Supplementary Figure 2a-l. Regional association plot of candidate SNPs identified in AA groups pages 17-22 Supplementary Figure 3a-f. Regional association plot of candidate SNPs identified in EA groups pages 23-25 Supplementary Figure 4. Molecular Interaction network of candidate genes for renal, cardiovascular and immunological diseases pages 26-27 Supplementary Figure 5. Molecular Interaction network of candidate genes for renal diseases pages 28-29 Supplementary Figure 6. ARRDC4 LD map page 30 Author list and affiliations page 31 1 Supplemental Materials Genotyping Genotyping was performed on a total of 3,635 CRIC participants who provided specific consent for investigations of inherited genetics (of a total of 3,939 CRIC participants). Genotyping was conducted at the Genetic Analysis Platform, Broad Institute of MIT and Harvard, using the Illumina HumanOmni1- Quad v1.0 microarray, which comprised of 1,140,419 SNPs. SNP genotypes were called using Illumina’s BeadStudio Genotyping Module (Illumina Inc, San Diego, CA, USA). Among the 3,635 participant samples, 81 were excluded based on quality control metrics including sex discordance (n=7), reduced or excess genotypic heterozygosity (heterozygosity rate ± 3 standard deviations from the mean; n=25), cryptic relatedness (π >0.185; n=64) or sample ID mismatches (n=28); 16 individuals were excluded because of more than one of these four metrics. All samples had a genotype call rate of at least 97%. This resulted in 3,527 samples that passed our quality control metrics. We then undertook principle component (PC) analyses to infer genetic ancestry on these remaining 3,527 samples. Using the cut points of PC 1 ≥ 0.0298 and PC 2 ≥ 0.0651, we identified a cluster of 1,581 participants of European ancestry (whites) and a cluster of 1,493 participants of African ancestry (blacks). Of the potential 1,140,419 markers on the genotype array, 184,789 were removed based on poor clustering or replication failures (i.e. assay performance). Then separately within the white and black subgroup populations, we completed SNP marker-level quality control on the remaining 955,630 SNP makers. In whites, we excluded a total of 214,437 markers (22.4%) because of minor allele frequency (MAF) < 0.03, 531 (0.06 %) because of deviation from Hardy-Weinberg equilibrium (P < 1 × 10-7), and 1,053 (0.11 %) because of genotype call rate < 0.95; 739,978 SNP markers remained in the discovery phase for white. Among blacks, we excluded a total of 111,478 markers (12.0%) because of MAF< 0.03, 1,047 (0.11 %) because of deviation from Hardy-Weinberg equilibrium (P < 1 × 10-7), and 1,109 (0.12 %) because of genotype call rate < 0.95; 839,205 SNP markers remained in the discovery phase for blacks. We completed genome-wide imputation based on the 1000 Genomes mixed race/ethnicity (ALL) genetic backbone (NCBI build 37, release date March 2012; n=1093) using IMPUTE2. Only those imputed SNP markers passing imputation quality control (info [a measure of r2] ≥ 0.80) were retained for analysis, resulting in a total of 8,191,067 and 11,726,769 variants with MAF >0.03 in our EA and AA cohorts, respectively. 2 Molecular pathway analysis information Of our 18 identified candidate SNPs, based on physical proximity, we were able to assign 16 SNP to specific genes. Interaction networks centered around the molecule products of these 16 genes were constructed through the use of Integrative pathway analysis (Ingenuity® Systems, www.ingenuity.com). An interaction network is a graphical representation of the molecular relationships between genes and gene products. Gene products are represented as nodes, and the biological relationship between two nodes is represented as a line that is supported by at least one reference from the literature. In the graphic network, type categories of the molecules are represented by various shapes, and the labels along with each line indicate varieties of interactions between the connected molecules. Supplementary Figures 4and 5 list the legends for the type categories of the molecules, and the nature of the interactions between molecules. Only experimentally verified interaction relationships, either direct physical contact relationship or indirect regulation relationship through intermediate molecules, were considered. Exogenous chemicals were excluded from the interaction network analysis. We found eight of our a-priori assigned gene products to be associated with renal or urological function (Supplementary Figure 5). Five of the SNPs were located within and three of the others SNPs or SNPs in LD with them were in close proximity (within 30kb) to the assigned genes (Supplementary Figures2- c,g,j,k and 3-b,c,d,e). In order to test for the statistical significance of our identified renal related molecular networks, the p-value was calculated based on Fisher's Exact Test to assess the degree of over- represented pathway compare to as expected by chance by accounting for: 1) the number of Functions/Pathways/Lists Eligible molecules that participate in that annotation as defined by the molecules in the selected Reference set; 2) the total number of molecules in the selected Reference set known to be associated with that function; 3) the total number of Functions/Pathways/Lists Eligible molecules in the selected Reference set; and 4) the total number of molecules in the Reference Set. 3 Replication Cohorts acknowledgments This study was supported by grants U01DK57292, U01DK57329, U01DK057300, U01DK057298, U01DK057249, U01DK57295, U01DK070657, U01DK057303, U01DK070657, U01DK57304 and DK57292-05 from the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) and, in part, by the Intramural Research Program of the NIDDK. This project has been funded in whole or in part with federal funds from the National FIND study (AA) Cancer Institute, National Institutes of Health (NIH), under contract N01-CO-12400 and the Intramural Research Program of the NIH, National Cancer Institute, Center for Cancer Research. This work was also supported by the National Center for Research Resources for the General Clinical Research Center grants: Case Western Reserve University, M01-RR-000080; Wake Forest University, M01-RR-07122; Harbor-University of California, Los Angeles Medical Center, M01-RR-00425; College of Medicine, University of California, Irvine, M01-RR-00827– 29; University of New Mexico, HSC M01-RR-00997; and Frederic C. Bartter, M01-RR-01346. Genotyping was performed by the Center for Inherited Disease Research, which is fully funded through a federal contract from the NIH to Johns Hopkins University (N01-HG-65403). The results of this analysis were obtained using the S.A.G.E. package of genetic epidemiology software, which is supported by a U.S. Public Health Service Resource Grant (RR03655) from the National Center for Research Resources. Wake Forest T2D- NIH R01 DK066358 (DWB), R01 DK053591 (DWB), R01 HL56266 (BIF), R01 DK070941 ESKD Study (BIF), General Clinical Research Center of Wake Forest University School of Medicine M01 RR07122. NIH R01 DK 070941 (BIF) and R01 DK53591 (DWB), and by the NIDDK Wake Forest non- and NCI Intramural Research Programs. MAB was supported by F32 diabetic ESKD Study DK080617 from the NIDDK. See funding sources of the CKDGen Cohorts in the manuscript with PMID 25493955. CKDGen CAB's work was supported by the Else Kröner-Fresenius-Stiftung. We thank all the FinnDiane researchers, as well as the physicians and nurses at each center participating in the recruitment of the patients and in the collection of samples and data (Sandholm et al. PloS Gen2012; Thorn LM et al., Diabetes Care 2005). The FinnDiane study was supported by grants from the Folkhälsan Research Foundation, Liv och Hälsa Foundation, the FinnDiane Willhelm and Else Stockmann Foundation, Helsinki University Central Hospital Research Funds (EVO), the Finnish Cultural Foundation, the Signe and Ane Gyllenberg Foundation, Finnish Medical Society (Finska läkaresällskapet), Academy of Finland, Novo Nordisk Foundation and Tekes. AASK U01 DK048689, M01 RR-00071 4 Supplementary Table 1. African American top hit region candidate gene function SNP to Gene Assignment (location or Gene name function and associations distance to gene) rs12057968 Chromosome 1 Open Reading Frame 100: Protein coding with unknown function. C1orf100 Expressed in kidney. Intronic Adenylosuccinate synthetase: nearby gene with overlapping SNPs in LD with C1orf100. ADSS catalyzes the first committed step in the conversion of inosine monophosphate to (ADSS-23kb) adenosine monophosphate. May affect energy metabolism. Associated with schizophrenia. Expressed in kidney. Bardet-Biedl Syndrome 9: protein coding gene associated with Bardet-Biedl syndrome, rs73690944 which includes renal malformations and CKD. The exact function of this gene has not yet BBS9 been determined. Thought to function as a coat complex required for sorting of specific Intronic membrane proteins to the primary cilia. Also involved in parathyroid hormone action in bones.
Recommended publications
  • Gene Expression Studies: from Case-Control to Multiple-Population-Based Studies
    From the Institute of Human Genetics, Helmholtz Zentrum Munchen,¨ Deutsches Forschungszentrum fur¨ Gesundheit und Umwelt (GmbH) Head: Prof. Dr. Thomas Meitinger Gene expression studies: From case-control to multiple-population-based studies Thesis Submitted for a Doctoral Degree in Natural Sciences at the Faculty of Medicine, Ludwig-Maximilians-Universitat¨ Munchen¨ Katharina Schramm Dachau, Germany 2016 With approval of the Faculty of Medicine Ludwig-Maximilians-Universit¨atM ¨unchen Supervisor/Examiner: Prof. Dr. Thomas Illig Co-Examiners: Prof. Dr. Roland Kappler Dean: Prof. Dr. med. dent. Reinhard Hickel Date of oral examination: 22.12.2016 II Dedicated to my family. III Abstract Recent technological developments allow genome-wide scans of gene expression levels. The reduction of costs and increasing parallelization of processing enable the quantification of 47,000 transcripts in up to twelve samples on a single microarray. Thereby the data collec- tion of large population-based studies was improved. During my PhD, I first developed a workflow for the statistical analyses of case-control stu- dies of up to 50 samples. With large population-based data sets generated I established a pipeline for quality control, data preprocessing and correction for confounders, which re- sulted in substantially improved data. In total, I processed more than 3,000 genome-wide expression profiles using the generated pipeline. With 993 whole blood samples from the population-based KORA (Cooperative Health Research in the Region of Augsburg) study we established one of the largest population-based resource. Using this data set we contributed to a number of transcriptome-wide association studies within national (MetaXpress) and international (CHARGE) consortia.
    [Show full text]
  • Differential Expression Profile Prioritization of Positional Candidate Glaucoma Genes the GLC1C Locus
    LABORATORY SCIENCES Differential Expression Profile Prioritization of Positional Candidate Glaucoma Genes The GLC1C Locus Frank W. Rozsa, PhD; Kathleen M. Scott, BS; Hemant Pawar, PhD; John R. Samples, MD; Mary K. Wirtz, PhD; Julia E. Richards, PhD Objectives: To develop and apply a model for priori- est because of moderate expression and changes in tization of candidate glaucoma genes. expression. Transcription factor ZBTB38 emerges as an interesting candidate gene because of the overall expres- Methods: This Affymetrix GeneChip (Affymetrix, Santa sion level, differential expression, and function. Clara, Calif) study of gene expression in primary cul- ture human trabecular meshwork cells uses a positional Conclusions: Only1geneintheGLC1C interval fits our differential expression profile model for prioritization of model for differential expression under multiple glau- candidate genes within the GLC1C genetic inclusion in- coma risk conditions. The use of multiple prioritization terval. models resulted in filtering 7 candidate genes of higher interest out of the 41 known genes in the region. Results: Sixteen genes were expressed under all condi- tions within the GLC1C interval. TMEM22 was the only Clinical Relevance: This study identified a small sub- gene within the interval with differential expression in set of genes that are most likely to harbor mutations that the same direction under both conditions tested. Two cause glaucoma linked to GLC1C. genes, ATP1B3 and COPB2, are of interest in the con- text of a protein-misfolding model for candidate selec- tion. SLC25A36, PCCB, and FNDC6 are of lesser inter- Arch Ophthalmol. 2007;125:117-127 IGH PREVALENCE AND PO- identification of additional GLC1C fami- tential for severe out- lies7,18-20 who provide optimal samples for come combine to make screening candidate genes for muta- adult-onset primary tions.7,18,20 The existence of 2 distinct open-angle glaucoma GLC1C haplotypes suggests that muta- (POAG) a significant public health prob- tions will not be limited to rare descen- H1 lem.
    [Show full text]
  • A Computational Approach for Defining a Signature of Β-Cell Golgi Stress in Diabetes Mellitus
    Page 1 of 781 Diabetes A Computational Approach for Defining a Signature of β-Cell Golgi Stress in Diabetes Mellitus Robert N. Bone1,6,7, Olufunmilola Oyebamiji2, Sayali Talware2, Sharmila Selvaraj2, Preethi Krishnan3,6, Farooq Syed1,6,7, Huanmei Wu2, Carmella Evans-Molina 1,3,4,5,6,7,8* Departments of 1Pediatrics, 3Medicine, 4Anatomy, Cell Biology & Physiology, 5Biochemistry & Molecular Biology, the 6Center for Diabetes & Metabolic Diseases, and the 7Herman B. Wells Center for Pediatric Research, Indiana University School of Medicine, Indianapolis, IN 46202; 2Department of BioHealth Informatics, Indiana University-Purdue University Indianapolis, Indianapolis, IN, 46202; 8Roudebush VA Medical Center, Indianapolis, IN 46202. *Corresponding Author(s): Carmella Evans-Molina, MD, PhD ([email protected]) Indiana University School of Medicine, 635 Barnhill Drive, MS 2031A, Indianapolis, IN 46202, Telephone: (317) 274-4145, Fax (317) 274-4107 Running Title: Golgi Stress Response in Diabetes Word Count: 4358 Number of Figures: 6 Keywords: Golgi apparatus stress, Islets, β cell, Type 1 diabetes, Type 2 diabetes 1 Diabetes Publish Ahead of Print, published online August 20, 2020 Diabetes Page 2 of 781 ABSTRACT The Golgi apparatus (GA) is an important site of insulin processing and granule maturation, but whether GA organelle dysfunction and GA stress are present in the diabetic β-cell has not been tested. We utilized an informatics-based approach to develop a transcriptional signature of β-cell GA stress using existing RNA sequencing and microarray datasets generated using human islets from donors with diabetes and islets where type 1(T1D) and type 2 diabetes (T2D) had been modeled ex vivo. To narrow our results to GA-specific genes, we applied a filter set of 1,030 genes accepted as GA associated.
    [Show full text]
  • Microrna Binding Site Polymorphisms Are Associated with the Development of Gastric Cancer
    Journal of Cell and Molecular Research (2017) 9 (2), 73-77 DOI: 10.22067/jcmr.v9i2.71487 Research Article MicroRNA Binding Site Polymorphisms are Associated with the Development of Gastric Cancer Jamal Asgarpour1, Jamshid Mehrzad* 2, Seyyed Abbas Tabatabaei Yazdi3 1Department of Genetic, Islamic Azad University, Neyshabur, Iran 2 Department of Biochemistry, Islamic Azad University, Neyshabur, Iran 3 Department of Pathology, Mashhad University of Medical Sciences, Mashhad, Iran Received 20 November 2017 Accepted 15 December 2017 Abstract MicroRNAs (miRNAs) binding to the 3'-untranslated regions (3'-UTRs) of messenger RNAs (mRNAs), affect several cellular mechanisms such as translation, differentiation, tumorigenesis, carcinogenesis and apoptosis. The occurrence of genetic polymorphisms in 3'-UTRs of target genes affect the binding affinity of miRNAs with the target genes resulting in their altered expression. The current case-control study of 100 samples (50% cancer patients and 50% health persons as control) was aimed to evaluate the genotyping of microRNA single nucleotide polymorphisms (SNPs) located at the binding site of the 3'-UTR of C14orf101 (rs4901706) and mir-124 (- rs531564) genes and their correlations with gastric cancer (GC) development. The Statistical analysis results indicated the significant association of AG development risk with the SNPs located in the 3'-UTR of C14orf101 and mir-124. It could be concluded from the results that these genes are associated with the gastric adenocarcinoma. Keywords: Carcinogenesis, C14orf101 polymorphism, Mir-124, Gastric cancer development Introduction The alternative name of C14orF101 is TMR260 Gastric cancer is regarded as one of the most (NCBI gene number; rs4901706) which is common malignancies in the world and is the responsible to translate a transmembrane second leading cause of death around the world.
    [Show full text]
  • Redefining the Specificity of Phosphoinositide-Binding by Human
    bioRxiv preprint doi: https://doi.org/10.1101/2020.06.20.163253; this version posted June 21, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC 4.0 International license. Redefining the specificity of phosphoinositide-binding by human PH domain-containing proteins Nilmani Singh1†, Adriana Reyes-Ordoñez1†, Michael A. Compagnone1, Jesus F. Moreno Castillo1, Benjamin J. Leslie2, Taekjip Ha2,3,4,5, Jie Chen1* 1Department of Cell & Developmental Biology, University of Illinois at Urbana-Champaign, Urbana, IL 61801; 2Department of Biophysics and Biophysical Chemistry, Johns Hopkins University School of Medicine, Baltimore, MD 21205; 3Department of Biophysics, Johns Hopkins University, Baltimore, MD 21218; 4Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21205; 5Howard Hughes Medical Institute, Baltimore, MD 21205, USA †These authors contributed equally to this work. *Correspondence: [email protected]. bioRxiv preprint doi: https://doi.org/10.1101/2020.06.20.163253; this version posted June 21, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC 4.0 International license. ABSTRACT Pleckstrin homology (PH) domains are presumed to bind phosphoinositides (PIPs), but specific interaction with and regulation by PIPs for most PH domain-containing proteins are unclear. Here we employed a single-molecule pulldown assay to study interactions of lipid vesicles with full-length proteins in mammalian whole cell lysates.
    [Show full text]
  • 1 Supporting Information for a Microrna Network Regulates
    Supporting Information for A microRNA Network Regulates Expression and Biosynthesis of CFTR and CFTR-ΔF508 Shyam Ramachandrana,b, Philip H. Karpc, Peng Jiangc, Lynda S. Ostedgaardc, Amy E. Walza, John T. Fishere, Shaf Keshavjeeh, Kim A. Lennoxi, Ashley M. Jacobii, Scott D. Rosei, Mark A. Behlkei, Michael J. Welshb,c,d,g, Yi Xingb,c,f, Paul B. McCray Jr.a,b,c Author Affiliations: Department of Pediatricsa, Interdisciplinary Program in Geneticsb, Departments of Internal Medicinec, Molecular Physiology and Biophysicsd, Anatomy and Cell Biologye, Biomedical Engineeringf, Howard Hughes Medical Instituteg, Carver College of Medicine, University of Iowa, Iowa City, IA-52242 Division of Thoracic Surgeryh, Toronto General Hospital, University Health Network, University of Toronto, Toronto, Canada-M5G 2C4 Integrated DNA Technologiesi, Coralville, IA-52241 To whom correspondence should be addressed: Email: [email protected] (M.J.W.); yi- [email protected] (Y.X.); Email: [email protected] (P.B.M.) This PDF file includes: Materials and Methods References Fig. S1. miR-138 regulates SIN3A in a dose-dependent and site-specific manner. Fig. S2. miR-138 regulates endogenous SIN3A protein expression. Fig. S3. miR-138 regulates endogenous CFTR protein expression in Calu-3 cells. Fig. S4. miR-138 regulates endogenous CFTR protein expression in primary human airway epithelia. Fig. S5. miR-138 regulates CFTR expression in HeLa cells. Fig. S6. miR-138 regulates CFTR expression in HEK293T cells. Fig. S7. HeLa cells exhibit CFTR channel activity. Fig. S8. miR-138 improves CFTR processing. Fig. S9. miR-138 improves CFTR-ΔF508 processing. Fig. S10. SIN3A inhibition yields partial rescue of Cl- transport in CF epithelia.
    [Show full text]
  • Predicting Gene Ontology Biological Process from Temporal Gene Expression Patterns Astrid Lægreid,1,4 Torgeir R
    Methods Predicting Gene Ontology Biological Process From Temporal Gene Expression Patterns Astrid Lægreid,1,4 Torgeir R. Hvidsten,2 Herman Midelfart,2 Jan Komorowski,2,3,4 and Arne K. Sandvik1 1Department of Cancer Research and Molecular Medicine, Norwegian University of Science and Technology, N-7489 Trondheim, Norway; 2Department of Information and Computer Science, Norwegian University of Science and Technology, N-7491 Trondheim, Norway; 3The Linnaeus Centre for Bioinformatics, Uppsala University, SE-751 24 Uppsala, Sweden The aim of the present study was to generate hypotheses on the involvement of uncharacterized genes in biological processes. To this end,supervised learning was used to analyz e microarray-derived time-series gene expression data. Our method was objectively evaluated on known genes using cross-validation and provided high-precision Gene Ontology biological process classifications for 211 of the 213 uncharacterized genes in the data set used. In addition,new roles in biological process were hypothesi zed for known genes. Our method uses biological knowledge expressed by Gene Ontology and generates a rule model associating this knowledge with minimal characteristic features of temporal gene expression profiles. This model allows learning and classification of multiple biological process roles for each gene and can predict participation of genes in a biological process even though the genes of this class exhibit a wide variety of gene expression profiles including inverse coregulation. A considerable number of the hypothesized new roles for known genes were confirmed by literature search. In addition,many biological process roles hypothesi zed for uncharacterized genes were found to agree with assumptions based on homology information.
    [Show full text]
  • Targeting PH Domain Proteins for Cancer Therapy
    The Texas Medical Center Library DigitalCommons@TMC The University of Texas MD Anderson Cancer Center UTHealth Graduate School of The University of Texas MD Anderson Cancer Biomedical Sciences Dissertations and Theses Center UTHealth Graduate School of (Open Access) Biomedical Sciences 12-2018 Targeting PH domain proteins for cancer therapy Zhi Tan Follow this and additional works at: https://digitalcommons.library.tmc.edu/utgsbs_dissertations Part of the Bioinformatics Commons, Medicinal Chemistry and Pharmaceutics Commons, Neoplasms Commons, and the Pharmacology Commons Recommended Citation Tan, Zhi, "Targeting PH domain proteins for cancer therapy" (2018). The University of Texas MD Anderson Cancer Center UTHealth Graduate School of Biomedical Sciences Dissertations and Theses (Open Access). 910. https://digitalcommons.library.tmc.edu/utgsbs_dissertations/910 This Dissertation (PhD) is brought to you for free and open access by the The University of Texas MD Anderson Cancer Center UTHealth Graduate School of Biomedical Sciences at DigitalCommons@TMC. It has been accepted for inclusion in The University of Texas MD Anderson Cancer Center UTHealth Graduate School of Biomedical Sciences Dissertations and Theses (Open Access) by an authorized administrator of DigitalCommons@TMC. For more information, please contact [email protected]. TARGETING PH DOMAIN PROTEINS FOR CANCER THERAPY by Zhi Tan Approval page APPROVED: _____________________________________________ Advisory Professor, Shuxing Zhang, Ph.D. _____________________________________________
    [Show full text]
  • SHROOM3 Is a Novel Candidate for Heterotaxy Identified by Whole Exome Sequencing Tariq Et Al
    SHROOM3 is a novel candidate for heterotaxy identified by whole exome sequencing Tariq et al. Tariq et al. Genome Biology 2011, 12:R91 http://genomebiology.com/2011/12/9/R91 (21 September 2011) Tariq et al. Genome Biology 2011, 12:R91 http://genomebiology.com/2011/12/9/R91 RESEARCH Open Access SHROOM3 is a novel candidate for heterotaxy identified by whole exome sequencing Muhammad Tariq1, John W Belmont2, Seema Lalani2, Teresa Smolarek3 and Stephanie M Ware1,3* Abstract Background: Heterotaxy-spectrum cardiovascular disorders are challenging for traditional genetic analyses because of clinical and genetic heterogeneity, variable expressivity, and non-penetrance. In this study, high-resolution SNP genotyping and exon-targeted array comparative genomic hybridization platforms were coupled to whole-exome sequencing to identify a novel disease candidate gene. Results: SNP genotyping identified absence-of-heterozygosity regions in the heterotaxy proband on chromosomes 1, 4, 7, 13, 15, 18, consistent with parental consanguinity. Subsequently, whole-exome sequencing of the proband identified 26,065 coding variants, including 18 non-synonymous homozygous changes not present in dbSNP132 or 1000 Genomes. Of these 18, only 4 - one each in CXCL2, SHROOM3, CTSO, RXFP1 - were mapped to the absence-of- heterozygosity regions, each of which was flanked by more than 50 homozygous SNPs, confirming recessive segregation of mutant alleles. Sanger sequencing confirmed the SHROOM3 homozygous missense mutation and it was predicted as pathogenic by four bioinformatic tools. SHROOM3 has been identified as a central regulator of morphogenetic cell shape changes necessary for organogenesis and can physically bind ROCK2, a rho kinase protein required for left-right patterning.
    [Show full text]
  • A Novel De Novo 20Q13.32&Ndash;Q13.33
    Journal of Human Genetics (2015) 60, 313–317 & 2015 The Japan Society of Human Genetics All rights reserved 1434-5161/15 www.nature.com/jhg ORIGINAL ARTICLE Anovelde novo 20q13.32–q13.33 deletion in a 2-year-old child with poor growth, feeding difficulties and low bone mass Meena Balasubramanian1, Edward Atack2, Kath Smith2 and Michael James Parker1 Interstitial deletions of the long arm of chromosome 20 are rarely reported in the literature. We report a 2-year-old child with a 2.6 Mb deletion of 20q13.32–q13.33, detected by microarray-based comparative genomic hybridization, who presented with poor growth, feeding difficulties, abnormal subcutaneous fat distribution with the lack of adipose tissue on clinical examination, facial dysmorphism and low bone mass. This report adds to rare publications describing constitutional aberrations of chromosome 20q, and adds further evidence to the fact that deletion of the GNAS complex may not always be associated with an Albright’s hereditary osteodystrophy phenotype as described previously. Journal of Human Genetics (2015) 60, 313–317; doi:10.1038/jhg.2015.22; published online 12 March 2015 INTRODUCTION resuscitation immediately after birth and Apgar scores were 9 and 9 at 1 and Reports of isolated subtelomeric deletions of the long arm of 10 min, respectively, of age. Birth parameters were: weight ~ 1.56 kg (0.4th–2nd chromosome 20 are rare, but a few cases have been reported in the centile), length ~ 40 cm (o0.4th centile) and head circumference ~ 28.2 cm o fi literature over the past 30 years.1–13 Traylor et al.12 provided an ( 0.4th centile).
    [Show full text]
  • Direct in Vivo Mapping of Functional Suppressors in Glioblastoma Genome
    bioRxiv preprint doi: https://doi.org/10.1101/153460; this version posted June 22, 2017. 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. Direct in vivo mapping of functional suppressors in glioblastoma genome Ryan D. Chow *,1,2,3, Christopher D. Guzman*,1,2,4,5,6, Guangchuan Wang *,1,2, Florian Schmidt *,7,8, Mark W. Youngblood **,1,3,9, Lupeng Ye **,1,2, Youssef Errami 1,2, Matthew B. Dong 1,2,3, Michael A. Martinez 1,2, Sensen Zhang 1,2, Paul Renauer, Kaya Bilguvar 1,10, Murat Gunel 1,3,9,10, Phillip A. Sharp 11,12, Feng Zhang 13,14, Randall J. Platt 7,8,#, and Sidi Chen 1,2,3,4,5,15,16# Affiliations 1. Department of Genetics 2. System Biology Institute 3. MD-PhD Program 4. Biological and Biomedical Sciences Program 5. Immunobiology Program 6. Department of Immunobiology Yale University School of Medicine 333 Cedar Street, SHM I-308, New Haven, CT 06520, USA 7. Department of Biosystems Science and Engineering, ETH Zurich Mattenstrasse 26, 4058 Basel, Switzerland 8. Department of Chemistry, University of Basel Petersplatz 1, 4003 Basel, Switzerland 9. Department of Neurosurgery 10. Yale Center for Genome Analysis Yale University School of Medicine 300 Cedar Street, New Haven, CT 06520-8043 11. Koch Institute for Integrative Cancer Research, MIT 12. Department of Biology, MIT 77 Massachusetts Avenue, Cambridge, MA 02139-4307, USA 13. Broad Institute of MIT and Harvard 14. Department of Biological Engineering, MIT 75 Ames Street, Cambridge, MA 02142, USA 15.
    [Show full text]
  • Mouse Osbpl8 Conditional Knockout Project (CRISPR/Cas9)
    https://www.alphaknockout.com Mouse Osbpl8 Conditional Knockout Project (CRISPR/Cas9) Objective: To create a Osbpl8 conditional knockout Mouse model (C57BL/6J) by CRISPR/Cas-mediated genome engineering. Strategy summary: The Osbpl8 gene (NCBI Reference Sequence: NM_175489 ; Ensembl: ENSMUSG00000020189 ) is located on Mouse chromosome 10. 24 exons are identified, with the ATG start codon in exon 2 and the TAG stop codon in exon 24 (Transcript: ENSMUST00000105275). Exon 5 will be selected as conditional knockout region (cKO region). Deletion of this region should result in the loss of function of the Mouse Osbpl8 gene. To engineer the targeting vector, homologous arms and cKO region will be generated by PCR using BAC clone RP23-456P17 as template. Cas9, gRNA and targeting vector will be co-injected into fertilized eggs for cKO Mouse production. The pups will be genotyped by PCR followed by sequencing analysis. Note: Mice homozygous for a gene trap allele exhibit elevated of HDL and gender-specific alterations in lipid metabolism. Exon 5 starts from about 8.17% of the coding region. The knockout of Exon 5 will result in frameshift of the gene. The size of intron 4 for 5'-loxP site insertion: 21997 bp, and the size of intron 5 for 3'-loxP site insertion: 7421 bp. The size of effective cKO region: ~571 bp. The cKO region does not have any other known gene. Page 1 of 7 https://www.alphaknockout.com Overview of the Targeting Strategy Wildtype allele gRNA region 5' gRNA region 3' 1 5 24 Targeting vector Targeted allele Constitutive KO allele (After Cre recombination) Legends Exon of mouse Osbpl8 Homology arm cKO region loxP site Page 2 of 7 https://www.alphaknockout.com Overview of the Dot Plot Window size: 10 bp Forward Reverse Complement Sequence 12 Note: The sequence of homologous arms and cKO region is aligned with itself to determine if there are tandem repeats.
    [Show full text]