Cytosine-5 DNA methyltransferases (C5-MTases) and methyl-CpG-binding-domain (MBD) genes can be co-expressed. They directly control target gene expression by enhancing their DNA methylation levels in humans; however, the presence of this kind of cooperative relationship in plants has not been determined. A popular garden plant worldwide, petunia (Petunia hybrida) is also a model plant in molecular biology. In this study, 9 PhC5-MTase and 11 PhMBD proteins were identified in petunia, and they were categorized into four and six subgroups, respectively, on the basis of phylogenetic analyses. An expression correlation analysis was performed to explore the co-expression relationships between PhC5-MTases and PhMBDs using RNA-seq data, and 11 PhC5-MTase/PhMBD pairs preferentially expressed in anthers were identified as having the most significant correlations (Pearson’s correlation coefficients > 0.9). Remarkably, the stability levels of the PhC5-MTase and PhMBD pairs significantly decreased in different tissues and organs compared with that in anthers, and most of the selected PhC5-MTases and PhMBDs responded to the abiotic and hormonal stresses. However, highly correlated expression relationships between most pairs were not observed under different stress conditions, indicating that anther developmental processes are preferentially influenced by the co-expression of PhC5-MTases and PhMBDs. Interestingly, the nuclear localization genes PhDRM2 and PhMBD2 still had higher correlations under GA treatment conditions, implying that they play important roles in the GA-mediated development of petunia. Collectively, our study suggests a regulatory role for DNA methylation by C5-MTase and MBD genes in petunia anther maturation processes and multi-stress responses, and it provides a framework for the functional characterization of C5-MTases and MBDs in the future.