Plasticity associated with adoption of social roles in clown anemonefish
- Department of Biology, Boston University, United States
- Marine Eco-Evo-Devo Unit, Okinawa Institute of Science and Technology, Japan
- Department of Ecology and Evolutionary Ecology, Cornell University, United States
Figures
Figure 1
Overview of experimental design and methods.
(A) Schematic of the experimental design. Replicate groups of juvenile clownfish (Amphiprion percula) housed as size-matched pairs (n=30 pairs) and size-matched solitary individuals (n=30 individuals). Standard length (SL) was measured at the beginning and end of the 5-week experiment, and individuals were denoted as dominant (P1), subordinate (P2), and solitary (S) at the end of the experiment based on their emergent social positions. (B) Feeding behavior scoring method for pairs and solitary fish. During feeding, pairs were separated by a divider to prevent competition; the divider was included in solitary tanks as a sham control, and individuals were provided with equal amounts of pre-portioned food daily. Feeding behavior was quantified weekly from video recordings by counting the number of bites taken by each fish over a 5 min period starting from the first bite. Initially, paired fish were denoted as P-max and P-min based on the total number of bites taken, because individuals could not be reliably identified due to their similar size and lack of distinguishing markings. From week 3 onwards, bites taken were scored separately for dominant (P1) and subordinate (P2) fish, in pairs where fish could be reliably identified based on clear size differences. (C) Agonistic behavior scoring method for pairs. Agonistic behaviors (aggression, submission) were quantified during a 5-min observation period from video recordings of paired interactions. Aggression was quantified by the number of chases, headbutts, and bites, while submissive behavior included flees and shakes. Initially, total aggression and submission were measured at the tank level, because fish were similar in size and lacked distinguishing markings. From week 3 onward, agonistic behaviors were scored separately for dominant (P1) and subordinate (P2) fish, in pairs where individuals could be reliably identified due to the emergence of clear size differences. In the accompanying illustrations, the focal fish is shown in orange, with its partner depicted in black and white. (D) Gene expression profiling for pairs and solitary fish. At the end of the experiment (week 5), whole-body samples were collected from a subset of individuals n=15 per social role: dominant (P1), subordinate (P2), solitary (S). To maximize detection of gene expression patterns associated with social roles, replicates with the largest size differences between dominant and subordinate fish were selected for transcriptomic profiling. The resulting RNA-seq data were used to compare molecular signatures across emergent social roles. In all panels, the timeline at the top indicates when (week 1 – week 5) and how (calipers, videos, tissue samples) data was collected.
Figure 2
Growth, feeding and agonistic behavior over time.
(A) Growth over time. Standard length (mm) of individuals of each social position (pair rank 1=P1; pair rank 2=P2; solitary = S; n=27 per social position) over time with mean and standard error of mean (SEM). Standard length differed by social position and time (GLMM with social positon and time as fixed effects and replicate ID as a random effect; time: DF = 1, F = 159.794, p < 0.0001; social position: DF = 2, F = 31.030, p < 0.0001; interaction: DF = 2, F = 10.534, p < 0.0001). (A.1) Overall growth (mm) over 5 weeks of each social position is shown with mean (dashed lines) and SEM. (B) Feeding behavior over time. Individual food intake, measured in number of bites of size-matched pairs (P-max and P-min) vs solitary (S) individuals over time with mean and SEM (n=27 pairs; n=27 solitaries in weeks 2-5). GLMM with social position and time as fixed effects and replicate ID as a random effect:time DF = 4, F = 27.6042, p < 0.0001; social position DF = 2, F = 58.3621, p < 0.0001; interaction DF = 7, F = 2.2745, p = 0.0288. (B.1) Food intake for a subset of replicates in which dominant (P1), subordinate (P2), and solitary (S) individuals could be reliably identified due to the emergence of clear size differences (week 3, n=5; week 4, n=17; week 5, n=19), shown over time with mean ± SEM. P1 ate more than P2 and S (social position DF = 2, F = 29.706, p < 0.0001; contrasts all p < 0.0001), which did not differ (p = 0.882) (C) Agonistic behavior of pairs over time. Aggression (the total number of chases, bites, and headbutts) and submission (based on the total number of flees and shakes) within pairs over time with mean and SEM (n=27 pairs). Agonistic behavior declined over time (aggression: DF = 4, F = 16.073; submission: DF = 4, F = 18.055; both p < 0.001) (C.1) Aggression and (C.2) Submission for a subset of pairs in which dominant (P1) and subordinate (P2) individuals could be reliably identified based on clear size differences (week 3, n=5; week 4, n=17; week 5, n=19), shown over time with mean and SEM. P1 aggression decreased (weeks 3–4, p = 0.0009; weeks 3–5, p = 0.0001) and P2 submission decreased (weeks 3–4, p = 0.0004; weeks 3–5, p < 0.0001). Data points shown jittered on all figures for visual aid, but collected at the same time.
Figure 3 with 3 supplements
Overall gene expression profiles showed no significant differences between P2 and S, compared to P1 individuals.
(A) Principal Component Analysis (PCA) of all genes across social positions, accounting for clutch ID (genetic background), with PERMANOVA results showing a significant effect of social position. Ellipses show the estimated confidence intervals from a 95% multivariate t-distribution. (B) Venn diagram illustrating differentially expressed genes (DEGs) from pairwise comparisons of social positions (FDR-adjusted p-value of 0.05). Orange upward arrows and corresponding numbers represent significantly upregulated DEGs in the first term, relative to the second, while blue downward arrows and numbers indicate downregulated DEGs (Appendix 2—table 3).
Figure 3—figure supplement 1
Principal component analysis (PCA) of gene expression across gene subsets and PC axis combinations.
Each panel shows a PCA scatter plot of individual samples (n=45) colored by social position (P1, P2, S), and shaped by genetic background (clutch ID). Rows represent gene subsets used for PCA: all genes (A–C), top 50% (D–F), top 25% (G–I), and top 5% most variable genes (J–L). Columns represent PC axis combinations: PC1 vs PC2 (left), PC2 vs PC3 (center), and PC1 vs PC3 (right). Ellipses show 95% confidence intervals estimated from a multivariate t-distribution for each social position. PERMANOVA p-values indicate the effect of social position on sample distribution, accounting for the genetic background of individuals (clutch ID). Panel B corresponds to Figure 3A in the main text.
Figure 3—figure supplement 2
Overall gene expression profiles showed no significant differences between P2 and S, compared to P1 individuals.
Principal Component Analysis (PCA) of all genes across social positions, accounting for clutch ID (genetic background), with PERMANOVA results showing a significant effect of social position. Samples are labeled, and replicates sharing the same numeric ID are connected by lines. Ellipses show the estimated confidence intervals from a 95% multivariate t-distribution.
Figure 3—figure supplement 3
Volcano plot of differential gene expression between pairwise comparisons of social positions (P1: dominant, P2: subordinate, S: solitary) of A. percula.
Pairwise comparisons of differentially expressed genes (DEGs) between (A) P1 and P2, (B) P1 and S, and (C) P2 and S individuals. Significant DEGs (adjusted p<0.05) are colored according to the social position in which they are upregulated. The top 10 most significant DEGs in each pairwise comparison are labeled with gene name or Ensembl gene ID.
Figure 4 with 2 supplements
Gene ontology (GO) enrichment of the brown module associated with growth and ossification.
An exhaustive list of GO enrichment results for genes assigned to the brown WGCNA module (Figure 5). Bubble size reflects numbers of genes associated with each term, color indicates GO category. Terms are ordered by adjusted p-value significance on the x-axis (–log₁₀ scale), with more significant terms on the right.
Figure 4—figure supplement 1
Modules of genes associated with phenotypes of interest.
(A) Heatmap (merging threshold of 0) of Weighted Gene Correlation Network Analysis (WGCNA) to correlate gene expression pattern with social position (P1, P2, S), growth, appetite, orange hue, and saturation of individuals. Only significant correlations are shown while non-significant correlations are marked with a dash. (B) The number of genes associated with each module found in WGCNA.
Figure 4—figure supplement 2
Gene Ontology (GO) terms identified in WGCNA modules showing correlation to phenotypes of interest.
Trees generated from GO enrichment analyses using Fisher’s Exact Tests (presence/absence in a module) to identify over-represented pathways within each module across GO divisions of Biological Processes (BP), Molecular Functions (MF), and Cellular Components (CC). Key GO terms, marked with stars, were explored in greater detail by extracting their associated genes and visualizing expression patterns of genes (FDR <0.05) using heatmaps. Terms with red stars indicated GO terms with clear patterns (Figure 5B), while yellow stars indicate terms which did not show clear patterns of DEGs dependent on social position (heatmaps not shown).
Figure 5
P1 individuals showed upregulation of genes associated with growth and ossification.
(A) Average eigengene values across social positions for the brown module (n=15 per social position). Points show individual samples, boxes show the median and interquartile range, and whiskers extend to 1.5 × the interquartile range. (B) Heatmap of significant DEGs (FDR P-value ≤0.05) of growth and ossification genes in the brown Biological Process (BP) and Molecular Function (MF) GO terms across the three social positions. Samples are grouped a priori by social position (dominant, P1; subordinate, P2; solitary, S) and ordered within each group by column dendrograms reflecting expression similarity within social positions. The row dendrograms show relatedness among genes based on expression profile similarity across all samples. Genes that are significantly differentially expressed between individuals of P1 vs P2 are denoted with an asterisk (*), those between P1 vs S with double asterisks (**), and three asterisks (***) represent both comparisons. Genes shown in bold and italics indicate significance at an FDR adjusted p-value <0.1.
Figure 6
P1 individuals showed downregulation of metabolic and appetite suppression genes.
Heatmap of significantly differentially expressed genes (DEGs) (FDR p-value ≤0.01) of genes associated with appetite regulation (APT), Krebs cycle (TCA), and glycolysis (GLY) pathways across the three social positions indicating genetic background of individuals (clutch ID). Samples (columns) are grouped a priori by social position (dominant, P1; subordinate, P2; solitary, S) and ordered within each group by the column dendrogram, which reflects expression similarity among samples within each social position. Genes (rows) are grouped a priori by pathway (APT, TCA, or GLY), and within each pathway group, the row dendrogram reflects relatedness among genes based on expression profile similarity across all samples. Genes that are significantly differentially expressed between individuals of P1 vs P2 are denoted with an asterisk (*), those between P1 vs S with double asterisks (**), and three asterisks (***) represent both comparisons. Genes shown in bold and italics indicate significance at an FDR adjusted p-value <0.1.
Appendix 1—figure 1
Clutch assignment validation via hierarchical clustering of microsatellite genotypes.
The dendrogram illustrates genetic relatedness among 36 individuals from three clutches, based on microsatellite genotypes (60 loci) processed using HotShot genotyping. (A) Initial assignment of individuals to their presumed clutches, with pairs incorrectly assigned highlighted in blue and orange. Branch lengths reflect genetic divergence; tightly clustered branches indicate sibling relationships within clutches. (B) Corrected clustering after reassignment of clutches for mismatched pairs.
Appendix 1—figure 2
Clutch assignment validation using TagSeq-derived SNPs and hierarchical clustering.
The dendrogram illustrates genetic relatedness among 45 individuals from three clutches, based on SNPs from TagSeq data. (A) Initial assignment of individuals to their presumed clutches, with pairs incorrectly assigned highlighted in colors. Branch lengths reflect genetic divergence; tightly clustered branches indicate sibling relationships within clutches. (B) Corrected clustering after reassignment of clutches for mismatched pairs.
Appendix 3—figure 1
P1 individuals showed upregulation of genes associated with growth and ossification.
Heatmap of significant DEGs (FDR p-value ≤0.05) of growth and ossification genes in the brown Biological Process (BP) and Molecular Function (MF) GO terms across the three social positions. The bottom panel shows the size ratio (P2 SL/P1 SL) at week 5 for each replicate pair. Black bars indicate pairs with a size ratio ≤0.8 (more dissimilar in size), and red bars indicate pairs with a size ratio >0.8 (more similar in size). The dashed line denotes a size ratio of 0.8. Samples (columns) are grouped a priori by social position (dominant, P1; subordinate, P2; solitary, S) and ordered within each group by the column dendrogram, which reflects expression similarity among samples within each social position. The gene (rows) dendrogram reflects relatedness among genes based on expression profile similarity across all samples. Genes that are significantly differentially expressed between individuals of P1 vs P2 are denoted with an asterisk (*), those between P1 vs S with a double asterisk (**), and three asterisks (***) represent both comparisons. Genes shown in bold and italics indicate significance at an FDR adjusted p-value <0.1.
Appendix 3—figure 2
P1 individuals showed downregulation of appetite suppression and metabolic genes.
Heatmap of significantly differentially expressed genes (DEGs) (FDR p-value ≤0.01) of genes associated with glycolysis (GLY), Krebs cycle (TCA), and appetite regulation (APT), TCA cycle, and glycolysis pathways across the three social positions. The bottom panel shows the size ratio at week 5 (P2 SL/P1 SL) for each replicate pair. Black bars indicate pairs with a size ratio ≤0.8 (greater size divergence), and red bars indicate pairs with a size ratio >0.8 (more similar in size). The dashed line denotes a size ratio of 0.8. Samples (columns) are grouped a priori by social position (dominant, P1; subordinate, P2; solitary, S) and ordered within each group by the column dendrogram, which reflects expression similarity among samples within each social position. Genes (rows) are grouped a priori by pathway (GLY, TCA, or APT), and within each pathway group, the row dendrogram reflects relatedness among genes based on expression profile similarity across all samples. Genes that are significantly differentially expressed between individuals of P1 vs P2 are denoted with an asterisk (*), those between P1 vs S with a double asterisks (**), and three asterisks (***) represent both comparisons. Genes shown in bold and italics indicate significance at an FDR adjusted p-value <0.1.
Appendix 3—figure 3
Principal component analysis (PCA) of gene expression in paired individuals (P1, P2) based on heatmaps of (A) growth, (B) appetite, and metabolic genes.
PCA was performed on specific gene sets (growth: Figure 5B; appetite and metabolism: Figure 6), (P1 and P2; n=30) using row Z-score scaled expression values, consistent with the heatmap scaling. Each point represents one individual. Ellipses represent 95% confidence intervals around each social position group. A PERMANOVA (adonis2; 999 restricted permutations; shown bottom right corners) found that size ratio was not a significant predictor of gene expression variation after accounting for social position and genetic background (A: psize_ratio = 0.858; B: psize_ratio = 0.844).
Appendix 4—figure 1
Orange hue and saturation over time.
(A) Hue of individuals of three social positions (pair rank 1=P1; pair rank 2=P2; solitary = S; n=27 per social position) over time with mean and SEM. Hue decreased in all social positions (LMM with social position and time as fixed effects and replicate ID as a random effect; time: F₁ = 639.88, p < 0.001; social position: F₂ = 12.14, p < 0.001). (B) Saturation of individuals of social positions (pair rank 1=P1; pair rank 2=P2; solitary = S; n=27 per social position) over time with mean and SEM. Saturation increased in all social positions (time: F₁ = 870.04, p < 0.01), with no significant effect of social position. Data points shown staggered for visual aid, but collected at the same time.
Author response image 1
Principal component analysis (PCA) of gene expression based on heatmaps of (A) growth, (B) appetite, and metabolism genes.
PCA was performed on gene sets of the main text (growth: Fig 5B and appetite and metabolism: Fig. 6), using row Z-score scaled expression values, consistent with the heatmap scaling. Each point represents one individual. Ellipses represent 95% confidence intervals around each social position group. PERMANOVA results (adonis2; shown in the bottom right corners).
Tables
Appendix 2—table 1
Post-hoc pairwise contrasts of total aggression (upper triangle) and total submission (lower triangle), displaying p-values of all comparisons.
Significant differences are denoted with *.
| Week 1 | Week 2 | Week 3 | Week 4 | Week 5 | |
|---|---|---|---|---|---|
| Week 1 | 0.688 | 0.026* | <0.001* | <0.001* | |
| Week 2 | 0.285 | 0.686 | <0.001* | <0.001* | |
| Week 3 | <0.001* | 0.1469 | 0.028* | 0.002* | |
| Week 4 | <0.001* | 0.002* | 0.4039 | 0.937 | |
| Week 5 | <0.001* | <0.001* | 0.1477 | 0.9845 |
Appendix 2—table 2
Mapped read count data.
A. percula replicate, social position (P1 – pair rank 1; P2 – pair rank 2; S – solitary) and clutch ID. Read counts of raw, filtered and trimmed sequences, and number of uniquely mapped reads against A. percula genome using STAR, as well as the corresponding mapping efficiency.
| Sample ID | Social position | Clutch ID | Raw read counts | Trimmed read counts | Uniquely mapped reads | Mapping efficiency (%) |
|---|---|---|---|---|---|---|
| 1 | P1 | L3 | 5,986,911 | 2,737,307 | 2,352,774 | 85.952 |
| 1 | P2 | L2 | 6,803,145 | 2,954,216 | 2,554,648 | 86.475 |
| 1 | S | L6 | 5,837,482 | 2,648,669 | 2,294,955 | 86.646 |
| 2 | P1 | L6 | 6,217,888 | 2,787,974 | 2,394,766 | 85.896 |
| 2 | P2 | L3 | 7,010,907 | 3,100,482 | 2,669,479 | 86.099 |
| 2 | S | L2 | 8,019,177 | 3,585,504 | 3,072,414 | 85.690 |
| 3 | P1 | L2 | 8,383,004 | 3,518,423 | 2,955,511 | 84.001 |
| 3 | P2 | L6 | 7,386,810 | 3,360,910 | 2,862,178 | 85.161 |
| 3 | S | L3 | 7,249,814 | 3,348,687 | 2,858,570 | 85.364 |
| 5 | P1 | L6 | 7,202,343 | 3,143,829 | 2,699,683 | 85.872 |
| 5 | P2 | L2 | 6,390,202 | 2,808,877 | 2,416,787 | 86.041 |
| 5 | S | L3 | 5,400,980 | 2,350,517 | 2,050,627 | 87.242 |
| 6 | P1 | L3 | 6,649,594 | 2,507,798 | 2,160,966 | 86.170 |
| 6 | P2 | L2 | 7,091,099 | 3,238,004 | 2,789,147 | 86.138 |
| 6 | S | L6 | 7,237,912 | 3,074,694 | 2,646,890 | 86.086 |
| 7 | P1 | L6 | 7,204,499 | 2,999,504 | 2,547,277 | 84.923 |
| 7 | P2 | L2 | 6,970,070 | 2,899,158 | 2,483,166 | 85.651 |
| 7 | S | L3 | 6,461,090 | 2,879,935 | 2,473,878 | 85.900 |
| 10 | P1 | L6 | 6,168,327 | 2,471,974 | 2,128,137 | 86.091 |
| 10 | P2 | L3 | 6,778,638 | 2,834,810 | 2,446,430 | 86.300 |
| 10 | S | L2 | 6,303,813 | 2,753,496 | 2,356,609 | 85.586 |
| 11 | P1 | L2 | 7,700,902 | 2,871,152 | 2,433,468 | 84.756 |
| 11 | P2 | L6 | 10,410,675 | 3,926,127 | 3,378,190 | 86.044 |
| 11 | S | L3 | 8,457,694 | 3,369,689 | 2,904,659 | 86.200 |
| 13 | P1 | L2 | 8,747,232 | 2,756,847 | 2,388,352 | 86.633 |
| 13 | P2 | L3 | 11,355,428 | 3,983,796 | 3,392,123 | 85.148 |
| 13 | S | L6 | 8,493,446 | 3,498,055 | 2,976,781 | 85.098 |
| 14 | P1 | L3 | 9,649,613 | 3,377,460 | 2,892,559 | 85.643 |
| 14 | P2 | L6 | 11,096,303 | 3,812,633 | 3,307,234 | 86.744 |
| 14 | S | L2 | 7,458,006 | 2,925,835 | 2,543,736 | 86.941 |
| 15 | P1 | L6 | 8,315,175 | 2,770,012 | 2,448,710 | 88.401 |
| 15 | P2 | L2 | 7,809,486 | 3,122,052 | 2,657,505 | 85.120 |
| 15 | S | L3 | 11,188,485 | 4,419,431 | 3,841,559 | 86.924 |
| 20 | P1 | L6 | 10,293,955 | 3,756,124 | 3,194,578 | 85.050 |
| 20 | P2 | L3 | 9,793,618 | 3,759,257 | 3,248,705 | 86.419 |
| 20 | S | L2 | 10,601,642 | 3,926,835 | 3,360,769 | 85.585 |
| 21 | P1 | L3 | 9,434,952 | 3,630,988 | 3,124,380 | 86.048 |
| 21 | P2 | L6 | 9,464,753 | 3,554,760 | 3,070,925 | 86.389 |
| 21 | S | L2 | 10,612,954 | 3,949,458 | 3,418,245 | 86.550 |
| 23 | P1 | L3 | 6,565,976 | 2,610,543 | 2,201,455 | 84.329 |
| 23 | P2 | L2 | 10,032,443 | 4,035,965 | 3,453,739 | 85.574 |
| 23 | S | L6 | 10,407,184 | 3,981,383 | 3,457,368 | 86.838 |
| 26 | P1 | L6 | 7,843,938 | 2,824,540 | 2,461,006 | 87.129 |
| 26 | P2 | L2 | 7,756,092 | 3,046,595 | 2,636,393 | 86.536 |
| 26 | S | L3 | 10,632,681 | 4,163,850 | 3,607,742 | 86.644 |
Appendix 2—table 3
List of significantly differentially expressed genes shared between the pairwise comparisons of P1 vs P2 and P1 vs S (main text Figure 3B).
A total of 25 genes were consistently downregulated in P1 (relative to both P2 and S), and 84 genes were consistently upregulated in P1. Gene names are shown where available; otherwise, Ensembl gene IDs are provided.
| Shared downregulated genes | Shared upregulated genes | |||
|---|---|---|---|---|
| c1qtnf5 | BRSK2 | IGFBP1 | SHMT2 | ENSAPEG00000005722 |
| COL10A1 | aldh6a1 | IQSEC3 | slc13a5a | ENSAPEG00000006596 |
| col1a1a | aldh8a1 | LAMP5 | slc16a9a | ENSAPEG00000015523 |
| col1a2 | aldob | LECT2 | slc6a1b | ENSAPEG00000015722 |
| dpt | ambp | LRTM2 | SLC7A5 | ENSAPEG00000006747 |
| hspa1b | ATCAY | mapk8a | sybu | ENSAPEG00000008011 |
| loxa | C8A | minpp1a | syt5b | ENSAPEG00000012490 |
| MRC2 | cabp1b | MPP4 | ucp1 | ENSAPEG00000016217 |
| mustn1b | cacna1da | OGDHL | ENSAPEG00000000371 | ENSAPEG00000020211 |
| mylpfa | camk1gb | oxct1a | ENSAPEG00000003493 | |
| myo1ca | CHIA | pax10 | ENSAPEG00000004162 | |
| PRKAB2 | cp | pax6b | ENSAPEG00000011721 | |
| pvalb4 | defbl2 | pck2 | ENSAPEG00000011919 | |
| runx3 | diabloa | PCYT2 | ENSAPEG00000013967 | |
| s100a1 | fam219aa | pde6ha | ENSAPEG00000018654 | |
| tcf3a | ftcd | PLPPR1 | ENSAPEG00000021138 | |
| tcf7 | GABRG2 | ppa1a | ENSAPEG00000022463 | |
| TNNI2 | gatd3l | prph2lb | ENSAPEG00000024362 | |
| VAT1 | glrbb | rbp2b | ENSAPEG00000006913 | |
| ENSAPEG00000006693 | glud1b | rcvrna | ENSAPEG00000016229 | |
| ENSAPEG00000008023 | got1 | rnf34b | ENSAPEG00000022341 | |
| ENSAPEG00000008236 | HABP2 | saga | ENSAPEG00000023250 | |
| ENSAPEG00000010334 | hlfa | SCG2 | ENSAPEG00000024133 | |
| ENSAPEG00000021032 | hpda | SERPINA4 | ENSAPEG00000024384 | |
| ENSAPEG00000023977 | hsd11b1la | SEZ6L | ENSAPEG00000001601 | |
Appendix 2—table 4
Complete list of Amphiprion percula genes associated with appetite regulation, glycolysis, and the TCA cycle, retrieved based on orthologs with A. ocellaris genes (Herrera et al., 2025).
Orthologs were identified using a reciprocal BLAST approach, and genes present in the whole-body gene expression dataset were retained. The table includes each gene and its associated p-values from pairwise contrasts between social positions (P1 – pair 1; P2 – pair 2; S – solitary). Genes with p-value ≤0.01, which were visualized in the heatmap shown in the main text (Figure 6), are indicated in bold and italics.
| Appetite regulation | Glycolysis | TCA cycle gene | ||||||
|---|---|---|---|---|---|---|---|---|
| gene ID | pval P1 vs P2 | pval P1 vs S | gene ID | pval P1 vs P2 | pval P1 vs S | gene ID | pval P1 vs P2 | pval P1 vs S |
| adipor1a | 0.6761 | 0.4666 | aldob | 0.0000 | 0.0000 | aco1 | 0.1639 | 0.2966 |
| adipor2 | 0.9723 | 0.7837 | aldocb | 0.3631 | 0.1511 | aco2a | 0.3443 | 0.5425 |
| adycap1a | 0.6028 | 0.6839 | eno1a | 0.0393 | 0.0653 | aco2b | 0.9936 | 0.1806 |
| adycap1b | 0.0184 | 0.0062 | eno1b | 0.3382 | 0.4826 | cs | 0.4455 | 0.4812 |
| agrp | 0.3415 | 0.2067 | eno2b | 0.0017 | 0.0001 | dld | 0.1028 | 0.2988 |
| cart1 | 0.6801 | 0.7373 | eno3 | 0.5279 | 0.5356 | dlsta | 0.8086 | 0.1922 |
| cart2 | 0.2357 | 0.8397 | eno4 | 0.7267 | 0.7755 | dlstb | 0.3865 | 0.8454 |
| cart3 | 0.0484 | 0.0634 | gapdha | 0.2781 | 0.5315 | fh | 0.8748 | 0.7770 |
| cart4 | 0.5422 | 0.1872 | gapdhb | 0.0304 | 0.0045 | idh1 | 0.0110 | 0.0604 |
| ccka | 0.9695 | 0.2773 | gapdhs | 0.5494 | 0.2076 | idh2 | 0.2936 | 0.1994 |
| cckar | 0.5921 | 0.7164 | gpia | 0.7274 | 0.8939 | idh3aa | 0.6119 | 0.3906 |
| cckb | 0.1367 | 0.0032 | gpib | 0.9137 | 0.5803 | idh3ab | 0.0753 | 0.0919 |
| cckbra | 0.4533 | 0.7310 | hk1a | 0.6127 | 0.6745 | idh3b | 0.2814 | 0.2014 |
| cckbrb | 0.1600 | 0.2561 | hk1b | 0.0982 | 0.0356 | idh3ga | 0.2682 | 0.3490 |
| cnr1 | 0.2586 | 0.1884 | hk2 | 0.9152 | 0.3750 | idh3gb | 0.9080 | 0.4846 |
| dbi | 0.1883 | 0.6002 | pfkl | 0.8873 | 0.3654 | mdh1aa | 0.0115 | 0.0003 |
| gal | 0.6902 | 0.9579 | pfkma | 0.9350 | 0.6929 | mdh1ab | 0.9899 | 0.8136 |
| gcga | 0.2682 | 0.5551 | pfkmb | 0.5924 | 0.4384 | ogdha | 0.4576 | 0.6432 |
| gcgb | 0.4859 | 0.9787 | pfkpa | 0.5099 | 0.1077 | ogdhb | 0.0053 | 0.0006 |
| ghrl | 0.0383 | 0.0022 | pfkpb | 0.1420 | 0.0234 | ogdhl | 0.0001 | 0.0011 |
| gnrhr4 | 0.9896 | 0.6204 | pgam1b | 0.6447 | 0.2204 | sdha | 0.1765 | 0.0034 |
| hcrt | 0.7410 | 0.2989 | pgam2 | 0.0213 | 0.0273 | sdhaf2 | 0.2860 | 0.0719 |
| ins | 0.4233 | 0.9990 | pgk1 | 0.0865 | 0.3382 | sdhaf3 | 0.2866 | 0.1430 |
| insra | 0.0108 | 0.0038 | pklr | 0.4707 | 0.2840 | sdhaf4 | 0.0480 | 0.2568 |
| insrb | 0.6936 | 0.5334 | pkma | 0.0773 | 0.1232 | sdhba | 0.6299 | 0.4958 |
| irs1 | 0.6567 | 0.4022 | pkmb | 0.5434 | 0.3802 | sdhbb | 0.2943 | 0.0387 |
| irs2a | 0.0098 | 0.0014 | tpi1a | 0.0004 | 0.0000 | sdhc | 0.0139 | 0.1278 |
| irs2b | 0.0278 | 0.0083 | tpi1b | 0.8185 | 0.8390 | sdhd | 0.5440 | 0.0946 |
| kiss1rb | 0.8019 | 0.8638 | sucla2 | 0.0136 | 0.1696 | |||
| lepa | 0.6358 | 0.2605 | suclg1 | 0.0047 | 0.0253 | |||
| lepr | 0.5288 | 0.7118 | suclg2 | 0.0071 | 0.0182 | |||
| mchr2 | 0.7329 | 0.5194 | ||||||
| nmu | 0.0486 | 0.4880 | ||||||
| npy | 0.0748 | 0.0179 | ||||||
| nr5a1a | 0.5679 | 0.2892 | ||||||
| nucb2a | 0.5676 | 0.3864 | ||||||
| nucb2b | 0.1028 | 0.4651 | ||||||
| oxt | 0.9410 | 0.5240 | ||||||
| pmch | 0.8590 | 0.5869 | ||||||
| pomcaa | 0.5722 | 0.5171 | ||||||
| pomcab | 0.0677 | 0.7574 | ||||||
| pomcb | 0.2758 | 0.7570 | ||||||
| pyy | 0.3126 | 0.0305 | ||||||
| scg2 | 0.0514 | 0.0330 | ||||||
| slc6a3 | 0.7674 | 0.7345 | ||||||
| spx | 0.4817 | 0.9035 | ||||||
| sst | 0.0617 | 0.0416 | ||||||
| TRH | 0.2166 | 0.1003 | ||||||
Additional files
Download links
A two-part list of links to download the article, or parts of the article, in various formats.
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Plasticity associated with adoption of social roles in clown anemonefish
eLife 15:RP109871.
https://doi.org/10.7554/eLife.109871.3
