See Also Figure 1
Total Page:16
File Type:pdf, Size:1020Kb
Figure S1. Box-and-whisker plots depicting the range of expression values per developmental stage, with DESeq normalization (A) or quantile normalization (B). See also Figure 1. Figure S2. Lv-Setmar expression has low variation over developmental time. A. A plot of Lv-setmar versus Lv-ubiquitin expression over time demonstrates that Lv-setmar exhibits less temporal variation than Lv-ubiquitin. B. A representative gel showing Lv-setmar qPCR products amplified from cDNAs representing each sequenced stage in this study, demonstrating comparable product levels and an absence of spurious amplification products. See also Figure 1E. Figure S3. LvEDGE database. Screen shots showing the home page (A), the search window (B), an example search with a temporal expression plot (C), and the numerical data reflected in the plot (D) for the LvEDGE public database, which hosts the data described herein. stage 1 2 3 4 5 6 7 8 9 10 11 Category Subcategory 2-cell 60-cell EB HB TVP MB EG MG LG EP LP meiotic Cell Division Cytokinesis Mitosis checkpoint cell division recombination cell cycle stem cell left-right cell left-right Development maintenance asymmetry morphogenesis asymmetry regulation of multicellular organismal process cell soma cell soma Gene Expression chromatin SWI/SNF Control Chromatin modification chromatin binding complex methylated histone Binding negative sequence- sequence- sequence- regulation of sequence- specific DNA specific DNA specific DNA transcription specific DNA sequence-specific DNA binding binding binding binding factor activity binding DNA binding nuclear mRNA nuclear mRNA nuclear mRNA nuclear mRNA nuclear mRNA nuclear mRNA splicing, via splicing, via splicing, via splicing, via splicing, via splicing, via Processing spliceosome spliceosome spliceosome spliceosome spliceosome spliceosome mRNA processing translation translation translation translation initiation factor initiation factor initiation factor initiation factor Translation activity activity activity activity NADH fatty acid dehydrogenase primary glycosaminoglycan ATP synthesis metabolic (ubiquinone) cell redox cell redox metabolic mitochondrial mitochondrial biosynthetic coupled proton Metabolic process activity homeostasis homeostasis process inner membrane matrix process transport oxidoreductase activity, acting on paired donors, with incorporation or reduction of positive regulation molecular monooxygenase monooxygenase oxidoreductase of cellular oxidoreductase oxygen activity activity activity metabolic process activity oxidoreductase activity, acting on the aldehyde transferase or oxo group of sulfuric ester activity, donors, NAD or hydrolase ATP catabolic transferring NADP as activity process hexosyl groups acceptor ATP synthesis coupled proton transport DNA Metabolic Process ciliary or microtubule Regulation cytoskeleton flagellar motility motor activity actin binding metallo- metallo- serine-type serine-type threonine threonine metallo- metallo- endopeptidase endopeptidase endopeptidase endopeptidase endopeptidase endopeptidase endopeptidase endopeptidase serine-type protease activity activity activity activity activity activity activity activity peptidase activity peptidase metallopeptidase metallopeptidase activity activity activity platelet platelet activation activation identical protein protein protein other binding ubiquitination polyubiquitination positive regulation regulation of of cellular cellular process metabolic process G-protein G-protein coupled coupled G-protein G-protein G-protein G-protein G-protein G-protein G-protein receptor receptor coupled receptor coupled coupled receptor coupled receptor coupled receptor coupled G-protein coupled G-protein coupled coupled Signaling activity activity activity receptor activity activity activity activity receptor activity receptor activity receptor activity receptor activity phosphoinositide cGMP phospholipase C neuropeptide Y neuropeptide Y neuropeptide Y biosynthetic regulation of activity receptor activity receptor activity receptor activity process apoptosis non-membrane spanning protein tyrosine kinase scavenger scavenger MAPKKK activity receptor activity receptor activity receptor binding cascade small GTPase small GTPase small GTPase growth factor mediated signal mediated signal mediated signal activity transduction transduction transduction positive regulation of GTPase activity axonemal dynein microtubule- Structural cytoskeleton complex based process dendritic shaft dendritic spine microtubule cytoskeletal part sarcomere organization nucleosome organelle assembly nucleosome nucleosome nucleosome nucleosome nuclear envelope integral to Golgi intracellular intracellular membrane organelle part organelle structural structural structural structural structural structural structural structural constituent of constituent of constituent of constituent of constituent of constituent of constituent of constituent of ribosome ribosome ribosome ribosome ribosome ribosome ribosome ribosome ribosome ribosome Proteinaceous protein extracellular Extracellular other cytoplasmic part polymerization matrix intracellular part intracellular part region part Extracellular region Figure S4. GO term enrichment is shown at each sequenced developmental stage. 96 GO terms identified as enriched in one or more of the sequenced stages were sorted into nine categories (various colors), and are shown here for each instance of enrichment relative to the adjacent timepoints. See also Fig. 1D and Table S3. Figure S5. GO term enrichment among expression phases. Only transcripts whose expression is limited to a specific expression phase were considered for this analysis. The results are shown as an enrichment heat map. See also Table S4. Figure S6. k-means cluster analysis of the non-ubiquitous transcriptome reveals complete adherence to the expression phases. K-means clustering (k=19) was performed on annotated sequences whose expression is ≤ 1% of the maximum expression level per gene for at least one timepoint (n = 4792 sequences). Left: A heat map depicts the clusters; bright red is maximal expression. Right: Plots of the average normalized expression values for each cluster. Expression phase boundaries are indicated by vertical dashed lines. Clusters a2, b1, c1, d3 and e2 are also shown in Figure 4. See also Table S5. Figure S7. Cluster analysis of the ubiquitous transcriptome reveals clusters that do not adhere to the expression phases. K-means clustering (k=37) was performed on annotated sequences whose expression is > 1% of the maximum expression level for all timepoints (n = 14,773 sequences). Left: A heat map depicts the clusters; bright red is maximal expression. Right: Plots of the average normalized expression values for each cluster. Expression phase boundaries are indicated by vertical dashed lines. Cluster f1 is also shown in Figure 4. See also Table S6. Figure S8. Cluster analysis of the transcription factors reveals that most are expressed in three or more sequential stages. K-means clustering (k = 23) was performed on all the identified transcription factors in the Lv transcriptome (n = 521). Left: A heat map depicts the clusters; bright red is maximum expression. Right: Plots of the average normalized expression value for each cluster. Phase boundaries are indicated by vertical dashed lines. See also Figure 4 and Table S7. Figure S9. Temporal metabolic network analysis. A. The S. purpuratus-specific network, produced by the built-in filter in iPath. B.-L. The metabolic network expressed in each indicated developmental stage. For these maps, zero expression is depicted in black; expression from 0 – 100 is in light blue, 100 – 1000 is in medium blue, and >1000 is in dark blue. Thin gray lines indicate edges expected to be present based on the filter, but for which expression was not detected in any developmental stage. See also Figure 6. Table S1. qPCR Primer Sequences Gene 5' Primer (5' -> 3') 3' Primer (5' -> 3') LvAlx1 CCGGTACCAGTACCACCAAC AGGATGGGAGTGTTGACTGC LvBlimp GACGATGAGGACGAGTTGGT CGGTTTCCATTGCATCTTCT LvBMP2/4 CCATCAGGAGCTATCATCACG ATCACTGTGTGTCTGTGGTG LvBMP5-8 TACCCTCCACAACGATACACCC GCCTTTAGCAACGGACCATACA LvBrachyury CAGGCTAGAGGACGTGGAAC CCACTCCCCATTGACGATCT LvChordin CCCACTCAAGGCAGAGATAC GAAGGTACACCCTGGAACAC LvDelta CTGCCAAAACAACATCAACG AAACCTCCCATCACATCTCG LvEndo16 CGGATCGGAATCAGAGTCAT AGCAGCATCCAACTCCAAAT LvErg GGACACTTGCTAACCCGGTA TTCTCCCCATCTTCTGGCTA LvEts1 CGGAGGCATTGAGAATCAAG TAGGCATGCGTTCATCAGAG LvEve CATCCCCTATCACGATCACC TGTGGAGGCCAAACTCTAGG LvFoxG GAGCGATACAGAGCCGGAGA CTGAAGGGAGGTTTTTCACC LvFoxO AGCAATAGTTCAGCGGGTTG ACTCTCCCTCGCTTCCTCTC LvFoxQ2 GAGAAGCGTCTATTGCTGTGTG ACGAATACTGTTCCTCCAGCTC LvGataC GGGCCGAGAATGTGTAAACT AAGCGGTCTGTTCTGTCCAT LvGataE TAAACCTCAGCGCCAATACC CATTCCGATTCTGCTTTGCT LvGcm TCCAAAGACACCTGTCTACCAG GGAGAAGCACTGGAACAAAGAC LvGoosecoid TCACCGAAGAACAACTGGAA TAAGCCTCACATGCCCTCTT LvHh GGGAATCAAACTTCGACCAA AGATCCTTGCCTGGTCACTG LvHnf6 GGGAGGAGATCAACACGAAA CTGCCGCTCTTGAGCTTACT LvHox11/13b ACCTCATCCACCACCTCATC GAAGTTGCGATTGATGACCA LvIrxA CCAGGAATCGATGCGGTGATGG CATCTCGATGATCACGGTGGTGG LvNodal CCATGTTCATGGTGGAATTG GAATGTCGGGAAGGAGGTCT LvOtp CCCAACTCAACGAACTGGAG GAGGCCGTGAGATGGTAGAA LvOtx GCCAAACGTTGAACTCATGC AGCTGAGCTCGGGTAAAAGTC LvPtc AGGTTGTCGTGGAAACCAGT GTGGGGGTACTGATGGTGAG LvSetmar1 GCCATCATGTCCTTGTCTCA CACATGAAGCTTGATCAGGTAGTC LvSix1/2 CCCACAACCCTTACCCTTCT TTCCCATTCATAGGCAGGTC LvSix3 GGCTCAAGCCACAGGACTAA GGGTTGGGACTGAGTAAGCA LvSM30 CCACAGCTTCTGCAACAACA AGGCTGGATGGCGTAGTATG LvSmo GATGCTCCTGGTATGGTTGG GTGCTCAATCCCGCTATCAC LvSnail CCGAGGTCTTTTCTCGTCAA ACTTGCCACTGACACAACACC