1
|
Liu J, Mosti F, Zhao HT, Sotelo-Fonseca JE, Escobar-Tomlienovich CF, Lollis D, Musso CM, Mao Y, Massri AJ, Doll HM, Sousa AM, Wray GA, Schmidt E, Silver DL. A human-specific enhancer fine-tunes radial glia potency and corticogenesis. bioRxiv 2024:2024.04.10.588953. [PMID: 38645099 PMCID: PMC11030412 DOI: 10.1101/2024.04.10.588953] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/23/2024]
Abstract
Humans evolved an extraordinarily expanded and complex cerebral cortex, associated with developmental and gene regulatory modifications 1-3 . Human accelerated regions (HARs) are highly conserved genomic sequences with human-specific nucleotide substitutions. Although there are thousands of annotated HARs, their functional contribution to human-specific cortical development is largely unknown 4,5 . HARE5 is a HAR transcriptional enhancer of the WNT signaling receptor Frizzled8 (FZD8) active during brain development 6 . Here, using genome-edited mouse and primate models, we demonstrate that human (Hs) HARE5 fine-tunes cortical development and connectivity by controlling the proliferative and neurogenic capacity of neural progenitor cells (NPCs). Hs-HARE5 knock-in mice have significantly enlarged neocortices containing more neurons. By measuring neural dynamics in vivo we show these anatomical features correlate with increased functional independence between cortical regions. To understand the underlying developmental mechanisms, we assess progenitor fate using live imaging, lineage analysis, and single-cell RNA sequencing. This reveals Hs-HARE5 modifies radial glial progenitor behavior, with increased self-renewal at early developmental stages followed by expanded neurogenic potential. We use genome-edited human and chimpanzee (Pt) NPCs and cortical organoids to assess the relative enhancer activity and function of Hs-HARE5 and Pt-HARE5. Using these orthogonal strategies we show four human-specific variants in HARE5 drive increased enhancer activity which promotes progenitor proliferation. These findings illustrate how small changes in regulatory DNA can directly impact critical signaling pathways and brain development. Our study uncovers new functions for HARs as key regulatory elements crucial for the expansion and complexity of the human cerebral cortex.
Collapse
|
2
|
Gartner V, Redelings BD, Gaither C, Parr JB, Kalonji A, Phanzu F, Brazeau NF, Juliano JJ, Wray GA. Genomic insights into Plasmodium vivax population structure and diversity in central Africa. Malar J 2024; 23:27. [PMID: 38238806 PMCID: PMC10797969 DOI: 10.1186/s12936-024-04852-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2023] [Accepted: 01/12/2024] [Indexed: 01/22/2024] Open
Abstract
BACKGROUND Though Plasmodium vivax is the second most common malaria species to infect humans, it has not traditionally been considered a major human health concern in central Africa given the high prevalence of the human Duffy-negative phenotype that is believed to prevent infection. Increasing reports of asymptomatic and symptomatic infections in Duffy-negative individuals throughout Africa raise the possibility that P. vivax is evolving to evade host resistance, but there are few parasite samples with genomic data available from this part of the world. METHODS Whole genome sequencing of one new P. vivax isolate from the Democratic Republic of the Congo (DRC) was performed and used in population genomics analyses to assess how this central African isolate fits into the global context of this species. RESULTS Plasmodium vivax from DRC is similar to other African populations and is not closely related to the non-human primate parasite P. vivax-like. Evidence is found for a duplication of the gene PvDBP and a single copy of PvDBP2. CONCLUSION These results suggest an endemic P. vivax population is present in central Africa. Intentional sampling of P. vivax across Africa would further contextualize this sample within African P. vivax diversity and shed light on the mechanisms of infection in Duffy negative individuals. These results are limited by the uncertainty of how representative this single sample is of the larger population of P. vivax in central Africa.
Collapse
Affiliation(s)
- Valerie Gartner
- Biology Department, Duke University, Durham, NC, 27708, USA
- University Program in Genetics and Genomics, Duke University, Durham, NC, 27708, USA
| | - Benjamin D Redelings
- Biology Department, Duke University, Durham, NC, 27708, USA
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS, 66045, USA
- Ronin Institute, Durham, NC, 27705, USA
| | | | | | - Albert Kalonji
- SANRU Asbl, 149 A/B, Boulevard du 30 Juin, Kinshasa, Gombe, Democratic Republic of Congo
| | - Fernandine Phanzu
- SANRU Asbl, 149 A/B, Boulevard du 30 Juin, Kinshasa, Gombe, Democratic Republic of Congo
| | | | | | - Gregory A Wray
- Biology Department, Duke University, Durham, NC, 27708, USA.
| |
Collapse
|
3
|
Devens HR, Davidson PL, Byrne M, Wray GA. Hybrid Epigenomes Reveal Extensive Local Genetic Changes to Chromatin Accessibility Contribute to Divergence in Embryonic Gene Expression Between Species. Mol Biol Evol 2023; 40:msad222. [PMID: 37823438 PMCID: PMC10638671 DOI: 10.1093/molbev/msad222] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2023] [Revised: 06/14/2023] [Accepted: 07/27/2023] [Indexed: 10/13/2023] Open
Abstract
Chromatin accessibility plays an important role in shaping gene expression, yet little is known about the genetic and molecular mechanisms that influence the evolution of chromatin configuration. Both local (cis) and distant (trans) genetic influences can in principle influence chromatin accessibility and are based on distinct molecular mechanisms. We, therefore, sought to characterize the role that each of these plays in altering chromatin accessibility in 2 closely related sea urchin species. Using hybrids of Heliocidaris erythrogramma and Heliocidaris tuberculata, and adapting a statistical framework previously developed for the analysis of cis and trans influences on the transcriptome, we examined how these mechanisms shape the regulatory landscape at 3 important developmental stages, and compared our results to similar analyses of the transcriptome. We found extensive cis- and trans-based influences on evolutionary changes in chromatin, with cis effects generally larger in effect. Evolutionary changes in accessibility and gene expression are correlated, especially when expression has a local genetic basis. Maternal influences appear to have more of an effect on chromatin accessibility than on gene expression, persisting well past the maternal-to-zygotic transition. Chromatin accessibility near gene regulatory network genes appears to be distinctly regulated, with trans factors appearing to play an outsized role in the configuration of chromatin near these genes. Together, our results represent the first attempt to quantify cis and trans influences on evolutionary divergence in chromatin configuration in an outbred natural study system and suggest that chromatin regulation is more genetically complex than was previously appreciated.
Collapse
Affiliation(s)
| | | | - Maria Byrne
- School of Medical Science, The University of Sydney, Sydney, New South Wales, Australia
- School of Life and Environmental Science, The University of Sydney, Sydney, New South Wales, Australia
| | - Gregory A Wray
- Department of Biology, Duke University, Durham, NC, USA
- Center for Genomic and Computational Biology, Duke University, Durham, NC, USA
| |
Collapse
|
4
|
Massri AJ, McDonald B, Wray GA, McClay DR. Feedback circuits are numerous in embryonic gene regulatory networks and offer a stabilizing influence on evolution of those networks. EvoDevo 2023; 14:10. [PMID: 37322563 DOI: 10.1186/s13227-023-00214-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2022] [Accepted: 06/06/2023] [Indexed: 06/17/2023] Open
Abstract
The developmental gene regulatory networks (dGRNs) of two sea urchin species, Lytechinus variegatus (Lv) and Strongylocentrotus purpuratus (Sp), have remained remarkably similar despite about 50 million years since a common ancestor. Hundreds of parallel experimental perturbations of transcription factors with similar outcomes support this conclusion. A recent scRNA-seq analysis suggested that the earliest expression of several genes within the dGRNs differs between Lv and Sp. Here, we present a careful reanalysis of the dGRNs in these two species, paying close attention to timing of first expression. We find that initial expression of genes critical for cell fate specification occurs during several compressed time periods in both species. Previously unrecognized feedback circuits are inferred from the temporally corrected dGRNs. Although many of these feedbacks differ in location within the respective GRNs, the overall number is similar between species. We identify several prominent differences in timing of first expression for key developmental regulatory genes; comparison with a third species indicates that these heterochronies likely originated in an unbiased manner with respect to embryonic cell lineage and evolutionary branch. Together, these results suggest that interactions can evolve even within highly conserved dGRNs and that feedback circuits may buffer the effects of heterochronies in the expression of key regulatory genes.
Collapse
Affiliation(s)
| | - Brennan McDonald
- Department of Biology, Duke University, Box 90338, Durham, NC, 27708, USA
| | - Gregory A Wray
- Department of Biology, Duke University, Box 90338, Durham, NC, 27708, USA
| | - David R McClay
- Department of Biology, Duke University, Box 90338, Durham, NC, 27708, USA.
| |
Collapse
|
5
|
Davidson PL, Lessios HA, Wray GA, McMillan WO, Prada C. Near-Chromosomal-Level Genome Assembly of the Sea Urchin Echinometra lucunter, a Model for Speciation in the Sea. Genome Biol Evol 2023; 15:7190441. [PMID: 37280750 DOI: 10.1093/gbe/evad093] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/19/2023] [Indexed: 06/08/2023] Open
Abstract
Echinometra lucunter, the rock-boring sea urchin, is a widely distributed echinoid and a model for ecological studies of reproduction, responses to climate change, and speciation. We present a near chromosome-level genome assembly of E. lucunter, including 21 scaffolds larger than 10 Mb predicted to represent each of the chromosomes of the species. The 760.4 Mb assembly includes a scaffold N50 of 30.0 Mb and BUSCO (benchmarking universal single-copy orthologue) single copy and a duplicated score of 95.8% and 1.4%, respectively. Ab-initio gene model prediction and annotation with transcriptomic data constructed 33,989 gene models composing 50.4% of the assembly, including 37,036 transcripts. Repetitive elements make up approximately 39.6% of the assembly, and unresolved gap sequences are estimated to be 0.65%. Whole genome alignment with Echinometra sp. EZ revealed high synteny and conservation between the two species, further bolstering Echinometra as an emerging genus for comparative genomics studies. This genome assembly represents a high-quality genomic resource for future evolutionary and developmental studies of this species and more broadly of echinoderms.
Collapse
Affiliation(s)
- Phillip L Davidson
- Department of Biology, Duke University, Durham, North Carolina, USA
- Department of Biology, Indiana University, Bloomington, Indiana, USA
| | | | - Gregory A Wray
- Department of Biology, Duke University, Durham, North Carolina, USA
| | | | - Carlos Prada
- Smithsonian Tropical Research Institute, Gamboa, Panama
- Department of Biological Sciences, University of Rhode Island, Kingston, Rhode Island, USA
| |
Collapse
|
6
|
Devens HR, Davidson PL, Byrne M, Wray GA. Hybrid epigenomes reveal extensive local genetic changes to chromatin accessibility contribute to divergence in embryonic gene expression between species. bioRxiv 2023:2023.01.04.522781. [PMID: 36711588 PMCID: PMC9881966 DOI: 10.1101/2023.01.04.522781] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Chromatin accessibility plays an important role in shaping gene expression patterns across development and evolution; however, little is known about the genetic and molecular mechanisms that influence chromatin configuration itself. Because cis and trans influences can both theoretically influence the accessibility of the epigenome, we sought to better characterize the role that both of these mechanisms play in altering chromatin accessibility in two closely related sea urchin species. Using hybrids of the two species, and adapting a statistical framework previously developed for the analysis of cis and trans influences on the transcriptome, we examined how these mechanisms shape the regulatory landscape at three important developmental stages, and compared our results to similar patterns in the transcriptome. We found extensive cis- and trans-based influences on evolutionary changes in chromatin, with cis effects slightly more numerous and larger in effect. Genetic mechanisms influencing gene expression and chromatin configuration are correlated, but differ in several important ways. Maternal influences also appear to have more of an effect on chromatin accessibility than on gene expression, persisting well past the maternal-to-zygotic transition. Furthermore, chromatin accessibility near GRN genes appears to be regulated differently than the rest of the epigenome, and indicates that trans factors may play an outsized role in the configuration of chromatin near these genes. Together, our results represent the first attempt to quantify cis and trans influences on evolutionary divergence in chromatin configuration in an outbred natural study system, and suggest that the regulation of chromatin is more genetically complex than was previously appreciated.
Collapse
Affiliation(s)
| | | | - Maria Byrne
- School of Medical Science, The University of Sydney, NSW 2006, Australia
- School of Life and Environmental Science, The University of Sydney, NSW 2006, Australia
| | - Gregory A. Wray
- Department of Biology, Duke University, Durham, NC 27708, USA
- Center for Genomic and Computational Biology, Duke University, Durham, NC 27708, USA
| |
Collapse
|
7
|
Davidson PL, Guo H, Swart JS, Massri AJ, Edgar A, Wang L, Berrio A, Devens HR, Koop D, Cisternas P, Zhang H, Zhang Y, Byrne M, Fan G, Wray GA. Recent reconfiguration of an ancient developmental gene regulatory network in Heliocidaris sea urchins. Nat Ecol Evol 2022; 6:1907-1920. [PMID: 36266460 DOI: 10.1038/s41559-022-01906-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2022] [Accepted: 09/12/2022] [Indexed: 12/15/2022]
Abstract
Changes in developmental gene regulatory networks (dGRNs) underlie much of the diversity of life, but the evolutionary mechanisms that operate on regulatory interactions remain poorly understood. Closely related species with extreme phenotypic divergence provide a valuable window into the genetic and molecular basis for changes in dGRNs and their relationship to adaptive changes in organismal traits. Here we analyse genomes, epigenomes and transcriptomes during early development in two Heliocidaris sea urchin species that exhibit highly divergent life histories and in an outgroup species. Positive selection and chromatin accessibility modifications within putative regulatory elements are enriched on the branch leading to the derived life history, particularly near dGRN genes. Single-cell transcriptomes reveal a dramatic delay in cell fate specification in the derived state, which also has far fewer open chromatin regions, especially near conserved cell fate specification genes. Experimentally perturbing key transcription factors reveals profound evolutionary changes to early embryonic patterning events, disrupting regulatory interactions previously conserved for ~225 million years. These results demonstrate that natural selection can rapidly reshape developmental gene expression on a broad scale when selective regimes abruptly change. More broadly, even highly conserved dGRNs and patterning mechanisms in the early embryo remain evolvable under appropriate ecological circumstances.
Collapse
Affiliation(s)
| | - Haobing Guo
- BGI-Qingdao, BGI-Shenzhen, Qingdao, China
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen, China
| | - Jane S Swart
- Department of Biology, Duke University, Durham, NC, USA
| | | | - Allison Edgar
- Department of Biology, Duke University, Durham, NC, USA
| | - Lingyu Wang
- Department of Biology, Duke University, Durham, NC, USA
| | | | | | - Demian Koop
- School of Life and Environmental Sciences, The University of Sydney, Sydney, New South Wales, Australia
| | - Paula Cisternas
- School of Life and Environmental Sciences, The University of Sydney, Sydney, New South Wales, Australia
| | - He Zhang
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen, China
| | - Yaolei Zhang
- BGI-Qingdao, BGI-Shenzhen, Qingdao, China
- Qingdao-Europe Advanced Institute for Life Sciences, BGI-Qingdao, BGI-Shenzhen, Qingdao, China
| | - Maria Byrne
- School of Life and Environmental Sciences, The University of Sydney, Sydney, New South Wales, Australia
| | - Guangyi Fan
- BGI-Qingdao, BGI-Shenzhen, Qingdao, China
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen, China
| | - Gregory A Wray
- Department of Biology, Duke University, Durham, NC, USA.
- Center for Genomic and Computational Biology, Duke University, Durham, NC, USA.
| |
Collapse
|
8
|
Ketchum RN, Davidson PL, Smith EG, Wray GA, Burt JA, Ryan JF, Reitzel AM. A chromosome-level genome assembly of the highly heterozygous sea urchin Echinometra sp. EZ reveals adaptation in the regulatory regions of stress response genes. Genome Biol Evol 2022; 14:6717576. [PMID: 36161313 PMCID: PMC9557091 DOI: 10.1093/gbe/evac144] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/20/2022] [Indexed: 11/14/2022] Open
Abstract
Echinometra is the most widespread genus of sea urchin and has been the focus of a wide range of studies in ecology, speciation, and reproduction. However, available genetic data for this genus are generally limited to a few select loci. Here, we present a chromosome-level genome assembly based on 10x Genomics, PacBio, and Hi-C sequencing for Echinometra sp. EZ from the Persian/Arabian Gulf. The genome is assembled into 210 scaffolds totaling 817.8 Mb with an N50 of 39.5 Mb. From this assembly, we determined that the E. sp. EZ genome consists of 2n = 42 chromosomes. BUSCO analysis showed that 95.3% of BUSCO genes were complete. Ab initio and transcript-informed gene modeling and annotation identified 29,405 genes, including a conserved Hox cluster. E. sp. EZ can be found in high-temperature and high-salinity environments, and we therefore compared E. sp. EZ gene families and transcription factors associated with environmental stress response (“defensome”) with other echinoid species with similar high-quality genomic resources. While the number of defensome genes was broadly similar for all species, we identified strong signatures of positive selection in E. sp. EZ noncoding elements near genes involved in environmental response pathways as well as losses of transcription factors important for environmental response. These data provide key insights into the biology of E. sp. EZ as well as the diversification of Echinometra more widely and will serve as a useful tool for the community to explore questions in this taxonomic group and beyond.
Collapse
Affiliation(s)
- Remi N Ketchum
- Department of Biological Sciences, University of North Carolina at Charlotte, Charlotte, NC, USA.,Whitney Laboratory for Marine Bioscience, University of Florida, FL, USA
| | | | - Edward G Smith
- Department of Biological Sciences, University of North Carolina at Charlotte, Charlotte, NC, USA
| | | | - John A Burt
- Water Research Center & Center for Genomics and Systems Biology, New York University Abu Dhabi, Abu Dhabi, UAE
| | - Joseph F Ryan
- Whitney Laboratory for Marine Bioscience, University of Florida, FL, USA
| | - Adam M Reitzel
- Department of Biological Sciences, University of North Carolina at Charlotte, Charlotte, NC, USA
| |
Collapse
|
9
|
Davidson PL, Byrne M, Wray GA. Evolutionary changes in the chromatin landscape contribute to reorganization of a developmental gene network during rapid life history evolution in sea urchins. Mol Biol Evol 2022; 39:6659243. [PMID: 35946348 PMCID: PMC9435058 DOI: 10.1093/molbev/msac172] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022] Open
Abstract
Chromatin configuration is highly dynamic during embryonic development in animals, exerting an important point of control in transcriptional regulation. Yet there exists remarkably little information about the role of evolutionary changes in chromatin configuration to the evolution of gene expression and organismal traits. Genome-wide assays of chromatin configuration, coupled with whole-genome alignments, can help address this gap in knowledge in several ways. In this study we present a comparative analysis of regulatory element sequences and accessibility throughout embryogenesis in three sea urchin species with divergent life histories: a lecithotroph Heliocidaris erythrogramma, a closely related planktotroph H. tuberculata, and a distantly related planktotroph Lytechinus variegatus. We identified distinct epigenetic and mutational signatures of evolutionary modifications to the function of putative cis-regulatory elements in H. erythrogramma that have accumulated nonuniformly throughout the genome, suggesting selection, rather than drift, underlies many modifications associated with the derived life history. Specifically, regulatory elements composing the sea urchin developmental gene regulatory network are enriched for signatures of positive selection and accessibility changes which may function to alter binding affinity and access of developmental transcription factors to these sites. Furthermore, regulatory element changes often correlate with divergent expression patterns of genes involved in cell type specification, morphogenesis, and development of other derived traits, suggesting these evolutionary modifications have been consequential for phenotypic evolution in H. erythrogramma. Collectively, our results demonstrate that selective pressures imposed by changes in developmental life history rapidly reshape the cis-regulatory landscape of core developmental genes to generate novel traits and embryonic programs.
Collapse
Affiliation(s)
| | - Maria Byrne
- School of Life and Environmental Sciences, University of Sydney, Sydney, NSW, Australia
| | | |
Collapse
|
10
|
Abstract
As analyses of developmental mechanisms extend to ever more species, it becomes important to understand not just what is conserved or altered during evolution, but why. Closely related species that exhibit extreme phenotypic divergence can be uniquely informative in this regard. A case in point is the sea urchin genus Heliocidaris, which contains species that recently evolved a life history involving nonfeeding larvae following nearly half a billion years of prior evolution with feeding larvae. The resulting shift in selective regimes produced rapid and surprisingly extensive changes in developmental mechanisms that are otherwise highly conserved among echinoderm species. The magnitude and extent of these changes challenges the notion that conservation of early development in echinoderms is largely due to internal constraints that prohibit modification and instead suggests that natural selection actively maintains stability of inherently malleable trait developmental mechanisms over immense time periods. Knowing how and why natural selection changed during the evolution of nonfeeding larvae can also reveal why developmental mechanisms do and do not change in particular ways.
Collapse
Affiliation(s)
- Gregory A Wray
- Department of Biology, Duke University, Durham, NC, United States.
| |
Collapse
|
11
|
Massri AJ, Greenstreet L, Afanassiev A, Berrio A, Wray GA, Schiebinger G, McClay DR. Developmental single-cell transcriptomics in the Lytechinus variegatus sea urchin embryo. Development 2021; 148:271986. [PMID: 34463740 DOI: 10.1242/dev.198614] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2020] [Accepted: 08/20/2021] [Indexed: 12/30/2022]
Abstract
Using scRNA-seq coupled with computational approaches, we studied transcriptional changes in cell states of sea urchin embryos during development to the larval stage. Eighteen closely spaced time points were taken during the first 24 h of development of Lytechinus variegatus (Lv). Developmental trajectories were constructed using Waddington-OT, a computational approach to 'stitch' together developmental time points. Skeletogenic and primordial germ cell trajectories diverged early in cleavage. Ectodermal progenitors were distinct from other lineages by the 6th cleavage, although a small percentage of ectoderm cells briefly co-expressed endoderm markers that indicated an early ecto-endoderm cell state, likely in cells originating from the equatorial region of the egg. Endomesoderm cells also originated at the 6th cleavage and this state persisted for more than two cleavages, then diverged into distinct endoderm and mesoderm fates asynchronously, with some cells retaining an intermediate specification status until gastrulation. Seventy-nine out of 80 genes (99%) examined, and included in published developmental gene regulatory networks (dGRNs), are present in the Lv-scRNA-seq dataset and are expressed in the correct lineages in which the dGRN circuits operate.
Collapse
Affiliation(s)
- Abdull J Massri
- Department of Biology, Duke University, Durham, NC 27708, USA
| | - Laura Greenstreet
- Department of Mathematics, University of British Columbia, 121-1984 Mathematics Road, Vancouver, BC V6T 1Z2, Canada
| | - Anton Afanassiev
- Department of Mathematics, University of British Columbia, 121-1984 Mathematics Road, Vancouver, BC V6T 1Z2, Canada
| | | | - Gregory A Wray
- Department of Biology, Duke University, Durham, NC 27708, USA
| | - Geoffrey Schiebinger
- Department of Mathematics, University of British Columbia, 121-1984 Mathematics Road, Vancouver, BC V6T 1Z2, Canada
| | - David R McClay
- Department of Biology, Duke University, Durham, NC 27708, USA
| |
Collapse
|
12
|
Davidson PL, Guo H, Wang L, Berrio A, Zhang H, Chang Y, Soborowski AL, McClay DR, Fan G, Wray GA. Chromosomal-Level Genome Assembly of the Sea Urchin Lytechinus variegatus Substantially Improves Functional Genomic Analyses. Genome Biol Evol 2021; 12:1080-1086. [PMID: 32433766 DOI: 10.1093/gbe/evaa101] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/12/2020] [Indexed: 11/13/2022] Open
Abstract
Lytechinus variegatus is a camarodont sea urchin found widely throughout the western Atlantic Ocean in a variety of shallow-water marine habitats. Its distribution, abundance, and amenability to developmental perturbation make it a popular model for ecologists and developmental biologists. Here, we present a chromosomal-level genome assembly of L. variegatus generated from a combination of PacBio long reads, 10× Genomics sequencing, and HiC chromatin interaction sequencing. We show L. variegatus has 19 chromosomes with an assembly size of 870.4 Mb. The contiguity and completeness of this assembly are reflected by a scaffold length N50 of 45.5 Mb and BUSCO completeness score of 95.5%. Ab initio and transcript-informed gene modeling and annotation identified 27,232 genes with an average gene length of 12.6 kb, comprising an estimated 39.5% of the genome. Repetitive regions, on the other hand, make up 45.4% of the genome. Physical mapping of well-studied developmental genes onto each chromosome reveals nonrandom spatial distribution of distinct genes and gene families, which provides insight into how certain gene families may have evolved and are transcriptionally regulated in this species. Lastly, aligning RNA-seq and ATAC-seq data onto this assembly demonstrates the value of highly contiguous, complete genome assemblies for functional genomics analyses that is unattainable with fragmented, incomplete assemblies. This genome will be of great value to the scientific community as a resource for genome evolution, developmental, and ecological studies of this species and the Echinodermata.
Collapse
Affiliation(s)
| | - Haobing Guo
- Beijing Genomics Institute-Qingdao, China.,Beijing Genomics Institute-Shenzhen, China
| | | | | | - He Zhang
- Beijing Genomics Institute-Qingdao, China.,Beijing Genomics Institute-Shenzhen, China
| | - Yue Chang
- Beijing Genomics Institute-Qingdao, China.,Beijing Genomics Institute-Shenzhen, China
| | - Andrew L Soborowski
- Program in Computational Biology and Bioinformatics, Duke University.,Center for Genomic and Computational Biology, Duke University
| | | | - Guangyi Fan
- Beijing Genomics Institute-Qingdao, China.,Beijing Genomics Institute-Shenzhen, China
| | - Gregory A Wray
- Department of Biology, Duke University.,Program in Computational Biology and Bioinformatics, Duke University.,Center for Genomic and Computational Biology, Duke University
| |
Collapse
|
13
|
Benito-Kwiecinski S, Giandomenico SL, Sutcliffe M, Riis ES, Freire-Pritchett P, Kelava I, Wunderlich S, Martin U, Wray GA, McDole K, Lancaster MA. An early cell shape transition drives evolutionary expansion of the human forebrain. Cell 2021; 184:2084-2102.e19. [PMID: 33765444 PMCID: PMC8054913 DOI: 10.1016/j.cell.2021.02.050] [Citation(s) in RCA: 103] [Impact Index Per Article: 34.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2020] [Revised: 12/10/2020] [Accepted: 02/22/2021] [Indexed: 12/12/2022]
Abstract
The human brain has undergone rapid expansion since humans diverged from other great apes, but the mechanism of this human-specific enlargement is still unknown. Here, we use cerebral organoids derived from human, gorilla, and chimpanzee cells to study developmental mechanisms driving evolutionary brain expansion. We find that neuroepithelial differentiation is a protracted process in apes, involving a previously unrecognized transition state characterized by a change in cell shape. Furthermore, we show that human organoids are larger due to a delay in this transition, associated with differences in interkinetic nuclear migration and cell cycle length. Comparative RNA sequencing (RNA-seq) reveals differences in expression dynamics of cell morphogenesis factors, including ZEB2, a known epithelial-mesenchymal transition regulator. We show that ZEB2 promotes neuroepithelial transition, and its manipulation and downstream signaling leads to acquisition of nonhuman ape architecture in the human context and vice versa, establishing an important role for neuroepithelial cell shape in human brain expansion. Human brain organoids are expanded relative to nonhuman apes prior to neurogenesis Ape neural progenitors go through a newly identified transition morphotype state Delayed morphological transition with shorter cell cycles underlie human expansion ZEB2 is as an evolutionary regulator of this transition
Collapse
Affiliation(s)
- Silvia Benito-Kwiecinski
- MRC Laboratory of Molecular Biology, Cambridge Biomedical Campus, Francis Crick Avenue, Cambridge CB2 0QH, UK
| | - Stefano L Giandomenico
- MRC Laboratory of Molecular Biology, Cambridge Biomedical Campus, Francis Crick Avenue, Cambridge CB2 0QH, UK
| | - Magdalena Sutcliffe
- MRC Laboratory of Molecular Biology, Cambridge Biomedical Campus, Francis Crick Avenue, Cambridge CB2 0QH, UK
| | - Erlend S Riis
- Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, UK
| | - Paula Freire-Pritchett
- MRC Laboratory of Molecular Biology, Cambridge Biomedical Campus, Francis Crick Avenue, Cambridge CB2 0QH, UK
| | - Iva Kelava
- MRC Laboratory of Molecular Biology, Cambridge Biomedical Campus, Francis Crick Avenue, Cambridge CB2 0QH, UK
| | - Stephanie Wunderlich
- Leibniz Research Laboratories for Biotechnology and Artificial Organs (LEBAO), REBIRTH-Research Center for Translational and Regenerative Medicine, Hannover Medical School, 30625 Hannover, Germany; Biomedical Research in Endstage and Obstructive Lung Disease (BREATH), Member of the German Center for Lung Research (DZL), Hannover Medical School, 30625 Hannover, Germany
| | - Ulrich Martin
- Leibniz Research Laboratories for Biotechnology and Artificial Organs (LEBAO), REBIRTH-Research Center for Translational and Regenerative Medicine, Hannover Medical School, 30625 Hannover, Germany; Biomedical Research in Endstage and Obstructive Lung Disease (BREATH), Member of the German Center for Lung Research (DZL), Hannover Medical School, 30625 Hannover, Germany
| | - Gregory A Wray
- Department of Biology, Duke University, Biological Sciences Building, 124 Science Drive, Durham, NC 27708, USA
| | - Kate McDole
- MRC Laboratory of Molecular Biology, Cambridge Biomedical Campus, Francis Crick Avenue, Cambridge CB2 0QH, UK
| | - Madeline A Lancaster
- MRC Laboratory of Molecular Biology, Cambridge Biomedical Campus, Francis Crick Avenue, Cambridge CB2 0QH, UK.
| |
Collapse
|
14
|
Song H, Guo X, Sun L, Wang Q, Han F, Wang H, Wray GA, Davidson P, Wang Q, Hu Z, Zhou C, Yu Z, Yang M, Feng J, Shi P, Zhou Y, Zhang L, Zhang T. The hard clam genome reveals massive expansion and diversification of inhibitors of apoptosis in Bivalvia. BMC Biol 2021; 19:15. [PMID: 33487168 PMCID: PMC7831173 DOI: 10.1186/s12915-020-00943-9] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2020] [Accepted: 12/17/2020] [Indexed: 12/24/2022] Open
Abstract
BACKGROUND Inhibitors of apoptosis (IAPs) are critical regulators of programmed cell death that are essential for development, oncogenesis, and immune and stress responses. However, available knowledge regarding IAP is largely biased toward humans and model species, while the distribution, function, and evolutionary novelties of this gene family remain poorly understood in many taxa, including Mollusca, the second most speciose phylum of Metazoa. RESULTS Here, we present a chromosome-level genome assembly of an economically significant bivalve, the hard clam Mercenaria mercenaria, which reveals an unexpected and dramatic expansion of the IAP gene family to 159 members, the largest IAP gene repertoire observed in any metazoan. Comparative genome analysis reveals that this massive expansion is characteristic of bivalves more generally. Reconstruction of the evolutionary history of molluscan IAP genes indicates that most originated in early metazoans and greatly expanded in Bivalvia through both lineage-specific tandem duplication and retroposition, with 37.1% of hard clam IAPs located on a single chromosome. The expanded IAPs have been subjected to frequent domain shuffling, which has in turn shaped their architectural diversity. Further, we observed that extant IAPs exhibit dynamic and orchestrated expression patterns among tissues and in response to different environmental stressors. CONCLUSIONS Our results suggest that sophisticated regulation of apoptosis enabled by the massive expansion and diversification of IAPs has been crucial for the evolutionary success of hard clam and other molluscan lineages, allowing them to cope with local environmental stresses. This study broadens our understanding of IAP proteins and expression diversity and provides novel resources for studying molluscan biology and IAP function and evolution.
Collapse
Affiliation(s)
- Hao Song
- CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
- Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao, 266237, China
- Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao, 266071, China
- CAS Engineering Laboratory for Marine Ranching, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
| | - Ximing Guo
- Haskin Shellfish Research Laboratory, Department of Marine and Coastal Sciences, Rutgers University, Port Norris, NJ, USA
| | - Lina Sun
- CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
- Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao, 266237, China
- Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao, 266071, China
- CAS Engineering Laboratory for Marine Ranching, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
- University of the Chinese Academy of Sciences, Beijing, 100049, China
| | - Qianghui Wang
- Novogene Bioinformatics Institute, Beijing, 100029, China
| | - Fengming Han
- Novogene Bioinformatics Institute, Beijing, 100029, China
| | - Haiyan Wang
- CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
- Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao, 266237, China
- Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao, 266071, China
- University of the Chinese Academy of Sciences, Beijing, 100049, China
| | | | | | - Qing Wang
- University of the Chinese Academy of Sciences, Beijing, 100049, China
- Research and Development Center for Efficient Utilization of Coastal Bioresources, Yantai Institute of Coastal Zone Research, Chinese Academy of Sciences, Yantai, 264003, China
| | - Zhi Hu
- CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
- Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao, 266237, China
- University of the Chinese Academy of Sciences, Beijing, 100049, China
| | - Cong Zhou
- CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
- Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao, 266237, China
- University of the Chinese Academy of Sciences, Beijing, 100049, China
| | - Zhenglin Yu
- CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
- Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao, 266237, China
- Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao, 266071, China
- CAS Engineering Laboratory for Marine Ranching, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
| | - Meijie Yang
- CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
- Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao, 266237, China
- University of the Chinese Academy of Sciences, Beijing, 100049, China
| | - Jie Feng
- CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
- Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao, 266237, China
- Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao, 266071, China
- CAS Engineering Laboratory for Marine Ranching, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
| | - Pu Shi
- CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
- Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao, 266237, China
- University of the Chinese Academy of Sciences, Beijing, 100049, China
| | - Yi Zhou
- CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
- Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao, 266237, China
- Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao, 266071, China
- CAS Engineering Laboratory for Marine Ranching, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
- University of the Chinese Academy of Sciences, Beijing, 100049, China
| | - Libin Zhang
- CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
- Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao, 266237, China
- Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao, 266071, China
- CAS Engineering Laboratory for Marine Ranching, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
- University of the Chinese Academy of Sciences, Beijing, 100049, China
| | - Tao Zhang
- CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China.
- Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao, 266237, China.
- Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao, 266071, China.
- CAS Engineering Laboratory for Marine Ranching, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China.
- University of the Chinese Academy of Sciences, Beijing, 100049, China.
| |
Collapse
|
15
|
Markunas AM, Manivannan PKR, Ezekian JE, Agarwal A, Eisner W, Alsina K, Allen HD, Wray GA, Kim JJ, Wehrens XHT, Landstrom AP. TBX5-encoded T-box transcription factor 5 variant T223M is associated with long QT syndrome and pediatric sudden cardiac death. Am J Med Genet A 2020; 185:923-929. [PMID: 33369127 DOI: 10.1002/ajmg.a.62037] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2020] [Revised: 11/27/2020] [Accepted: 12/01/2020] [Indexed: 12/14/2022]
Abstract
Long QT syndrome (LQTS) is a genetic disease resulting in a prolonged QT interval on a resting electrocardiogram, predisposing affected individuals to polymorphic ventricular tachycardia and sudden death. Although a number of genes have been implicated in this disease, nearly one in four individuals exhibiting the LQTS phenotype are genotype-negative. Whole-exome sequencing identified a missense T223M variant in TBX5 that cosegregates with prolonged QT interval in a family with otherwise genotype-negative LQTS and sudden death. The TBX5-T223M variant was absent among large ostensibly healthy populations (gnomAD) and predicted to be pathogenic by in silico modeling based on Panther, PolyPhen-2, Provean, SIFT, SNAP2, and PredictSNP prediction tools. The variant was located in a highly conserved region of TBX5 predicted to be part of the DNA-binding interface. A luciferase assay identified a 57.5% reduction in the ability of TBX5-T223M to drive expression at the atrial natriuretic factor promotor compared to wildtype TBX5 in vitro. We conclude that the variant is pathogenic in this family, and we put TBX5 forward as a disease susceptibility allele for genotype-negative LQTS. The identification of this familial variant may serve as a basis for the identification of previously unknown mechanisms of LQTS with broader implications for cardiac electrophysiology.
Collapse
Affiliation(s)
- Alexandra M Markunas
- Department of Pediatrics, Division of Cardiology, Duke University School of Medicine, Durham, North Carolina, USA
| | - Perathu K R Manivannan
- Department of Pediatrics, Division of Cardiology, Duke University School of Medicine, Durham, North Carolina, USA
| | - Jordan E Ezekian
- Department of Pediatrics, Division of Cardiology, Duke University School of Medicine, Durham, North Carolina, USA
| | - Agnim Agarwal
- Department of Biology, Duke University, Durham, North Carolina, USA
| | - William Eisner
- Department of Pediatrics, Division of Cardiology, Duke University School of Medicine, Durham, North Carolina, USA
| | - Katherina Alsina
- Departments of Molecular Physiology & Biophysics and Medicine (Cardiology), Baylor College of Medicine, Houston, Texas, USA
| | - Hugh D Allen
- Department of Pediatrics, Section of Cardiology, Baylor College of Medicine, Houston, Texas, USA
| | - Gregory A Wray
- Department of Biology, Duke University, Durham, North Carolina, USA
| | - Jeffrey J Kim
- Department of Pediatrics, Section of Cardiology, Baylor College of Medicine, Houston, Texas, USA
| | - Xander H T Wehrens
- Departments of Molecular Physiology & Biophysics and Medicine (Cardiology), Baylor College of Medicine, Houston, Texas, USA
| | - Andrew P Landstrom
- Department of Pediatrics, Division of Cardiology, Duke University School of Medicine, Durham, North Carolina, USA.,Department of Cell Biology, Duke University School of Medicine, Durham, North Carolina, USA
| |
Collapse
|
16
|
Devens HR, Davidson PL, Deaker DJ, Smith KE, Wray GA, Byrne M. Ocean acidification induces distinct transcriptomic responses across life history stages of the sea urchin Heliocidaris erythrogramma. Mol Ecol 2020; 29:4618-4636. [PMID: 33002253 PMCID: PMC8994206 DOI: 10.1111/mec.15664] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2020] [Revised: 08/23/2020] [Accepted: 08/25/2020] [Indexed: 09/01/2023]
Abstract
Ocean acidification (OA) from seawater uptake of rising carbon dioxide emissions impairs development in marine invertebrates, particularly in calcifying species. Plasticity in gene expression is thought to mediate many of these physiological effects, but how these responses change across life history stages remains unclear. The abbreviated lecithotrophic development of the sea urchin Heliocidaris erythrogramma provides a valuable opportunity to analyse gene expression responses across a wide range of life history stages, including the benthic, post-metamorphic juvenile. We measured the transcriptional response to OA in H. erythrogramma at three stages of the life cycle (embryo, larva, and juvenile) in a controlled breeding design. The results reveal a broad range of strikingly stage-specific impacts of OA on transcription, including changes in the number and identity of affected genes; the magnitude, sign, and variance of their expression response; and the developmental trajectory of expression. The impact of OA on transcription was notably modest in relation to gene expression changes during unperturbed development and much smaller than genetic contributions from parentage. The latter result suggests that natural populations may provide an extensive genetic reservoir of resilience to OA. Taken together, these results highlight the complexity of the molecular response to OA, its substantial life history stage specificity, and the importance of contextualizing the transcriptional response to pH stress in light of normal development and standing genetic variation to better understand the capacity for marine invertebrates to adapt to OA.
Collapse
Affiliation(s)
| | | | - Dione J Deaker
- School of Life and Environmental Science, The University of Sydney, Sydney, NSW, Australia
| | - Kathryn E Smith
- The Laboratory, The Marine Biological Association, Plymouth, UK
| | - Gregory A Wray
- Department of Biology, Duke University, Durham, NC, USA
- Center for Genomic and Computational Biology, Duke University, Durham, NC, USA
| | - Maria Byrne
- School of Life and Environmental Science, The University of Sydney, Sydney, NSW, Australia
| |
Collapse
|
17
|
Berrio A, Gartner V, Wray GA. Positive selection within the genomes of SARS-CoV-2 and other Coronaviruses independent of impact on protein function. PeerJ 2020; 8:e10234. [PMID: 33088633 PMCID: PMC7571416 DOI: 10.7717/peerj.10234] [Citation(s) in RCA: 36] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2020] [Accepted: 10/04/2020] [Indexed: 12/20/2022] Open
Abstract
BACKGROUND The emergence of a novel coronavirus (SARS-CoV-2) associated with severe acute respiratory disease (COVID-19) has prompted efforts to understand the genetic basis for its unique characteristics and its jump from non-primate hosts to humans. Tests for positive selection can identify apparently nonrandom patterns of mutation accumulation within genomes, highlighting regions where molecular function may have changed during the origin of a species. Several recent studies of the SARS-CoV-2 genome have identified signals of conservation and positive selection within the gene encoding Spike protein based on the ratio of synonymous to nonsynonymous substitution. Such tests cannot, however, detect changes in the function of RNA molecules. METHODS Here we apply a test for branch-specific oversubstitution of mutations within narrow windows of the genome without reference to the genetic code. RESULTS We recapitulate the finding that the gene encoding Spike protein has been a target of both purifying and positive selection. In addition, we find other likely targets of positive selection within the genome of SARS-CoV-2, specifically within the genes encoding Nsp4 and Nsp16. Homology-directed modeling indicates no change in either Nsp4 or Nsp16 protein structure relative to the most recent common ancestor. These SARS-CoV-2-specific mutations may affect molecular processes mediated by the positive or negative RNA molecules, including transcription, translation, RNA stability, and evasion of the host innate immune system. Our results highlight the importance of considering mutations in viral genomes not only from the perspective of their impact on protein structure, but also how they may impact other molecular processes critical to the viral life cycle.
Collapse
Affiliation(s)
| | - Valerie Gartner
- Department of Biology, Duke University, Durham, NC, USA
- University Program in Genetics and Genomics, Duke University, Durham, NC, USA
| | - Gregory A. Wray
- Department of Biology, Duke University, Durham, NC, USA
- Center for Genomic and Computational Biology, Duke University, Durham, NC, USA
| |
Collapse
|
18
|
Abstract
An epithelial-mesenchymal transition (EMT) occurs in almost every metazoan embryo at the time mesoderm begins to differentiate. Several embryos have a long record as models for studying an EMT given that a known population of cells enters the EMT at a known time thereby enabling a detailed study of the process. Often, however, it is difficult to learn the molecular details of these model EMT systems because the transitioning cells are a minority of the population of cells in the embryo and in most cases there is an inability to isolate that population. Here we provide a method that enables an examination of genes expressed before, during, and after the EMT with a focus on just the cells that undergo the transition. Single cell RNA-seq (scRNA-seq) has advanced as a technology making it feasible to study the trajectory of gene expression specifically in the cells of interest, in vivo, and without the background noise of other cell populations. The sea urchin skeletogenic cells constitute only 5% of the total number of cells in the embryo yet with scRNA-seq it is possible to study the genes expressed by these cells without background noise. This approach, though not perfect, adds a new tool for uncovering the mechanism of EMT in this cell type.
Collapse
Affiliation(s)
| | | | | | - Lingyu Wang
- Department of Biology, Duke University, Durham, NC, USA
| | | | - David R McClay
- Department of Biology, Duke University, Durham, NC, USA.
| |
Collapse
|
19
|
Berrio A, Haygood R, Wray GA. Identifying branch-specific positive selection throughout the regulatory genome using an appropriate proxy neutral. BMC Genomics 2020; 21:359. [PMID: 32404186 PMCID: PMC7222330 DOI: 10.1186/s12864-020-6752-4] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2019] [Accepted: 04/21/2020] [Indexed: 01/09/2023] Open
Abstract
BACKGROUND Adaptive changes in cis-regulatory elements are an essential component of evolution by natural selection. Identifying adaptive and functional noncoding DNA elements throughout the genome is therefore crucial for understanding the relationship between phenotype and genotype. RESULTS We used ENCODE annotations to identify appropriate proxy neutral sequences and demonstrate that the conservativeness of the test can be modulated during the filtration of reference alignments. We applied the method to noncoding Human Accelerated Elements as well as open chromatin elements previously identified in 125 human tissues and cell lines to demonstrate its utility. Then, we evaluated the impact of query region length, proxy neutral sequence length, and branch count on test sensitivity and specificity. We found that the length of the query alignment can vary between 150 bp and 1 kb without affecting the estimation of selection, while for the reference alignment, we found that a length of 3 kb is adequate for proper testing. We also simulated sequence alignments under different classes of evolution and validated our ability to distinguish positive selection from relaxation of constraint and neutral evolution. Finally, we re-confirmed that a quarter of all non-coding Human Accelerated Elements are evolving by positive selection. CONCLUSION Here, we introduce a method we called adaptiPhy, which adds significant improvements to our earlier method that tests for branch-specific directional selection in noncoding sequences. The motivation for these improvements is to provide a more sensitive and better targeted characterization of directional selection and neutral evolution across the genome.
Collapse
Affiliation(s)
- Alejandro Berrio
- Department of Biology, Duke University, Biological Sciences Building, 124 Science Drive, Durham, NC, 27708, USA.
| | - Ralph Haygood
- Ronin Institute for Independent Scholarship, 127 Haddon Pl., Montclair, NJ, 07043, USA
| | - Gregory A Wray
- Department of Biology, Duke University, Biological Sciences Building, 124 Science Drive, Durham, NC, 27708, USA
| |
Collapse
|
20
|
Wang L, Israel JW, Edgar A, Raff RA, Raff EC, Byrne M, Wray GA. Genetic basis for divergence in developmental gene expression in two closely related sea urchins. Nat Ecol Evol 2020; 4:831-840. [PMID: 32284581 DOI: 10.1038/s41559-020-1165-y] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2019] [Accepted: 03/03/2020] [Indexed: 12/13/2022]
Abstract
The genetic basis for divergence in developmental gene expression among species is poorly understood, despite growing evidence that such changes underlie many interesting traits. Here we quantify transcription in hybrids of Heliocidaris tuberculata and Heliocidaris erythrogramma, two closely related sea urchins with highly divergent developmental gene expression and life histories. We find that most expression differences between species result from genetic influences that affect one stage of development, indicating limited pleiotropic consequences for most mutations that contribute to divergence in gene expression. Activation of zygotic transcription is broadly delayed in H. erythrogramma, the species with the derived life history, despite its overall faster premetamorphic development. Altered expression of several terminal differentiation genes associated with the derived larval morphology of H. erythrogramma is based largely on differences in the expression or function of their upstream regulators, providing insights into the genetic basis for the evolution of key life history traits.
Collapse
Affiliation(s)
- Lingyu Wang
- Department of Biology, Duke University, Durham, NC, USA
| | | | - Allison Edgar
- Department of Biology, Duke University, Durham, NC, USA
| | - Rudolf A Raff
- Department of Biology, Indiana University, Bloomington, IN, USA
| | | | - Maria Byrne
- School of Medical Science, The University of Sydney, Sydney, New South Wales, Australia.,School of Life and Environmental Sciences, The University of Sydney, Sydney, New South Wales, Australia
| | - Gregory A Wray
- Department of Biology, Duke University, Durham, NC, USA. .,Center for Genomic and Computational Biology, Duke University, Durham, NC, USA.
| |
Collapse
|
21
|
Swain-Lenz D, Berrio A, Safi A, Crawford GE, Wray GA. Comparative Analyses of Chromatin Landscape in White Adipose Tissue Suggest Humans May Have Less Beigeing Potential than Other Primates. Genome Biol Evol 2020; 11:1997-2008. [PMID: 31233101 PMCID: PMC6648876 DOI: 10.1093/gbe/evz134] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/19/2019] [Indexed: 12/20/2022] Open
Abstract
Humans carry a much larger percentage of body fat than other primates. Despite the central role of adipose tissue in metabolism, little is known about the evolution of white adipose tissue in primates. Phenotypic divergence is often caused by genetic divergence in cis-regulatory regions. We examined the cis-regulatory landscape of fat during human origins by performing comparative analyses of chromatin accessibility in human and chimpanzee adipose tissue using rhesus macaque as an outgroup. We find that many regions that have decreased accessibility in humans are enriched for promoter and enhancer sequences, are depleted for signatures of negative selection, are located near genes involved with lipid metabolism, and contain a short sequence motif involved in the beigeing of fat, the process in which lipid-storing white adipocytes are transdifferentiated into thermogenic beige adipocytes. The collective closing of many putative regulatory regions associated with beigeing of fat suggests a mechanism that increases body fat in humans.
Collapse
Affiliation(s)
| | | | - Alexias Safi
- Duke Center for Genomic and Computational Biology, Duke University
| | - Gregory E Crawford
- Duke Center for Genomic and Computational Biology, Duke University.,Division of Medical Genetics, Department of Pediatrics, Duke University
| | - Gregory A Wray
- Biology Department, Duke University.,Duke Center for Genomic and Computational Biology, Duke University
| |
Collapse
|
22
|
Edgar A, Byrne M, Wray GA. Embryo microinjection of the lecithotrophic sea urchin Heliocidaris erythrogramma. J Biol Methods 2019; 6:e119. [PMID: 31772951 PMCID: PMC6875645 DOI: 10.14440/jbm.2019.292] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2019] [Revised: 08/06/2019] [Accepted: 08/06/2019] [Indexed: 01/20/2023] Open
Abstract
Microinjection is a common embryological technique used for many types of experiments, including lineage tracing, manipulating gene expression, or genome editing. Injectable reagents include mRNA overexpression, mis-expression, or dominant-negative experiments to examine a gene of interest, a morpholino antisense oligo to prevent translation of an mRNA or spliceoform of interest and CRISPR-Cas9 reagents. Thus, the technique is broadly useful for basic embryological studies, constructing gene regulatory networks, and directly testing hypotheses about cis-regulatory and coding sequence changes underlying the evolution of development. However, the methods for microinjection in typical planktotrophic marine invertebrates may not work well in the highly modified eggs and embryos of lecithotrophic species. This protocol is optimized for the lecithotrophic sea urchin Heliocidaris erythrogramma.
Collapse
Affiliation(s)
- Allison Edgar
- Department of Biology, Duke University, Durham, NC 27710, USA
| | - Maria Byrne
- School of Medical Science and Bosch Institute, Department of Anatomy and Histology, The University of Sydney, Sydney, NSW 2006, Australia.,School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW 2006, Australia
| | - Gregory A Wray
- Department of Biology, Duke University, Durham, NC 27710, USA.,Center for Genomic and Computational Biology, Duke University, Durham, NC 27708, USA
| |
Collapse
|
23
|
Edsall LE, Berrio A, Majoros WH, Swain-Lenz D, Morrow S, Shibata Y, Safi A, Wray GA, Crawford GE, Allen AS. Evaluating Chromatin Accessibility Differences Across Multiple Primate Species Using a Joint Modeling Approach. Genome Biol Evol 2019; 11:3035-3053. [PMID: 31599933 PMCID: PMC6821351 DOI: 10.1093/gbe/evz218] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/07/2019] [Indexed: 12/13/2022] Open
Abstract
Changes in transcriptional regulation are thought to be a major contributor to the evolution of phenotypic traits, but the contribution of changes in chromatin accessibility to the evolution of gene expression remains almost entirely unknown. To address this important gap in knowledge, we developed a new method to identify DNase I Hypersensitive (DHS) sites with differential chromatin accessibility between species using a joint modeling approach. Our method overcomes several limitations inherent to conventional threshold-based pairwise comparisons that become increasingly apparent as the number of species analyzed rises. Our approach employs a single quantitative test which is more sensitive than existing pairwise methods. To illustrate, we applied our joint approach to DHS sites in fibroblast cells from five primates (human, chimpanzee, gorilla, orangutan, and rhesus macaque). We identified 89,744 DHS sites, of which 41% are identified as differential between species using the joint model compared with 33% using the conventional pairwise approach. The joint model provides a principled approach to distinguishing single from multiple chromatin accessibility changes among species. We found that nondifferential DHS sites are enriched for nucleotide conservation. Differential DHS sites with decreased chromatin accessibility relative to rhesus macaque occur more commonly near transcription start sites (TSS), while those with increased chromatin accessibility occur more commonly distal to TSS. Further, differential DHS sites near TSS are less cell type-specific than more distal regulatory elements. Taken together, these results point to distinct classes of DHS sites, each with distinct characteristics of selection, genomic location, and cell type specificity.
Collapse
Affiliation(s)
- Lee E Edsall
- Center for Genomic and Computational Biology, Duke University
- Division of Medical Genetics, Department of Pediatrics, Duke University
- University Program in Genetics and Genomics, Duke University
| | | | | | | | | | - Yoichiro Shibata
- Center for Genomic and Computational Biology, Duke University
- Division of Medical Genetics, Department of Pediatrics, Duke University
| | - Alexias Safi
- Center for Genomic and Computational Biology, Duke University
- Division of Medical Genetics, Department of Pediatrics, Duke University
| | - Gregory A Wray
- Center for Genomic and Computational Biology, Duke University
- Department of Biology, Duke University
| | - Gregory E Crawford
- Center for Genomic and Computational Biology, Duke University
- Division of Medical Genetics, Department of Pediatrics, Duke University
| | - Andrew S Allen
- Department of Biostatistics and Bioinformatics, Duke University
| |
Collapse
|
24
|
Davidson PL, Thompson JW, Foster MW, Moseley MA, Byrne M, Wray GA. A comparative analysis of egg provisioning using mass spectrometry during rapid life history evolution in sea urchins. Evol Dev 2019; 21:188-204. [PMID: 31102332 PMCID: PMC7232848 DOI: 10.1111/ede.12289] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2018] [Revised: 12/20/2018] [Accepted: 02/27/2019] [Indexed: 01/20/2023]
Abstract
A dramatic life history switch that has evolved numerous times in marine invertebrates is the transition from planktotrophic (feeding) to lecithotrophic (nonfeeding) larval development-an evolutionary tradeoff with many important developmental and ecological consequences. To attain a more comprehensive understanding of the molecular basis for this switch, we performed untargeted lipidomic and proteomic liquid chromatography-tandem mass spectrometry on eggs and larvae from three sea urchin species: the lecithotroph Heliocidaris erythrogramma, the closely related planktotroph Heliocidaris tuberculata, and the distantly related planktotroph Lytechinus variegatus. We identify numerous molecular-level changes possibly associated with the evolution of lecithotrophy in H. erythrogramma. We find the massive lipid stores of H. erythrogramma eggs are largely composed of low-density, diacylglycerol ether lipids that, contrary to expectations, appear to support postmetamorphic development and survivorship. Rapid premetamorphic development in this species may instead be powered by upregulated carbohydrate metabolism or triacylglycerol metabolism. We also find proteins involved in oxidative stress regulation are upregulated in H. erythrogramma eggs, and apoB-like lipid transfer proteins may be important for echinoid oogenic nutrient provisioning. These results demonstrate how mass spectrometry can enrich our understanding of life history evolution and organismal diversity by identifying specific molecules associated with distinct life history strategies and prompt new hypotheses about how and why these adaptations evolve.
Collapse
Affiliation(s)
| | - J. Will Thompson
- Center for Genomic and Computational Biology, Duke University, Durham, North Carolina
- Department of Pharmacology and Cancer Biology, Duke University, Durham, North Carolina
- Proteomics and Metabolomics Shared Resource, Duke University, Durham, North Carolina
| | - Matthew W. Foster
- Center for Genomic and Computational Biology, Duke University, Durham, North Carolina
- Department of Medicine, Duke University, Durham, North Carolina
- Proteomics and Metabolomics Shared Resource, Duke University, Durham, North Carolina
| | - M. Arthur Moseley
- Center for Genomic and Computational Biology, Duke University, Durham, North Carolina
- Department of Medicine, Duke University, Durham, North Carolina
- Proteomics and Metabolomics Shared Resource, Duke University, Durham, North Carolina
| | - Maria Byrne
- School of Medical Sciences, University of Sydney, Sydney, New South Wales, Australia
| | - Gregory A. Wray
- Department of Biology, Duke University, Durham, North Carolina
- Center for Genomic and Computational Biology, Duke University, Durham, North Carolina
| |
Collapse
|
25
|
Pizzollo J, Nielsen WJ, Shibata Y, Safi A, Crawford GE, Wray GA, Babbitt CC. Comparative Serum Challenges Show Divergent Patterns of Gene Expression and Open Chromatin in Human and Chimpanzee. Genome Biol Evol 2018; 10:826-839. [PMID: 29608722 PMCID: PMC5848805 DOI: 10.1093/gbe/evy041] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/19/2018] [Indexed: 12/13/2022] Open
Abstract
Humans experience higher rates of age-associated diseases than our closest living evolutionary relatives, chimpanzees. Environmental factors can explain many of these increases in disease risk, but species-specific genetic changes can also play a role. Alleles that confer increased disease susceptibility later in life can persist in a population in the absence of selective pressure if those changes confer positive adaptation early in life. One age-associated disease that disproportionately affects humans compared with chimpanzees is epithelial cancer. Here, we explored genetic differences between humans and chimpanzees in a well-defined experimental assay that mimics gene expression changes that happen during cancer progression: A fibroblast serum challenge. We used this assay with fibroblasts isolated from humans and chimpanzees to explore species-specific differences in gene expression and chromatin state with RNA-Seq and DNase-Seq. Our data reveal that human fibroblasts increase expression of genes associated with wound healing and cancer pathways; in contrast, chimpanzee gene expression changes are not concentrated around particular functional categories. Chromatin accessibility dramatically increases in human fibroblasts, yet decreases in chimpanzee cells during the serum response. Many regions of opening and closing chromatin are in close proximity to genes encoding transcription factors or genes involved in wound healing processes, further supporting the link between changes in activity of regulatory elements and changes in gene expression. Together, these expression and open chromatin data show that humans and chimpanzees have dramatically different responses to the same physiological stressor, and how a core physiological process can evolve quickly over relatively short evolutionary time scales.
Collapse
Affiliation(s)
- Jason Pizzollo
- Molecular and Cellular Biology Graduate Program, University of Massachusetts Amherst.,Department of Biology, University of Massachusetts Amherst
| | | | - Yoichiro Shibata
- Division of Medical Genetics, Department of Pediatrics, Duke University
| | - Alexias Safi
- Division of Medical Genetics, Department of Pediatrics, Duke University
| | - Gregory E Crawford
- Division of Medical Genetics, Department of Pediatrics, Duke University.,Center for Genomic and Computational Biology, Duke University
| | - Gregory A Wray
- Department of Biology, Duke University.,Center for Genomic and Computational Biology, Duke University.,Department of Evolutionary Anthropology, Duke University
| | | |
Collapse
|
26
|
Bryois J, Garrett ME, Song L, Safi A, Giusti-Rodriguez P, Johnson GD, Shieh AW, Buil A, Fullard JF, Roussos P, Sklar P, Akbarian S, Haroutunian V, Stockmeier CA, Wray GA, White KP, Liu C, Reddy TE, Ashley-Koch A, Sullivan PF, Crawford GE. Evaluation of chromatin accessibility in prefrontal cortex of individuals with schizophrenia. Nat Commun 2018; 9:3121. [PMID: 30087329 PMCID: PMC6081462 DOI: 10.1038/s41467-018-05379-y] [Citation(s) in RCA: 99] [Impact Index Per Article: 16.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2018] [Accepted: 06/28/2018] [Indexed: 01/19/2023] Open
Abstract
Schizophrenia genome-wide association studies have identified >150 regions of the genome associated with disease risk, yet there is little evidence that coding mutations contribute to this disorder. To explore the mechanism of non-coding regulatory elements in schizophrenia, we performed ATAC-seq on adult prefrontal cortex brain samples from 135 individuals with schizophrenia and 137 controls, and identified 118,152 ATAC-seq peaks. These accessible chromatin regions in the brain are highly enriched for schizophrenia SNP heritability. Accessible chromatin regions that overlap evolutionarily conserved regions exhibit an even higher heritability enrichment, indicating that sequence conservation can further refine functional risk variants. We identify few differences in chromatin accessibility between cases and controls, in contrast to thousands of age-related differential accessible chromatin regions. Altogether, we characterize chromatin accessibility in the human prefrontal cortex, the effect of schizophrenia and age on chromatin accessibility, and provide evidence that our dataset will allow for fine mapping of risk variants.
Collapse
Affiliation(s)
- Julien Bryois
- Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, SE-17177, Stockholm, Sweden
| | | | - Lingyun Song
- Center for Genomic and Computational Biology, Duke University, Durham, NC, 27708, USA
| | - Alexias Safi
- Center for Genomic and Computational Biology, Duke University, Durham, NC, 27708, USA
| | | | - Graham D Johnson
- Center for Genomic and Computational Biology, Duke University, Durham, NC, 27708, USA
| | - Annie W Shieh
- Department of Psychiatry, SUNY Upstate Medical University, Syracuse, NY, 13210, USA
| | - Alfonso Buil
- Research Institute of Biological Psychiatry, Mental Health Center Sct. Hans, Roskilde, 4000, Denmark
| | - John F Fullard
- Department of Psychiatry and Neuroscience, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Panos Roussos
- Department of Psychiatry and Neuroscience, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Department of Genetics and Genomic Sciences and Institute for Genomics and Multiscale Biology, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Mental Illness Research Education and Clinical Center (MIRECC), James J. Peters VA Medical Center, Bronx, NY, 10468, USA
| | - Pamela Sklar
- Department of Psychiatry and Neuroscience, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Schahram Akbarian
- Department of Psychiatry and Neuroscience, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Vahram Haroutunian
- Department of Psychiatry and Neuroscience, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- MIRECC, JJ Peters VA Medical Center, Bronx, NY, 10468, USA
| | - Craig A Stockmeier
- Department of Psychiatry and Human Behavior, Center for Psychiatric Neuroscience, University of Mississippi Medical Center, Jackson, MS, 39216, USA
| | - Gregory A Wray
- Center for Genomic and Computational Biology, Duke University, Durham, NC, 27708, USA
- Department of Biology, Duke University, Durham, NC, 27708, USA
| | - Kevin P White
- Department of Human Genetics, University of Chicago, Chicago, IL, 60637, USA
| | - Chunyu Liu
- Department of Psychiatry, SUNY Upstate Medical University, Syracuse, NY, 13210, USA
| | - Timothy E Reddy
- Center for Genomic and Computational Biology, Duke University, Durham, NC, 27708, USA
- Department of Biostatistics and Bioinformatics, Duke University, Durham, NC, 27708, USA
| | - Allison Ashley-Koch
- Duke Molecular Physiology Institute, Durham, NC, 27701, USA
- Department of Medicine, Duke University, Durham, NC, 27708, USA
| | - Patrick F Sullivan
- Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, SE-17177, Stockholm, Sweden.
- Department of Genetics, University of North Carolina, Chapel Hill, NC, 27599-7264, USA.
- Department of Psychiatry, University of North Carolina, Chapel Hill, NC, 27599-7264, USA.
| | - Gregory E Crawford
- Center for Genomic and Computational Biology, Duke University, Durham, NC, 27708, USA.
- Department of Pediatrics, Division of Medical Genetics, Duke University, Durham, NC, 27708, USA.
| |
Collapse
|
27
|
Byrne M, Koop D, Morris VB, Chui J, Wray GA, Cisternas P. Expression of genes and proteins of the pax-six-eya-dach network in the metamorphic sea urchin: Insights into development of the enigmatic echinoderm body plan and sensory structures. Dev Dyn 2017; 247:239-249. [PMID: 28850769 DOI: 10.1002/dvdy.24584] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2017] [Revised: 07/25/2017] [Accepted: 07/31/2017] [Indexed: 11/08/2022] Open
Abstract
BACKGROUND Photoreception-associated genes of the Pax-Six-Eya-Dach network (PSEDN) are deployed for many roles in addition to photoreception development. In this first study of PSEDN genes during development of the pentameral body in sea urchins, we investigated their spatial expression in Heliocidaris erythrogramma. RESULTS Expression of PSEDN genes in the hydrocoele of early (Dach, Eya, Six1/2) and/or late (Pax6, Six3/6) larvae, and the five hydrocoele lobes, the first morphological expression of pentamery, supports a role in body plan development. Pax6, Six1/2, and Six3/6 were localized to the primary and/or secondary podia and putative sensory/neuronal cells. Six1/2 and Six3/6 were expressed in the neuropil region in the terminal disc of the podia. Dach was localized to spines. Sequential up-regulation of gene expression as new podia and spines formed was evident. Rhabdomeric opsin and pax6 protein were localized to cells in the primary podia and spines. CONCLUSIONS Our results support roles for PSEDN genes in development of the pentameral body plan, contributing to our understanding of how the most unusual body plan in the Bilateria may have evolved. Development of sensory cells within the Pax-Six expression field is consistent with the role of these genes in sensory cell development in diverse species. Developmental Dynamics 247:239-249, 2018. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Maria Byrne
- School of Medical Sciences, The University of Sydney, NSW, Australia.,School of Life and Environmental Sciences, The University of Sydney, NSW, Australia
| | - Demian Koop
- School of Medical Sciences, The University of Sydney, NSW, Australia
| | - Valerie B Morris
- School of Life and Environmental Sciences, The University of Sydney, NSW, Australia
| | - Juanita Chui
- School of Medical Sciences, The University of Sydney, NSW, Australia
| | - Gregory A Wray
- Department of Biology and Center for Genomic and Computational Biology, Duke University, Durham, North Carolina
| | - Paula Cisternas
- School of Medical Sciences, The University of Sydney, NSW, Australia
| |
Collapse
|
28
|
Linchangco GV, Foltz DW, Reid R, Williams J, Nodzak C, Kerr AM, Miller AK, Hunter R, Wilson NG, Nielsen WJ, Mah CL, Rouse GW, Wray GA, Janies DA. The phylogeny of extant starfish (Asteroidea: Echinodermata) including Xyloplax, based on comparative transcriptomics. Mol Phylogenet Evol 2017; 115:161-170. [PMID: 28757447 DOI: 10.1016/j.ympev.2017.07.022] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2017] [Revised: 06/05/2017] [Accepted: 07/25/2017] [Indexed: 11/28/2022]
Abstract
Multi-locus phylogenetic studies of echinoderms based on Sanger and RNA-seq technologies and the fossil record have provided evidence for the Asterozoa-Echinozoa hypothesis. This hypothesis posits a sister relationship between asterozoan classes (Asteroidea and Ophiuroidea) and a similar relationship between echinozoan classes (Echinoidea and Holothuroidea). Despite this consensus around Asterozoa-Echinozoa, phylogenetic relationships within the class Asteroidea (sea stars or starfish) have been controversial for over a century. Open questions include relationships within asteroids and the status of the enigmatic taxon Xyloplax. Xyloplax is thought by some to represent a newly discovered sixth class of echinoderms - and by others to be an asteroid. To address these questions, we applied a novel workflow to a large RNA-seq dataset that encompassed a broad taxonomic and genomic sample. This study included 15 species sampled from all extant orders and 13 families, plus four ophiuroid species as an outgroup. To expand the taxonomic coverage, the study also incorporated five previously published transcriptomes and one previously published expressed sequence tags (EST) dataset. We developed and applied methods that used a range of alignment parameters with increasing permissiveness in terms of gap characters present within an alignment. This procedure facilitated the selection of phylogenomic data subsets from large amounts of transcriptome data. The results included 19 nested data subsets that ranged from 37 to 4,281loci. Tree searches on all data subsets reconstructed Xyloplax as a velatid asteroid rather than a new class. This result implies that asteroid morphology remains labile well beyond the establishment of the body plan of the group. In the phylogenetic tree with the highest average asteroid nodal support several monophyletic groups were recovered. In this tree, Forcipulatida and Velatida are monophyletic and form a clade that includes Brisingida as sister to Forcipulatida. Xyloplax is consistently recovered as sister to Pteraster. Paxillosida and Spinulosida are each monophyletic, with Notomyotida as sister to the Paxillosida. Valvatida is recovered as paraphyletic. The results from other data subsets are largely consistent with these results. Our results support the hypothesis that the earliest divergence event among extant asteroids separated Velatida and Forcipulatacea from Valvatacea and Spinulosida.
Collapse
Affiliation(s)
- Gregorio V Linchangco
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, NC, USA.
| | - David W Foltz
- Department of Biological Sciences, Louisiana State University, Baton Rouge, LA, USA
| | - Rob Reid
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, NC, USA
| | - John Williams
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, NC, USA
| | - Conor Nodzak
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, NC, USA
| | | | | | - Rebecca Hunter
- Department of Biology, Abilene Christian University, Abilene, TX, USA
| | - Nerida G Wilson
- Scripps Institution of Oceanography, University of California San Diego, La Jolla, CA, USA; Western Australian Museum, Locked Bag 49, Welshpool DC, Western Australia 6986, Australia
| | - William J Nielsen
- Department of Biology and Center for Genomic and Computational Biology, Duke University, Durham, NC, USA
| | - Christopher L Mah
- Department of Invertebrate Zoology, Smithsonian Institution, Washington, District of Columbia, USA
| | - Greg W Rouse
- Scripps Institution of Oceanography, University of California San Diego, La Jolla, CA, USA
| | - Gregory A Wray
- Department of Biology and Center for Genomic and Computational Biology, Duke University, Durham, NC, USA
| | - Daniel A Janies
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, NC, USA.
| |
Collapse
|
29
|
Babbitt CC, Haygood R, Nielsen WJ, Wray GA. Gene expression and adaptive noncoding changes during human evolution. BMC Genomics 2017; 18:435. [PMID: 28583075 PMCID: PMC5460488 DOI: 10.1186/s12864-017-3831-2] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2016] [Accepted: 05/31/2017] [Indexed: 01/14/2023] Open
Abstract
Background Despite evidence for adaptive changes in both gene expression and non-protein-coding, putatively regulatory regions of the genome during human evolution, the relationship between gene expression and adaptive changes in cis-regulatory regions remains unclear. Results Here we present new measurements of gene expression in five tissues of humans and chimpanzees, and use them to assess this relationship. We then compare our results with previous studies of adaptive noncoding changes, analyzing correlations at the level of gene ontology groups, in order to gain statistical power to detect correlations. Conclusions Consistent with previous studies, we find little correlation between gene expression and adaptive noncoding changes at the level of individual genes; however, we do find significant correlations at the level of biological function ontology groups. The types of function include processes regulated by specific transcription factors, responses to genetic or chemical perturbations, and differentiation of cell types within the immune system. Among functional categories co-enriched with both differential expression and noncoding adaptation, prominent themes include cancer, particularly epithelial cancers, and neural development and function. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3831-2) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Courtney C Babbitt
- Department of Biology, Duke University, Durham, NC, 27708, USA. .,Institute for Genome Sciences & Policy, Duke University, Durham, NC, 27708, USA. .,Present Address: Department of Biology, University of Massachusetts Amherst, Amherst, MA, 01003, USA.
| | | | | | - Gregory A Wray
- Department of Biology, Duke University, Durham, NC, 27708, USA.,Institute for Genome Sciences & Policy, Duke University, Durham, NC, 27708, USA.,Department of Evolutionary Anthropology, Duke University, Durham, NC, 27708, USA
| |
Collapse
|
30
|
|
31
|
Wray GA, Raff RA. RAPID EVOLUTION OF GASTRULATION MECHANISMS IN A SEA URCHIN WITH LECITHOTROPHIC LARVAE. Evolution 2017; 45:1741-1750. [DOI: 10.1111/j.1558-5646.1991.tb02684.x] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/1990] [Accepted: 04/08/1991] [Indexed: 11/29/2022]
Affiliation(s)
- Gregory A. Wray
- Department of Biology; Indiana University; Bloomington IN 47405 USA
| | - Rudolf A. Raff
- Department of Biology; Indiana University; Bloomington IN 47405 USA
| |
Collapse
|
32
|
Runcie DE, Dorey N, Garfield DA, Stumpp M, Dupont S, Wray GA. Genomic Characterization of the Evolutionary Potential of the Sea Urchin Strongylocentrotus droebachiensis Facing Ocean Acidification. Genome Biol Evol 2017; 8:3672-3684. [PMID: 28082601 PMCID: PMC5521728 DOI: 10.1093/gbe/evw272] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/10/2016] [Indexed: 12/19/2022] Open
Abstract
Ocean acidification (OA) is increasing due to anthropogenic CO2 emissions and poses a threat to marine species and communities worldwide. To better project the effects of acidification on organisms’ health and persistence, an understanding is needed of the 1) mechanisms underlying developmental and physiological tolerance and 2) potential populations have for rapid evolutionary adaptation. This is especially challenging in nonmodel species where targeted assays of metabolism and stress physiology may not be available or economical for large-scale assessments of genetic constraints. We used mRNA sequencing and a quantitative genetics breeding design to study mechanisms underlying genetic variability and tolerance to decreased seawater pH (-0.4 pH units) in larvae of the sea urchin Strongylocentrotus droebachiensis. We used a gene ontology-based approach to integrate expression profiles into indirect measures of cellular and biochemical traits underlying variation in larval performance (i.e., growth rates). Molecular responses to OA were complex, involving changes to several functions such as growth rates, cell division, metabolism, and immune activities. Surprisingly, the magnitude of pH effects on molecular traits tended to be small relative to variation attributable to segregating functional genetic variation in this species. We discuss how the application of transcriptomics and quantitative genetics approaches across diverse species can enrich our understanding of the biological impacts of climate change.
Collapse
Affiliation(s)
- Daniel E Runcie
- Department of Biology, Duke University, Durham, NC, USA.,Department of Plant Sciences, University of California, Davis, USA
| | - Narimane Dorey
- Department of Biological and Environmental Sciences, University of Gothenburg, Fiskebäckskil, Sweden
| | - David A Garfield
- Department of Biology, Duke University, Durham, NC, USA.,Integrative Research Institute for the Life Sciences, Humboldt University, Berlin, Germany
| | - Meike Stumpp
- Department of Biological and Environmental Sciences, University of Gothenburg, Fiskebäckskil, Sweden.,Helmholtz Centre for Ocean Sciences (GEOMAR), Kiel, Germany
| | - Sam Dupont
- Department of Biological and Environmental Sciences, University of Gothenburg, Fiskebäckskil, Sweden
| | - Gregory A Wray
- Department of Biology, Duke University, Durham, NC, USA.,Center for Genomic and Computational Biology, Duke University, Durham, NC, USA
| |
Collapse
|
33
|
Abstract
The timing of early animal evolution remains poorly resolved, yet remains critical for understanding nervous system evolution. Methods for estimating divergence times from sequence data have improved considerably, providing a more refined understanding of key divergences. The best molecular estimates point to the origin of metazoans and bilaterians tens to hundreds of millions of years earlier than their first appearances in the fossil record. Both the molecular and fossil records are compatible, however, with the possibility of tiny, unskeletonized, low energy budget animals during the Proterozoic that had planktonic, benthic, or meiofaunal lifestyles. Such animals would likely have had relatively simple nervous systems equipped primarily to detect food, avoid inhospitable environments and locate mates. The appearance of the first macropredators during the Cambrian would have changed the selective landscape dramatically, likely driving the evolution of complex sense organs, sophisticated sensory processing systems, and diverse effector systems involved in capturing prey and avoiding predation.
Collapse
Affiliation(s)
- Gregory A Wray
- Department of Biology and Center for Genomic and Computational Biology, Duke University, Durham, NC 27708, USA
| |
Collapse
|
34
|
Akbarian S, Liu C, Knowles JA, Vaccarino FM, Farnham PJ, Crawford GE, Jaffe AE, Pinto D, Dracheva S, Geschwind DH, Mill J, Nairn AC, Abyzov A, Pochareddy S, Prabhakar S, Weissman S, Sullivan PF, State MW, Weng Z, Peters MA, White KP, Gerstein MB, Amiri A, Armoskus C, Ashley-Koch AE, Bae T, Beckel-Mitchener A, Berman BP, Coetzee GA, Coppola G, Francoeur N, Fromer M, Gao R, Grennan K, Herstein J, Kavanagh DH, Ivanov NA, Jiang Y, Kitchen RR, Kozlenkov A, Kundakovic M, Li M, Li Z, Liu S, Mangravite LM, Mattei E, Markenscoff-Papadimitriou E, Navarro FCP, North N, Omberg L, Panchision D, Parikshak N, Poschmann J, Price AJ, Purcaro M, Reddy TE, Roussos P, Schreiner S, Scuderi S, Sebra R, Shibata M, Shieh AW, Skarica M, Sun W, Swarup V, Thomas A, Tsuji J, van Bakel H, Wang D, Wang Y, Wang K, Werling DM, Willsey AJ, Witt H, Won H, Wong CCY, Wray GA, Wu EY, Xu X, Yao L, Senthil G, Lehner T, Sklar P, Sestan N. The PsychENCODE project. Nat Neurosci 2016; 18:1707-12. [PMID: 26605881 DOI: 10.1038/nn.4156] [Citation(s) in RCA: 272] [Impact Index Per Article: 34.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
| | | | - Chunyu Liu
- University of Illinois at Chicago, Chicago, Illinois, USA
| | - James A Knowles
- University of Southern California, Los Angeles, California, USA
| | | | - Peggy J Farnham
- University of Southern California, Los Angeles, California, USA
| | | | - Andrew E Jaffe
- Lieber Institute for Brain Development, Baltimore, Maryland, USA
| | - Dalila Pinto
- Icahn School of Medicine at Mount Sinai, New York, New York, USA
| | - Stella Dracheva
- Icahn School of Medicine at Mount Sinai, New York, New York, USA
| | | | - Jonathan Mill
- King's College London, London, and University of Exeter, Devon, UK
| | | | | | | | | | | | | | | | - Zhiping Weng
- University of Massachusetts Medical School, Worcester, Massachusetts, USA
| | | | | | | | | | - Chris Armoskus
- University of Southern California, Los Angeles, California, USA
| | | | | | | | | | | | | | - Nancy Francoeur
- Icahn School of Medicine at Mount Sinai, New York, New York, USA
| | - Menachem Fromer
- Icahn School of Medicine at Mount Sinai, New York, New York, USA
| | - Robert Gao
- University of Southern California, Los Angeles, California, USA
| | - Kay Grennan
- University of Illinois at Chicago, Chicago, Illinois, USA
| | | | - David H Kavanagh
- Icahn School of Medicine at Mount Sinai, New York, New York, USA
| | - Nikolay A Ivanov
- Lieber Institute for Brain Development, Baltimore, Maryland, USA
| | - Yan Jiang
- Icahn School of Medicine at Mount Sinai, New York, New York, USA
| | | | - Alexey Kozlenkov
- Icahn School of Medicine at Mount Sinai, New York, New York, USA
| | | | - Mingfeng Li
- Yale University, New Haven, Connecticut, USA
| | - Zhen Li
- Yale University, New Haven, Connecticut, USA
| | - Shuang Liu
- Yale University, New Haven, Connecticut, USA
| | | | - Eugenio Mattei
- University of Massachusetts Medical School, Worcester, Massachusetts, USA
| | | | | | - Nicole North
- National Institute of Mental Health, Bethesda, Maryland, USA
| | | | | | | | | | - Amanda J Price
- Lieber Institute for Brain Development, Baltimore, Maryland, USA
| | - Michael Purcaro
- University of Massachusetts Medical School, Worcester, Massachusetts, USA
| | | | - Panos Roussos
- Icahn School of Medicine at Mount Sinai, New York, New York, USA
| | | | | | - Robert Sebra
- Icahn School of Medicine at Mount Sinai, New York, New York, USA
| | | | - Annie W Shieh
- University of Illinois at Chicago, Chicago, Illinois, USA
| | | | - Wenjie Sun
- Genome Institute of Singapore, Singapore
| | - Vivek Swarup
- University of California, Los Angeles, California, USA
| | - Amber Thomas
- The University of Chicago, Chicago, Illinois, USA
| | - Junko Tsuji
- University of Massachusetts Medical School, Worcester, Massachusetts, USA
| | - Harm van Bakel
- Icahn School of Medicine at Mount Sinai, New York, New York, USA
| | | | - Yongjun Wang
- University of Illinois at Chicago, Chicago, Illinois, USA
| | - Kai Wang
- University of Southern California, Los Angeles, California, USA
| | | | | | - Heather Witt
- University of Southern California, Los Angeles, California, USA
| | - Hyejung Won
- University of California, Los Angeles, California, USA
| | - Chloe C Y Wong
- King's College London, London, and University of Exeter, Devon, UK
| | | | - Emily Y Wu
- University of California, Los Angeles, California, USA
| | - Xuming Xu
- Yale University, New Haven, Connecticut, USA
| | - Lijing Yao
- University of Southern California, Los Angeles, California, USA
| | - Geetha Senthil
- National Institute of Mental Health, Bethesda, Maryland, USA
| | - Thomas Lehner
- National Institute of Mental Health, Bethesda, Maryland, USA
| | - Pamela Sklar
- Icahn School of Medicine at Mount Sinai, New York, New York, USA
| | | |
Collapse
|
35
|
Israel JW, Martik ML, Byrne M, Raff EC, Raff RA, McClay DR, Wray GA. Comparative Developmental Transcriptomics Reveals Rewiring of a Highly Conserved Gene Regulatory Network during a Major Life History Switch in the Sea Urchin Genus Heliocidaris. PLoS Biol 2016; 14:e1002391. [PMID: 26943850 PMCID: PMC4778923 DOI: 10.1371/journal.pbio.1002391] [Citation(s) in RCA: 60] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2015] [Accepted: 01/26/2016] [Indexed: 01/08/2023] Open
Abstract
The ecologically significant shift in developmental strategy from planktotrophic (feeding) to lecithotrophic (nonfeeding) development in the sea urchin genus Heliocidaris is one of the most comprehensively studied life history transitions in any animal. Although the evolution of lecithotrophy involved substantial changes to larval development and morphology, it is not known to what extent changes in gene expression underlie the developmental differences between species, nor do we understand how these changes evolved within the context of the well-defined gene regulatory network (GRN) underlying sea urchin development. To address these questions, we used RNA-seq to measure expression dynamics across development in three species: the lecithotroph Heliocidaris erythrogramma, the closely related planktotroph H. tuberculata, and an outgroup planktotroph Lytechinus variegatus. Using well-established statistical methods, we developed a novel framework for identifying, quantifying, and polarizing evolutionary changes in gene expression profiles across the transcriptome and within the GRN. We found that major changes in gene expression profiles were more numerous during the evolution of lecithotrophy than during the persistence of planktotrophy, and that genes with derived expression profiles in the lecithotroph displayed specific characteristics as a group that are consistent with the dramatically altered developmental program in this species. Compared to the transcriptome, changes in gene expression profiles within the GRN were even more pronounced in the lecithotroph. We found evidence for conservation and likely divergence of particular GRN regulatory interactions in the lecithotroph, as well as significant changes in the expression of genes with known roles in larval skeletogenesis. We further use coexpression analysis to identify genes of unknown function that may contribute to both conserved and derived developmental traits between species. Collectively, our results indicate that distinct evolutionary processes operate on gene expression during periods of life history conservation and periods of life history divergence, and that this contrast is even more pronounced within the GRN than across the transcriptome as a whole.
Collapse
Affiliation(s)
- Jennifer W. Israel
- Department of Biology, Duke University, Durham, North Carolina, United States of America
| | - Megan L. Martik
- Department of Biology, Duke University, Durham, North Carolina, United States of America
| | - Maria Byrne
- Schools of Medical and Biological Sciences, The University of Sydney, Sydney, New South Wales, Australia
| | - Elizabeth C. Raff
- Department of Biology, Indiana University, Bloomington, Indiana, United States of America
| | - Rudolf A. Raff
- Department of Biology, Indiana University, Bloomington, Indiana, United States of America
| | - David R. McClay
- Department of Biology, Duke University, Durham, North Carolina, United States of America
| | - Gregory A. Wray
- Department of Biology, Duke University, Durham, North Carolina, United States of America
- Center for Genomic and Computational Biology, Duke University, Durham, North Carolina, United States of America
| |
Collapse
|
36
|
Schierwater B, Holland PWH, Miller DJ, Stadler PF, Wiegmann BM, Wörheide G, Wray GA, DeSalle R. Never Ending Analysis of a Century Old Evolutionary Debate: “Unringing” the Urmetazoon Bell. Front Ecol Evol 2016. [DOI: 10.3389/fevo.2016.00005] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
|
37
|
Janies DA, Witter Z, Linchangco GV, Foltz DW, Miller AK, Kerr AM, Jay J, Reid RW, Wray GA. EchinoDB, an application for comparative transcriptomics of deeply-sampled clades of echinoderms. BMC Bioinformatics 2016; 17:48. [PMID: 26800861 PMCID: PMC4724074 DOI: 10.1186/s12859-016-0883-2] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2015] [Accepted: 01/07/2016] [Indexed: 02/06/2023] Open
Abstract
Background One of our goals for the echinoderm tree of life project (http://echinotol.org) is to identify orthologs suitable for phylogenetic analysis from next-generation transcriptome data. The current dataset is the largest assembled for echinoderm phylogeny and transcriptomics. We used RNA-Seq to profile adult tissues from 42 echinoderm specimens from 24 orders and 37 families. In order to achieve sampling members of clades that span key evolutionary divergence, many of our exemplars were collected from deep and polar seas. Description A small fraction of the transcriptome data we produced is being used for phylogenetic reconstruction. Thus to make a larger dataset available to researchers with a wide variety of interests, we made a web-based application, EchinoDB (http://echinodb.uncc.edu). EchinoDB is a repository of orthologous transcripts from echinoderms that is searchable via keywords and sequence similarity. Conclusions From transcripts we identified 749,397 clusters of orthologous loci. We have developed the information technology to manage and search the loci their annotations with respect to the Sea Urchin (Strongylocentrotus purpuratus) genome. Several users have already taken advantage of these data for spin-off projects in developmental biology, gene family studies, and neuroscience. We hope others will search EchinoDB to discover datasets relevant to a variety of additional questions in comparative biology. Electronic supplementary material The online version of this article (doi:10.1186/s12859-016-0883-2) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Daniel A Janies
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, 9201 University City Blvd, Charlotte, NC, 28223-0001, USA.
| | - Zach Witter
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, 9201 University City Blvd, Charlotte, NC, 28223-0001, USA
| | - Gregorio V Linchangco
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, 9201 University City Blvd, Charlotte, NC, 28223-0001, USA
| | - David W Foltz
- Department of Biological Sciences, Louisiana State University, Baton Rouge, LA, 70803, USA
| | - Allison K Miller
- Marine Laboratory, University of Guam, University Dr, Mangilao, 96923, Guam
| | - Alexander M Kerr
- Marine Laboratory, University of Guam, University Dr, Mangilao, 96923, Guam
| | - Jeremy Jay
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, 9201 University City Blvd, Charlotte, NC, 28223-0001, USA
| | - Robert W Reid
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, 9201 University City Blvd, Charlotte, NC, 28223-0001, USA
| | - Gregory A Wray
- Department of Biology, Duke University, Durham, NC, 27708, USA
| |
Collapse
|
38
|
Byrne M, Koop D, Cisternas P, Strbenac D, Yang JYH, Wray GA. Transcriptomic analysis of Nodal- and BMP-associated genes during juvenile development of the sea urchin Heliocidaris erythrogramma. Mar Genomics 2015; 24 Pt 1:41-5. [DOI: 10.1016/j.margen.2015.05.019] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2015] [Revised: 05/30/2015] [Accepted: 05/30/2015] [Indexed: 10/23/2022]
|
39
|
Bauernfeind AL, Reyzer ML, Caprioli RM, Ely JJ, Babbitt CC, Wray GA, Hof PR, Sherwood CC. High spatial resolution proteomic comparison of the brain in humans and chimpanzees. J Comp Neurol 2015; 523:2043-61. [PMID: 25779868 DOI: 10.1002/cne.23777] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2014] [Revised: 02/03/2015] [Accepted: 03/11/2015] [Indexed: 12/30/2022]
Abstract
We performed high-throughput mass spectrometry at high spatial resolution from individual regions (anterior cingulate and primary motor, somatosensory, and visual cortices) and layers of the neocortex (layers III, IV, and V) and cerebellum (granule cell layer), as well as the caudate nucleus in humans and chimpanzees. A total of 39 mass spectrometry peaks were matched with probable protein identifications in both species, allowing for comparison in expression. We explored how the pattern of protein expression varies across regions and cortical layers to provide insights into the differences in molecular phenotype of these neural structures between species. The expression of proteins differed principally in a region- and layer-specific pattern, with more subtle differences between species. Specifically, human and chimpanzee brains were similar in their distribution of proteins related to the regulation of transcription and enzyme activity but differed in their expression of proteins supporting aerobic metabolism. Whereas most work assessing molecular expression differences in the brains of primates has been performed on gene transcripts, this dataset extends current understanding of the differential molecular expression that may underlie human cognitive specializations.
Collapse
Affiliation(s)
- Amy L Bauernfeind
- Department of Anatomy and Neurobiology, Washington University School of Medicine, St. Louis, Missouri, 63110.,Department of Anthropology, Washington University in St. Louis, St. Louis, Missouri, 63130.,Department of Anthropology, The George Washington University, Washington, DC, 20052
| | - Michelle L Reyzer
- Mass Spectrometry Research Center, Vanderbilt University Medical Center, Nashville, Tennessee, 37232.,Department of Biochemistry, Vanderbilt University Medical Center, Nashville, Tennessee, 37232
| | - Richard M Caprioli
- Mass Spectrometry Research Center, Vanderbilt University Medical Center, Nashville, Tennessee, 37232.,Department of Biochemistry, Vanderbilt University Medical Center, Nashville, Tennessee, 37232
| | - John J Ely
- MAEBIOS-TM, Alamogordo, New Mexico, 88310
| | - Courtney C Babbitt
- Department of Biology, University of Massachusetts Amherst, Amherst, Massachusetts 01003
| | - Gregory A Wray
- Institute for Genome Sciences & Policy, Duke University, Durham, North Carolina, 27708.,Department of Biology, Duke University, Durham, North Carolina, 27708.,Department of Evolutionary Anthropology, Duke University, Durham, North Carolina, 27708
| | - Patrick R Hof
- Fishberg Department of Neuroscience and Friedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York, New York, 10029.,New York Consortium in Evolutionary Primatology, New York, New York
| | - Chet C Sherwood
- Department of Anthropology, The George Washington University, Washington, DC, 20052
| |
Collapse
|
40
|
Bauernfeind AL, Soderblom EJ, Turner ME, Moseley MA, Ely JJ, Hof PR, Sherwood CC, Wray GA, Babbitt CC. Evolutionary Divergence of Gene and Protein Expression in the Brains of Humans and Chimpanzees. Genome Biol Evol 2015; 7:2276-88. [PMID: 26163674 PMCID: PMC4558850 DOI: 10.1093/gbe/evv132] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Although transcriptomic profiling has become the standard approach for exploring molecular differences in the primate brain, very little is known about how the expression levels of gene transcripts relate to downstream protein abundance. Moreover, it is unknown whether the relationship changes depending on the brain region or species under investigation. We performed high-throughput transcriptomic (RNA-Seq) and proteomic (liquid chromatography coupled with tandem mass spectrometry) analyses on two regions of the human and chimpanzee brain: The anterior cingulate cortex and caudate nucleus. In both brain regions, we found a lower correlation between mRNA and protein expression levels in humans and chimpanzees than has been reported for other tissues and cell types, suggesting that the brain may engage extensive tissue-specific regulation affecting protein abundance. In both species, only a few categories of biological function exhibited strong correlations between mRNA and protein expression levels. These categories included oxidative metabolism and protein synthesis and modification, indicating that the expression levels of mRNA transcripts supporting these biological functions are more predictive of protein expression compared with other functional categories. More generally, however, the two measures of molecular expression provided strikingly divergent perspectives into differential expression between human and chimpanzee brains: mRNA comparisons revealed significant differences in neuronal communication, ion transport, and regulatory processes, whereas protein comparisons indicated differences in perception and cognition, metabolic processes, and organization of the cytoskeleton. Our results highlight the importance of examining protein expression in evolutionary analyses and call for a more thorough understanding of tissue-specific protein expression levels.
Collapse
Affiliation(s)
- Amy L Bauernfeind
- Department of Anatomy and Neurobiology, Washington University Medical School Department of Anthropology, Washington University in St. Louis Department of Anthropology and Center for the Advanced Study of Human Paleobiology, The George Washington University
| | - Erik J Soderblom
- Proteomics and Metabolomics Shared Resource, Duke University School of Medicine Center for Genomic and Computational Biology, Duke University
| | - Meredith E Turner
- Proteomics and Metabolomics Shared Resource, Duke University School of Medicine Center for Genomic and Computational Biology, Duke University
| | - M Arthur Moseley
- Proteomics and Metabolomics Shared Resource, Duke University School of Medicine Center for Genomic and Computational Biology, Duke University
| | | | - Patrick R Hof
- Fishberg Department of Neuroscience and Friedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York, New York New York Consortium in Evolutionary Primatology, New York, New York
| | - Chet C Sherwood
- Department of Anthropology and Center for the Advanced Study of Human Paleobiology, The George Washington University
| | - Gregory A Wray
- Center for Genomic and Computational Biology, Duke University Department of Biology, Duke University Department of Evolutionary Anthropology, Duke University
| | | |
Collapse
|
41
|
Boyd JL, Skove SL, Rouanet JP, Pilaz LJ, Bepler T, Gordân R, Wray GA, Silver DL. Human-chimpanzee differences in a FZD8 enhancer alter cell-cycle dynamics in the developing neocortex. Curr Biol 2015; 25:772-779. [PMID: 25702574 DOI: 10.1016/j.cub.2015.01.041] [Citation(s) in RCA: 170] [Impact Index Per Article: 18.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2014] [Revised: 08/28/2014] [Accepted: 01/16/2015] [Indexed: 12/14/2022]
Abstract
The human neocortex differs from that of other great apes in several notable regards, including altered cell cycle, prolonged corticogenesis, and increased size [1-5]. Although these evolutionary changes most likely contributed to the origin of distinctively human cognitive faculties, their genetic basis remains almost entirely unknown. Highly conserved non-coding regions showing rapid sequence changes along the human lineage are candidate loci for the development and evolution of uniquely human traits. Several studies have identified human-accelerated enhancers [6-14], but none have linked an expression difference to a specific organismal trait. Here we report the discovery of a human-accelerated regulatory enhancer (HARE5) of FZD8, a receptor of the Wnt pathway implicated in brain development and size [15, 16]. Using transgenic mice, we demonstrate dramatic differences in human and chimpanzee HARE5 activity, with human HARE5 driving early and robust expression at the onset of corticogenesis. Similar to HARE5 activity, FZD8 is expressed in neural progenitors of the developing neocortex [17-19]. Chromosome conformation capture assays reveal that HARE5 physically and specifically contacts the core Fzd8 promoter in the mouse embryonic neocortex. To assess the phenotypic consequences of HARE5 activity, we generated transgenic mice in which Fzd8 expression is under control of orthologous enhancers (Pt-HARE5::Fzd8 and Hs-HARE5::Fzd8). In comparison to Pt-HARE5::Fzd8, Hs-HARE5::Fzd8 mice showed marked acceleration of neural progenitor cell cycle and increased brain size. Changes in HARE5 function unique to humans thus alter the cell-cycle dynamics of a critical population of stem cells during corticogenesis and may underlie some distinctive anatomical features of the human brain.
Collapse
Affiliation(s)
- J Lomax Boyd
- Department of Molecular Genetics and Microbiology, Duke University Medical Center, Durham, NC 27710, USA
| | - Stephanie L Skove
- Department of Molecular Genetics and Microbiology, Duke University Medical Center, Durham, NC 27710, USA
| | - Jeremy P Rouanet
- Department of Molecular Genetics and Microbiology, Duke University Medical Center, Durham, NC 27710, USA
| | - Louis-Jan Pilaz
- Department of Molecular Genetics and Microbiology, Duke University Medical Center, Durham, NC 27710, USA
| | - Tristan Bepler
- Center for Genomic and Computational Biology, Duke University, Durham, NC 27708, USA
| | - Raluca Gordân
- Center for Genomic and Computational Biology, Duke University, Durham, NC 27708, USA; Department of Biostatistics and Bioinformatics, Duke University Medical Center, Durham, NC 27710, USA
| | - Gregory A Wray
- Center for Genomic and Computational Biology, Duke University, Durham, NC 27708, USA; Department of Biology, Duke University, Durham, NC 27708, USA; Department of Evolutionary Anthropology, Duke University, Durham, NC 27708, USA
| | - Debra L Silver
- Department of Molecular Genetics and Microbiology, Duke University Medical Center, Durham, NC 27710, USA; Department of Cell Biology, Duke University Medical Center, Durham, NC 27710, USA; Department of Neurobiology, Duke University Medical Center, Durham, NC 27710, USA; Duke Institute for Brain Sciences, Durham, NC 27710 USA.
| |
Collapse
|
42
|
Abstract
With their complex life cycle and highly derived body plan, echinoderms are unique among bilaterians. Although early development has been intensively studied, the molecular mechanisms underlying development of the adult echinoderm and its unusual radial body plan are largely unknown. To investigate the evolution of developmental changes in gene expression underlying radial body plan development and metamorphosis, we assembled a reference transcriptome de novo and used RNA-seq to measure gene expression profiles across larval, metamorphic, and postmetamorphic life cycle phases in the sea urchin Heliocidaris erythrogramma. Our results present a high-resolution view of gene expression dynamics during the complex transition from pre- to postmetamorphic development and suggest that distinct sets of regulatory and effector proteins are used during different life history phases. These analyses provide an important foundation for more detailed analyses of the evolution of the radial adult body of echinoderms.
Collapse
|
43
|
Pfefferle LW, Wray GA. Insights from a chimpanzee adipose stromal cell population: opportunities for adult stem cells to expand primate functional genomics. Genome Biol Evol 2014; 5:1995-2005. [PMID: 24092797 PMCID: PMC3814206 DOI: 10.1093/gbe/evt148] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Comparisons between humans and chimpanzees are essential for understanding traits unique to each species. However, linking important phenotypic differences to underlying molecular changes is often challenging. The ability to generate, differentiate, and profile adult stem cells provides a powerful but underutilized opportunity to investigate the molecular basis for trait differences between species within specific cell types and in a controlled environment. Here, we characterize adipose stromal cells (ASCs) from Clint, the chimpanzee whose genome was first sequenced. Using imaging and RNA-Seq, we compare the chimpanzee ASCs with three comparable human cell lines. Consistent with previous studies on ASCs in humans, the chimpanzee cells have fibroblast-like morphology and express genes encoding components of the extracellular matrix at high levels. Differentially expressed genes are enriched for distinct functional classes between species: immunity and protein processing are higher in chimpanzees, whereas cell cycle and DNA processing are higher in humans. Although hesitant to draw definitive conclusions from these data given the limited sample size, we wish to stress the opportunities that adult stem cells offer for studying primate evolution. In particular, adult stem cells provide a powerful means to investigate the profound disease susceptibilities unique to humans and a promising tool for conservation efforts with nonhuman primates. By allowing for experimental perturbations in relevant cell types, adult stem cells promise to complement classic comparative primate genomics based on in vivo sampling.
Collapse
|
44
|
Horvath JE, Ramachandran GL, Fedrigo O, Nielsen WJ, Babbitt CC, St Clair EM, Pfefferle LW, Jernvall J, Wray GA, Wall CE. Genetic comparisons yield insight into the evolution of enamel thickness during human evolution. J Hum Evol 2014; 73:75-87. [PMID: 24810709 DOI: 10.1016/j.jhevol.2014.01.005] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2013] [Revised: 10/29/2013] [Accepted: 01/09/2014] [Indexed: 12/29/2022]
Abstract
Enamel thickness varies substantially among extant hominoids and is a key trait with significance for interpreting dietary adaptation, life history trajectory, and phylogenetic relationships. There is a strong link in humans between enamel formation and mutations in the exons of the four genes that code for the enamel matrix proteins and the associated protease. The evolution of thick enamel in humans may have included changes in the regulation of these genes during tooth development. The cis-regulatory region in the 5' flank (upstream non-coding region) of MMP20, which codes for enamelysin, the predominant protease active during enamel secretion, has previously been shown to be under strong positive selection in the lineages leading to both humans and chimpanzees. Here we examine evidence for positive selection in the 5' flank and 3' flank of AMELX, AMBN, ENAM, and MMP20. We contrast the human sequence changes with other hominoids (chimpanzees, gorillas, orangutans, gibbons) and rhesus macaques (outgroup), a sample comprising a range of enamel thickness. We find no evidence for positive selection in the protein-coding regions of any of these genes. In contrast, we find strong evidence for positive selection in the 5' flank region of MMP20 and ENAM along the lineage leading to humans, and in both the 5' flank and 3' flank regions of MMP20 along the lineage leading to chimpanzees. We also identify putative transcription factor binding sites overlapping some of the species-specific nucleotide sites and we refine which sections of the up- and downstream putative regulatory regions are most likely to harbor important changes. These non-coding changes and their potential for differential regulation by transcription factors known to regulate tooth development may offer insight into the mechanisms that allow for rapid evolutionary changes in enamel thickness across closely-related species, and contribute to our understanding of the enamel phenotype in hominoids.
Collapse
Affiliation(s)
- Julie E Horvath
- North Carolina Museum of Natural Sciences, Nature Research Center, Raleigh, NC 27601, USA; Department of Biology, North Carolina Central University, Durham, NC 27707, USA; Department of Evolutionary Anthropology, Duke University, Durham, NC 27708, USA; Duke Institute for Genome Sciences and Policy, Duke University, Durham, NC 27708, USA
| | | | - Olivier Fedrigo
- Duke Institute for Genome Sciences and Policy, Duke University, Durham, NC 27708, USA
| | | | - Courtney C Babbitt
- Duke Institute for Genome Sciences and Policy, Duke University, Durham, NC 27708, USA; Department of Biology, Duke University, Durham, NC 27708, USA
| | | | | | - Jukka Jernvall
- Institute for Biotechnology, University of Helsinki, Helsinki, Finland
| | - Gregory A Wray
- Department of Evolutionary Anthropology, Duke University, Durham, NC 27708, USA; Duke Institute for Genome Sciences and Policy, Duke University, Durham, NC 27708, USA; Department of Biology, Duke University, Durham, NC 27708, USA
| | - Christine E Wall
- Department of Evolutionary Anthropology, Duke University, Durham, NC 27708, USA.
| |
Collapse
|
45
|
Muntané G, Horvath JE, Hof PR, Ely JJ, Hopkins WD, Raghanti MA, Lewandowski AH, Wray GA, Sherwood CC. Analysis of synaptic gene expression in the neocortex of primates reveals evolutionary changes in glutamatergic neurotransmission. ACTA ACUST UNITED AC 2014; 25:1596-607. [PMID: 24408959 DOI: 10.1093/cercor/bht354] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Increased relative brain size characterizes the evolution of primates, suggesting that enhanced cognition plays an important part in the behavioral adaptations of this mammalian order. In addition to changes in brain anatomy, cognition can also be regulated by molecular changes that alter synaptic function, but little is known about modifications of synapses in primate brain evolution. The aim of the current study was to investigate the expression patterns and evolution of 20 synaptic genes from the prefrontal cortex of 12 primate species. The genes investigated included glutamate receptors, scaffolding proteins, synaptic vesicle components, as well as factors involved in synaptic vesicle release and structural components of the nervous system. Our analyses revealed that there have been significant changes during primate brain evolution in the components of the glutamatergic signaling pathway in terms of gene expression, protein expression, and promoter sequence changes. These results could entail functional modifications in the regulation of specific genes related to processes underlying learning and memory.
Collapse
Affiliation(s)
- Gerard Muntané
- Department of Anthropology, The George Washington University, Washington, DC 20052, USA
| | - Julie E Horvath
- Department of Evolutionary Anthropology, Duke University, Durham, NC 27707, USA North Carolina Museum of Natural Sciences, Nature Research Center, Raleigh, NC 27601, USA Department of Biology, NC Central University, Durham, NC 27707, USA
| | - Patrick R Hof
- Fishberg Department of Neuroscience and Friedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York, NY 10029 New York Consortium for Primate Evolution, New York, NY 10024, USA
| | - John J Ely
- Alamogordo Primate Facility, Holloman Air Force Base, NM 88330, USA
| | - William D Hopkins
- Neuroscience Institute and Language Research Center, Georgia State University, Atlanta, GA 30302, USA
| | - Mary Ann Raghanti
- Department of Anthropology and School of Biomedical Sciences, Kent State University, Kent, OH 44242, USA
| | | | - Gregory A Wray
- Department of Evolutionary Anthropology, Duke University, Durham, NC 27707, USA Biology Department, Duke University, Durham, NC, USA
| | - Chet C Sherwood
- Department of Anthropology, The George Washington University, Washington, DC 20052, USA
| |
Collapse
|
46
|
Abstract
Evolutionary genetics has entered an unprecedented era of discovery, catalyzed in large part by the development of technologies that provide information about genome sequence and function. An important benefit is the ability to move beyond a handful of model organisms in lab settings to identify the genetic basis for evolutionarily interesting traits in many organisms in natural settings. Other benefits are the abilities to identify causal mutations and validate their phenotypic consequences more readily and in many more species. Genomic technologies have reinvigorated interest in some of the most fundamental and persistent questions in evolutionary genetics, revealed previously unsuspected evolutionary phenomena, and opened the door to a wide range of new questions.
Collapse
Affiliation(s)
- Gregory A. Wray
- Department of Biology and Institute for Genome Sciences & Policy, Duke University, Durham, North Carolina 27701
| |
Collapse
|
47
|
Garfield DA, Runcie DE, Babbitt CC, Haygood R, Nielsen WJ, Wray GA. The impact of gene expression variation on the robustness and evolvability of a developmental gene regulatory network. PLoS Biol 2013; 11:e1001696. [PMID: 24204211 PMCID: PMC3812118 DOI: 10.1371/journal.pbio.1001696] [Citation(s) in RCA: 64] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2013] [Accepted: 09/16/2013] [Indexed: 11/18/2022] Open
Abstract
Regulatory interactions buffer development against genetic and environmental perturbations, but adaptation requires phenotypes to change. We investigated the relationship between robustness and evolvability within the gene regulatory network underlying development of the larval skeleton in the sea urchin Strongylocentrotus purpuratus. We find extensive variation in gene expression in this network throughout development in a natural population, some of which has a heritable genetic basis. Switch-like regulatory interactions predominate during early development, buffer expression variation, and may promote the accumulation of cryptic genetic variation affecting early stages. Regulatory interactions during later development are typically more sensitive (linear), allowing variation in expression to affect downstream target genes. Variation in skeletal morphology is associated primarily with expression variation of a few, primarily structural, genes at terminal positions within the network. These results indicate that the position and properties of gene interactions within a network can have important evolutionary consequences independent of their immediate regulatory role.
Collapse
Affiliation(s)
- David A. Garfield
- Department of Biology, Duke University, Durham, North Carolina, United States of America
| | - Daniel E. Runcie
- Department of Biology, Duke University, Durham, North Carolina, United States of America
| | - Courtney C. Babbitt
- Department of Biology, Duke University, Durham, North Carolina, United States of America
- Institute for Genome Sciences & Policy, Duke University, Durham, North Carolina, United States of America
| | - Ralph Haygood
- Department of Biology, Duke University, Durham, North Carolina, United States of America
- Center for Systems Biology, Duke University, Durham, North Carolina, United States of America
| | - William J. Nielsen
- Department of Biology, Duke University, Durham, North Carolina, United States of America
- Institute for Genome Sciences & Policy, Duke University, Durham, North Carolina, United States of America
| | - Gregory A. Wray
- Department of Biology, Duke University, Durham, North Carolina, United States of America
- Institute for Genome Sciences & Policy, Duke University, Durham, North Carolina, United States of America
- Center for Systems Biology, Duke University, Durham, North Carolina, United States of America
- * E-mail:
| |
Collapse
|
48
|
Runcie DE, Wiedmann RT, Archie EA, Altmann J, Wray GA, Alberts SC, Tung J. Social environment influences the relationship between genotype and gene expression in wild baboons. Philos Trans R Soc Lond B Biol Sci 2013; 368:20120345. [PMID: 23569293 DOI: 10.1098/rstb.2012.0345] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Variation in the social environment can have profound effects on survival and reproduction in wild social mammals. However, we know little about the degree to which these effects are influenced by genetic differences among individuals, and conversely, the degree to which social environmental variation mediates genetic reaction norms. To better understand these relationships, we investigated the potential for dominance rank, social connectedness and group size to modify the effects of genetic variation on gene expression in the wild baboons of the Amboseli basin. We found evidence for a number of gene-environment interactions (GEIs) associated with variation in the social environment, encompassing social environments experienced in adulthood as well as persistent effects of early life social environment. Social connectedness, maternal dominance rank and group size all interacted with genotype to influence gene expression in at least one sex, and either in early life or in adulthood. These results suggest that social and behavioural variation, akin to other factors such as age and sex, can impact the genotype-phenotype relationship. We conclude that GEIs mediated by the social environment are important in the evolution and maintenance of individual differences in wild social mammals, including individual differences in responses to social stressors.
Collapse
Affiliation(s)
- Daniel E Runcie
- Department of Biology, Duke University, Durham, NC 27708, USA
| | | | | | | | | | | | | |
Collapse
|
49
|
Garfield D, Haygood R, Nielsen WJ, Wray GA. Population genetics of cis-regulatory sequences that operate during embryonic development in the sea urchin Strongylocentrotus purpuratus. Evol Dev 2013; 14:152-67. [PMID: 23017024 DOI: 10.1111/j.1525-142x.2012.00532.x] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Despite the fact that noncoding sequences comprise a substantial fraction of functional sites within all genomes, the evolutionary mechanisms that operate on genetic variation within regulatory elements remain poorly understood. In this study, we examine the population genetics of the core, upstream cis-regulatory regions of eight genes (AN, CyIIa, CyIIIa, Endo16, FoxB, HE, SM30 a, and SM50) that function during the early development of the purple sea urchin, Strongylocentrotus purpuratus. Quantitative and qualitative measures of segregating variation are not conspicuously different between cis-regulatory and closely linked "proxy neutral" noncoding regions containing no known functional sites. Length and compound mutations are common in noncoding sequences; conventional descriptive statistics ignore such mutations, under-representing true genetic variation by approximately 28% for these loci in this population. Patterns of variation in the cis-regulatory regions of six of the genes examined (CyIIa, CyIIIa, Endo16, FoxB, AN, and HE) are consistent with directional selection. Genetic variation within annotated transcription factor binding sites is comparable to, and frequently greater than, that of surrounding sequences. Comparisons of two paralog pairs (CyIIa/CyIIIa and AN/HE) suggest that distinct evolutionary processes have operated on their cis-regulatory regions following gene duplication. Together, these analyses provide a detailed view of the evolutionary mechanisms operating on noncoding sequences within a natural population, and underscore how little is known about how these processes operate on cis-regulatory sequences.
Collapse
Affiliation(s)
- David Garfield
- Department of Biology and Institute for Genome Sciences & Policy, Duke University, Box 90338, Durham, NC 27708, USA
| | | | | | | |
Collapse
|
50
|
Wray GA. Adaptation and Gene Networks. Bioscience 2012. [DOI: 10.1525/bio.2012.62.12.10] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
|