1
|
Wu EY, Singh NP, Choi K, Zakeri M, Vincent M, Churchill GA, Ackert-Bicknell CL, Patro R, Love MI. SEESAW: detecting isoform-level allelic imbalance accounting for inferential uncertainty. Genome Biol 2023; 24:165. [PMID: 37438847 PMCID: PMC10337143 DOI: 10.1186/s13059-023-03003-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2022] [Accepted: 06/29/2023] [Indexed: 07/14/2023] Open
Abstract
Detecting allelic imbalance at the isoform level requires accounting for inferential uncertainty, caused by multi-mapping of RNA-seq reads. Our proposed method, SEESAW, uses Salmon and Swish to offer analysis at various levels of resolution, including gene, isoform, and aggregating isoforms to groups by transcription start site. The aggregation strategies strengthen the signal for transcripts with high uncertainty. The SEESAW suite of methods is shown to have higher power than other allelic imbalance methods when there is isoform-level allelic imbalance. We also introduce a new test for detecting imbalance that varies across a covariate, such as time.
Collapse
Affiliation(s)
- Euphy Y Wu
- Department of Biostatistics, University of North Carolina-Chapel Hill, Chapel Hill, NC, USA
| | - Noor P Singh
- Department of Computer Science, University of Maryland, College Park, MD, USA
| | | | - Mohsen Zakeri
- Department of Computer Science, University of Maryland, College Park, MD, USA
| | | | | | - Cheryl L Ackert-Bicknell
- Department of Orthopedics, School of Medicine, University of Colorado, Anschutz Campus, Aurora, CO, USA
| | - Rob Patro
- Department of Computer Science, University of Maryland, College Park, MD, USA
| | - Michael I Love
- Department of Biostatistics, University of North Carolina-Chapel Hill, Chapel Hill, NC, USA.
- Department of Genetics, University of North Carolina-Chapel Hill, Chapel Hill, NC, USA.
| |
Collapse
|
2
|
Fraser HB. Existing methods are effective at measuring natural selection on gene expression. Nat Ecol Evol 2022; 6:1836-1837. [PMID: 36344679 DOI: 10.1038/s41559-022-01889-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2022] [Accepted: 08/17/2022] [Indexed: 11/09/2022]
Affiliation(s)
- Hunter B Fraser
- Department of Biology, Stanford University, Stanford, CA, USA.
| |
Collapse
|
3
|
Conrad T, Altmüller J. Single cell- and spatial 'Omics revolutionize physiology. Acta Physiol (Oxf) 2022; 235:e13848. [PMID: 35656634 DOI: 10.1111/apha.13848] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2022] [Revised: 04/24/2022] [Accepted: 05/27/2022] [Indexed: 11/29/2022]
Abstract
Single cell multi- 'Omics and Spatial Transcriptomics are prominent technological highlights of recent years, and both fields still witness a ceaseless firework of novel approaches for high resolution profiling of additional omics layers. As all life processes in organs and organisms are based on the functions of their fundamental building blocks, the individual cells and their interactions, these methods are of utmost worth for the study of physiology in health and disease. Recent discoveries on embryonic development, tumor immunology, the detailed cellular composition and function of complex tissues like for example the kidney or the brain, different roles of the same cell type in different organs, the oncogenic program of individual tumor entities, or the architecture of immunopathology in infected tissue are based on single cell and spatial transcriptomics experiments. In this review, we will give a broad overview of technological concepts for single cell and spatial analysis, showing both advantages and limitations, and illustrate their impact with some particularly impressive case studies.
Collapse
Affiliation(s)
- Thomas Conrad
- Genomics Technology Platform Max Delbrück Center for Molecular Medicine in the Helmholtz Association (MDC) Berlin Germany
| | - Janine Altmüller
- Genomics Technology Platform Max Delbrück Center for Molecular Medicine in the Helmholtz Association (MDC) Berlin Germany
- Core Facility Genomics Berlin Institute of Health at Charité ‐ Universitätsmedizin Berlin Berlin Germany
- Center for Molecular Medicine Cologne (CMMC) Cologne Germany
| |
Collapse
|
4
|
Mu W, Sarkar H, Srivastava A, Choi K, Patro R, Love MI. Airpart: interpretable statistical models for analyzing allelic imbalance in single-cell datasets. Bioinformatics 2022; 38:2773-2780. [PMID: 35561168 PMCID: PMC9113279 DOI: 10.1093/bioinformatics/btac212] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2021] [Revised: 03/05/2022] [Accepted: 04/05/2022] [Indexed: 02/03/2023] Open
Abstract
MOTIVATION Allelic expression analysis aids in detection of cis-regulatory mechanisms of genetic variation, which produce allelic imbalance (AI) in heterozygotes. Measuring AI in bulk data lacking time or spatial resolution has the limitation that cell-type-specific (CTS), spatial- or time-dependent AI signals may be dampened or not detected. RESULTS We introduce a statistical method airpart for identifying differential CTS AI from single-cell RNA-sequencing data, or dynamics AI from other spatially or time-resolved datasets. airpart outputs discrete partitions of data, pointing to groups of genes and cells under common mechanisms of cis-genetic regulation. In order to account for low counts in single-cell data, our method uses a Generalized Fused Lasso with Binomial likelihood for partitioning groups of cells by AI signal, and a hierarchical Bayesian model for AI statistical inference. In simulation, airpart accurately detected partitions of cell types by their AI and had lower Root Mean Square Error (RMSE) of allelic ratio estimates than existing methods. In real data, airpart identified differential allelic imbalance patterns across cell states and could be used to define trends of AI signal over spatial or time axes. AVAILABILITY AND IMPLEMENTATION The airpart package is available as an R/Bioconductor package at https://bioconductor.org/packages/airpart. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Wancen Mu
- To whom correspondence should be addressed. or
| | - Hirak Sarkar
- Department of Biomedical Informatics, Harvard Medical School, Boston, MA 02115, USA
| | | | | | - Rob Patro
- Department of Computer Science, University of Maryland, College Park, MD 20742, USA
| | | |
Collapse
|
5
|
Kopania EEK, Larson EL, Callahan C, Keeble S, Good JM. Molecular Evolution across Mouse Spermatogenesis. Mol Biol Evol 2022; 39:6517785. [PMID: 35099536 PMCID: PMC8844503 DOI: 10.1093/molbev/msac023] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Genes involved in spermatogenesis tend to evolve rapidly, but we lack a clear understanding of how protein sequences and patterns of gene expression evolve across this complex developmental process. We used fluorescence-activated cell sorting (FACS) to generate expression data for early (meiotic) and late (postmeiotic) cell types across 13 inbred strains of mice (Mus) spanning ∼7 My of evolution. We used these comparative developmental data to investigate the evolution of lineage-specific expression, protein-coding sequences, and expression levels. We found increased lineage specificity and more rapid protein-coding and expression divergence during late spermatogenesis, suggesting that signatures of rapid testis molecular evolution are punctuated across sperm development. Despite strong overall developmental parallels in these components of molecular evolution, protein and expression divergences were only weakly correlated across genes. We detected more rapid protein evolution on the X chromosome relative to the autosomes, whereas X-linked gene expression tended to be relatively more conserved likely reflecting chromosome-specific regulatory constraints. Using allele-specific FACS expression data from crosses between four strains, we found that the relative contributions of different regulatory mechanisms also differed between cell types. Genes showing cis-regulatory changes were more common late in spermatogenesis, and tended to be associated with larger differences in expression levels and greater expression divergence between species. In contrast, genes with trans-acting changes were more common early and tended to be more conserved across species. Our findings advance understanding of gene evolution across spermatogenesis and underscore the fundamental importance of developmental context in molecular evolutionary studies.
Collapse
Affiliation(s)
- Emily E K Kopania
- Division of Biological Sciences, University of Montana, Missoula, MT, 59812, USA
| | - Erica L Larson
- Department of Biological Sciences, University of Denver, Denver, CO, 80208, USA
| | - Colin Callahan
- Division of Biological Sciences, University of Montana, Missoula, MT, 59812, USA
| | - Sara Keeble
- Division of Biological Sciences, University of Montana, Missoula, MT, 59812, USA
| | - Jeffrey M Good
- Division of Biological Sciences, University of Montana, Missoula, MT, 59812, USA
| |
Collapse
|
6
|
Benowitz KM, Coleman JM, Allan CW, Matzkin LM. Contributions of cis- and trans-Regulatory Evolution to Transcriptomic Divergence across Populations in the Drosophila mojavensis Larval Brain. Genome Biol Evol 2021; 12:1407-1418. [PMID: 32653899 PMCID: PMC7495911 DOI: 10.1093/gbe/evaa145] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/06/2020] [Indexed: 12/22/2022] Open
Abstract
Natural selection on gene expression was originally predicted to result primarily in cis- rather than trans-regulatory evolution, due to the expectation of reduced pleiotropy. Despite this, numerous studies have ascribed recent evolutionary divergence in gene expression predominantly to trans-regulation. Performing RNA-seq on single isofemale lines from genetically distinct populations of the cactophilic fly Drosophila mojavensis and their F1 hybrids, we recapitulated this pattern in both larval brains and whole bodies. However, we demonstrate that improving the measurement of brain expression divergence between populations by using seven additional genotypes considerably reduces the estimate of trans-regulatory contributions to expression evolution. We argue that the finding of trans-regulatory predominance can result from biases due to environmental variation in expression or other sources of noise, and that cis-regulation is likely a greater contributor to transcriptional evolution across D. mojavensis populations. Lastly, we merge these lines of data to identify several previously hypothesized and intriguing novel candidate genes, and suggest that the integration of regulatory and population-level transcriptomic data can provide useful filters for the identification of potentially adaptive genes.
Collapse
Affiliation(s)
| | - Joshua M Coleman
- Department of Entomology, University of Arizona.,Department of Biological Sciences, University of Alabama in Huntsville
| | | | - Luciano M Matzkin
- Department of Entomology, University of Arizona.,Department of Ecology and Evolutionary Biology, University of Arizona.,BIO5 Institute, University of Arizona
| |
Collapse
|
7
|
Agoglia RM, Sun D, Birey F, Yoon SJ, Miura Y, Sabatini K, Pașca SP, Fraser HB. Primate cell fusion disentangles gene regulatory divergence in neurodevelopment. Nature 2021; 592:421-427. [PMID: 33731928 PMCID: PMC8719633 DOI: 10.1038/s41586-021-03343-3] [Citation(s) in RCA: 39] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2020] [Accepted: 02/08/2021] [Indexed: 12/31/2022]
Abstract
Among primates, humans display a unique trajectory of development that is responsible for the many traits specific to our species. However, the inaccessibility of primary human and chimpanzee tissues has limited our ability to study human evolution. Comparative in vitro approaches using primate-derived induced pluripotent stem cells have begun to reveal species differences on the cellular and molecular levels1,2. In particular, brain organoids have emerged as a promising platform to study primate neural development in vitro3-5, although cross-species comparisons of organoids are complicated by differences in developmental timing and variability of differentiation6,7. Here we develop a new platform to address these limitations by fusing human and chimpanzee induced pluripotent stem cells to generate a panel of tetraploid hybrid stem cells. We applied this approach to study species divergence in cerebral cortical development by differentiating these cells into neural organoids. We found that hybrid organoids provide a controlled system for disentangling cis- and trans-acting gene-expression divergence across cell types and developmental stages, revealing a signature of selection on astrocyte-related genes. In addition, we identified an upregulation of the human somatostatin receptor 2 gene (SSTR2), which regulates neuronal calcium signalling and is associated with neuropsychiatric disorders8,9. We reveal a human-specific response to modulation of SSTR2 function in cortical neurons, underscoring the potential of this platform for elucidating the molecular basis of human evolution.
Collapse
Affiliation(s)
- Rachel M Agoglia
- Department of Genetics, Stanford University School of Medicine, Stanford, CA, USA
| | - Danqiong Sun
- Department of Biology, Stanford University, Stanford, CA, USA
| | - Fikri Birey
- Department of Psychiatry and Behavioral Sciences, Stanford University School of Medicine, Stanford, CA, USA
| | - Se-Jin Yoon
- Department of Psychiatry and Behavioral Sciences, Stanford University School of Medicine, Stanford, CA, USA
| | - Yuki Miura
- Department of Psychiatry and Behavioral Sciences, Stanford University School of Medicine, Stanford, CA, USA
| | - Karen Sabatini
- Department of Psychiatry and Behavioral Sciences, Stanford University School of Medicine, Stanford, CA, USA
| | - Sergiu P Pașca
- Department of Psychiatry and Behavioral Sciences, Stanford University School of Medicine, Stanford, CA, USA.
| | - Hunter B Fraser
- Department of Biology, Stanford University, Stanford, CA, USA.
| |
Collapse
|
8
|
Liao J, Lu X, Shao X, Zhu L, Fan X. Uncovering an Organ's Molecular Architecture at Single-Cell Resolution by Spatially Resolved Transcriptomics. Trends Biotechnol 2020; 39:43-58. [PMID: 32505359 DOI: 10.1016/j.tibtech.2020.05.006] [Citation(s) in RCA: 119] [Impact Index Per Article: 29.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2019] [Revised: 05/11/2020] [Accepted: 05/12/2020] [Indexed: 01/17/2023]
Abstract
Revealing fine-scale cellular heterogeneity among spatial context and the functional and structural foundations of tissue architecture is fundamental within biological research and pharmacology. Unlike traditional approaches involving single molecules or bulk omics, cutting-edge, spatially resolved transcriptomics techniques offer near-single-cell or even subcellular resolution within tissues. Massive information across higher dimensions along with position-coordinating labels can better map the whole 3D transcriptional landscape of tissues. In this review, we focus on developments and strategies in spatially resolved transcriptomics, compare the cell and gene throughput and spatial resolution in detail for existing methods, and highlight the enormous potential in biomedical research.
Collapse
Affiliation(s)
- Jie Liao
- Pharmaceutical Informatics Institute, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Xiaoyan Lu
- Pharmaceutical Informatics Institute, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Xin Shao
- Pharmaceutical Informatics Institute, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Ling Zhu
- The Save Sight Institute, Faculty of Medicine and Health, the University of Sydney, Sydney, NSW 2000, Australia
| | - Xiaohui Fan
- Pharmaceutical Informatics Institute, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China; The Save Sight Institute, Faculty of Medicine and Health, the University of Sydney, Sydney, NSW 2000, Australia.
| |
Collapse
|
9
|
The 3D Genome Shapes the Regulatory Code of Developmental Genes. J Mol Biol 2020; 432:712-723. [DOI: 10.1016/j.jmb.2019.10.017] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2019] [Revised: 10/11/2019] [Accepted: 10/24/2019] [Indexed: 02/06/2023]
|
10
|
Kalay G, Atallah J, Sierra NC, Tang AM, Crofton AE, Murugesan MK, Wykoff-Clary S, Lott SE. Evolution of larval segment position across 12 Drosophila species. Evolution 2019; 74:1409-1422. [PMID: 31886902 PMCID: PMC7496318 DOI: 10.1111/evo.13911] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2019] [Accepted: 12/12/2019] [Indexed: 12/25/2022]
Abstract
Many developmental traits that are critical to the survival of the organism are also robust. These robust traits are resistant to phenotypic change in the face of variation. This presents a challenge to evolution. In this article, we asked whether and how a well‐established robust trait, Drosophila segment patterning, changed over the evolutionary history of the genus. We compared segment position scaled to body length at the first‐instar larval stage among 12 Drosophila species. We found that relative segment position has changed many times across the phylogeny. Changes were frequent, but primarily small in magnitude. Phylogenetic analysis demonstrated that rates of change in segment position are variable along the Drosophila phylogenetic tree, and that these changes can occur in short evolutionary timescales. Correlation between position shifts of segments decreased as the distance between two segments increased, suggesting local control of segment position. The posterior‐most abdominal segment showed the highest magnitude of change on average, had the highest rate of evolution between species, and appeared to be evolving more independently as compared to the rest of the segments. This segment was exceptionally elongated in the cactophilic species in our dataset, raising questions as to whether this change may be adaptive.
Collapse
Affiliation(s)
- Gizem Kalay
- Department of Evolution and Ecology, University of California, Davis, One Shields Avenue, Davis, California, 95616
| | - Joel Atallah
- Department of Evolution and Ecology, University of California, Davis, One Shields Avenue, Davis, California, 95616.,current address: Department of Biological Sciences, University of New Orleans, 2000 Lakeshore Drive, New Orleans, LA, 70148
| | - Noemie C Sierra
- Earth and Planetary Sciences Department, University of California, Davis, One Shields Avenue, Davis, California, 95616
| | - Austin M Tang
- Department of Evolution and Ecology, University of California, Davis, One Shields Avenue, Davis, California, 95616
| | - Amanda E Crofton
- Department of Evolution and Ecology, University of California, Davis, One Shields Avenue, Davis, California, 95616
| | - Mohan K Murugesan
- Department of Evolution and Ecology, University of California, Davis, One Shields Avenue, Davis, California, 95616
| | - Sherri Wykoff-Clary
- Department of Evolution and Ecology, University of California, Davis, One Shields Avenue, Davis, California, 95616
| | - Susan E Lott
- Department of Evolution and Ecology, University of California, Davis, One Shields Avenue, Davis, California, 95616
| |
Collapse
|
11
|
Colgan W, Leanza A, Hwang A, DeBiasse MB, Llosa I, Rodrigues D, Adhikari H, Barreto Corona G, Bock S, Carillo-Perez A, Currie M, Darkoa-Larbi S, Dellal D, Gutow H, Hokama P, Kibby E, Linhart N, Moody S, Naganuma A, Nguyen D, Stanton R, Stark S, Tumey C, Velleca A, Ryan JF, Davidson B. Variable levels of drift in tunicate cardiopharyngeal gene regulatory elements. EvoDevo 2019; 10:24. [PMID: 31632631 PMCID: PMC6790052 DOI: 10.1186/s13227-019-0137-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2019] [Accepted: 09/13/2019] [Indexed: 12/29/2022] Open
Abstract
BACKGROUND Mutations in gene regulatory networks often lead to genetic divergence without impacting gene expression or developmental patterning. The rules governing this process of developmental systems drift, including the variable impact of selective constraints on different nodes in a gene regulatory network, remain poorly delineated. RESULTS Here we examine developmental systems drift within the cardiopharyngeal gene regulatory networks of two tunicate species, Corella inflata and Ciona robusta. Cross-species analysis of regulatory elements suggests that trans-regulatory architecture is largely conserved between these highly divergent species. In contrast, cis-regulatory elements within this network exhibit distinct levels of conservation. In particular, while most of the regulatory elements we analyzed showed extensive rearrangements of functional binding sites, the enhancer for the cardiopharyngeal transcription factor FoxF is remarkably well-conserved. Even minor alterations in spacing between binding sites lead to loss of FoxF enhancer function, suggesting that bound trans-factors form position-dependent complexes. CONCLUSIONS Our findings reveal heterogeneous levels of divergence across cardiopharyngeal cis-regulatory elements. These distinct levels of divergence presumably reflect constraints that are not clearly associated with gene function or position within the regulatory network. Thus, levels of cis-regulatory divergence or drift appear to be governed by distinct structural constraints that will be difficult to predict based on network architecture.
Collapse
Affiliation(s)
| | - Alexis Leanza
- Thomas Jefferson University Sidney Kimmel Medical College, Philadelphia, USA
| | - Ariel Hwang
- University of North Carolina, Chapel Hill, USA
| | | | | | | | | | | | | | | | | | | | - Daniel Dellal
- Icahn School of Medicine at Mount Sinai, New York, USA
| | | | | | - Emily Kibby
- University of Colorado Boulder, Boulder, USA
| | | | | | | | | | | | - Sierra Stark
- University of California San Francisco, San Francisco, USA
| | | | | | - Joseph F. Ryan
- Whitney Laboratory for Marine Bioscience, St. Augustine, USA
| | | |
Collapse
|
12
|
Affiliation(s)
- Sofie Y. N. Delbare
- Molecular Biology and Genetics, Cornell University, Ithaca, New York, United States of America
| | - Andrew G. Clark
- Molecular Biology and Genetics, Cornell University, Ithaca, New York, United States of America
- * E-mail:
| |
Collapse
|