Next Article in Journal
Gene-Based Resistance to Erysiphe Species Causing Powdery Mildew Disease in Peas (Pisum sativum L.)
Next Article in Special Issue
Epigenetic and Physiological Responses to Varying Root-Zone Temperatures in Greenhouse Rocket
Previous Article in Journal
GUCY2D-Related Retinal Dystrophy with Autosomal Dominant Inheritance—A Multicenter Case Series and Review of Reported Data
Previous Article in Special Issue
Identification of Genomic Regions Controlling Chalkiness and Grain Characteristics in a Recombinant Inbred Line Rice Population Based on High-Throughput SNP Markers
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 13
Issue 2
10.3390/genes13020314
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploration of the Potential Transcriptional Regulatory Mechanisms of DNA Methyltransferases and MBD Genes in Petunia Anther Development and Multi-Stress Responses

by
Lisha Shi
1,2,
Huimin Shen
1,2,
Jiawei Liu
1,2,
Hongmin Hu
1,2,
Hongyan Tan
1,2,
Xiulian Yang
1,2,
Lianggui Wang
1,2 and
Yuanzheng Yue
1,2,*
1
Key Laboratory of Landscape Architecture, College of Landscape Architecture, Nanjing Forestry University, Nanjing 210037, China
2
Co-Innovation Center for Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing 210037, China
*
Author to whom correspondence should be addressed.
Genes 2022, 13(2), 314; https://doi.org/10.3390/genes13020314
Submission received: 16 December 2021 / Revised: 28 January 2022 / Accepted: 1 February 2022 / Published: 8 February 2022
(This article belongs to the Special Issue Plant Genomics and Epigenomics in Breeding for Yield, Quality, and Sustainability)
Browse Figures
Review Reports Versions Notes

Abstract

:
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.

1. Introduction

DNA methylation is an important epigenetic regulatory feature, and it plays critical roles in multiple plant biology processes by regulating the expression levels of target genes. Additionally, it can be stably inherited by offspring without altering primary DNA sequences [1]. DNA methylation mainly exists in three forms: 5-methylcytosine, 6-methyadenine and 4-methylcytosine, among which the latter two are the predominant forms in prokaryotes, and the former is the prevalent form in plants [2,3]. DNA methylation (5-methylcytosine) usually occurs in three cytosine contexts, CpG, CpHpG and CpHpH, and it is mainly regulated by cytosine-5 DNA methyltransferases (C5-MTases), which can transfer the methyl group from S-adenosyl-l-methionine to the 5 position of the pyrimidine [4]. DOMAINS REARRANGED METHYLASE 2 (DRM2), METHYLTRANSFERASE 1 (MET1), and CHROMOMETHYLASE 3 (CMT3) are three important members of the C5-MTase family, and the establishment of DNA methylation depends on DRM2, whereas MET1 and CMT3 are mainly responsible for maintaining DNA methylation [5].
Methyl-CpG-binding-domain (MBD) proteins can specifically bind to the symmetrically methylated CpG sites of the up-stream regions of the target genes and recruit C5-MTases to silence gene expression by increasing DNA methylation [6]. MBD family members contain a conserved MBD domain, which is composed of approximately 70 aa [7]. The MBD protein family is divided into six subclasses in plants, and it exhibits two evolutionary branches, dicot and monocot, among which subclasses VI and IV are only present in the former, indicating dicot specificity [8,9]. Studies on the MBD family in plants have mainly focused on model plants, such as Arabidopsis, rice, maize and populous, which have between 13 and 18 MBD members [9,10,11]. In comparison, knowledge on MBD proteins in ornamental crops is still limited.
C5-MTase and MBD proteins can be cooperatively expressed and regulate many important biological activities in mammals. For example, DNMT1, which is a member of the C5-MTases, cooperates with MBD3 in human to down-regulate TREM2 expression and cause neuronal injury by enhancing the DNA methylation level of the TREM2 cis-acting element region [12]. MeCP2 is a member of the human MBD family, and it interacts directly with DNMT1. The resulting MeCP2-DNMT1 complex plays a critical role in genome-wide DNA methylation and affects the differentiation of human dental pulp cells. A similar synergy between MeCP2 and DNMT1 has been confirmed in rat cancer cells [13,14]. However, we find that there is no potential interaction relationship between C5-MTase and MBD proteins in plants by searching in BioGRID [15].
DNA methylation plays an important regulatory role in plant developmental processes, such as leaf growth [16], fruit ripening [17], seed development [18] and vernalization [19]. Remarkably, DNA methylation also participates in anther developmental processes, including tapetum differentiation and tapetum programmed cell death (PCD) [20,21,22]. In Brassica napus and Nicotiana tabacum, DNA methylation may promote the progress of tapetum PCD. When tapetum cells are active, the DNA methylation level is low; however, as the tapetum enters the degradation stage, the DNA methylation level begins to increase, and reaches its highest level during the tapetum PCD stage [23]. In Malus, the level of DNA methylation was found to be increased during the anther developmental process, and highly modified DNA methylation is a characteristic of tapetum differentiation [24]. The genes NtMET1 and BnMET1 are members of the C5-MTase family, and they participate in the developmental processes from microspore to early pollen. They are responsible for the increase in DNA methylation during tapetum PCD processes, indicating that MET1 plays an essential role in the maintenance and establishment of DNA methylation during anther development [23]. Thus, DNA methylation plays critical roles in plant growth and anther development.
In plants, DNA methylation is closely associated with responses to various abiotic stresses, and the expression levels of some C5-MTase and MBD proteins are induced by drought and cold stresses [17,25]. The major plant hormones, such as IAA and GA, also play important roles in defense-related stress responses, but studies on the plant hormone-based regulation of DNA methylation levels are limited [26]. Particularly, GA is also involved in plant growth and differentiation, especially in anther development, which is closely related to male fertility and fertilization [20,21]. The C5-MTase and MBD proteins have been identified in Arabidopsis thaliana [11,27], Zea mays [25,28], Solanum lycopersicum [10,17] and Triticum aestivum [29,30]. However, to the best of our knowledge, no research has yet focused on the identification of C5-MTase and MBD proteins in Petunia.
Petunia, which is one of the world’s most famous garden plants, has a long-standing history as a model for molecular biological studies [31,32]. Here, we identified and characterized the C5-MTase and MBD proteins in the whole genome of petunia, including phylogenetic, protein-conserved motif, gene structure and expression pattern analyses. The expression correlations of PhC5-MTases and PhMBDs in anther developmental processes and different tissues, along with the expression correlations under multiple abiotic and hormonal stresses, were analyzed, and one PhC5-MTase and PhMBD gene pair that can be co-expressed was selected and subjected to subcellular localization experiments. Our study provides valuable clues to the functional identification and evolutionary relationships of C5-MTase and MBD proteins, and it indicates that they have a synergistic effect in plants.

2. Materials and Methods

2.1. Identification of C5-MTase and MBD Genes

