Objectives. Although myocardial fibrosis is a common pathophysiological process associated with many heart diseases, the molecular mechanisms regulating the development of fibrosis have not been fully determined. Recently, single cell RNA sequencing (scRNA-seq) analysis has been used to examine cellular fate and function during cellular differentiation and has contributed to elucidating the mechanisms of various diseases. The main purpose of this study was to characterize the fate of cardiac fibroblasts (CFs) and the dynamic gene expression patterns in a model of cardiac pressure overload using scRNA-seq analysis. Methods. The public scRNA-seq dataset of the transverse aortic coarctation (TAC) model in mice was downloaded from the GEO database, GSE155882. First, we performed quality control, dimensionality reduction, clustering, and annotation of the data through the Seurat R package (v4.0.5). Then, we constructed the pseudotime trajectory of cell development and identified key regulatory genes using the Monocle R package (v2.22.0). Different cell fates and groups were fully characterized by Gene Set Enrichment Analysis (GSEA) analysis and Transcription factor (TF) activity analysis. Finally, we used Cytoscape (3.9.1) to extensively examine the gene regulatory network related to cell fate. Results. Pseudotime analysis showed that CFs differentiated into two distinct cell fates, one of which produced activated myofibroblasts, and the other which produced protective cells that were associated with reduced fibrosis levels, increased antioxidative stress responses, and the ability to promote angiogenesis. In the TAC model, activated CFs were significantly upregulated, while protective cells were downregulated. Treatment with the bromodomain inhibitor JQ1 reversed this change and improved fibrosis. Analysis of dynamic gene expression revealed that Gpx3 was significantly upregulated during cell differentiation into protective cells. Gpx3 expression was affected by JQ1 treatment. Furthermore, Gpx3 expression levels were negatively correlated with the different levels of fibrosis observed in the various treatment groups. Finally, we found that transcription factors Jun, Fos, Atf3, and Egr1 were upregulated in protective cells, especially Egr1 was predicted to be involved in the regulation of genes related to antioxidant stress and angiogenesis, suggesting a role in promoting differentiation into this cell phenotype. Conclusions. The scRNA-seq analysis was used to characterize the dynamic changes associated with fibroblast differentiation and identified Gpx3 as a factor that might be involved in the regulation of myocardial fibrosis under cardiac pressure overload. These findings will help to further understanding of the mechanism of fibrosis and provide potential intervention targets.