Click
here to close Hello! We notice that
you are using Internet Explorer, which is not supported by Echinobase
and may cause the site to display incorrectly. We suggest using a
current version of Chrome,
FireFox,
or Safari.
Sci Rep
2019 Nov 07;91:16201. doi: 10.1038/s41598-019-52577-9.
Show Gene links
Show Anatomy links
Developmental transcriptomes of the sea star, Patiria miniata, illuminate how gene expression changes with evolutionary distance.
Gildor T
,
Lalzar M
,
Hinman VF
,
Ben-Tabou de-Leon S
.
Abstract
Understanding how changes in developmental gene expression alter morphogenesis is a fundamental problem in development and evolution. A promising approach to address this problem is to compare the developmental transcriptomes between related species. The echinoderm phylum consists of several model species that have significantly contributed to the understanding of gene regulation and evolution. Particularly, the regulatory networks of the sea star, Patiria miniata (P. miniata), have been extensively studied, however developmental transcriptomes for this species were lacking. Here we generated developmental transcriptomes of P. miniata and compared these with those of two sea urchins species. We demonstrate that the conservation of gene expression depends on gene function, cell type and evolutionary distance. With increasing evolutionary distance the interspecies correlations in gene expression decreases. The reduction is more severe in the correlations between morphologically equivalent stages (diagonal elements) than in the correlation between morphologically distinct stages (off-diagonal elements). This could reflect a decrease in the morphological constraints compared to other constraints that shape gene expression at large evolutionary divergence. Within this trend, the interspecies correlations of developmental control genes maintain their diagonality at large evolutionary distance, and peak at the onset of gastrulation, supporting the hourglass model of phylotypic stage conservation.
MCB 1715721 National Science Foundation (NSF), 2304/15 Israel Science Foundation (ISF), 41/14 Israel Science Foundation (ISF), 2015031 United States - Israel Binational Science Foundation (BSF), P41 HD071837 NICHD NIH HHS
Figure 1. Developmental time points studied and examples for gene expression profiles on the three species. (A) Images of adult and larval stage of P. miniata, P. lividus and S. purpuratus. Arrows point to the sea urchin skeletogenic rods and arrow heads point to the sea urchin pigments. (B) Images of P. miniata, P. lividus and S. purpuratus embryos at the developmental stages that were studied in this work. Time point 6 hpf in S. purpuratus does not have RNA-seq data. (C) Relative gene expression in the three species measured in the current paper by RNA-seq for P. miniata (orange curves), by RNA-seq for P. lividus19 (purple curves) and by nanostring for S. purpuratus28 (black curves). Error bars indicate standard deviation. To obtain relative expression levels for each species we divide the level at each time point in the maximal mRNA level measured for this species in this time interval; so 1 is the maximal expression in this time interval.
Figure 2. Venn diagram and NMDS analysis of 1:1:1 homologous genes. (A) Venn diagram showing the number of 1:1:1 homologous genes expressed in all three species, in two of the species or only in one species. (B) First two principal components of expression variation (NMDS) between different developmental time points in P. miniata (orange), P. lividus (purple) and S. purpuratus (black).
Figure 3. Interspecies Pearson correlations for different GO terms, ordered by the level of matrix diagonality (MD) of Pl-Sp matrices. In each panel, from (AâJ) we present the Pearson correlation of the expression levels of genes with specific GO term between different developmental stages in two species. Upper matrix in each panel shows the Pearson correlation between the two sea urchins (P. lividus and S. purpuratus) and the bottom matrix is the Pearson correlation between the sea star, P. miniata and the sea urchin P. lividus. These matrices include the seven developmental points that have RNA-seq data in all species (Fig. 1B, excluding 9 hpf in P. miniata and 4 hpf in P. lividus). In each panel we indicate the GO term tested, the number of genes in each set, the average correlation strength in the diagonal (AC) and the matrix diagonality (MD), see text for explanation. Linear color scale of Pearson correlations 0â1 is identical for all graphs and given at the middle of the figure. (F) Shows the Pearson interspecies correlation for all 1:1:1 genes.
Figure 4. The average correlation strength and the matrix diagonality are independent parameters that reflect different properties of expression conservation. (A) The average correlation strength (AC) between Pl-Sp (black bars) and Pl-Pm (red bars) in receding order of Pl-Pm correlation strength. Error bars indicate standard deviation of the correlation strength along the matrix diagonal. (B) Matrix diagonality (MD) of the interspecies correlations between Pl-Sp (black bars) and Pl-Pm (cyan bars) in receding order of Pl-Pm matrix diagonality. In (A,B) the different gene sets are colored by the following key: developmental GO terms in green, housekeeping GO terms in blue, lineage specific genes in red and all the other sets in black. (C) The matrix diagonality changes independently of the average correlation for both Pl-Sp (black dots) and Pl-Pm (orange dots). (D) The interspecies average correlation between S. purpuratus and P. lividus corresponds to the interspecies average correlation between P. lividus and P. miniata (R pearsonâ=â0.68). (E) The matrix diagonality of the interspecies correlations between S. purpuratus and P. lividus relates to the matrix diagonality between P. lividus and P. miniata (R Pearsonâ=â0.78). In (D,E) we didnât include the skeletogenic lineage data point as it is an outlier.
Figure 5. The correlation matrix diagonality, MD, reflects how the dominance between cellular and developmental constraints changes with evolutionary distance. Illustration of typical interspecies correlation matrices of developmental control genes and housekeeping genes between closely related and further diverged species (upper and lower panels, respectively). With increasing evolutionary distance, that is, between the sea urchin and the sea star, the average correlation and the diagonality decrease for all gene sets but the diagonality of developmental control genes is least affected and they still maintain the hourglass pattern (left panels). On the other hand, the interspecies correlation matrices of housekeeping genes are strong and non-diagonal even between the two sea urchins and remain non-diagonal between the sea urchin and the sea star (right panels).
Anders,
HTSeq--a Python framework to work with high-throughput sequencing data.
2015, Pubmed
Anders,
HTSeq--a Python framework to work with high-throughput sequencing data.
2015,
Pubmed
Aubry,
Deep evolutionary comparison of gene expression identifies parallel recruitment of trans-factors in two independent origins of C4 photosynthesis.
2014,
Pubmed
Barsi,
Genome-wide assessment of differential effector gene use in embryogenesis.
2015,
Pubmed
,
Echinobase
Ben-Tabou de-Leon,
Gene regulatory control in the sea urchin aboral ectoderm: spatial initiation, signaling inputs, and cell fate lockdown.
2013,
Pubmed
,
Echinobase
Bogdanović,
Active DNA demethylation at enhancers during the vertebrate phylotypic period.
2016,
Pubmed
Bolger,
Trimmomatic: a flexible trimmer for Illumina sequence data.
2014,
Pubmed
Bray,
Near-optimal probabilistic RNA-seq quantification.
2016,
Pubmed
Cameron,
SpBase: the sea urchin genome database and web site.
2009,
Pubmed
,
Echinobase
Cary,
Echinoderm development and evolution in the post-genomic era.
2017,
Pubmed
,
Echinobase
Dobin,
STAR: ultrafast universal RNA-seq aligner.
2013,
Pubmed
Dylus,
Developmental transcriptomics of the brittle star Amphiura filiformis reveals gene regulatory network rewiring in echinoderm larval skeleton evolution.
2018,
Pubmed
,
Echinobase
Flores,
The skeletal proteome of the sea star Patiria miniata and evolution of biomineralization in echinoderms.
2017,
Pubmed
,
Echinobase
Gildor,
Comparative Study of Regulatory Circuits in Two Sea Urchin Species Reveals Tight Control of Timing and High Conservation of Expression Dynamics.
2015,
Pubmed
,
Echinobase
Gildor,
Regulatory heterochronies and loose temporal scaling between sea star and sea urchin regulatory circuits.
2017,
Pubmed
,
Echinobase
Gildor,
Comparative Studies of Gene Expression Kinetics: Methodologies and Insights on Development and Evolution.
2018,
Pubmed
Grabherr,
Full-length transcriptome assembly from RNA-Seq data without a reference genome.
2011,
Pubmed
Haas,
De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis.
2013,
Pubmed
Hinman,
Developmental gene regulatory network architecture across 500 million years of echinoderm evolution.
2003,
Pubmed
,
Echinobase
Irie,
Comparative transcriptome analysis reveals vertebrate phylotypic period during organogenesis.
2011,
Pubmed
Irie,
The developmental hourglass model: a predictor of the basic body plan?
2014,
Pubmed
Israel,
Comparative Developmental Transcriptomics Reveals Rewiring of a Highly Conserved Gene Regulatory Network during a Major Life History Switch in the Sea Urchin Genus Heliocidaris.
2016,
Pubmed
,
Echinobase
Kalinka,
Gene expression divergence recapitulates the developmental hourglass model.
2010,
Pubmed
Levin,
The mid-developmental transition and the evolution of animal body plans.
2016,
Pubmed
Li,
Direct and indirect control of oral ectoderm regulatory gene expression by Nodal signaling in the sea urchin embryo.
2012,
Pubmed
,
Echinobase
Love,
Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2.
2014,
Pubmed
Malik,
Parallel embryonic transcriptional programs evolve under distinct constraints and may enable morphological conservation amidst adaptation.
2017,
Pubmed
,
Echinobase
Materna,
High accuracy, high-resolution prevalence measurement for the majority of locally expressed regulatory genes in early sea urchin development.
2010,
Pubmed
,
Echinobase
McCauley,
A conserved gene regulatory network subcircuit drives different developmental fates in the vegetal pole of highly divergent echinoderm embryos.
2010,
Pubmed
,
Echinobase
McCauley,
Development of an embryonic skeletogenic mesenchyme lineage in a sea cucumber reveals the trajectory of change for the evolution of novel structures in echinoderms.
2012,
Pubmed
,
Echinobase
McCauley,
Expression of wnt and frizzled genes during early sea star development.
2013,
Pubmed
,
Echinobase
Morgulis,
Possible cooption of a VEGF-driven tubulogenesis program for biomineralization in echinoderms.
2019,
Pubmed
,
Echinobase
Oliveri,
Global regulatory logic for specification of an embryonic cell lineage.
2008,
Pubmed
,
Echinobase
Piasecka,
The hourglass and the early conservation models--co-existing patterns of developmental constraints in vertebrates.
2013,
Pubmed
Raff,
The active evolutionary lives of echinoderm larvae.
2006,
Pubmed
,
Echinobase
Robinson,
edgeR: a Bioconductor package for differential expression analysis of digital gene expression data.
2010,
Pubmed
She,
Definition, conservation and epigenetics of housekeeping and tissue-enriched genes.
2009,
Pubmed
Solek,
An ancient role for Gata-1/2/3 and Scl transcription factor homologs in the development of immunocytes.
2013,
Pubmed
,
Echinobase
Soneson,
Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences.
2015,
Pubmed
Tu,
Quantitative developmental transcriptomes of the sea urchin Strongylocentrotus purpuratus.
2014,
Pubmed
,
Echinobase
Tu,
Gene structure in the sea urchin Strongylocentrotus purpuratus based on transcriptome analysis.
2012,
Pubmed
,
Echinobase
Waterhouse,
BUSCO Applications from Quality Assessments to Gene Prediction and Phylogenomics.
2018,
Pubmed
Yaguchi,
A Wnt-FoxQ2-nodal pathway links primary and secondary axis specification in sea urchin embryos.
2008,
Pubmed
,
Echinobase
Yanai,
Mapping gene expression in two Xenopus species: evolutionary constraints and developmental flexibility.
2011,
Pubmed
Yankura,
Gene regulatory network for neurogenesis in a sea star embryo connects broad neural specification and localized patterning.
2013,
Pubmed
,
Echinobase
Yankura,
Uncoupling of complex regulatory patterning during evolution of larval development in echinoderms.
2010,
Pubmed
,
Echinobase
Young,
Gene ontology analysis for RNA-seq: accounting for selection bias.
2010,
Pubmed
Zabidi,
Enhancer-core-promoter specificity separates developmental and housekeeping gene regulation.
2015,
Pubmed