The members of C5-MTase and MBD families of petunia were identified from the SGN repository database (https://solgenomics.net/organism/Petunia_axillaris/genome, accessed on 8 January 2021). The Hidden Markov model (HMM) profiles of DNA methylase domain (PF00145) and MBD domains (PF01429) were downloaded from the PFAM database (https://pfam.xfam.org/, accessed on 2 January 2021), and HMMER software (version 3.0) was used to scan all the C5-MTase and MBD protein sequences (E-value ≤ 1 × 10−10) [33]. The web tools Batch CD-search [34] and SMART [35] were used to detect the integrity of the conserved domains of the candidate genes, and genes with incomplete conserved domains were manually removed. The molecular weights and theoretical isoelectric point values of C5-MTase and MBD proteins were calculated using ExPASy (https://web.expasy.org/compute_pi/, accessed on 7 February 2021). The database BioGRID was used to identify interactions between C5-MTase and MBD proteins [15].

2.2. Phylogenetic and Gene Structure Analysis

The neighbor-joining phylogenetic trees of C5-MTase and MBD proteins were constructed using MEGA 11 software (bootstrap value = 1000), and the protein sequences of 17 species used in the evolutionary tree were obtained from Phytozome (http://phytozome.jgi.doe.gov/pz/portal.html, accessed on 5 May 2021); they are listed in Tables S1 and S2. The sequences of C5-MTase and MBD proteins were visualized and analyzed using Jalview (version 2.0) [36]. MEME was used to identify the preserved motifs [37], and the structures of C5-MTase and MBD genes were visualized using TBtools software [38]. The gene structures were constructed by IBS 1.0 [39].

2.3. Plant Materials and Stress Treatments

Petunia ‘W115’ seedlings were grown under the following conditions in a controlled growth chamber: light/dark of 16/8 h, day/night temperature of 23 °C/22 °C, relative air humidity of 60% and light intensity of 200 μmol m−2s−1. Anthers were collected when longitudinal lengths of buds reached 0.2, 0.3, 0.5, 1.2 and 3.5 cm. Leaves, roots, stem leaves, petals, and the four whirls of the floral organs (pistil, anther, petal and sepal) were collected at the 0.3 cm flower bud stage and used in the RNA-seq analysis. Petunia ‘W115’ seedlings were selected at 25 days of age for stress treatments, including 10% PEG 6000, cold (4 °C), IAA (100 μmol/L) and GA (100 μmol/L), which were performed in accordance with methods described by a previous study [40]. For the cold treatment, the petunia seedlings were transported to a 4 °C growth chamber. For drought stress, the petunia seedlings were irrigated with a 10% PEG 6000 solution. For hormone treatments, IAA and GA were applied using a spray application method, in which they were sprayed onto the leaf surface until it was sufficiently wet but without droplet condensation [41]. Each stress treatment had three biological replicates. The above-ground tissues of the seedlings were collected at 0, 3 and 24 h. All the samples were immediately frozen in liquid nitrogen and stored at −80 °C for subsequent experiments.

2.4. Isolation of RNA, Real-Time PCR, and Gene Expression Analyses

An EASYspin Plus Plant RNA Kit (Aidlab, Beijing, China) was used to extract total RNA from petunia. Total RNA quantity was assayed with the spectrophotometer method (Denovix, DS-11 FX+, Wilmington, DE, USA) for absorbance ratio of 260/280 (1.8 to 2.0). Total RNA was reverse-transcribed using the TransScript One-Step gDNA Remover and cDNA Synthesis SuperMix (Transgen, Beijing, China), following the manufacturer’s protocol at 42 °C for 30 min and 85 °C for 5 s, and was carried out on a Veriti 96 Well Thermal Cycler (Applied Biosystems™, Thermo Fisher Scientific Inc, Waltham, MA, USA). Approximately 5 μg of RNA was used for cDNA was utilized to synthesize the first-strand cDNA. The cDNA was diluted 1:20 with ddH2O. Next, 1 μL of the dilution was used in qRT-PCR analyses [42]. The qRT-PCR reaction reagent was prepared using TB Green™ Premix Ex Taq™ (TaKaRa, Kusatsu, Japan). Each reaction (10 μL) consisted of 5 μL TB Green Premix Ex Taq II, 0.4 μL of 10 μM forward and reserve primers, 1 μL of diluted cDNA, 3 μL ddH2O, and 0.2 μL ROX reference dye. The qRT-PCR reactions were carried out on a on Step One Plus real-time PCR (Applied Biosystems™, Thermo Fisher Scientific Inc.). The conditions were 95 °C for 30 s, followed by 40 cycles of 95 °C for 5 s and 60 °C for 30 s. Each qRT-PCR analysis had three biological and technical replicates. On the basis of the Ct values of PhC5-MTase and PhMBD genes, the raw data were quantified using the 2−ΔΔCt method by Excel 2019 [43,44]. The level of β-actin transcripts was selected as a reference [45]. All the real-time PCR reactions were carried out on an Applied Biosystems 7500 Real-Time PCR System, according to the manufacture’s instructions. The PCR conditions were 94 °C for 4 min, followed by 32–35 cycles of 94 °C for 30 s, 55–60 °C for 30 s, 72 °C for 1–3 min, and a final extension of 72 °C for 10 min. The qRT-PCR primers were designed using Primer 5.0 and are listed in Table S3 [46]. A heatmap was constructed using the TBtools software with the RPKM data, which was retrieved from the transcriptome data (Table S4).

2.5. Subcellular Localization

The coding sequences of PhDRM2 and PhMBD2 without the stop codons were cloned independently into the pCAMBIA Super 1300-GFP vector at the HindIII and KpnI restriction sites to construct 35S::GFP-PhDRM2/PhMBD2. For transient expression investigation, tobacco (Nicotiana benthamiana) leaves were infiltrated with Agrobacterium tumefaciens (GV3101) harboring 35S::GFP-PhDRM2/PhMBD2 and a control vector. The infiltrated plants were grown for 3 days before being stained with the DAPI and observed using a LSM710 microscope (Zeiss, Germany) for the green fluorescent protein (GFP) fluorescence signal.

3. Results

3.1. Identification of PhC5-MTase and PhMBD Genes

To identify PhC5-MTase and PhMBD genes across the petunia genome, the HMM profiles and the amino acid sequences from Arabidopsis were used as simultaneous queries to search for homologous petunia genes (E-values < 0.001). In this study, nine PhC5-MTases were identified; each contained a complete DNA methylase domain. They were divided into four subfamilies: three CMT, one MET, four DRM and one DNMT (Table S5). The coding polypeptides of PhC5-MTases ranged from 397 (PhDNMT2) to 1,557 (PhMET1) aa, the theoretical isoelectric points varied from 4.80 (PhDRM2) to 6.32 (PhCMT2) and the protein molecular weights ranged from 45.22 kDa (PhDNMT2) to 174.84 kDa (PhMET1) (Table S5). All the PhC5-MTase genes contained conserved DNA methylase domains at the C-termini, but there were different conserved domain combinations at the N-termini among the subfamilies. The MET group members contained two bromo adjacent homology (BAH) domains, but only one BAH domain was observed in CMT group members. The members of the DRM group possessed two UBA domains. Furthermore, C5-MTases that lacked a conserved region at the N-terminal belonged to the DNMT2 group (Figure S1).
In this study, the protein-coding genes containing a conserved MBD domain were considered to be members of MBD family in petunia (Figure S2). After deleting the mispredicted genes, 11 MBD proteins were identified. The molecular weights of the PhMBDs varied from 19.70 kDa to 213.81 kDa, and the theoretical isoelectric points ranged from 4.43 kDa to 9.92 kDa. The lengths varied from 167 to 1948 aa, and they were designated as PhMBD1-11 (Table S6). The petunia genome has not been assembled at the chromosomal level; therefore, the PhC5-MTases and PhMBDs were named in accordance with the scaffold order.

3.2. Phylogenetic Analysis and Classification

To understand the evolutionary relationships and history of C5-MTase and MBD genes in plants, the full-length protein sequences of 15 (A. thaliana, Brachpodium distachyon, Cynara cardunculus, Fragaria ananassa, Glycine max, Oryza sativa, P. hybrida, Populus trichocarpa, Ricinus communis, Sorghum bicolor, S. lycopersicum, Salvia miltiorrhiza, Solanum melanogena, Solanum tuberosum and Z. mays) and 7 (A. thaliana, Capsicum annuum, O. sativa, P. hybrida, S. lycopersicum, S. tuberosum and T. aestivum) species were used to construct the rootless evolutionary trees of the C5-MTase and MBD genes, respectively. On the basis of the similarities in protein structure and sequence, the evolutionary tree of C5-MTase was divided into four groups, corresponding to the MET, CMT, DNMT2 and DRM subfamilies (Figure 1a). The MBD family was divided into six subclasses (I, II, III, IV, V and VI) on the basis of the sequence similarity within the MBD motif. Each C5-MTase and MBD subfamily was further separated into two groups, dicot and monocot, indicating that the evolution of the C5-MTase and MBD genes may differ between dicotyledons and monocotyledons (Figure 1b). In terms of the phylogenetic relationships, PhC5-MTases and PhMBDs were closely related to their orthologs in eggplant, potato, tomato and pepper, suggesting that the C5-MTase and MBD families may be conserved within Solanaceae.

3.3. Multiple-Sequence Alignment and Motif Analysis

To understand the conservation and divergence of the C5-MTase and MBD proteins, we performed a multiple-sequence alignment and analyzed their gene structures and motif compositions using MEME. In total, 15 common motifs of between 28 and 100 aa were predicted in the C5-MTase family (Table S7). The motif positions of the C5-MTase family members were conserved within each class. The C-terminal DNA methyltransferase domain contained motifs 1, 2, 3, 5, 6, 9 and 13, and there also six highly conserved motifs (I, IV, VI, VIII, IX and X). These six highly conserved sequences were present in all the members of petunia C5-MTase family, which can serve as guides for the identification of potential C5-MTases, and their order of distribution in the four subfamilies differed (Figure S3). Among them, the motifs I and X can bind to AdoMet, and the motif IV, which contains a proline-cysteine doublet, is a functional catalytic site for C5-MTases [47]. Motif VI can provide a glutamate, which is critical for the binding of the target cytosine [48]. Motif VIII can make non-specific contact with cytosine and motif IX is associated with the organization of the target recognition domain (TRD) [28]. The N-terminal domain contained the remaining eight motifs, for the MET subfamily, motifs 10 and 11 represented two RFD domains, that targeting proteins to replication foci in a non-catalytic context [49]. Motifs 2 and 8 represented two BAH domains in the MET subfamily, and motif 14 corresponded to the BAH domain in the CMT family, which links DNA methylation, replication, and transcription [50]. The DRM subfamily contained the unique motifs 3, 7 and 12 (Figure 2a). The exon–intron organization of PhC5-MTases was examined further to demonstrate the structural diversity of PhC5-MTases. The number of introns in PhC5-MTases varied greatly, from 8 to 20, with the majority of PhC5-MTases having eight introns (Figure 2a). To confirm the structural properties of C5-MTase and MBD proteins, the conserved domains of PhC5-MTases and PhMBDs were analyzed. In this study, 15 common motifs were predicted in the MBD proteins (Table S8). Each MBD protein contained three to six motifs, and the lengths of the motifs were between 18 and 46 aa. The motif positions of PhMBDs revealed that all of the MBD proteins contained motif 1, and seven MBD proteins contained motif 5. Motifs 1 and 5 corresponded to the conserved MBD domain that has the function of binding to symmetrically methylated CpG dinucleotides [51] (Figure 2b). The exon–intron patterns of the PhMBD family were studied to learn more about the petunia MBD family’s evolution. Except for PhMBD4, which has 15 introns, most PhMBDs contained one to five introns (Figure 2b).

3.4. Expression Analysis of PhC5-MTases and PhMBDs

To investigate the transcriptional abundance levels of all the selected PhC5-MTase and PhMBD genes, their expression profiles were analyzed using RNA-seq data. The clustered heatmap showed that most of the genes were expressed in all the petunia tissues and showed any specificity in different tissues/stages (Figure 3a,b). Some genes expressed preferentially in anthers, such as PhDRM2 and PhMBD2 (Figure 3a,b). Five genes (PhDRM4, PhDNMT2, PhDRM1, PhMBD7 and PhMBD8) had the lowest expression levels in all the tissues. In addition, correlations between the expression levels of PhC5-MTase and PhMBD genes among the five anther developmental stages and different tissues were explored. The expression levels of the PhC5-MTase and PhMBD genes were significantly positively correlated in anthers, and the 11 pairs of genes (PhMBD1 and PhCMT1/CMT2/DRM3, PhDRM2 and PhMBD2/MBD3/MBD1 and PhMBD10 and PhCMT3/MET1/CMT1/CMT2/DRM3) with the highest Pearson’s correlation coefficient (PCC) were selected (Figure 3d) (PCCs > 0.9, p < 0.05). However, the significance levels of PhC5-MTase and PhMBD generally decreased in different organs and tissues, with the pairs PhMBD1 and PhDRM3 (0.91) /PhCMT2 (0.89) /PhDRM2 (0.85), PhMBD10 and PhDRM3 (0.90) /PhCMT3 (0.89) /PhMET1 (0.85), PhMBD2 and PhDRM2 (0.87) still maintaining significant positive correlations (PCCs > 0.85, p < 0.05) (Figure 4).

3.5. Expression Analysis of Phc5-Mtases and Phmbds in Response to Abiotic and Hormonal Treatments

It has been reported that the C5-MTase and MBD proteins are induced by drought and cold stresses in plants, and the plant hormones (such as IAA and GA) often play important roles in defense-related stress responses [17,25,26].
To further determine whether the expression levels of the PhC5-MTase and PhMBD genes were influenced by abiotic and hormonal stresses, the 11 PhC5-MTase and PhMBD gene pairs, which consisted of five PhC5-MTases and five PhMBDs, having the highest Pearson’s correlation coefficients in anthers, were examined under four abiotic and hormonal stresses: cold, drought (PEG), IAA and GA treatments. We used qRT-PCR to analyze the gene expression patterns in the petunia seedlings, and the expression levels of all the selected PhC5-MTase and PhMBD genes were induced by a variety of abiotic and hormonal stresses (Figure 5a and Figure 6a). The expression levels of all the PhC5-MTase and PhMBD genes were regulated to different degrees under various conditions, and most PhC5-MTase and PhMBD genes were up-regulated under all four stresses, except PhCMT2. The correlations of the PhC5-MTase and PhMBD gene pairs under multiple stresses were further analyzed, and the correlations between the members of the 11 PhC5-MTase and PhMBD gene pairs generally decreased under abiotic stresses and hormone treatments (Figure 5b and Figure 6b). Remarkably, some gene pairs maintained significant positive correlations under one or more stresses. PhMBD1 and PhCMT1 maintained a positive correlation (0.92) under cold stress but not under drought stress. PhMBD1 and PhDRM3 maintained a significant positive correlation (0.85) under IAA-treatment conditions, whereas significant positive correlations, such as PhDRM2 and PhMBD2 (0.85) were observed under GA treatment conditions.

3.6. Subcellular Localization of PhDRM2 and PhMBD2

PhDRM2 and PhMBD2 are preferentially expressed in anthers, and they are significantly positively correlated in anthers (0.92) and all tissues (0.87), as well as under GA-treatment conditions (0.85). To explore the potential functions of the PhDRM2 and PhMBD2 genes, Agrobacterium cultures containing 35S::GFP-PhDRM2 and 35S::GFP-PhMBD2 constructs were injected into the tobacco leaves for subcellular localization studies. As shown in Figure 7, similar to the empty 35S::GFP, the GFP fluorescence of 35S::GFP-PhDRM2 and 35S::GFP-PhMBD2 were both found in the nucleus.

4. Discussion

DNA methylation is a common phenomenon in plants, which can lead to gene and transposon silencing, and also participates in plant biological processes, which affects many aspects of plant growth and development [1]. As the main regulators of DNA methylation and the interpreter of DNA methylation sites, C5-MTase and MBD proteins, respectively, affect the development of plants by regulating the methylation level and participating in responses to multiple stresses [8,52], and have been identified in maize [25,28] and tomato [10,17]. Petunia is one of the world’s best known ornamental plants and a model for molecular biology research [31]. However, C5-Mtase and MBD proteins have not been identified in petunia.
In the present study, 9 PhC5-MTases and 11 PhMBDs were found in the petunia genome, and the multiple-sequence alignment showed that each of the PhC5-MTase contained a complete DNA methylase domain at the C-terminus (Figure S1), and all the PhMBD genes contained an MBD domain (Figure S2), indicating that they had the typical characteristics of C5-MTase and MBD family members, respectively. Interestingly, despite the genome size of petunia (1.4 GB) being almost 1.5 times greater than that of tomato (900 MB), 1.6 times greater than that of potato (844 MB), 3 times greater than that of rice (466 MB) and 11 times greater than that of Arabidopsis (125 MB), the number of C5-MTase members in petunia is very similar, ranging from 7 to 11, indicating that there is no strong correlation between the number of C5-MTase members and genome size. However, unlike the PhC5-MTases, the number of MBD genes in petunia was less than in tomato (18), potato (15), rice (17) and Arabidopsis (14), and the number of DNA methylation-related genes is associated with genome-wide duplication and gene loss events [10,25,53]. Consequently, we suspect that gene duplication and loss events during evolution may have contributed to the difference in the number of MBD genes, which is consistent with the research on maize [25].
On the basis of the phylogenetic tree, both C5-MTase and MBD families were classified into two groups, dicot and monocot (Figure 1), indicating that the evolution of C5-MTase and MBD genes may differ between dicotyledons and monocotyledons, which is similar to previous research [9]. The PhC5-MTases are classified into four subfamilies (Figure 1a), MET, DRM, CMT and DNMT2, which is consistent with the results obtained for Arabidopsis, tomato, eggplant and strawberry [17,54,55]. However, there are four MET subfamily genes in Arabidopsis, but only one in petunia, suggesting that the MET subfamily of genes were lost during petunia evolution. In addition, most of the C5-MTase and MBD genes of petunia, potato and tomato in the phylogenetic tree tended to cluster on one branch (Figure 1), indicating that the C5-MTase and MBD genes are highly homologous and conserved in Solanaceae.
DNA methylation plays important roles in many aspects of plant developmental processes, such as leaf growth [16], fruit ripening [17] and seed development [18]. DRM1 and DRM2 are members of the C5-MTase family, and Arabidopsis thaliana drm1drm2 mutants display retarded growth and late flowering phenotypes [56]. MBD11 is a member of the MBD family, and Arabidopsis plants with inhibited AtMBD11 display a phenotype of reduced fertility and abnormal leaf and flower development [11]. Remarkably, DNA methylation also participates in the processes of anther development in Malus and tobacco, including tapetum differentiation and tapetum PCD [23,24]. The expression levels of C5-MTase and MBD genes showed obvious tissue-specificity in tomato and strawberry, respectively, and play different functions in plant growth and development [10,55]. Similarly, we investigated the expression levels of C5-MTase and MBD families in different tissues of petunia, and found that most of the genes were expressed in all the tissues (Figure 3a,b), suggesting that C5-MTase and MBD genes participate in plant biological processes. Some of the C5-MTase and MBD genes exhibited floral organ-specific expression, especially in anthers. Furthermore, the expression levels of most PhC5-MTase and PhMBD genes were significantly positively correlated in anthers, and the 11 pairs of genes with the highest PCCs were selected (Figure 3d) (PCCs > 0.9, p < 0.05). However, the correlation coefficients generally decreased in different organs and tissues, with only two gene pairs maintaining significant positive correlations (Figure 4), implying that the co-expression of PhC5-MTases and PhMBDs preferably influenced anther development.
DNA methylation is regulated by various abiotic stresses, such as cold and drought [57,58]. As the key factors in DNA methylation, C5-MTase and MBD genes are sensitive to various abiotic stresses in many plants, and their expression levels are greatly altered by cold and drought stresses [10,59]. Here, we evaluated the expression levels of the 11 PhC5-MTase/PhMBD pairs with the highest PCCs in anthers in response to abiotic and hormonal treatments (Figure 5a and Figure 6a). Most of the selected PhC5-MTase and PhMBD genes were induced by multiple stresses, and their expression levels increased under cold and drought stresses (Figure 5a), which is similar to in other plants [10,60]. In addition, plant hormones play important roles in defense responses, but studies on DNA methylation induced by plant hormones are limited [26]. In this study, in response to hormonal treatments, most of the PhC5-MTase and PhMBD genes were up-regulated under GA-treatment conditions, but some were not significantly changed under IAA-treatment conditions (Figure 6a), implying that the PhC5-MTase and PhMBD genes were preferentially influenced by GA. Furthermore, the correlations between the members of the 11 pairs of PhC5-MTase and PhMBD genes generally decreased under abiotic stresses and hormone treatments, with only a few pairs maintaining significant positive correlations (Figure 5b and Figure 6b), implying that the highly correlated expression levels in anthers were unstable under multiple stresses. Remarkably, the two nuclear localization genes (Figure 7), PhDRM2 and PhMBD2, were preferentially expressed during anther development and had a significant positive correlation in anthers (0.92), which was maintained under GA-treatment conditions (0.85) (Figure 6b), implying that they were involved in the regulation of anther development and play important roles in GA-mediated development in petunia.

5. Conclusions

In this study, we performed a systematic genome-wide analysis of the C5-MTase and MBD gene families in petunia. A total of 9 PhC5-MTase and 11 PhMBD genes were identified and categorized into four and six subgroups, respectively, on the basis of the phylogenetic analyses. The motifs and gene structures of the PhC5-MTases and PhMBDs within the subfamilies were similar, which validated the predicted classifications. Using the RNA-seq data, 11 pairs of PhC5-MTases and PhMBDs with the highest PCCs in anthers were identified. Interestingly, although most of the selected pairs had positive correlative relationships in different tissues and organs, the significance levels generally decreased under multi-stress conditions. The qRT-PCR results indicated that the expression of the PhC5-MTase and PhMBD gene pairs could be induced by multiple abiotic and hormonal stresses; however, the correlative relationships were unstable, with only a few pairs maintaining positive correlations. Remarkably, the PhDRM2 and PhMBD2 pair had high correlation values in both anthers and under GA-treatment conditions. They also both localized to the nucleus, implying that they play important roles in the GA-mediated development in petunia. Thus, our study indicates a potential epigenetic regulatory mechanism in petunia that functions through the co-operative expression of PhC5-MTases and PhMBDs, and the results provide a basis for further research on the biological functions and evolution of C5-MTase and MBD genes.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes13020314/s1. Figure S1: Schematic structures of PhC5-MTase genes. Figure S2: Multiple-sequence alignments of MBD genes. Figure S3: Multiple-sequence alignments of C5-MTase and MBD genes. Table S1: Sequences of C5-MTase proteins in the 15 tested species. Table S2: Sequences of MBD proteins in the 7 tested species. Table S3: Primers of PhC5-MTase and PhMBD genes used in the experiments. Table S4: The RPKM data of PhC5-MTase and PhMBD genes. Table S5: Nomenclature and physicochemical properties of PhC5-Mtases. Table S6: Nomenclature and physicochemical properties of PhMBD genes. Table S7: Fifteen different motifs commonly observed in PhC5-MTase genes. Table S8: Fifteen different motifs commonly observed in PhMBD genes.

Author Contributions

Methodology, Y.Y.; software, J.L.; validation, H.S.; formal analysis, H.H. and H.T.; writing—original draft preparation, L.S.; writing—review and editing, L.S. and Y.Y.; funding acquisition, Y.Y, X.Y. and L.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant Nos. 31601785 and 32071828), the Natural Science Fund for Colleges and Universities in Jiangsu Province (Grant No. 20KJB220005), and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hirsch, S.; Baumberger, R.; Grossniklaus, U. Epigenetic variation, inheritance, and selection in plant populations. Cold Spring Harbor Symp. Quant. Biol. 2012, 77, 97–104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Zhang, H.; Lang, Z.; Zhu, J.-K. Dynamics and function of DNA methylation in plants. Nat. Rev. Mol. Cell Biol. 2018, 19, 489–506. [Google Scholar] [CrossRef] [PubMed]
  3. Ye, P.; Luan, Y.; Chen, K.; Liu, Y.; Xiao, C.; Xie, Z. MethSMRT: An integrative database for DNA N6-methyladenine and N4-methylcytosine generated by single-molecular real-time sequencing. Nucleic Acids Res. 2017, 45, D85–D89. [Google Scholar] [CrossRef] [Green Version]
  4. Bender, J. DNA methylation and epigenetics. Ann. Rev. Plant Biol. 2004, 55, 41–68. [Google Scholar] [CrossRef] [PubMed]
  5. Wang, P.; Gao, C.; Bian, X.; Zhao, S.; Zhao, C.; Xia, H.; Song, H.; Hou, L.; Wan, S.; Wang, X. Genome-wide identification and comparative analysis of cytosine-5 DNA methyltransferase and demethylase families in wild and cultivated Peanut. Front. Plant Sci. 2016, 7, 7. [Google Scholar] [CrossRef] [Green Version]
  6. Zemach, A.; Grafi, G. Methyl-CpG-binding domain proteins in plants: Interpreters of DNA methylation. Trends Plant Sci. 2007, 12, 80–85. [Google Scholar] [CrossRef] [PubMed]
  7. Song, X.; Zhang, Y.; Zhong, Q.; Zhan, K.; Bi, J.; Tang, J.; Xie, J.; Li, B. Identification and functional characterization of methyl-CpG binding domain protein from Tribolium Castaneum. Genomics 2020, 112, 2223–2232. [Google Scholar] [CrossRef]
  8. Grafi, G.; Zemach, A.; Pitto, L. Methyl-CpG-binding domain (MBD) proteins in plants. Biochim. Biophys. Acta Gene Struct. Expr. 2007, 1769, 287–294. [Google Scholar] [CrossRef]
  9. Springer, N.M.; Kaeppler, S.M. Evolutionary divergence of monocot and dicot Methyl-CpG-dinding domain proteins. Plant Physiol. 2005, 138, 92–104. [Google Scholar] [CrossRef] [Green Version]
  10. Parida, A.P.; Raghuvanshi, U.; Pareek, A.; Singh, V.; Kumar, R.; Sharma, A.K. Genome-wide analysis of genes encoding MBD domain-containing proteins from tomato suggest their role in fruit development and abiotic stress responses. Mol. Biol. Rep. 2018, 45, 2653–2669. [Google Scholar] [CrossRef]
  11. Berg, A.; Meza, T.J.; Mahić, M.; Thorstensen, T.; Kristiansen, K.; Aalen, R.B. cgene family encoding methyl-CpG-binding domain proteins are transcriptionally active and at least one, AtMBD11, is crucial for normal development. Nucleic Acids Res. 2003, 31, 5291–5304. [Google Scholar] [CrossRef] [Green Version]
  12. Wang, H.; Feng, Y.; Sun, J.; Zhang, W.; Han, Z.; Yu, S.; Gu, Y.; Cheng, X.; Lin, Z.; Na, M. Methyl-CpG-binding domain protein 3 promotes seizures by recruiting Methyltransferase DNMT1 to enhance TREM2 methylation. Neurochem. Res. 2021, 46, 2451–2462. [Google Scholar] [CrossRef]
  13. Li, J.; Deng, Q.; Fan, W.; Zeng, Q.; He, H.; Huang, F. Melatonin-induced suppression of DNA methylation promotes odontogenic differentiation in human dental pulp cells. Bioengineered 2020, 11, 829–840. [Google Scholar] [CrossRef] [PubMed]
  14. Guan, Y.; Guan, X.; An, H.; Baihetiya, A.; Wang, W.; Shao, W.; Yang, H.; Wang, Y. Epigenetic silencing of miR-137 induces resistance to bicalutamide by targeting TRIM24 in prostate cancer cells. Am. J. Transl. Res. 2019, 11, 3226–3237. [Google Scholar] [PubMed]
  15. Chatraryamontri, A.; Oughtred, R.; Boucher, L.; Rust, J.; Chang, C.; Kolas, N.K.; O’Donnell, L.; Oster, S.; Theesfeld, C.; Sellam, A.; et al. The BioGRID interaction database: 2017 update. Nucleic Acids Res. 2016, 45, D369–D379. [Google Scholar] [CrossRef]
  16. Candaele, J.; Demuynck, K.; Mosoti, D.; Beemster, G.T.S.; Inzé, D.; Nelissen, H. Differential methylation during maize leaf growth targets developmentally regulated genes. Plant Physiol. 2014, 164, 1350–1364. [Google Scholar] [CrossRef] [Green Version]
  17. Cao, D.; Ju, Z.; Gao, C.; Mei, X.; Fu, D.; Zhu, H.; Luo, Y.; Zhu, B. Genome-wide identification of cytosine-5 DNA methyltransferases and demethylases in Solanum Lycopersicum. Gene 2014, 550, 230–237. [Google Scholar] [CrossRef] [PubMed]
  18. Xing, M.; Zhang, Y.; Zhou, S.; Hu, W.; Wu, X.; Ye, Y.; Wu, X.; Xiao, Y.; Li, X.; Xue, H. Global analysis reveals the crucial roles of DNA methylation during rice seed development. Plant Physiol. 2015, 168, 1417–1432. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Finnegan, E.J.; Kovac, K.A.; Jaligot, E.; Sheldon, C.C.; James Peacock, W.; Dennis, E.S. The downregulation of FLOWERING LOCUS C (FLC) expression in plants with low levels of DNA methylation and by vernalization occurs by distinct mechanisms. Plant J. Cell Mol. Biol. 2005, 44, 420–432. [Google Scholar] [CrossRef]
  20. Goto, N.; Pharis, R.P. Role of gibberellins in the development of floral organs of the gibberellin-deficient mutant, ga1-1, of Arabidopsis Thaliana. Can. J. Bot. 1999, 77, 944–954. [Google Scholar]
  21. Hirano, K.; Aya, K.; Hobo, T.; Sakakibara, H.; Kojima, M.; Shim, R.A.; Hasegawa, Y.; Ueguchi-Tanaka, M.; Matsuoka, M. Comprehensive transcriptome analysis of phytohormone biosynthesis and signaling genes in microspore/pollen and tapetum of rice. Plant Cell Physiol. 2008, 49, 1429–1450. [Google Scholar] [CrossRef] [PubMed]
  22. Vanyushin, B. Enzymatic DNA methylation is an epigenetic control for genetic functions of the cell. Biochemistry 2005, 70, 488–499. [Google Scholar] [CrossRef] [PubMed]
  23. Solís, M.-T.; Chakrabarti, N.; Corredor, E.; Cortés-Eslava, J.; Rodríguez-Serrano, M.; Biggiogera, M.; Risueño, M.C.; Testillano, P.S. Epigenetic changes accompany developmental programmed cell death in tapetum cells. Plant Cell Physiol. 2014, 55, 16–29. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Forino, L.M.C.; Andreucci, A.C.; Tredici, I.D.; Felici, C.; Giraldi, E.; Tagliasacchi, A.M. DNA methylation of tapetum cells during microsporogenesis in Malus domestica Borkh. Isr. J. Plant Sci. 2003, 51, 91–100. [Google Scholar] [CrossRef]
  25. Qian, Y.; Ren, Q.; Jiang, L.; Zhang, J.; Chen, C.; Chen, L. Genome-wide analysis of maize MBD gene family and expression profiling under abiotic stress treatment at the seedling stage. Plant Biotechnol. Rep. 2020, 14, 323–338. [Google Scholar] [CrossRef]
  26. Verma, V.; Ravindran, P.; Kumar, P.P. Plant hormone-mediated regulation of stress responses. BMC Plant Biol. 2016, 16, 86. [Google Scholar] [CrossRef] [Green Version]
  27. Ogneva, Z.V.; Dubrovina, A.S.; Kiselev, K.V. Age-associated alterations in DNA methylation and expression of methyltransferase and demethylase genes in Arabidopsis Thaliana. Biol. Plant 2016, 60, 628–634. [Google Scholar] [CrossRef]
  28. Qian, Y.; Xi, Y.; Cheng, B.; Zhu, S. Genome-wide identification and expression profiling of DNA methyltransferase gene family in maize. Plant Cell Rep. 2014, 33, 1661–1672. [Google Scholar] [CrossRef]
  29. Gahlaut, V.; Samtani, H.; Khurana, P. Genome-wide identification and expression profiling of cytosine-5 DNA methyltransferases during drought and heat stress in wheat (Triticum aestivum). Genomics 2020, 112, 4796–4807. [Google Scholar] [CrossRef]
  30. Li, Y.; Meng, F.; Yin, J.; Liu, H.; Si, Z.; Ni, Z.; Sun, Q.; Ren, J.; Niu, H. Isolation and comparative expression analysis of six MBD genes in wheat. Biochim. Biophys. Acta Gene Regul. Mech 2008, 1779, 90–98. [Google Scholar] [CrossRef]
  31. Bombarely, A.; Moser, M.; Amrad, A.; Bao, M.; Bapaume, L.; Barry, C.S.; Bliek, M.; Boersma, M.R.; Borghi, L.; Bruggmann, R.; et al. Insight into the evolution of the Solanaceae from the parental genomes of Petunia Hybrida. Nat. Plants 2016, 2, 16074. [Google Scholar] [CrossRef] [Green Version]
  32. Vandenbussche, M.; Chambrier, P.; Rodrigues Bento, S.; Morel, P. Petunia, your next supermodel? Front. Plant Sci. 2016, 7, 72. [Google Scholar] [CrossRef]
  33. Johnson, L.S.; Eddy, S.R.; Portugaly, E. Hidden markov model speed heuristic and iterative HMM search procedure. BMC Bioinf. 2010, 11, 431. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Marchler-Bauer, A.; Bryant, S.H. CD-Search: Protein domain annotations on the fly. Nucleic Acids Res. 2004, 32, W327–W331. [Google Scholar] [CrossRef] [PubMed]
  35. Letunic, I.; Doerks, T.; Bork, P. SMART: Recent updates, new developments and status in 2015. Nucleic Acids Res. 2014, 43, D257–D260. [Google Scholar] [CrossRef] [PubMed]
  36. Waterhouse, A.M.; Procter, J.B.; Martin, D.M.; Clamp, M.; Barton, G.J. Jalview Version 2—A multiple sequence alignment editor and analysis workbench. Bioinformatics 2009, 25, 1189–1191. [Google Scholar] [CrossRef] [Green Version]
  37. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. MEME Suite: Tools for motif discovery and searching. Nucleic Acids Res. 2009, 37, W202–W208. [Google Scholar] [CrossRef]
  38. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  39. Liu, W.; Xie, Y.; Ma, J.; Luo, X.; Nie, P.; Zuo, Z.; Lahrmann, U.; Zhao, Q.; Zheng, Y.; Zhao, Y.; et al. IBS: An illustrator for the presentation and visualization of biological sequences. Bioinformatics 2015, 31, 3359–3361. [Google Scholar] [CrossRef] [Green Version]
  40. Ali, F.; Qanmber, G.; Wei, Z.; Yu, D.; Li, Y.h.; Gan, L.; Li, F.; Wang, Z. Genome-wide characterization and expression analysis of geranylgeranyl diphosphate synthase genes in cotton (Gossypium spp.) in plant development and abiotic stresses. BMC Genom. 2020, 21, 561. [Google Scholar] [CrossRef]
  41. Liu, H.; Wang, L.; Jing, X.; Chen, Y.; Hu, F. Functional analysis of CgWRKY57 from Cymbidium goeringii in ABA response. PeerJ 2021, 9, e10982. [Google Scholar] [CrossRef] [PubMed]
  42. Yue, Y.; Yin, C.; Guo, R.; Peng, H.; Yang, Z.; Liu, G.; Bao, M.; Hu, H. An anther-specific gene PhGRP is regulated by PhMYC2 and causes male sterility when overexpressed in petunia anthers. Plant Cell Rep. 2017, 36, 1401–1415. [Google Scholar] [CrossRef] [PubMed]
  43. Pfaffl, M.W. A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acids Res. 2001, 29, e45. [Google Scholar] [CrossRef] [PubMed]
  44. Yang, X.; Zhou, T.; Wang, M.; Li, T.; Wang, G.; Fu, F.-F.; Cao, F. Systematic investigation and expression profiles of the GbR2R3-MYB transcription factor family in ginkgo (Ginkgo biloba L.). Int. J. Biol. Macromol. 2021, 172, 250–262. [Google Scholar] [CrossRef] [PubMed]
  45. Yue, Y.; Du, J.; Li, Y.; Thomas, H.R.; Frank, M.H.; Wang, L.; Hu, H. Insight into the petunia Dof transcription factor family reveals a new regulator of male-sterility. Ind. Crops Prod. 2021, 161, 113196. [Google Scholar] [CrossRef]
  46. Lalitha, S. Primer premier 5. Biotech. Softw. Internet Rep. 2000, 1, 270–272. [Google Scholar] [CrossRef]
  47. Jeltsch, A. Molecular enzymology of mammalian DNA methyltransferases. In DNA Methylation: Basic Mechanisms; Doerfler, W., Böhm, P., Eds.; Springer: Berlin/Heidelberg, Germany, 2006; pp. 203–225. [Google Scholar]
  48. Pavlopoulou, A.; Kossida, S. Plant cytosine-5 DNA methyltransferases: Structure, function, and molecular evolution. Genomics 2007, 90, 530–541. [Google Scholar] [CrossRef] [Green Version]
  49. Rountree, M.R.; Bachman, K.E.; Baylin, S.B. DNMT1 binds HDAC2 and a new co-repressor, DMAP1, to form a complex at replication foci. Nat. Genet. 2000, 25, 269–277. [Google Scholar] [CrossRef]
  50. Callebaut, I.; Courvalin, J.-C.; Mornon, J.-P. The BAH (bromo-adjacent homology) domain: A link between DNA methylation, replication and transcriptional regulation. FEBS Lett. 1999, 446, 189–193. [Google Scholar] [CrossRef]
  51. Du, Q.; Luu, P.-L.; Stirzaker, C.; Clark, S.J. Methyl-CpG-binding domain proteins: Readers of the epigenome. Epigenomics 2015, 7, 1051–1073. [Google Scholar] [CrossRef]
  52. Liu, G.; Xia, Y.; Liu, T.; Dai, S.; Hou, X. The DNA methylome and association of differentially methylated regions with differential gene expression during heat stress in Brassica rapa. Int. J. Mol. Sci. 2018, 19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Pei, L.; Zhang, L.; Li, J.; Shen, C.; Qiu, P.; Tu, L.; Zhang, X.; Wang, M. Tracing the origin and evolution history of methylation-related genes in plants. BMC Plant Biol. 2019, 19, 307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Moglia, A.; Gianoglio, S.; Acquadro, A.; Valentino, D.; Milani, A.M.; Lanteri, S.; Comino, C. Identification of DNA methyltransferases and demethylases in Solanum melongena L., and their transcription dynamics during fruit development and after salt and drought stresses. PLoS ONE 2019, 14, e0223581. [Google Scholar] [CrossRef] [Green Version]
  55. Gu, T.; Ren, S.; Wang, Y.; Han, Y.; Li, Y. Characterization of DNA methyltransferase and demethylase genes in Fragaria vesca. Mol. Genet. Genom. 2016, 291, 1333–1345. [Google Scholar] [CrossRef] [PubMed]
  56. Cao, X.; Jacobsen, S.E. Role of the Arabidopsis DRM methyltransferases in de novo DNA methylation and gene silencing. Curr. Biol. 2002, 12, 1138–1144. [Google Scholar] [CrossRef] [Green Version]
  57. Ganguly, D.R.; Crisp, P.A.; Eichten, S.R.; Pogson, B.J. The Arabidopsis DNA methylome is stable under transgenerational drought stress. Plant Physiol. 2017, 175, 1893–1912. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Steward, N.; Ito, M.; Yamaguchi, Y.; Koizumi, N.; Sano, H. Periodic DNA methylation in Maize nucleosomes and demethylation by environmental stress. J. Biol. Chem. 2002, 277, 37741–37746. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Zhu, C.; Zhang, S.; Zhou, C.; Chen, L.; Fu, H.; Li, X.; Lin, Y.; Lai, Z.; Guo, Y. Genome-wide investigation and transcriptional analysis of cytosine-5 DNA methyltransferase and DNA demethylase gene families in tea plant (Camellia sinensis) under abiotic stress and withering processing. PeerJ 2020, 8, e8432. [Google Scholar] [CrossRef] [Green Version]
  60. Kumar, R.; Chauhan, P.K.; Khurana, A. Identification and expression profiling of DNA methyltransferases during development and stress conditions in Solanaceae. Funct. Integr. Genom. 2016, 16, 513–528. [Google Scholar] [CrossRef]
Genes 13 00314 g001 550
Figure 1. Phylogenetic analysis of the C5-MTase (a) and MBD (b) genes. At, A. thaliana; Bd, B. distachyon; Ca, C. annuum; Cc, C. cardunculus; Fa, F. ananassa; Gm, G. max; Os, O. sativa; Ph, P. hybrida (blue); Pt, P. trichocarpa; Rc, R. communis; Sb, S. bicolor; Sl, S. lycopersicum; Sm, S. miltiorrhiza; Smel, S. melanogena; St, S. tuberosum; Ta, T. aestivum; Zm, Z. mays.
Figure 1. Phylogenetic analysis of the C5-MTase (a) and MBD (b) genes. At, A. thaliana; Bd, B. distachyon; Ca, C. annuum; Cc, C. cardunculus; Fa, F. ananassa; Gm, G. max; Os, O. sativa; Ph, P. hybrida (blue); Pt, P. trichocarpa; Rc, R. communis; Sb, S. bicolor; Sl, S. lycopersicum; Sm, S. miltiorrhiza; Smel, S. melanogena; St, S. tuberosum; Ta, T. aestivum; Zm, Z. mays.
Genes 13 00314 g001
Genes 13 00314 g002 550
Figure 2. Gene structure and motif analyses of petunia C5-MTase (a) and MBD genes (b). The full-length sequences of nine PhC5-MTases and 11 PhMBDs were used to construct the phylogenetic tree with MEGA 5 software. Lines represent introns, and filled yellow and green boxes represent UTRs and coding sequences, respectively. Different colored boxes indicate the motifs numbered 1–15.
Figure 2. Gene structure and motif analyses of petunia C5-MTase (a) and MBD genes (b). The full-length sequences of nine PhC5-MTases and 11 PhMBDs were used to construct the phylogenetic tree with MEGA 5 software. Lines represent introns, and filled yellow and green boxes represent UTRs and coding sequences, respectively. Different colored boxes indicate the motifs numbered 1–15.
Genes 13 00314 g002
Genes 13 00314 g003 550
Figure 3. Expression profiles of the identified PhC5-MTase and PhMBD genes, and the correlations between their expression patterns in anthers. (a) The expression patterns of PhC5-MTases; (b) the expression patterns of PhMBD genes; (c) a schematic diagram of the five stages of the anther developmental process; (d) the correlations between PhC5-MTase and PhMBD genes in anthers (RPKM > 10). The values in the boxes represent Pearson’s correlation coefficients (PCCs, p < 0.05), with red and blue indicating positive and negative correlations, respectively. The labels an1, an2, an3, an4 and an5 represent anthers from buds that reach 0.2, 0.3, 0.5, 1.2 and 3.5 cm in longitudinal length, respectively. The labels pe2, pe4 and ope represent the petals reach 0.3, 1.2 and 3.5 cm in longitudinal length, respectively. The labels pi2 and sep2 represent the pistils and sepals from buds reach 0.3 cm in longitudinal length, respectively.
Figure 3. Expression profiles of the identified PhC5-MTase and PhMBD genes, and the correlations between their expression patterns in anthers. (a) The expression patterns of PhC5-MTases; (b) the expression patterns of PhMBD genes; (c) a schematic diagram of the five stages of the anther developmental process; (d) the correlations between PhC5-MTase and PhMBD genes in anthers (RPKM > 10). The values in the boxes represent Pearson’s correlation coefficients (PCCs, p < 0.05), with red and blue indicating positive and negative correlations, respectively. The labels an1, an2, an3, an4 and an5 represent anthers from buds that reach 0.2, 0.3, 0.5, 1.2 and 3.5 cm in longitudinal length, respectively. The labels pe2, pe4 and ope represent the petals reach 0.3, 1.2 and 3.5 cm in longitudinal length, respectively. The labels pi2 and sep2 represent the pistils and sepals from buds reach 0.3 cm in longitudinal length, respectively.
Genes 13 00314 g003
Genes 13 00314 g004 550
Figure 4. The correlations between PhC5-MTase and PhMBD gene expression patterns in all the tissues.
Figure 4. The correlations between PhC5-MTase and PhMBD gene expression patterns in all the tissues.
Genes 13 00314 g004
Genes 13 00314 g005 550
Figure 5. Gene expression levels of six PhC5-MTase and four PhMBD genes under abiotic stresses (Cold and PEG). (a) Expression profiles of the selected PhC5-MTase and PhMBD genes under cold- and drought-stress conditions; (b) the correlations between the selected PhC5-MTase and PhMBD genes under cold- and drought-stress conditions. The values in the boxes represent Pearson’s correlation coefficients (PCCs, p < 0.05), with red and blue indicating positive and negative correlations, respectively. Each qRT-PCR analysis had three biological and technical replicates. Significant differences obtained by Duncan’s test (p < 0.05) are indicated by different letters above the error bars.
Figure 5. Gene expression levels of six PhC5-MTase and four PhMBD genes under abiotic stresses (Cold and PEG). (a) Expression profiles of the selected PhC5-MTase and PhMBD genes under cold- and drought-stress conditions; (b) the correlations between the selected PhC5-MTase and PhMBD genes under cold- and drought-stress conditions. The values in the boxes represent Pearson’s correlation coefficients (PCCs, p < 0.05), with red and blue indicating positive and negative correlations, respectively. Each qRT-PCR analysis had three biological and technical replicates. Significant differences obtained by Duncan’s test (p < 0.05) are indicated by different letters above the error bars.
Genes 13 00314 g005
Genes 13 00314 g006 550
Figure 6. Gene expression levels of six PhC5-MTase and four PhMBD genes under hormonal stresses (IAA and GA treatments) at the seedling stage. (a) Expression profiles of PhC5-MTase and PhMBD genes under IAA- and GA-treatment conditions; (b) the correlations between the selected PhC5-MTase and PhMBD genes under IAA- and GA-treatment conditions. The values in the boxes represent Pearson’s correlation coefficients (PCCs, p < 0.05), with red and blue indicating positive and negative correlations, respectively. Each qRT-PCR analysis had three biological and technical replicates. Significant differences obtained by Duncan’s test (p < 0.05) are indicated by different letters above the error bars.
Figure 6. Gene expression levels of six PhC5-MTase and four PhMBD genes under hormonal stresses (IAA and GA treatments) at the seedling stage. (a) Expression profiles of PhC5-MTase and PhMBD genes under IAA- and GA-treatment conditions; (b) the correlations between the selected PhC5-MTase and PhMBD genes under IAA- and GA-treatment conditions. The values in the boxes represent Pearson’s correlation coefficients (PCCs, p < 0.05), with red and blue indicating positive and negative correlations, respectively. Each qRT-PCR analysis had three biological and technical replicates. Significant differences obtained by Duncan’s test (p < 0.05) are indicated by different letters above the error bars.
Genes 13 00314 g006
Genes 13 00314 g007 550
Figure 7. Subcellular localization of PhDRM2 and PhMBD2. Tobacco leaves transformed with 35S::GFP were used as a control. The green represents the fluorescent signal of the green fluorescent protein (GFP), and blue represents the nucleus visualized by DAPI staining.
Figure 7. Subcellular localization of PhDRM2 and PhMBD2. Tobacco leaves transformed with 35S::GFP were used as a control. The green represents the fluorescent signal of the green fluorescent protein (GFP), and blue represents the nucleus visualized by DAPI staining.
Genes 13 00314 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shi, L.; Shen, H.; Liu, J.; Hu, H.; Tan, H.; Yang, X.; Wang, L.; Yue, Y. Exploration of the Potential Transcriptional Regulatory Mechanisms of DNA Methyltransferases and MBD Genes in Petunia Anther Development and Multi-Stress Responses. Genes 2022, 13, 314. https://doi.org/10.3390/genes13020314

AMA Style

Shi L, Shen H, Liu J, Hu H, Tan H, Yang X, Wang L, Yue Y. Exploration of the Potential Transcriptional Regulatory Mechanisms of DNA Methyltransferases and MBD Genes in Petunia Anther Development and Multi-Stress Responses. Genes. 2022; 13(2):314. https://doi.org/10.3390/genes13020314

Chicago/Turabian Style

Shi, Lisha, Huimin Shen, Jiawei Liu, Hongmin Hu, Hongyan Tan, Xiulian Yang, Lianggui Wang, and Yuanzheng Yue. 2022. "Exploration of the Potential Transcriptional Regulatory Mechanisms of DNA Methyltransferases and MBD Genes in Petunia Anther Development and Multi-Stress Responses" Genes 13, no. 2: 314. https://doi.org/10.3390/genes13020314

APA Style

Shi, L., Shen, H., Liu, J., Hu, H., Tan, H., Yang, X., Wang, L., & Yue, Y. (2022). Exploration of the Potential Transcriptional Regulatory Mechanisms of DNA Methyltransferases and MBD Genes in Petunia Anther Development and Multi-Stress Responses. Genes, 13(2), 314. https://doi.org/10.3390/genes13020314

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop