Fig. 1

Expression profile of lung tissue from male P. californicus differing in bonding history. a. Unsupervised hierarchical clustering of whole transcriptome RNAseq data. Siblings are indicated by the same color coding. For example, B2, BD2, and V2 are siblings. Squared boxes indicate cousins. For example, B3 and BD3 are siblings and V3 is their cousin. b. Differential gene expression in whole transcriptome data. The number of differentially expressed genes is indicated. Gene IDs, fold change and adjusted P values are indicated in Supplementary Table 1. FDR cutoff: 0.1 and fold change ≥ 2. c. Among the four groups derived from K-means clustering, pathway enrichment analysis (GO analysis for biological processes) revealed consistent functions in group 1 (Supplementary Table 3). d. Bubble plot of top 10 enriched pathways. B, bonded; BD, bond-disrupted; V, Virgins. FPKM values were used, and the 1,000 most variable genes were considered