A Worldwide Map of Swine Short Tandem Repeats and Their

Total Page:16

File Type:pdf, Size:1020Kb

A Worldwide Map of Swine Short Tandem Repeats and Their Wu et al. Genet Sel Evol (2021) 53:39 https://doi.org/10.1186/s12711-021-00631-4 Genetics Selection Evolution RESEARCH ARTICLE Open Access A worldwide map of swine short tandem repeats and their associations with evolutionary and environmental adaptations Zhongzi Wu1, Huanfa Gong1, Mingpeng Zhang1, Xinkai Tong1, Huashui Ai1, Shijun Xiao1, Miguel Perez‑Enciso2,3, Bin Yang1* and Lusheng Huang1* Abstract Background: Short tandem repeats (STRs) are genetic markers with a greater mutation rate than single nucleotide polymorphisms (SNPs) and are widely used in genetic studies and forensics. However, most studies in pigs have focused only on SNPs or on a limited number of STRs. Results: This study screened 394 deep‑sequenced genomes from 22 domesticated pig breeds/populations world‑ wide, wild boars from both Europe and Asia, and numerous outgroup Suidaes, and identifed a set of 878,967 poly‑ morphic STRs (pSTRs), which represents the largest repository of pSTRs in pigs to date. We found multiple lines of evidence that pSTRs in coding regions were afected by purifying selection. The enrichment of trinucleotide pSTRs in coding sequences (CDS), 5′UTR and H3K4me3 regions suggests that trinucleotide STRs serve as important com‑ ponents in the exons and promoters of the corresponding genes. We demonstrated that, compared to SNPs, pSTRs provide comparable or even greater accuracy in determining the breed identity of individuals. We identifed pSTRs that showed signifcant population diferentiation between domestic pigs and wild boars in Asia and Europe. We also observed that some pSTRs were signifcantly associated with environmental variables, such as average annual tem‑ perature or altitude of the originating sites of Chinese indigenous breeds, among which we identifed loss‑of‑function and/or expanded STRs overlapping with genes such as AHR, LAS1L and PDK1. Finally, our results revealed that several pSTRs show stronger signals in domestic pig—wild boar diferentiation or association with the analysed environmen‑ tal variables than the fanking SNPs within a 100‑kb window. Conclusions: This study provides a genome‑wide high‑density map of pSTRs in diverse pig populations based on genome sequencing data, enabling a more comprehensive characterization of their roles in evolutionary and environ‑ mental adaptation. Background animals [1], and they account for 3 and 3.4% of the human Short tandem repeats (STRs), also known as microsat- and mouse genomes, respectively [2]. STRs stem mainly ellite DNA loci, are defned as tandem repetitive DNA from DNA polymerase slippage events during DNA rep- elements with core unit lengths of one to six base pairs. lication [3]. In humans, the mutation rate of STRs ranges STRs are largely present in the genomes of plants and from 10–8 to 10–2 per locus per generation [4], which is much higher than the mutation rate of single nucleotide polymorphisms (SNPs) (10–9 to 10–8) [5]. Owing to their *Correspondence: [email protected]; [email protected] high genotype diversity and broad distribution across the 1 State Key Laboratory for Pig Genetic Improvement and Production Technology, Jiangxi Agricultural University, Nanchang, China genome, STRs have been widely used in genetic mapping Full list of author information is available at the end of the article of complex diseases [6], forensics [7], population genetics © The Author(s) 2021. 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. Wu et al. Genet Sel Evol (2021) 53:39 Page 2 of 16 [8], linkage analysis [9], and cell line identifcation [10]. pSTRs substantially augment the pSTRs landscape in In humans, over 40 Mendelian diseases are caused by pigs. We surveyed the pSTRs for their genomic distri- STRs, many of which are trinucleotide repeats. Tese butions by diferent unit lengths, genetic diversities in diseases include the fragile X syndrome, which is associ- diferent populations, and their accuracy in classifying ated with amplifcation of the CGG motif in the 5′ UTR individuals by breed identity. We annotated loss-of-func- of the FMR1 gene [11], and Huntington’s disease, which tion and expanded pSTRs. We also identifed a number is linked to CAG repeat expansion in the frst exon of the of pSTRs that showed signifcant diferentiation between HTT gene [12]. Studies have also shown important roles domesticated pigs and wild boars in Asia and Europe for STRs in the regulation of gene expression traits [13, and pSTRs that are associated with environmental vari- 14]. ables such as average annual temperature or altitude in Population genetics studies using dozens to hundreds Chinese indigenous pigs, and compared the signals of of STRs have been carried out since the early 1990s [15]. population diferentiation or environmental adaptation- Recently developed software tools such as lobSTR [16], related pSTRs with those of SNPs called from the same STR-FM [17] and HipSTR [18] have enabled accurate dataset. genotyping of STRs from genome sequence data and their characterization in terms of their association with Methods evolutionary and environmental adaptations and with Samples complex traits in cohorts of various species including Sequence data of 394 pigs or other Sus species were used humans [19], macaques [20], cattle [21] and ducks [22]. to genotype the genome-wide STRs, of which 241 were For example, Xu et al. [21] generated the frst high-den- generated by our laboratory, and the rest were down- sity bovine STR map and identifed STRs associated with loaded from public databases (see Additional fle 1: dairy traits. Fan and his colleagues [22] found that STRs Table S1 and Fig. 1). Te 394 individuals comprise 24 that diferentiated Pekin ducks from mallards are associ- Asian wild boars from South China and Korea; 157 ated with energy metabolism and nervous system genes. Chinese domesticated pigs representing 15 indigenous Forman et al. [23] observed that a (GAA)n repeat expan- breeds, including Bamaxiang (BMX), Luchuan (LC), sion in intron 35 of the ITPR1 gene potentially causes Wuzhishan (WZS), Laiwu (LWH), Hetao (HT), Min spinocerebellar ataxia in Italian spinone dogs. (MIN), Bamei (BM), Baoshan (BS), Neijiang (NJ), Jinhua Pigs are among the most important domestic animals. (JH), Erhualian (EHL), Yunnan Tibetan (YT), Sichuan To date, most of the population genetics studies that Tibetan (ST), Gansu Tibetan (GT) and Tibet Tibetan have been performed in pigs have been based on SNPs (TT); 21 European wild boars; 43 European domesti- called from high-density SNP arrays [24] or second- cated pigs; 130 international commercial pigs, including generation sequencing data [25, 26], or on a small num- Duroc (DU), Landrace (LR), Pietrain (PI), Hampshire ber (from dozens to hundreds) of STRs genotyped with (HP) and Large White (LW) breeds, and 26 outgroup the Sanger sequencing technology [27]. A recent study species, including seven Sus cebifrons, two Sus verru- identifed approximately 16,000 high-quality polymor- cosus, one Sus celebensis, one Sus barbatus, six Porcula phic STRs in sequence data from 102 pigs by referring salvania, one Phacochoerus africanus and one Babyrousa to the Sscrofa10.2 genome assembly and characterized babyrussa. their genetic diversity among diferent breeds [28]. How- ever, this study included a limited number of samples Annotation of repeat sequences in the reference genome and breeds, which restricted its ability to investigate We downloaded the Sscrofa11.1 reference genome the diversity of STRs and their evolutionary roles in sequences from ENSEMBL (ftp:// ftp. ensem bl. org/ pub/ pigs. Moreover, the reference genome of pigs has been relea se- 95/ fasta/ sus_ scrofa/). Te Tandem Repeat Finder updated to Sscrofa11.1, which is a much better assem- (TRF v4.09) software [30] was used to search for poten- bly than Sscrofa10.2 [29]. Terefore, characterizing the tial STR loci in the reference genome with the following STRs in larger cohorts and using the latest porcine refer- command: “${trf} genome.fa 2 7 7 80 10 20 6 -d–h”. All ence genome would be helpful to better understand the overlapping and equivocal STRs were excluded. Repeat- genetic diversity, evolutionary and environmental adap- Masker and RepeatModeler were used to detect other tations of STRs in pigs. repeat elements, including short interspersed nucleo- In this study, we investigated polymorphic STRs tide elements (SINE), long interspersed nucleotide ele- (pSTRs) in 330 domestic pigs from 22 indigenous pig ments (LINE), long terminal repeat retrotransposons populations, 24 Asian and 21 European wild boars, and (LTR), and DNA transposable elements (DNA_TE). 19 individuals from other Suidae species based on the Briefy, RepeatMasker (v4.0.7) [31] was used to search for Sscrofa11.1 genome assembly (Fig. 1). Te identifed repeat sequences by BLAST search of the existing repeat Wu et al. Genet Sel Evol (2021) 53:39 Page 3 of 16 Fig. 1 Geographical distributions and STR‑based neighbor joining tree on 394 samples investigated in this study.
Recommended publications
  • Figure S1. Representative Report Generated by the Ion Torrent System Server for Each of the KCC71 Panel Analysis and Pcafusion Analysis
    Figure S1. Representative report generated by the Ion Torrent system server for each of the KCC71 panel analysis and PCaFusion analysis. (A) Details of the run summary report followed by the alignment summary report for the KCC71 panel analysis sequencing. (B) Details of the run summary report for the PCaFusion panel analysis. A Figure S1. Continued. Representative report generated by the Ion Torrent system server for each of the KCC71 panel analysis and PCaFusion analysis. (A) Details of the run summary report followed by the alignment summary report for the KCC71 panel analysis sequencing. (B) Details of the run summary report for the PCaFusion panel analysis. B Figure S2. Comparative analysis of the variant frequency found by the KCC71 panel and calculated from publicly available cBioPortal datasets. For each of the 71 genes in the KCC71 panel, the frequency of variants was calculated as the variant number found in the examined cases. Datasets marked with different colors and sample numbers of prostate cancer are presented in the upper right. *Significantly high in the present study. Figure S3. Seven subnetworks extracted from each of seven public prostate cancer gene networks in TCNG (Table SVI). Blue dots represent genes that include initial seed genes (parent nodes), and parent‑child and child‑grandchild genes in the network. Graphical representation of node‑to‑node associations and subnetwork structures that differed among and were unique to each of the seven subnetworks. TCNG, The Cancer Network Galaxy. Figure S4. REVIGO tree map showing the predicted biological processes of prostate cancer in the Japanese. Each rectangle represents a biological function in terms of a Gene Ontology (GO) term, with the size adjusted to represent the P‑value of the GO term in the underlying GO term database.
    [Show full text]
  • Downloaded from [266]
    Patterns of DNA methylation on the human X chromosome and use in analyzing X-chromosome inactivation by Allison Marie Cotton B.Sc., The University of Guelph, 2005 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY in The Faculty of Graduate Studies (Medical Genetics) THE UNIVERSITY OF BRITISH COLUMBIA (Vancouver) January 2012 © Allison Marie Cotton, 2012 Abstract The process of X-chromosome inactivation achieves dosage compensation between mammalian males and females. In females one X chromosome is transcriptionally silenced through a variety of epigenetic modifications including DNA methylation. Most X-linked genes are subject to X-chromosome inactivation and only expressed from the active X chromosome. On the inactive X chromosome, the CpG island promoters of genes subject to X-chromosome inactivation are methylated in their promoter regions, while genes which escape from X- chromosome inactivation have unmethylated CpG island promoters on both the active and inactive X chromosomes. The first objective of this thesis was to determine if the DNA methylation of CpG island promoters could be used to accurately predict X chromosome inactivation status. The second objective was to use DNA methylation to predict X-chromosome inactivation status in a variety of tissues. A comparison of blood, muscle, kidney and neural tissues revealed tissue-specific X-chromosome inactivation, in which 12% of genes escaped from X-chromosome inactivation in some, but not all, tissues. X-linked DNA methylation analysis of placental tissues predicted four times higher escape from X-chromosome inactivation than in any other tissue. Despite the hypomethylation of repetitive elements on both the X chromosome and the autosomes, no changes were detected in the frequency or intensity of placental Cot-1 holes.
    [Show full text]
  • Cellular and Molecular Signatures in the Disease Tissue of Early
    Cellular and Molecular Signatures in the Disease Tissue of Early Rheumatoid Arthritis Stratify Clinical Response to csDMARD-Therapy and Predict Radiographic Progression Frances Humby1,* Myles Lewis1,* Nandhini Ramamoorthi2, Jason Hackney3, Michael Barnes1, Michele Bombardieri1, Francesca Setiadi2, Stephen Kelly1, Fabiola Bene1, Maria di Cicco1, Sudeh Riahi1, Vidalba Rocher-Ros1, Nora Ng1, Ilias Lazorou1, Rebecca E. Hands1, Desiree van der Heijde4, Robert Landewé5, Annette van der Helm-van Mil4, Alberto Cauli6, Iain B. McInnes7, Christopher D. Buckley8, Ernest Choy9, Peter Taylor10, Michael J. Townsend2 & Costantino Pitzalis1 1Centre for Experimental Medicine and Rheumatology, William Harvey Research Institute, Barts and The London School of Medicine and Dentistry, Queen Mary University of London, Charterhouse Square, London EC1M 6BQ, UK. Departments of 2Biomarker Discovery OMNI, 3Bioinformatics and Computational Biology, Genentech Research and Early Development, South San Francisco, California 94080 USA 4Department of Rheumatology, Leiden University Medical Center, The Netherlands 5Department of Clinical Immunology & Rheumatology, Amsterdam Rheumatology & Immunology Center, Amsterdam, The Netherlands 6Rheumatology Unit, Department of Medical Sciences, Policlinico of the University of Cagliari, Cagliari, Italy 7Institute of Infection, Immunity and Inflammation, University of Glasgow, Glasgow G12 8TA, UK 8Rheumatology Research Group, Institute of Inflammation and Ageing (IIA), University of Birmingham, Birmingham B15 2WB, UK 9Institute of
    [Show full text]
  • A Database of Gene Regulatory Network Models Across Human Conditions
    bioRxiv preprint doi: https://doi.org/10.1101/2021.06.18.448997; this version posted June 18, 2021. 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 4.0 International license. GRAND: A database of gene regulatory network models across human conditions Marouen Ben Guebila1†, Camila M Lopes-Ramos1†, Deborah Weighill1,7, Abhijeet Rajendra Sonawane2, Rebekka Burkholz1, Behrouz Shamsaei3, John Platig4, Kimberly Glass1,4, Marieke L Kuijjer5,6, John Quackenbush1,4,*. 1Department of Biostatistics, Harvard School of Public Health, Boston, MA, USA. 2Center for Interdisciplinary Cardiovascular Sciences, Cardiovascular Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts 02115, USA. 3Division of Biostatistics and Bioinformatics, Department of Environmental and Public Health Sciences, University of Cincinnati College of Medicine, Cincinnati, Ohio, USA. 4Channing Division of Network Medicine, Department of Medicine, Harvard Medical School and Brigham and Women’s Hospital 5Center for Molecular Medicine Norway, Faculty of Medicine, University of Oslo, Norway 6Leiden University Medical Center, Leiden, the Netherlands 7Current address: UNC Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA † Joint Authors. * To whom correspondence should be addressed. Email: [email protected]. 1 bioRxiv preprint doi: https://doi.org/10.1101/2021.06.18.448997; this version posted June 18, 2021. 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.
    [Show full text]
  • Supplementary Table S4. FGA Co-Expressed Gene List in LUAD
    Supplementary Table S4. FGA co-expressed gene list in LUAD tumors Symbol R Locus Description FGG 0.919 4q28 fibrinogen gamma chain FGL1 0.635 8p22 fibrinogen-like 1 SLC7A2 0.536 8p22 solute carrier family 7 (cationic amino acid transporter, y+ system), member 2 DUSP4 0.521 8p12-p11 dual specificity phosphatase 4 HAL 0.51 12q22-q24.1histidine ammonia-lyase PDE4D 0.499 5q12 phosphodiesterase 4D, cAMP-specific FURIN 0.497 15q26.1 furin (paired basic amino acid cleaving enzyme) CPS1 0.49 2q35 carbamoyl-phosphate synthase 1, mitochondrial TESC 0.478 12q24.22 tescalcin INHA 0.465 2q35 inhibin, alpha S100P 0.461 4p16 S100 calcium binding protein P VPS37A 0.447 8p22 vacuolar protein sorting 37 homolog A (S. cerevisiae) SLC16A14 0.447 2q36.3 solute carrier family 16, member 14 PPARGC1A 0.443 4p15.1 peroxisome proliferator-activated receptor gamma, coactivator 1 alpha SIK1 0.435 21q22.3 salt-inducible kinase 1 IRS2 0.434 13q34 insulin receptor substrate 2 RND1 0.433 12q12 Rho family GTPase 1 HGD 0.433 3q13.33 homogentisate 1,2-dioxygenase PTP4A1 0.432 6q12 protein tyrosine phosphatase type IVA, member 1 C8orf4 0.428 8p11.2 chromosome 8 open reading frame 4 DDC 0.427 7p12.2 dopa decarboxylase (aromatic L-amino acid decarboxylase) TACC2 0.427 10q26 transforming, acidic coiled-coil containing protein 2 MUC13 0.422 3q21.2 mucin 13, cell surface associated C5 0.412 9q33-q34 complement component 5 NR4A2 0.412 2q22-q23 nuclear receptor subfamily 4, group A, member 2 EYS 0.411 6q12 eyes shut homolog (Drosophila) GPX2 0.406 14q24.1 glutathione peroxidase
    [Show full text]
  • Target Gene Gene Description Validation Diana Miranda
    Supplemental Table S1. Mmu-miR-183-5p in silico predicted targets. TARGET GENE GENE DESCRIPTION VALIDATION DIANA MIRANDA MIRBRIDGE PICTAR PITA RNA22 TARGETSCAN TOTAL_HIT AP3M1 adaptor-related protein complex 3, mu 1 subunit V V V V V V V 7 BTG1 B-cell translocation gene 1, anti-proliferative V V V V V V V 7 CLCN3 chloride channel, voltage-sensitive 3 V V V V V V V 7 CTDSPL CTD (carboxy-terminal domain, RNA polymerase II, polypeptide A) small phosphatase-like V V V V V V V 7 DUSP10 dual specificity phosphatase 10 V V V V V V V 7 MAP3K4 mitogen-activated protein kinase kinase kinase 4 V V V V V V V 7 PDCD4 programmed cell death 4 (neoplastic transformation inhibitor) V V V V V V V 7 PPP2R5C protein phosphatase 2, regulatory subunit B', gamma V V V V V V V 7 PTPN4 protein tyrosine phosphatase, non-receptor type 4 (megakaryocyte) V V V V V V V 7 EZR ezrin V V V V V V 6 FOXO1 forkhead box O1 V V V V V V 6 ANKRD13C ankyrin repeat domain 13C V V V V V V 6 ARHGAP6 Rho GTPase activating protein 6 V V V V V V 6 BACH2 BTB and CNC homology 1, basic leucine zipper transcription factor 2 V V V V V V 6 BNIP3L BCL2/adenovirus E1B 19kDa interacting protein 3-like V V V V V V 6 BRMS1L breast cancer metastasis-suppressor 1-like V V V V V V 6 CDK5R1 cyclin-dependent kinase 5, regulatory subunit 1 (p35) V V V V V V 6 CTDSP1 CTD (carboxy-terminal domain, RNA polymerase II, polypeptide A) small phosphatase 1 V V V V V V 6 DCX doublecortin V V V V V V 6 ENAH enabled homolog (Drosophila) V V V V V V 6 EPHA4 EPH receptor A4 V V V V V V 6 FOXP1 forkhead box P1 V
    [Show full text]
  • A High-Throughput Approach to Uncover Novel Roles of APOBEC2, a Functional Orphan of the AID/APOBEC Family
    Rockefeller University Digital Commons @ RU Student Theses and Dissertations 2018 A High-Throughput Approach to Uncover Novel Roles of APOBEC2, a Functional Orphan of the AID/APOBEC Family Linda Molla Follow this and additional works at: https://digitalcommons.rockefeller.edu/ student_theses_and_dissertations Part of the Life Sciences Commons A HIGH-THROUGHPUT APPROACH TO UNCOVER NOVEL ROLES OF APOBEC2, A FUNCTIONAL ORPHAN OF THE AID/APOBEC FAMILY A Thesis Presented to the Faculty of The Rockefeller University in Partial Fulfillment of the Requirements for the degree of Doctor of Philosophy by Linda Molla June 2018 © Copyright by Linda Molla 2018 A HIGH-THROUGHPUT APPROACH TO UNCOVER NOVEL ROLES OF APOBEC2, A FUNCTIONAL ORPHAN OF THE AID/APOBEC FAMILY Linda Molla, Ph.D. The Rockefeller University 2018 APOBEC2 is a member of the AID/APOBEC cytidine deaminase family of proteins. Unlike most of AID/APOBEC, however, APOBEC2’s function remains elusive. Previous research has implicated APOBEC2 in diverse organisms and cellular processes such as muscle biology (in Mus musculus), regeneration (in Danio rerio), and development (in Xenopus laevis). APOBEC2 has also been implicated in cancer. However the enzymatic activity, substrate or physiological target(s) of APOBEC2 are unknown. For this thesis, I have combined Next Generation Sequencing (NGS) techniques with state-of-the-art molecular biology to determine the physiological targets of APOBEC2. Using a cell culture muscle differentiation system, and RNA sequencing (RNA-Seq) by polyA capture, I demonstrated that unlike the AID/APOBEC family member APOBEC1, APOBEC2 is not an RNA editor. Using the same system combined with enhanced Reduced Representation Bisulfite Sequencing (eRRBS) analyses I showed that, unlike the AID/APOBEC family member AID, APOBEC2 does not act as a 5-methyl-C deaminase.
    [Show full text]
  • Genome Architecture and Transcription Data Reveal Allelic Bias During the Cell Cycle
    bioRxiv preprint doi: https://doi.org/10.1101/2020.03.15.992164; this version posted April 1, 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. Genome architecture and transcription data reveal allelic bias during the cell cycle Stephen Lindsly1, Wenlong Jia2, Haiming Chen1, Sijia Liu3, Scott Ronquist1, Can Chen4;7, Xingzhao Wen5, Gilbert Omenn1, Shuai Cheng Li2, Max Wicha1;6, Alnawaz Rehemtulla6, Lindsey Muir1, Indika Rajapakse1;7;∗ 1Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor 2Department of Computer Science, City University of Hong Kong 3MIT-IBM Watson AI Lab, IBM Research 4Department of Electrical Engineering and Computer Science, University of Michigan, Ann Arbor 5Department of Bioengineering, University of California San Diego, La Jolla 6Department of Hematology/Oncology, University of Michigan, Ann Arbor 7Department of Mathematics, University of Michigan, Ann Arbor ∗To whom correspondence should be addressed; E-mail: [email protected]. bioRxiv preprint doi: https://doi.org/10.1101/2020.03.15.992164; this version posted April 1, 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. Highlights • Structural and functional differences between the maternal and paternal genomes, includ- ing similar allelic bias in related subsets of genes. • Coupling between the dynamics of gene expression and genome architecture for specific alleles illuminates an allele-specific 4D Nucleome. • Introduction of a novel allele-specific phasing algorithm for genome architecture and a quantitative framework for integration of gene expression and genome architecture.
    [Show full text]
  • Engineered Type 1 Regulatory T Cells Designed for Clinical Use Kill Primary
    ARTICLE Acute Myeloid Leukemia Engineered type 1 regulatory T cells designed Ferrata Storti Foundation for clinical use kill primary pediatric acute myeloid leukemia cells Brandon Cieniewicz,1* Molly Javier Uyeda,1,2* Ping (Pauline) Chen,1 Ece Canan Sayitoglu,1 Jeffrey Mao-Hwa Liu,1 Grazia Andolfi,3 Katharine Greenthal,1 Alice Bertaina,1,4 Silvia Gregori,3 Rosa Bacchetta,1,4 Norman James Lacayo,1 Alma-Martina Cepika1,4# and Maria Grazia Roncarolo1,2,4# Haematologica 2021 Volume 106(10):2588-2597 1Department of Pediatrics, Division of Stem Cell Transplantation and Regenerative Medicine, Stanford School of Medicine, Stanford, CA, USA; 2Stanford Institute for Stem Cell Biology and Regenerative Medicine, Stanford School of Medicine, Stanford, CA, USA; 3San Raffaele Telethon Institute for Gene Therapy, Milan, Italy and 4Center for Definitive and Curative Medicine, Stanford School of Medicine, Stanford, CA, USA *BC and MJU contributed equally as co-first authors #AMC and MGR contributed equally as co-senior authors ABSTRACT ype 1 regulatory (Tr1) T cells induced by enforced expression of interleukin-10 (LV-10) are being developed as a novel treatment for Tchemotherapy-resistant myeloid leukemias. In vivo, LV-10 cells do not cause graft-versus-host disease while mediating graft-versus-leukemia effect against adult acute myeloid leukemia (AML). Since pediatric AML (pAML) and adult AML are different on a genetic and epigenetic level, we investigate herein whether LV-10 cells also efficiently kill pAML cells. We show that the majority of primary pAML are killed by LV-10 cells, with different levels of sensitivity to killing. Transcriptionally, pAML sensitive to LV-10 killing expressed a myeloid maturation signature.
    [Show full text]
  • Structure of Voltage-Modulated Sodium-Selective NALCN-FAM155A Channel Complex
    bioRxiv preprint doi: https://doi.org/10.1101/2020.07.26.221747; this version posted July 26, 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-ND 4.0 International license. Structure of voltage-modulated sodium-selective NALCN-FAM155A channel complex Yunlu Kang 1, Jing-Xiang Wu 1,2,3, and Lei Chen 1,2,3* 1 State Key Laboratory of Membrane Biology, Institute of Molecular Medicine, Peking University, Beijing Key Laboratory of Cardiometabolic Molecular Medicine, Beijing 100871, China 2 Peking-Tsinghua Center for Life Sciences, Peking University, Beijing 100871, China 3Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China *Correspondence: Lei Chen, [email protected] bioRxiv preprint doi: https://doi.org/10.1101/2020.07.26.221747; this version posted July 26, 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-ND 4.0 International license. Abstract Resting membrane potential determines the excitability of the cell and is essential for the cellular electrical activities. NALCN channel mediates sodium leak currents, which positively tune the resting membrane potential towards depolarization. NALCN channel is involved in many important neurological processes and is implicated in a spectrum of neurodevelopmental diseases. Despite its functional importance, the mechanisms of ion permeation and voltage-modulation for NALCN channel remain elusive.
    [Show full text]
  • A Temporally Controlled Sequence of X-Chromosome Inactivation and Reactivation Defines Female Mouse in Vitro Germ Cells with Meiotic Potential
    bioRxiv preprint doi: https://doi.org/10.1101/2021.08.11.455976; this version posted August 11, 2021. 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. A temporally controlled sequence of X-chromosome inactivation and reactivation defines female mouse in vitro germ cells with meiotic potential Jacqueline Severino1†, Moritz Bauer1,9†, Tom Mattimoe1, Niccolò Arecco1, Luca Cozzuto1, Patricia Lorden2, Norio Hamada3, Yoshiaki Nosaka4,5,6, So Nagaoka4,5,6, Holger Heyn2, Katsuhiko Hayashi7, Mitinori Saitou4,5,6 and Bernhard Payer1,8* Abstract The early mammalian germ cell lineage is characterized by extensive epigenetic reprogramming, which is required for the maturation into functional eggs and sperm. In particular, the epigenome needs to be reset before parental marks can be established and then transmitted to the next generation. In the female germ line, reactivation of the inactive X- chromosome is one of the most prominent epigenetic reprogramming events, and despite its scale involving an entire chromosome affecting hundreds of genes, very little is known about its kinetics and biological function. Here we investigate X-chromosome inactivation and reactivation dynamics by employing a tailor-made in vitro system to visualize the X-status during differentiation of primordial germ cell-like cells (PGCLCs) from female mouse embryonic stem cells (ESCs). We find that the degree of X-inactivation in PGCLCs is moderate when compared to somatic cells and characterized by a large number of genes escaping full inactivation.
    [Show full text]
  • Structure of Voltage-Modulated Sodium-Selective NALCN-FAM155A Channel Complex ✉ Yunlu Kang 1, Jing-Xiang Wu1,2,3 & Lei Chen 1,2,3
    ARTICLE https://doi.org/10.1038/s41467-020-20002-9 OPEN Structure of voltage-modulated sodium-selective NALCN-FAM155A channel complex ✉ Yunlu Kang 1, Jing-Xiang Wu1,2,3 & Lei Chen 1,2,3 Resting membrane potential determines the excitability of the cell and is essential for the cellular electrical activities. The NALCN channel mediates sodium leak currents, which positively adjust resting membrane potential towards depolarization. The NALCN channel is 1234567890():,; involved in several neurological processes and has been implicated in a spectrum of neu- rodevelopmental diseases. Here, we report the cryo-EM structure of rat NALCN and mouse FAM155A complex to 2.7 Å resolution. The structure reveals detailed interactions between NALCN and the extracellular cysteine-rich domain of FAM155A. We find that the non- canonical architecture of NALCN selectivity filter dictates its sodium selectivity and calcium block, and that the asymmetric arrangement of two functional voltage sensors confers the modulation by membrane potential. Moreover, mutations associated with human diseases map to the domain-domain interfaces or the pore domain of NALCN, intuitively suggesting their pathological mechanisms. 1 State Key Laboratory of Membrane Biology, College of Future Technology, Institute of Molecular Medicine, Beijing Key Laboratory of Cardiometabolic Molecular Medicine, Peking University, 100871 Beijing, China. 2 Peking-Tsinghua Center for Life Sciences, Peking University, 100871 Beijing, China. ✉ 3 Academy for Advanced Interdisciplinary Studies, Peking University, 100871 Beijing, China. email: [email protected] NATURE COMMUNICATIONS | (2020) 11:6199 | https://doi.org/10.1038/s41467-020-20002-9 | www.nature.com/naturecommunications 1 ARTICLE NATURE COMMUNICATIONS | https://doi.org/10.1038/s41467-020-20002-9 ALCN channel is a voltage-modulated sodium-selective Therefore, we started with the cryo-EM studies of the NALCN– ion channel1.
    [Show full text]