The role of heterochronic gene expression and regulatory architecture in early developmental divergence
Figures
![](https://iiif.elifesciences.org/lax:93062%2Felife-93062-fig1-v1.tif/full/617,/0/default.jpg)
Embryology and development of S. benedicti.
(A) Timeline of early development of the two morphs. The ‘swimming’ stage is the onset of swimming ability in both morphs. The ‘1-week’ old stage (7 days post fertilization) occurs approximately 1–4 days following larval release from maternal brood pouches under natural conditions. Planktotrophic (PP) larvae at this stage were not fed. Stages used in this study are marked with a red dot. (B) Adult male. The adults of both morphs are indistinguishable outside of reproductive traits. (C) Table of larval traits. Lecithotrophic (LL) offspring are facultatively feeding, meaning that feeding is not obligatory, but they are capable of feeding (Pernet and McArthur, 2006). Gut formation from Pernet and McHugh, 2010, where the three-band stage is our swimming stage and the five-band stage is shortly before our 1-week stage.
![](https://iiif.elifesciences.org/lax:93062%2Felife-93062-fig2-v1.tif/full/617,/0/default.jpg)
Differentially expressed (DE) genes across six developmental stages in the two developmental morphs.
(A) Sequenced libraries. Colors indicate morphs: green = PP (planktotrophic), orange = LL (lecithotrophic), purple = (F1) PL, pink = (F1) LP. Numbers in boxes are the sample’s sequencing depth in million reads. Total number of replicates for each time point are in bold. Samples for which all six stages were collected from the same cross are indicated as ‘complete’ with a check mark. (B) Principal component analysis of PP and LL samples. PC1 represents developmental time, whereas PC2 separates pre-gastrulation from post-gastrulation (C) DE between morphs over development. Black line is the percentage of total expressed genes that are DE at each stage, ribbon displays the relative mean and median log2 fold-change between morphs for those genes that are DE.
![](https://iiif.elifesciences.org/lax:93062%2Felife-93062-fig2-figsupp1-v1.tif/full/617,/0/default.jpg)
Read mapping rate between samples.
Mapping rates in all pairwise comparisons do not differ significantly. p-Values are shown above the plot indicating no significant differences between the mapping rates of any two groups.
![](https://iiif.elifesciences.org/lax:93062%2Felife-93062-fig2-figsupp2-v1.tif/full/617,/0/default.jpg)
Number of genes expressed at each stage in planktotrophic (PP) and lecithotrophic (LL) samples.
![](https://iiif.elifesciences.org/lax:93062%2Felife-93062-fig3-v1.tif/full/617,/0/default.jpg)
Classifications with examples of gene expression categories.
Functional heterochrony is a change in the developmental stage at which a particular pattern of gene expression appears instead of a difference in absolute time.
![](https://iiif.elifesciences.org/lax:93062%2Felife-93062-fig4-v1.tif/full/617,/0/default.jpg)
Six gene expression profiles representing the patterns of expression.
The number of conserved genes and heteromorphic genes that match each cluster is listed to summarize overall trends in expression patterns that appear mostly similar in both morphs. Clusters 2 (early expression) and 5 (late expression) represent the most genes.
![](https://iiif.elifesciences.org/lax:93062%2Felife-93062-fig4-figsupp1-v1.tif/full/617,/0/default.jpg)
Cluster number optimization.
The correlation between the two most similar gene expression clusters (profiles) as cluster number increases. Threshold (solid line) at r=0.85.
![](https://iiif.elifesciences.org/lax:93062%2Felife-93062-fig4-figsupp2-v1.tif/full/617,/0/default.jpg)
Clusters with all genes shown plotted on top of cluster core genes.
Red lines indicate genes with high cluster membership scores while blue lines indicate genes with lower membership scores. Genes with no cluster membership score greater than 0.5 were not assigned to a cluster and are not presented in this study’s results.
![](https://iiif.elifesciences.org/lax:93062%2Felife-93062-fig5-v1.tif/full/617,/0/default.jpg)
Magnitude of expression profile changes.
(A) Number of switches between each cluster. Clusters are arranged on the axes in order from high early developmental expression (left/bottom) to high late developmental gene expression (right/top). The diagonal represents genes which are assigned to the same cluster in both planktotrophic (PP) and lecithotrophic (LL) samples, boxes near the diagonal are switches between similar cluster profiles, and boxes furthest from the diagonal are extreme cluster switches. Genes expressed earlier in PP are above diagonal and earlier in LL are below. (B) Density plot of genes correlations between PP and LL. Heterochronic genes are yellow and heteromorphic genes in blue. Most heterochronies are minor changes in gene expression timing but the left-tailed distribution shows that few genes have large differences in timing.
![](https://iiif.elifesciences.org/lax:93062%2Felife-93062-fig6-v1.tif/full/617,/0/default.jpg)
Morph-specific gene expression by stage.
These genes are never expressed in the opposite morph, but may be expressed at one to six stages within a morph. Planktotrophic (PP) is green and lecithotrophic (LL) is orange.
![](https://iiif.elifesciences.org/lax:93062%2Felife-93062-fig7-v1.tif/full/617,/0/default.jpg)
Distribution of the 772 genes with differentiated expression between morphs.
Heteromorphy makes up nearly half of all differentiated genes, while heterochrony and morph-specific genes split the remainder nearly equally.
![](https://iiif.elifesciences.org/lax:93062%2Felife-93062-fig8-v1.tif/full/617,/0/default.jpg)
Differential gene expression including F1 offspring.
(A) Principal component analysis (PCA) of the top 500 most variable genes, including F1 (PL or LP) offspring. F1s have more misexpression in early stages but converge after the swimming stage and become intermediate to planktotrophic (PP) and lecithotrophic (LL) samples. (B) Plot relating F1 expression back to parental types. Dots indicate the number of DEGs between PP and LL samples at each stage. Purple bars represent the number of those DEGs for which the F1 gene expression is more like its maternal parent, and gray bars represent the number of those DEGs for which the F1 gene expression is more similar to its paternal parent.
![](https://iiif.elifesciences.org/lax:93062%2Felife-93062-fig8-figsupp1-v1.tif/full/617,/0/default.jpg)
Principal component analysis (PCA) of all gene expression labeled by family (crossID).
![](https://iiif.elifesciences.org/lax:93062%2Felife-93062-fig8-figsupp2-v1.tif/full/617,/0/default.jpg)
Genes with parental effects.
Most genes match maternal expression patterns.
![](https://iiif.elifesciences.org/lax:93062%2Felife-93062-fig9-v1.tif/full/617,/0/default.jpg)
Distribution of regulatory mode over time.
The majority of regulatory changes act in trans, and compensatory changes emerge after gastrulation.
Tables
Logical table showing criterion for morph-specific expression classification of genes based on number of samples in each group that have gene expression greater than their expression level threshold as assessed by data-adaptive flag method (DAFS).
An allowance for a single mismatched replicate is made in this criterion. This is because we do not believe that a single incidence within a morph (across all developmental stages) is a functionally significant and biologically relevant level of gene expression in this case.
PP samples above expression threshold | LL samples above expression threshold | Classification |
---|---|---|
≤1 | ≥3 | Lecithotrophic-specific expression |
≥3 | ≤1 | Planktotrophic-specific expression |
Comparison (Read Mapping Rates) | p-value of two sided T-test |
---|---|
PP versus LL | 0.1118 |
PP versus PL | 0.453 |
PL versus LP | 0.418 |
LP versus LL | 0.908 |
Additional files
-
Supplementary file 1
Gene ontology (GO) enrichment results spreadsheet.
- https://cdn.elifesciences.org/articles/93062/elife-93062-supp1-v1.xlsx
-
Supplementary file 2
Mode of regulatory change classifications.
- https://cdn.elifesciences.org/articles/93062/elife-93062-supp2-v1.docx
-
MDAR checklist
- https://cdn.elifesciences.org/articles/93062/elife-93062-mdarchecklist1-v1.docx