Dynamic responses during early development of the sea urchin Strongylocentrotus intermedius to CO2-driven ocean acidification: A microRNA-mRNA integrated analysis.

To explore the dynamic molecular responses to CO2-driven ocean acidification (OA) during the early developmental stages of sea urchins, gametes of Strongylocentrotus intermedius were fertilized and developed to the four-armed larva stage in either natural seawater (as a control; pHNBS = 7.99 ± 0.01) or acidified conditions (ΔpHNBS = -0.3, -0.4, and - 0.5 units) according to the prediction for ocean pH by the end of this century. Specimens from five developmental stages (fertilization, cleavage, blastula, prism, and four-armed larva) were collected and comparative microRNA (miRNA) and mRNA transcriptome analyses were performed. The results showed that 1) a total of 22,224 differentially expressed genes (DEGs) and 51 differentially expressed miRNAs (DEMs) were identified in the OA-treated groups compared with the control group. 2) The numbers of both DEGs and DEMs were the largest at the blastula stage, indicating dramatic changes in gene expression. 3) Five "miR-1/DEG" modules were identified as potential biomarkers reflecting the response of sea urchins to OA during the early developmental period. 4) The PI3K/Akt signaling pathway was a key pathway involved in the response of S. intermedius to OA in its early developmental stages. This study deepens our understanding of the dynamic molecular regulatory mechanisms underlying sea urchin responses to CO2-driven OA.