Open AccessArticle

Characterization, Evolution, Expression and Functional Divergence of the DMP Gene Family in Plants

Zeeshan Ahmad

^1,2

Dingyan Tian

^1,2,

Yan Li

^1,2,

Isah Mansur Aminu

^1,2,

Javaria Tabusam

^1,2,

Yongshan Zhang

^1,2,*

and

Shouhong Zhu

^1,2,*

National Key Laboratory of Cotton Bio-Breeding and Integrated Utilization, Institute of Cotton Research, Chinese Academy of Agricultural Sciences, Anyang 455000, China

Western Agricultural Research Center, Chinese Academy of Agricultural Sciences, Changji 831100, China

Authors to whom correspondence should be addressed.

Int. J. Mol. Sci. 2024, 25(19), 10435; https://doi.org/10.3390/ijms251910435

Submission received: 3 September 2024 / Revised: 22 September 2024 / Accepted: 24 September 2024 / Published: 27 September 2024

(This article belongs to the Section Molecular Plant Sciences)

Download

Browse Figures

Figure 1
Phylogenetic relationship of 240 DMP genes from 24 species, including Arabidopsis and cotton, showing a classification of DMPs in various groups. (A) MEGA 11 software was used to construct phylogenetic trees using the neighbor-joining (NJ) method with 1000 bootstrap replications. Each color shows a different group (I–V). Species are denoted with different colors and symbols in the phylogenetic tree. (B) Classification and arrangement of DMP genes in various groups based on a phylogenetic tree. Plant species classes are highlighted in various colors on the left. Five groups on the right side are made based on the phylogenetic tree containing DMP genes. "> Figure 1 Cont.
Phylogenetic relationship of 240 DMP genes from 24 species, including Arabidopsis and cotton, showing a classification of DMPs in various groups. (A) MEGA 11 software was used to construct phylogenetic trees using the neighbor-joining (NJ) method with 1000 bootstrap replications. Each color shows a different group (I–V). Species are denoted with different colors and symbols in the phylogenetic tree. (B) Classification and arrangement of DMP genes in various groups based on a phylogenetic tree. Plant species classes are highlighted in various colors on the left. Five groups on the right side are made based on the phylogenetic tree containing DMP genes. "> Figure 2
Comparison of the gene structure, conserved motifs and domains in DMP genes between Arabidopsis, cotton and 22 other species. (A) The NJ phylogenetic tree was constructed in MEGA 11 software using DMP sequences from 24 species. (B) The conserved protein motif composition of all DMPs and motif 1–10 is shown in various colored boxes. (C) Representation of the domain of Arabidopsis and other species. Yellow boxes represent the DMP domain in the sequences. (D) Exon and Intron structure of DMP genes of all Selected species. The red box represents exon or coding region, and the grey line represents intron with the scale at the bottom for measurement length of sequences. "> Figure 3
Promoter region analysis of DMP sequences from two monocot and two dicot species along with phylogenetic tree of the DMPs for each species separately. (A,B) = Dicots (Arabidopsis thaliana and Gossypium raimondi). (C,D) = Monocots (Sorghum bicolor and Zea mays). "> Figure 4
Evolutionary events in the DMP gene family in plants. Numbers of gene losses are shown in red color after “+” and duplications in blue color after “−”in the phylogenetic tree on the left side. "> Figure 5
Expression profiles of DMP genes in various organ tissue from two monocot and two dicot species. Raw RNA-seq data were procured from relevant databases. (A) Arabidopsis thaliana from TAIR database (TAIR, v.11, <a href="http://www.arabidopsis.org" target="_blank">http://www.arabidopsis.org</a>). (B) Avena sativa and (C) Glycine max from Phytozome (<a href="https://phytozome.jgi.doe.gov/pz/portal.html/" target="_blank">https://phytozome.jgi.doe.gov/pz/portal.html/</a>, accessed on 5 March 2024). (D) Gossypium hirsutum from COTTONGEN (<a href="http://www.cottongen.org/" target="_blank">http://www.cottongen.org/</a>, accessed on 7 March 2024) and the Cotton Functional Genomics Database (CottonFGD, <a href="https://cottonfgd.org/" target="_blank">https://cottonfgd.org/</a>). DMP gene expression levels are represented with various colors, with blue color indicating the least expression and red indicates the highest expression. Eudicot species (A,C,D); monocot species (B). "> Figure 6
Expression levels of DMP sequences in various tissues of cotton. GhActin was used as a reference. Error bars denote standard deviation measured from three experiments. ">

Versions Notes

Abstract

The DMP (DOMAIN OF UNKNOWN FUNCTION 679 membrane protein) domain, containing a family of membrane proteins specific to green plants, is involved in numerous biological functions including physiological processes, reproductive development and senescence in Arabidopsis, but their evolutionary relationship and biological function in most crops remains unknown. In this study, we scrutinized phylogenetic relationships, gene structure, conserved domains and motifs, promoter regions, gene loss/duplication events and expression patterns. Overall, 240 DMPs were identified and analyzed in 24 plant species selected from lower plants to angiosperms. Comprehensive evolutionary analysis revealed that these DMPs underwent purifying selection and could be divided into five groups (I–V). DMP gene structure showed that it may have undergone an intron loss event during evolution. The five DMP groups had the same domains, which were distinct from each other in terms of the number of DMPs; group III was the largest, closely followed by group V. The DMP promotor region with various cis-regulatory elements was predicted to have a potential role in development, hormone induction and abiotic stresses. Based on transcriptomic data, expression profiling revealed that DMPs were primarily expressed in reproductive organs and were moderately expressed in other tissues. Evolutionary analysis suggested that gene loss events occurred more frequently than gene duplication events among all groups. Overall, this genome-wide study elucidates the potential function of the DMP gene family in selected plant species, but further research is needed in many crops to validate their biological roles.

Keywords:

DMP genes; genome-wide; cotton; characterization; evolution; gene expression

1. Introduction

Membrane proteins are essential for many biological and physiological functions, such as material transport, signal transmission and recognition, cell division and proliferation, and stress responses [1]. For example, AtHKT1 proteins facilitate the transport of ions and metabolites across cell membranes, which is essential for maintaining cellular homeostasis and overall growth [2] and development by mediating hormone signaling pathways and environmental interactions [3]. Another protein, AtEC1, regulates sperm cell recognition with egg cells, and the AtGEX2 protein plays a vital role in gamete attachment for successful fertilization in Arabidopsis [4].

A class of unique membrane proteins present only in green plants are known as DMPs (Domain of unknown function 679 membrane proteins), a family that typically contains four to five transmembrane domains with cytoplasmic amino and carboxyl termini [5]. In Arabidopsis thaliana, a total of ten DMPs have been identified and are predicted to have vital roles in various physiological processes, particularly in senescence, reproductive development and various stress responses [6]. Studies have reported that these DMPs exhibit different expression patterns in various plant tissues, where AtDMP1 and AtDMP2 are involved in developmental senescence and darkness-induced senescence and engage in membrane fission and fusion during senescing and root vacuole biogenesis [7]. Similarly, AtDMP3 and AtDMP4 are upregulated in senescent leaves and siliques, demonstrating overlapping roles during senescence. Moreover, AtDMP4 is co-expressed with developmental programmed cell death genes such as RIBONUCLEASE3 (RNS3), BIFUNCTIONAL NUCLEASE1 (BFN1) [8] and EXITUS1 (EXI1), which promote senescence, differentiation-induced and age-induced developmental-programmed cell death [9]. AtDMP3 and AtDMP4 are expressed in flowers and roots, while AtDMP2, AtDMP6 and AtDMP7 are expressed in all plant organs [6]. Moreover, AtDMP8 and AtDMP9 are involved in gamete fusion during double fertilization in flowering plants. The expression levels of DMPs strongly suggests their involvement in senescence, abscission, vacuole biogenesis, reproductive development and stresses, but some DMP promotor regions do not have the cis elements for these activities, which indicates their potential involvement in other cellular activities [5].

DMPs have been studied, analyzed and characterized in some of the few crop species for their roles and functions. For example, loss of function mutations in AtDMPs and mutations in qhir8 of ZmDMPs led to the development of an effective haploid induction system in Arabidopsis thaliana and Zea mays, respectively [10]. Similarly, a haploid induction system was also established in in allotetraploid tobacco [11], polyploid Brassica napus [12] and tomato [13] using DMP mutants. AsDMP1 and AsDMP19 have been found to regulate seed anti-aging and seed longevity in oats (Avena sativa) [14], while Glyma.09G237500 and Glyma.18G098300 in soybean have been found to control female organ-mediated haploid induction [14]. Furthermore, a systematic study of the DMP family in four cotton species identified 58 DMPs, and their expression patterns unveiled their probable involvement in vital biological process like reproductive development, plant aging and stress responses [5]. However, few studies have been carried out on the investigation and characterization of the functions of DMPs in plant crops. In this regard, we conducted a study on 24 plant species selected from different classes (mosses to angiosperms) in which we identified 240 DMP genes. Furthermore, we used bioinformatics tools to comprehensively analyze the physiochemical properties, gene structures, conserved motifs, phylogenetic relationships and promoter regions of the DMPs, and have characterized their expression patterns in various crops.

2. Results

2.1. Identification of the DMP Proteins in Various Species

To identify the DMP gene family, BLASTP searches were performed across animal, plant, protist, fungi, and prokaryote protein databases using the ten known DMP amino acid sequences as query sequences from the model plant Arabidopsis (Table S7). According to BLAST searches for DMP-like proteins, no related amino acid motifs were discovered in mammals, fungi, prokaryotes, or archaea, whereas all higher plants, mosses, and Chlamydomonas reinhardtii contain DMP homologues. Overall, the BLAST results revealed a total of 1 DMP candidate gene in Chlamydomonas reinhardtii (Cre), 2 in Physcomitrella patens (Pp), 2 in Selaginella moellendorffii (Sm), 5 in Amborella trichopoda (Atr), 14 in Aquilegia coerulea (Aqco), 6 in Solanum lycopersicum (Sl), 7 in Vitis vinifera (Vv), 10 in Arabidopsis thaliana (At), 7 in Carica papaya (Cp), 11 in Gossypium raimondii (Gr), 8 in Theobroma cacao (Tc), 8 in Citrus sinensis (Csi), 12 in Populus trichocarpa (Ptr), 13 in Glycine max (Gm), 13 in Acorus americanus (Aa), 12 in Musa acuminata (Ma), 7 in Ananas comosus (Aco), 11 in Brachypodium distachyon (Bd), 19 in Oryza sativa (Os), 14 in Zea mays (Zm), 15 in Sorghum bicolor (Sb), 14 in Panicum hallii (Ph), 13 in Setaria italica (Si) and 15 in Setaria viridis (Sv).

2.2. Phylogenetic Analysis of the DMP Gene Family

To explore the evolutionary history of the DMP gene family across plant species, an unrooted phylogenetic tree was constructed in MEGA 11 using the neighbor-joining (NJ) method to analyze 240 orthologous and paralogous DMP protein sequences from the 24 species (Figure 1A). The phylogenetic tree of the retrieved protein sequences suggested that the plant DMP gene family clustered into five different groups, namely I, II, III, IV and V (Figure 1B). Most DMP genes (62) were found in Group III, closely followed by Group V (61), while Group II had the fewest DMPs (23). The most primordial green plant, Chlamydomonas reinhardtii, a single-celled green alga, has one DMP gene. The bryophytes (P. patens) and Pteridophyta (S. moellendorffii) each have two DMP homologue proteins. From a basal angiosperm, Amborella trichopoda, the DMPs began to diverge into five groups, and each group contained at least a single DMP homologue from there on (Figure 1B). The ten eudicot species in the study have several DMP homologues, ranging from 6 (Solanum lycopersicum) to 14 (Aquilegia coerulea), whereas the number of DMPs in ten monocot species ranges from seven (Ananas cosmos) to 19 (Oryza sativa) homologues. Interestingly, in eudicots, the normal frequency of DMP gene distribution was one to four, which was similar in most of the groups, except for group I, which had higher number of DMP genes in Aquilegia coerulea (10), suggesting occurrence of duplication of the genes in this species (Figure 1B). In contrast, no DMP was detected in Solanum lycopersicum, Vitis vinefera and Glycine max of group V, suggesting that deletion of the DMP genes may have occurred from these species. Likewise in monocots, the normal gene distribution frequency ranged from one to six, which was similar within the groups, except for group V, which had higher number of DMP genes in Oryza sativa (11), suggesting the occurrence of duplication of the genes in this species. On the other hand, group II of the monocots have normal gene distribution frequency 1, where no DMP genes were detected in Acorus americanus and Brachypodiun distachyon, suggesting loss of DMP genes from these species.

Interestingly, in eudicots, the number of DMP genes is almost the same in all groups, except group I have higher number of DMP genes in Aquilegia coerulea, suggesting duplication of DMP genes occurred in this species. On the other hand, in some species of group V, we did not find DMP genes, indicating deletion of the genes from this group. Likewise in monocots, group V deviates from the normal gene distribution pattern by presenting a higher number of DMPs, providing evidence of gene duplication in the species in this group, whereas no or lower distribution pattern of group II reveals deletion of DMP genes from some species in that group.

2.3. Gene Structure, Motifs and Domain Analysis of DMPs

Gene function is related to motifs, domains and structure of the gene. The protein sequences of DMPs were used to analyze the gene structure, motifs and conserved domains of all 24 species (Figure 2). Among all species, most of the DMPs have no intron structure and have an exon region for a few genes in various groups. Surprisingly, the DMP genes have a striking similar intron/exon structure within the same group compared to other groups, i.e., four genes (CsDMP7, CpDMP1, AqDmp1, AaDMP13) in Group II harbors two exons and one intron, each of various lengths (Figure 2D). CsDMP7 exhibits a very large intron region of more than 90% which is unique as compared to other DMP genes, suggesting it might have special biological functions. Group III shows the same pattern, with two exons and one intron gene structure observed in VvDMP2, BdDMP6 and AaDMP2. Similarly, in Group IV, ZmDMP5, AtrDMP1 and CsDMP5 each harbor two exons and one intron. In contrast, Group I have TcDMP5, AqDMP14 and VvDMP1, each comprising two introns and one exon of highly conserved lengths, in comparison with the gene structure of CpDMP4, which contains two very small exons and two very large exons in the same group. Likewise, Group V has a total of seven genes having both exon-/intron-like gene structure. Among them, the gene structure of MaDMP7 and SBDMP4 has three exons and two introns of conserved length, while AaDMP3, AaDMP4, SbDMP5, SiDMP10 and PtDMP9 bear two exons and one intron.

The DMP sequences were run through MEME Suite (v5.5.7) motif analysis software, and the results revealed that motifs of the DMP family were distinct from each other; however, motifs were conserved within the groups (Figure 2B). All the genes were found to have a single DMP domain. In total, 10 motifs were identified, whereby each gene had three to ten motifs. Motifs 1 and 2 were commonly present in all groups, and members within groups shared similar motifs, as expected, indicated by the conserved protein architecture within the groups which suggests function similarity. However, functions of these conserved motifs are yet to be elucidated. In contrast, some groups lacking some of the specific motifs such as motif 10 was not present in all members of Groups I and V and most members of Groups II and IV. Motif 9 was missing from most of the members of Groups II and V. However, further research on whether specific motifs have unique functional roles in the plants is needed. These results fully support our phylogenic analysis, and such diversity in motif composition within subfamilies might be explained by their functional diversification.

The protein sequences of the 24 selected species from different families were subjected to domain analysis to better understand the divergence among five groups of DMPs. Results of logo sequences revealed that all sequences of DMP’s share only one typical conserved DUF679/DMP domain (Figure 2C) and having two to five trans-membrane domains like Arabidopsis (Table S6; Figures S2–S6). All trans-membrane domains were within the DUF/DMP domain, and the DMP domain locus within the same group was highly conserved, i.e., the DMP proteins in Groups I and II are most alike (Figure 2C, Figures S2 and S3). These results suggest that all DMP sequences among the species with the same DUF/DMP domain may have the same functions in all the species. However, identification of the specific functions of these genes requires further research.

2.4. Promoter Region Analysis in Dicots and Monocots

Cis elements in the upstream promoter sequence often influence gene expression. Cis elements in the non-coding DNA upstream of the gene transcription start site control gene expression in various settings and may be responsive to stress receptors or differentially regulated in tissues. Analyzing cis elements involved in DMP gene regulation can assist in understanding the gene’s regulatory mechanism and its probable functions. To identify the Cis-acting elements in the promoter region of DMPs, a 5′ terminal 2000 bp upstream region of the start codon was selected and analyzed using the PLANTCARE database [15]. Results disclosed that the promoter regions of DMPs of both monocot (Sorghum bicolor and Zea mays) and dicot (Arabidopsis thaliana and Gossypium raimondii) species had similar cis-acting elements, including responsive elements, and development, phytohormones and stress response related elements from 59 DMP genes (Figure 3).

Among the cis elements involved in growth and development of plants, the light-responsive elements accounted for more than 50% and were distributed along the promoter regions in both monocots and dicots (Table S3). These elements include Sp1, L-box, LS7, AE-box, I-box, G-box, Gap box, Box II, Box IV, GA motif, ATC motif, AAAC motif, GTI-motif and TCT motif, which are present in the promoter regions of AtDMP5, GrDMP4, GrDMP7, GrDMP9, ATDMP1, AtDMP2 and AtDMP3 (Table S8a). Cis elements responsible for metabolism regulation, plant growth and development include circadian, O₂ cite and HD-Zip1 (present in the promoter region of AtDMp8, AtDMP9, AtDMP10, GrDMP1, GrDMP9, ZmDMP6, ZmDMP7, SbDMP1 and SbDMP8), meristem expression-related CAT box (in the promoter region of AtDMP2, AtDMP9, GrDmp5, GrDMP6, GrDMP8, ZmDMP3, ZmDMP4, SbDMP3 and SbDMP7), endospermic expression-related GCN4 (in the promoter region of ATDMP4, GrDMP9 and ZmDMP9), and seed-specific regulation-related RY elements (in the promoter region of ATDMP10, GrDMP3 and SbDMP10) (Table S8a–c). The presence of GCN4 and RY-elements in the above genes suggests that DMPs may have vital roles in plant reproductive development.

Phytohormone-related elements, including TGA box, TGA element and AuxRR-core (in promoters of AtDMP1, GrDMP10, SbDMP2 and ZmDMP3) involved in auxin responsiveness, TCA element responsible for salicylic acid responsiveness (Table S3), GARE-motif and TATC-box (in promoters of AtDMP7, GrDMP5, SbDMP7 and ZmDMP10) involved in gibberellin responsiveness, and CGTCA- and TGACG-motifs involved in methyl jasmonic acid responsiveness signaling, were identified in the promoter regions of DMPs (Table S8a–d). Most of the cis-acting elements such as ABRE and AAGAA that are responsible for abscisic acid responsiveness were present in promoter regions of AtDMP5, AtDMP8, GrDMP1 and GrDMP2. However, most important cytokinin related cis elements were not found in any of the promoter regions.

In the promoter region of DMPs, we also identified cis elements related to stress responses; in particular, elements related to drought stress were enriched, such as dehydration-related MYC, Myb, Myc and Myb binding cites (GrDMP2, SbDMP15 and ZmDMP9), MBS involved in drought inducibility, TC-rich repeats for defense and stress responsiveness, and drought, salt and low temperature stress-responsive DRE cores (Figure 3). Cis-acting elements such as LTR and STRE (in promoter regions of AtDMP4, AtDMP8, ZmDMP1, ZmDMP2, SbDMP2, SbDMP3) that are involved in low temperature responsiveness were also identified (Table S3). Other cis elements like ARE (in promoter regions of AtDMP6, GrDMP11, SbDMP14 and ZmDMP2) for anaerobic induction, WUN-motif, WRE3 and W-box (in promoters of AtDMP9 and GrDMP6) involved in wound responsiveness were not present in monocot species (Table S8c,d).

2.5. Duplication and Loss of DMP Genes in Plants during Evolution

We performed gene duplication and loss analysis in plants using NOTUNG (v2.9.1.5) software to better understand the evolutionary events that occurred among these DMP genes [16]. Getting DMP genes from phylogenetic analyses at various phases of evolution, we further determined if internal nodes within each group were related to gene duplication, gene loss or lineage divergence events (Figure 4).

In Chlamydomonas reinhardtii, no gene duplication and gene loss event were detected. Evolving from chlorophytes to mosses, one gene duplication event occurred; however, no gene loss event was found. During the emergence of lycophytes, one gene duplication and one gene loss event took place. Selaginella moellendorffii in this taxonomic class also experienced one gene duplication and gene loss event. However, in the lower angiosperms, 19 genes were duplicated, and no gene was lost in the event. From Amborella trichopoda, being a lower angiosperm, the DMPs began to diverge into five groups, and each group contained at least a single DMP homologue within the group, showing more evolution, gene duplication and gene loss events (Figure S1). However, this important species lost 16 DMP genes during the evolution process. In eudicots, including ten important species, four gene duplication events occurred while no gene was lost. Arabidopsis, an important model species of this class, went through three gene duplication and six gene lost events during the evolution process. In monocots, containing ten species, three genes were duplicated while 16 genes were lost from this class. Oryza sativa, an important model species in the group, experienced five gene duplications and ten gene losses. However, Setaria virdis displayed no gene duplication or gene loss event.

2.6. Tissue Expression Levels of DMPs in Monocots and Eudicots

We performed the expression analysis of DMPs in various tissues across two monocot and two eudicot species using the transcriptomic data from various databases (Table S4). The results showed that different DMPs expressed in most of the organ tissues of plants, particularly in the reproductive organs as well as various growth and developmental stages of plants, suggests that these genes may have different functions and roles in various crop species (Figure 5).

DMP genes of Group I were highly expressed in roots, leaves, flowers and petals in monocot species, while in eudicots, expression of DMPs was limited to just roots and seeds. In Group II, the monocot species (Avena sativa) AsDMPs such as AsDMP12, AsDMP11, AsDMP17 and AsDMP26 are greatly expressed just in roots (Figure 5B). However, eudicots followed a different pattern, where AtDMPs (Arabidopsis thaliana) such as AtDMP1 and AtDMP2 are highly expressed in roots, GmDMPs (Glycine max) such as GmDMP8 and GmDMP5 are highly expressed in flowers, and GrDMs (Gossypium raimondii) are highly expressed in petals, pistils and anthers (Figure 5A,C,D). The Group III monocot AsDMP1, AsDMP2 and AsDMP3 were rarely expressed in roots, whereas in dicots, AtDMP4 and AtDMP6 of Arabidopsis and GrDMP4 and GrDMP6 of cotton were highly expressed in flowers and reproductive organs like sepals, petals and seeds. However, GmDMPs were not expressed in any of the tissues (Figure 5C). In Group IV, the monocot AsDMP10 and AsDMP22 were highly expressed in roots and leaves. In contrast, in dicot species, AtDMP9, GmDMP1 and GrDMP9A were greatly expressed in flower, sepal and petal tissues. The other gene in this group was expressed in various tissues. For Group V, AsDMP11, AsDMP12 and AsDMP24 of the monocot species Avena sativa were predominantly expressed in roots. However, dicot DMPs like AtDMP10, GrDMP10A and GrDMP10B were expressed in the reproductive tissues as well as in vegetative tissues of the plants.

2.7. Tissues Expression of DMPs in Cotton

We selected five DMP genes (GhDMP1A, GhDMP3A, GhDMP6A, GhDMP8A and GhDMP10A) based on RNA-Seq data to confirm the expression profiles of DMP genes in various tissues of Gossypium hirsutum by qRT-PCR. The gene expression patterns from qRT-PCR analysis showed a similar trend to the pattern detected using the RNA-seq data, as shown in Figure 6. GhDMP1A and GhDMP3A dominantly expressed in petals with the least expression in vegetative parts, whereas GhDMP6A, GhDMP8A and GhDMP10A expressed in reproductive as well as vegetative tissues of the plant (Table S5). Based on these results, DMP genes may have been involved in reproductive functions rather than playing a role in vegetative growth and development of plants.

3. Discussion

The DMP sequence resembled known membrane proteins, transporters or channels in any kingdom [5], suggesting functions in processes that also occur in plants. These protein sequences are limited to green plants solely, evolved and expanded particularly in flowering plant lineages [5,17]. In this study, we used various bioinformatics tools and publicly available transcriptomic datasets to mine genomes of different crops for DUF679 family genes. Various numbers of DMP genes were identified in different crop species from distinct plant classes. Polyploid and complex genomes, and gene duplication and loss events during the evolutionary process in plants may increase or decrease the number of DMPs. To construct a phylogenetic tree, we analyzed Arabidopsis DMP sequences along with 23 other species from various classes of plants and found that all species’ DMP members could be clustered into five groups, the same as in Arabidopsis (Figure 1A). Groupwise, the gene frequency distribution pattern was the same within groups, with some exceptions; for example, in dicot Group I, all members had almost equal numbers of genes except Aquilegia coerulea, which had eight times more genes than other species. Likewise, in monocot Group V, the number of genes in Oryza sativa was double that of other species in the group, which might be due to duplication of genes in their larger genome over time in comparison to smaller genomes (Figure 1B). Structure analysis of DMP sequences revealed that most protein sequences in all species had no introns, with some exceptions having single or double intron/exon structure (Figure 2D), indicating potential special functions of the single or double intron/exon genes in various plants. These results support previous studies on DMPs conducted in Arabidopsis and cotton [6,7]. Similar to our studies, previous findings on genetic sequencing also suggested that initially, eukaryotic predecessors had intron-rich genes individually, and most of them experienced losses of introns during the process of evolution [18]. Motif analysis results disclosed certain differences in conserved motifs, suggesting that DMPs take part in various biological functions of growth and development. Three motifs, namely, motif1, motif2 and motif3, were conserved in all sequences of DMPs particularly with motif1 and motif2 of the DUF679 domain among the species, while other were restricted to specific genes (Figure 2C). Moreover, the motif pattern was the same within the same group while various groups showed different combination of motifs. However, most of the amino acids of DMPs contain a single DMP domain with four to five transmembrane domains (Table S6; Figures S2–S6), which indicates that DMPs have a common ancestor and evolved independently and may have the same special functions. Some earlier studies proposed that a single amino acid mutation in the initial transmembrane of DMP gene in qhir8 in maize can induce haploids [19,20]. Furthermore, based on transcriptomic data, the expression pattern of DMP genes among four monocot and dicot species was almost similar (Table S4). Most of the DMP genes widely expressed in reproductive organs, some in vegetative parts, showing that DMPs may play a reproductive role rather than playing a part in vegetative growth in general [5,14]. Consistent with this, the qRT PCR analysis showed that each representative DMP gene from the group was expressed widely in floral organs and to some extent in vegetative tissues (Figure 6). These results signify previous findings that DMP may have roles in fertility and reproductivity [5]. Looking briefly at the functions of DMPs, DMP4 promotes senescence and interacts with DMP1 [21], AtDMP2 and AtDMP3 involved in programmed cell death [14,22]. Also, other studies reported that mutation in DMP8 and DMP9 in Arabidopsis resulted in decreased fertility and seed setting rates [22,23]. Equally, a recent study in soybean also demonstrated the probable involvement of DMP8 and DMP9 in pollination and fertilization. Previously, it was determined that DMP1 is involved in membrane fission and fusion during leaf senescence and vacuole biogenesis in roots [6,24]. Additionally, research in oat has found the involvement of DMP1 in anti-aging and seed longevity [25]. It can be concluded that DMPs play a vital role in development and fertility of monocot and dicot species. However, the functions of all the DMP genes are still unknown in many crop species, and further research is needed to investigate these thoroughly.

Looking into gene loss and duplication events, primordial plants contain a maximum of two DMPs per species, with no gene loss or duplication events recorded. DMP genes began to diverge in the basal angiosperms by having at least five copies of DMPs in a species (Figure 4). Thereafter, the greatest number of genes was duplicated in monocot species, while high numbers of DMP genes were frequently lost from dicot species, signifying that rapid birth and death events occurred in DMPs during that time (Figure S1). A single copy of genes in any genome is enough for essential viability, and loss or duplications leads to change of functions. For example, Cter-SUN proteins in Caenorhabditis elegans and Schizosaccharomyces pombe were responsible for nuclear migration and duplication of SBPs (spindle pole bodies) [26,27] However, when Cter-SUN proteins duplicated into SUN1 and SUN2, they started to have vital roles in chromosome movement [28]. These results suggest that DMP functions may have changed over time due to frequent loss or duplication events.

We studied the distribution and frequency of cis-acting regulatory elements (CAREs) in four crop species (monocots and dicots) to explain the possible cis-regulatory roles of DMPs in growth and stress responses. Interestingly, all the studied species from monocots and dicots had almost the same cis elements and promoter regions with minor exceptions (Figure 3; Table S8a–d). Regarding growth, light responsive elements (LREs) were found to be predominant in DMP promoter regions. G-Box, Box4, GTI-motif, LS7 and I-box were the most common cis elements and are considered critical for regulating light-mediated transcriptional activity [29,30]. The most common G-box (CACGTG) is reported to be involved in regulation of chlorophyll biosynthesis in Arabidopsis [31], while the GATA motif has a vital role in tissue specificity and development of tissues [32]. G-box and I-box present in the promoter regions of DMPs are considered to be photosynthetic-responsive elements, and combined, both are essential for activation of phytochrome, cryptochrome and plastid signals [33,34]. Interestingly, Arabidopsis, maize, oat and G. raimondii had a site for an sp1 (GGGCGG) element (Table S8a–d), a mammalian promoter element associated with the regulation of a large number of housekeeping and tissue-specific genes [35]. Previous studies concluded that it was absent from the rice genome, lacking significant functions [36]. Relating growth and development of the plant O₂-site, GCN4_motif and activation sequence 1 (TGACG) were frequently found in all of the DMP sequences (Table S3; Table S8a,c). Activation sequence 1 element has been found to play a role in auxin- or salicylic acid-dependent enhanced expression in leaves [37]. Furthermore, we identified seed-specific regulation related to the RY element in promoter regions of numerous DMPs (Figure 3; Table S8a–d). This element interacts and binds with ABI3 protein and FUS3 protein during development of mature seeds in Arabidopsis [38].

Phytohormones, i.e., plant hormones, are chemical substances that regulate diverse biological processes in plants, involving growth, development and responses to environmental stresses [39]. Key phytohormones controlling regulatory elements including methyl jasmonic acid responsive elements, abscisic acid responsive, auxins and gibberellin responsive elements were present in promoters of the studied DMPs. Among them, (abscisic acid responsive elements (ABREs) have higher binding sites in most of the DMP promoter regions (Figure 3; Table S8a–d), acting as regulators of biological processes like stomatal closure, seed and bud dormancy, and plant responses to drought, cold, heat and saline stresses [40]. Methyl jasmonic acid responsive elements (CGTCA-motif and TGACG-motif) were found in several promoter regions (Table S3). They are considered vital in the regulation of TF-mediated gene [41]. For instance, the TGACG motif from the rice promoter region regulates 12-oxophytodienoic acid reductase-1 and plays a crucial role in defense responses [42]. Moreover, some promoter regions have auxin responsive elements (Figure 3; Table S3), which induces the expression of WOX11 and its homologue WOX12 gene, thus stimulating the cell fate transition [43].

DMP promoter regions also contained cis-regulatory elements accountable for modulating gene expression in response to various stresses. For stresses, DRE cores, MBS, MYC, MBSI, Myb, Myc and Myb-binding sites were identified in most promoter regions (Figure 3; Table S3). Previous reports suggested they play inevitable roles in drought-inducible expression, indicating that DMP expression is also related to abiotic stress [44]. Additionally, the presence of low temperature-responsive (LTR) elements in some of the DMP sequence promoters shows its involvement in cold stress [45]. Some of the sequences contained promoter regions for anaerobic-responsive elements (AREs), which are known for having a role in response to low oxygen and dehydration [46]. The WUN-motif and W-box (TTGACC) were present in promoter regions of AtDMP9 and GrDMP6 (Figure 3; Table S8c,d), and is believed to help plants in healing during stress conditions [47]. Furthermore, the W-box element interacts with WRKY transcription factors to regulate the expression of defense related genes and helps plant seeds during dormancy and senescence [48]. The presence of diverse stress-responsive cis-regulator elements and varied expression patterns of DMPs using public data indicates their potential roles in improving stress tolerance and survival (Figure 3). However, this research only preliminarily analyzed phylogenetic relations and characterized basic functions of DMPs, and further research is needed on the functional validation of genes to understand their relevance during stress and various biological processes in various monocot and dicot crop species.

4. Materials and Methods

4.1. Identification and Retrieval of DMP Gene Family

The published protein sequences of Arabidopsis thaliana (AtDMPs) were retrieved from the Arabidopsis Information Resource database (TAIR, v.11, http://www.arabidopsis.org/, accessed on 2 March 2024), while the genome sequences of Chlamydomonas reinhardtii (v5.6), Physcomitrella patens (v3.3), Selaginella moellendorffii (v1.0), Amborella trichopoda (v1.0), Aquilegia coerulea (v3.1), Solanum lycopersicum (ITAG4.0), Vitis vinifera (v2.1), Carica papaya (v0.4), Theobroma cacao (v2.1), Citrus sinensis (v1.1), Populus trichorpa (v3.1), Glycine max (Wm82.a4.v1), Acorus americanus (v1.1), Musa acuminata (v1.0), Ananas comosus v3.0, Brachypodium distachyon v3.1, Oryza sativa (v7.0), Zea mays (RefGen_V4), Sorghum bicolor (v3.1.1), Panicum hallii (v3.1), Setaria italica (v2.2) and Setaria viridis (v2.0) were downloaded from phytozome (https://phytozome.jgi.doe.gov/pz/portal.html/, accessed on 5 March 2024). The genome sequences of Gossypium raimondii (CRIv2.1) were acquired from the COTTONGEN (http://www.cottongen.org/, accessed on 7 March 2024) [49] and reconfirmed in the Cotton Functional Genomics Database (CottonFGD) (https://cottonfgd.org/, accessed on 11 March 2024) [50].

Arabidopsis protein sequences were used as query sequences (Table S7), and the BLASTP program was used to search the candidate DMP sequences in the above species protein databases. Subsequently, the acquired protein sequences were verified for the DMP domain (IPR007770) in the candidate sequences using SMART (http://smart.embl.de/, accessed on 18 March 2024) [51] and Interproscan 05 (https://www.ebi.ac.uk/jdispatcher/pfa/iprscan5, accessed on 22 March 2024) [52], and eventually the DMP sequences were identified. Later, Deep TMHMM server v.2.0 (https://services.healthtech.dtu.dk/services/TMHMM-2.0/, accessed on 25 March 2024) was used to predict the number of transmembrane domains in the targeted sequences [53]. This information, along with gene numbers, is presented in Supplementary Table S1.

4.2. Amino-Acid Sequence Alignment and Phylogenetic Analysis

The full-length protein sequences of DMPs were multiple aligned using DNAMan v2.0 for multiple sequence logo analysis, the transmembrane domain sequences of DMPs from all the selected species were aligned, and the multiple alignment result was used for generating sequence logos through an online tool, WEBLOGO (https://weblogo.berkeley.edu/logo.cgi, accessed on 27 March 2024) [54]. The extracted DMP sequences from 24 plant species were aligned using MEGA Suite software (v.11.0) followed by a neighbor-joining (NJ) tree method to construct a phylogenetic tree [55]. The phylogenetic analysis was performed using bootstrap method with 1000 replicates and the substitution was then evaluated by a Poisson model with default parameters.

4.3. Conserved Motifs, Domain and Gene Structure Analysis

The conserved motifs of DMP sequences were predicted using Multiple Expectation Maximization for Motif Elicitation (MEME suite v.5.5.5) online website (https://meme-suite.org/meme/info/status?service=MEME&id=appMEME_5.5.51716255002353-1793066522, accessed on 29 May 2024) for motif identification [56]. For motif identification, the p-value was kept lower than 10⁻⁵, mode was set to “classic mode”, the site distributions were “zero or once per sequence”, the number of motifs was set to “10”, and the width of motifs was set to “between 6 and 50”. The DMP gene exon/intron structures, conserved motifs and domains along with phylogenetic tree was drawn and visualized with TBtools software (v.2.102) [57], and the Gene Structure Display Server (GSDS v2.0, https://gsds.gao-lab.org/index.php, accessed on 9 June 2024) by using NWK file from the phylogenetic tree, MAST file from the MEME website, CDS and genome files of all selected species [17].

4.4. Promoter Analysis of Monocot and Eudicot Species

Among the studied species, two monocot (S. bicolor and Z. mays) and two dicot (A. thaliana and G. raimondii) species were selected, and their promoter regions were analyzed to predict the cis-acting elements and compared with each other to find any differences among them. The DNA sequences of the 2000 bp upstream region of the start codon for DMPs were retrieved from Phytozome v.13 (https://phytozome-next.jgi.doe.gov/, accessed on 7 April 2024) website for selected species. The obtained sequences were then submitted to the PlantCARE website (https://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 17 June 2024) to predict and analyze the cis-regulator elements related to phytohormones, plant growth and development, stresses and light responses in the promoter region of DMPs [15]. TBtool (v2.119) software was used to draw and visualize the graphics for cis-acting elements.

4.5. Gene Loss and Duplication Event Analysis

The plant species tree was constructed using an online website (Time Tree v.5.0, https://www.timetree.org/, accessed on 3 April 2024) [58] in NWK format, while the gene tree was adapted from MEGA (v.11) software [59] in NWK file format using protein sequences. Sequences were first aligned in MEGA and then a NWK tree was constructed for the protein sequences. Similarly, a separate NWK species tree was constructed for all species from which the sequences were extracted. The gene tree was then reconciled with the plant species tree by Notung-DM software (v.2.9.1.5) with default parameters to obtain the trees’ gene loss and duplication events [16].

4.6. Expression Profiles of DMPs in Monocots and Eudicots

RNA-seq data were downloaded from the TAIR database (v.11, http://www.arabidopsis.org, 10 May 2024) for Arabidopsis thaliana, from the Cotton Omics database (http://cotton.zju.edu.cn/10.rnasearch.html, accessed on 13 May 2024) for cotton [60] and from the Ensemble Plants (v.59) database (https://plants.ensembl.org/, accessed on 16 May 2024) for Avena sativa and Glycine max. Individual experimental design and sampling methods used for RNA extraction in different monocot and dicot species are explained below:

For Arabidopsis, total RNA was isolated from various tissues such as young leaves (20, 34 and 50 days after sowing), green and ageing cauline leaves, siliques, root tips, plant callus, sepals, petals, young and old stems. For young stems, the upper main part was used, whereas old stems included only the brown part below the first node. Roots and root tips samples were collected from lower parts of plants covered in soil [6].

Oat seeds of the Baiyan 2 variety were utilized in the experiment. They were naturally stored for 0, 1, 2 and 3 years at an average temperature of 10 °C. The harvested seeds were artificially aged by lowering their moisture content to 10–14% and storing them for 24, 48, 72 and 96 h at 45 °C and 95% humidity. Afterwards, the seeds were dried and kept at 4 °C. Unaged seeds were used as control treatment [25].

Soybean seeds of the Williams 82 variety were cultivated at 25 °C with 16 h of light and 8 h of darkness. After 60 days, RNA was collected from leaves, stems and flowers. In total, three fresh leaves were cut off from the top, and stems were cut 3 cm below the top part of the plant. Plant floral parts such as petals, sepals, pistils and pollen were taken from 100 open flowers on a total of four plants [14].

The cotton (Gossipium hirsutum) cultivar Zhongmiansuo 100 was used to extract total RNA from stems, roots, sepals, leaves, petals, bracts, stigma, anthers (flower buds of different sizes: <3 mm, 3–5 mm, 5–8 mm and >8 mm), pollen and fiber (6, 12, 18 and 24 days post anthesis). The RN38-EASYspin-Plus Plant RNA kit was used for RNA extraction [5].

Different isolation kits were used for DNA extraction in the above experiments, followed by agarose gel electrophoresis to assess the quality and Nanodrop or Nanopro spectrophotometry to determine the concentration and purity of the samples. Thereafter, RNA was reverse transcribed with the prime Script^TM RT reagent kit and gDNA Eraser. Specific primers were used in every experiment according to their materials. PCR conditions were consistent across all tests, beginning with denaturation at 95 °C for 15 min, followed by amplification (40 cycles) at 95 °C for 10 s and annealing at 58 °C for 30 s. The melt curve analysis was set to default, which includes 1 s at 95 °C, 15 s at 65 °C, and 1 s at 95 °C. All experiments were replicated thrice, and the 2^−ΔΔCt protocol was applied to obtain the expression levels of genes. Heat map charts along with phylogenetic trees were then generated through TBtool (v2.119) software according to gene expression values (fragments per kilobase million, FPKM).

4.7. RNA Extraction and qRT-PCR Analysis

Total RNA was extracted from various tissues of the cotton (Gossypium hirsutum) cultivar, including root, leaf, stem, and petal, using an RN38-EASYspin-Plus Plant RNA Kit (Aid lab Co., Ltd., Beijing, China). A 60-day-old plant was selected to take samples of newly formed leaves from the top and a stem sample was taken 3 cm below the shoot apex. Five newly opened flowers were selected from each plant (with a total of 10 plants), and fresh roots along with tips were sampled. Samples were taken from the cotton cultivars and immediately stored in liquid nitrogen for subsequent RNA extraction. The pure RNA was then reverse transcribed to cDNA using a Prime Script™ RT reagent kit (Takara Biomedical Technology Co., Ltd., Beijing, China) following the manufacturer’s protocols. Bio-Rad 7500 fast fluorescence quantitative PCR system with SYBR^® Premix Ex Taq™ (Tli RNaseH Plus, Takara Biomedical Technology Co., Ltd., Beijing, China) was used for qRT-PCR amplifications according to the manufacturer’s instructions. The PCR machine was set for primary denaturation at 95 °C for 30 s followed by 40 cycles of amplifications at 95 °C for 5 s and at 65 °C for 30 s. The The qRT-PCR experiments compromised biological replicates and three technical replicates with the gene GhActin serving as the internal reference, whereas the relative gene expression levels of GhDMPs were calculate by the 2^−ΔΔCt method [61]. As DMP genes are divided into five groups, an individual representative gene was selected from each group and then their expression levels were compared. The specific primers of selected genes for qRT-PCR are listed in Supplementary Table S2.

5. Conclusions

This study provides a comprehensive genome-wide analysis of the DMP gene family across 24 plant species, highlighting its evolutionary history, gene structure and potential biological functions. Our findings suggest that DMP genes have undergone purifying selection and significant evolutionary changes, including intron loss and gene loss events, contributing to their current diversity. The conserved domain structure across five distinct groups indicates functional conservation, while promoter analysis and expression profiling suggest a key role for DMP genes in reproductive development, hormone response and abiotic stress adaptation. Despite these insights, further functional studies are required, particularly in major crop species, to fully understand the biological roles of DMP genes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms251910435/s1.

Author Contributions

Conceptualization, Z.A. and Y.Z.; data curation, Z.A. and D.T.; formal analysis, Z.A.; funding acquisition, Y.Z.; investigation, D.T.; methodology, Z.A.; project administration, S.Z.; resources, Y.L.; software Z.A.; supervision, Y.Z. and S.Z.; validation, Z.A.; writing—original draft, Z.A.; writing—review & editing, Y.Z., S.Z., I.M.A. and J.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Biological Breeding-Major Projects in National Science and Technology 2023ZD04038.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article or Supplementary Materials.

Acknowledgments

We are extremely thankful to the administrative and technical support staff of the Cotton research Institute and our lab students for their support throughout the experiment.

Conflicts of Interest

The authors declare no conflicts of interests.

References

Xicluna, J.; Lacombe, B.; Dreyer, I.; Alcon, C.; Jeanguenin, L.; Sentenac, H.; Thibaud, J.-B.; Chérel, I. Increased functional diversity of plant K+ channels by preferential heteromerization of the shaker-like subunits AKT2 and KAT2. J. Biol. Chem. 2007, 282, 486–494. [Google Scholar] [CrossRef] [PubMed]
Mäser, P.; Thomine, S.; Schroeder, J.I.; Ward, J.M.; Hirschi, K.; Sze, H.; Talke, I.N.; Amtmann, A.; Maathuis, F.J.; Sanders, D. Phylogenetic relationships within cation transporter families of Arabidopsis. Plant Physiol. 2001, 126, 1646–1667. [Google Scholar] [CrossRef] [PubMed]
Kang JooHyun, K.J.; Hwang JaeUng, H.J.; Lee MiYoung, L.M.; Kim YuYoung, K.Y.; Assmann, S.; Martinoia, E.; Lee YoungSook, L.Y. PDR-type ABC transporter mediates cellular uptake of the phytohormone abscisic acid. Proc. Natl. Acad. Sci. USA 2010, 107, 2355–2360. [Google Scholar] [CrossRef] [PubMed]
Mori, T.; Igawa, T.; Tamiya, G.; Miyagishima, S.-Y.; Berger, F. Gamete attachment requires GEX2 for successful fertilization in Arabidopsis. Curr. Biol. 2014, 24, 170–175. [Google Scholar] [CrossRef]
Zhu, S.; Wang, X.; Chen, W.; Yao, J.; Li, Y.; Fang, S.; Lv, Y.; Li, X.; Pan, J.; Liu, C. Cotton DMP gene family: Characterization, evolution, and expression profiles during development and stress. Int. J. Biol. Macromol. 2021, 183, 1257–1269. [Google Scholar] [CrossRef]
Kasaras, A.; Kunze, R. Expression, localisation and phylogeny of a novel family of plant-specific membrane proteins. Plant Biol. 2010, 12, 140–152. [Google Scholar] [CrossRef]
Kasaras, A.; Melzer, M.; Kunze, R. Arabidopsis senescence-associated protein DMP1 is involved in membrane remodeling of the ER and tonoplast. BMC Plant Biol. 2012, 12, 54. [Google Scholar] [CrossRef]
Olvera-Carrillo, Y.; Van Bel, M.; Van Hautegem, T.; Fendrych, M.; Huysmans, M.; Simaskova, M.; van Durme, M.; Buscaill, P.; Rivas, S.; Coll, N.S.; et al. A Conserved Core of Programmed Cell Death Indicator Genes Discriminates Developmentally and Environmentally Induced Programmed Cell Death in Plants. Plant Physiol. 2015, 169, 2684–2699. [Google Scholar] [CrossRef]
Gao, Z.; Daneva, A.; Salanenka, Y.; Van Durme, M.; Huysmans, M.; Lin, Z.; De Winter, F.; Vanneste, S.; Karimi, M.; Van de Velde, J.; et al. KIRA1 and ORESARA1 terminate flower receptivity by promoting cell death in the stigma of Arabidopsis. Nat. Plants 2018, 4, 365–375. [Google Scholar] [CrossRef]
Prigge, V.; Xu, X.; Li, L.; Babu, R.; Chen, S.; Atlin, G.N.; Melchinger, A.E. New Insights into the Genetics of in Vivo Induction of Maternal Haploids, the Backbone of Doubled Haploid Technology in Maize. Genetics 2012, 190, 781–793. [Google Scholar] [CrossRef]
Zhang, X.; Zhang, L.; Zhang, J.; Jia, M.; Cao, L.; Yu, J.; Zhao, D. Haploid induction in allotetraploid tobacco using DMPs mutation. Planta 2022, 255, 98. [Google Scholar] [CrossRef] [PubMed]
Zhong, Y.; Wang, Y.; Chen, B.; Liu, J.; Wang, D.; Li, M.; Qi, X.; Liu, C.; Boutilier, K.; Chen, S. Establishment of a DMP based maternal haploid induction system for polyploid Brassica napus and Nicotiana tabacum. J. Integr. Plant Biol. 2022, 64, 1281–1294. [Google Scholar] [CrossRef] [PubMed]
Zhong, Y.; Chen, B.; Wang, D.; Zhu, X.; Li, M.; Zhang, J.; Chen, M.; Wang, M.; Riksen, T.; Liu, J.; et al. In vivo maternal haploid induction in tomato. Plant Biotechnol. J. 2022, 20, 250–252. [Google Scholar] [CrossRef] [PubMed]
Nawade, B.; Bosamia, T.C.; Lee, J.H.; Jang, J.H.; Lee, O.R. Genome-wide characterization of the soybean DOMAIN OF UNKNOWN FUNCTION 679 membrane protein gene family highlights their potential involvement in growth and stress response. Front. Plant Sci. 2023, 14, 1216082. [Google Scholar] [CrossRef]
Lescot, M.; Déhais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Van de Peer, Y.; Rouzé, P.; Rombauts, S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30, 325–327. [Google Scholar] [CrossRef]
Chen, K.; Durand, D.; Farach-Colton, M. NOTUNG: A program for dating gene duplications and optimizing gene family trees. J. Comput. Biol. 2000, 7, 429–447. [Google Scholar] [CrossRef]
Guo, A.Y.; Zhu, Q.H.; Chen, X.; Luo, J.C. GSDS: A gene structure display server. Yi Chuan 2007, 29, 1023–1026. [Google Scholar] [CrossRef]
Roy, S.W.; Gilbert, W. The evolution of spliceosomal introns: Patterns, puzzles and progress. Nat. Rev. Genet. 2006, 7, 211–221. [Google Scholar] [CrossRef]
Zhong, Y.; Liu, C.; Qi, X.; Jiao, Y.; Wang, D.; Wang, Y.; Liu, Z.; Chen, C.; Chen, B.; Tian, X.; et al. Mutation of ZmDMP enhances haploid induction in maize. Nat. Plants 2019, 5, 575–580. [Google Scholar] [CrossRef]
Galluzzi, L.; Baehrecke, E.H.; Ballabio, A.; Boya, P.; Bravo-San Pedro, J.M.; Cecconi, F.; Choi, A.M.; Chu, C.T.; Codogno, P.; Colombo, M.I.; et al. Molecular definitions of autophagy and related processes. Embo J. 2017, 36, 1811–1836. [Google Scholar] [CrossRef]
Kasaras, A. Characterization of the Senescence-Associated Membrane Protein DMP1 and the DMP Family in Arabidopsis thaliana. Ph.D. Thesis, University of Berlin, Berlin, Germany, 2013. [Google Scholar] [CrossRef]
Takahashi, T.; Mori, T.; Ueda, K.; Yamada, L.; Nagahara, S.; Higashiyama, T.; Sawada, H.; Igawa, T. The male gamete membrane protein DMP9/DAU2 is required for double fertilization in flowering plants. Development 2018, 145, dev170076. [Google Scholar] [CrossRef] [PubMed]
Zhong, Y.; Chen, B.; Li, M.; Wang, D.; Jiao, Y.; Qi, X.; Wang, M.; Liu, Z.; Chen, C.; Wang, Y.; et al. A DMP-triggered in vivo maternal haploid induction system in the dicotyledonous Arabidopsis. Nat. Plants 2020, 6, 466–472. [Google Scholar] [CrossRef]
Kasaras, A.; Kunze, R. Dual-targeting of Arabidopsis DMP1 isoforms to the tonoplast and the plasma membrane. PLoS ONE 2017, 12, e0174062. [Google Scholar] [CrossRef] [PubMed]
Ma, Y.; Liu, H.; Wang, J.; Zhao, G.; Niu, K.; Zhou, X.; Zhang, R.; Yao, R. Genome-wide identification of the oat DMP gene family and its expression analysis in response to seed aging. BMC Genom. 2024, 25, 863. [Google Scholar] [CrossRef]
Hagan, I.; Yanagida, M. The product of the spindle formation gene sad1+ associates with the fission yeast spindle pole body and is essential for viability. J. Cell Biol. 1995, 129, 1033–1047. [Google Scholar] [CrossRef]
Malone, C.J.; Fixsen, W.D.; Horvitz, H.R.; Han, M. UNC-84 localizes to the nuclear envelope and is required for nuclear migration and anchoring during C. elegans development. Development 1999, 126, 3171–3181. [Google Scholar] [CrossRef]
Yuan, L.; Pan, J.; Zhu, S.; Li, Y.; Yao, J.; Li, Q.; Fang, S.; Liu, C.; Wang, X.; Li, B.; et al. Evolution and Functional Divergence of SUN Genes in Plants. Front. Plant Sci. 2021, 12, 646622. [Google Scholar] [CrossRef]
Ezer, D.; Shepherd, S.J.K.; Brestovitsky, A.; Dickinson, P.; Cortijo, S.; Charoensawan, V.; Box, M.S.; Biswas, S.; Jaeger, K.E.; Wigge, P.A. The G-Box Transcriptional Regulatory Code in Arabidopsis. Plant Physiol. 2017, 175, 628–640. [Google Scholar] [CrossRef]
Gangappa, S.N.; Maurya, J.P.; Yadav, V.; Chattopadhyay, S. The regulation of the Z-and G-box containing promoters by light signaling components, SPA1 and MYC2, in Arabidopsis. PLoS ONE 2013, 8, e62194. [Google Scholar] [CrossRef]
Kobayashi, K.; Obayashi, T.; Masuda, T. Role of the G-box element in regulation of chlorophyll biosynthesis in Arabidopsis roots. Plant Signal. Behav. 2012, 7, 922–926. [Google Scholar] [CrossRef]
Trishla, V.S.; Marriboina, S.; Boyidi, P.; Kirti, P.B. GUS-reporter based analysis of the promoter activity of Gossypium hirsutum NAC transcription factor, GhNAC4 that is induced by phytohormones and environmental stresses. Plant Cell Tissue Organ Cult. (PCTOC) 2020, 141, 643–654. [Google Scholar] [CrossRef]
Martínez-Hernández, A.; López-Ochoa, L.; Argüello-Astorga, G.; Herrera-Estrella, L. Functional properties and regulatory complexity of a minimal RBCS light-responsive unit activated by phytochrome, cryptochrome, and plastid signals. Plant Physiol. 2002, 128, 1223–1233. [Google Scholar] [CrossRef] [PubMed]
López-Ochoa, L.; Acevedo-Hernández, G.; Martínez-Hernández, A.; Argüello-Astorga, G.; Herrera-Estrella, L. Structural relationships between diverse cis-acting elements are critical for the functional properties of a rbcS minimal light regulatory unit. J. Exp. Bot. 2007, 58, 4397–4406. [Google Scholar] [CrossRef]
Hagen, G.; Guilfoyle, T. Auxin-responsive gene expression: Genes, promoters and regulatory factors. Plant Mol. Biol. 2002, 49, 373–385. [Google Scholar] [CrossRef]
Yamamoto, Y.Y.; Ichida, H.; Abe, T.; Suzuki, Y.; Sugano, S.; Obokata, J. Differentiation of core promoter architecture between plants and mammals revealed by LDSS analysis. Nucleic Acids Res. 2007, 35, 6219–6226. [Google Scholar] [CrossRef]
Niggeweg, R.; Thurow, C.; Kegler, C.; Gatz, C. Tobacco transcription factor TGA2.2 is the main component of as-1-binding factor ASF-1 and is involved in salicylic acid- and auxin-inducible expression of as-1-containing target promoters. J. Biol. Chem. 2000, 275, 19897–19905. [Google Scholar] [CrossRef]
Mönke, G.; Altschmied, L.; Tewes, A.; Reidt, W.; Mock, H.P.; Bäumlein, H.; Conrad, U. Seed-specific transcription factors ABI3 and FUS3: Molecular interaction with DNA. Planta 2004, 219, 158–166. [Google Scholar] [CrossRef]
Huang, J.; Hai, Z.; Wang, R.; Yu, Y.; Chen, X.; Liang, W.; Wang, H. Genome-wide analysis of HSP20 gene family and expression patterns under heat stress in cucumber (Cucumis sativus L.). Front. Plant Sci. 2022, 13, 968418. [Google Scholar] [CrossRef]
Yoshida, T.; Fujita, Y.; Maruyama, K.; Mogami, J.; Todaka, D.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Four Arabidopsis AREB/ABF transcription factors function predominantly in gene expression downstream of SnRK2 kinases in abscisic acid signalling in response to osmotic stress. Plant Cell Environ. 2015, 38, 35–49. [Google Scholar] [CrossRef]
Rouster, J.; Leah, R.; Mundy, J.; Cameron-Mills, V. Identification of a methyl jasmonate-responsive region in the promoter of a lipoxygenase 1 gene expressed in barley grain. Plant J. 1997, 11, 513–523. [Google Scholar] [CrossRef]
Sobajima, H.; Tani, T.; Chujo, T.; Okada, K.; Suzuki, K.; Mori, S.; Minami, E.; Nishiyama, M.; Nojiri, H.; Yamane, H. Identification of a Jasmonic Acid-Responsive Region in the Promoter of the Rice 12-Oxophytodienoic Acid Reductase 1 Gene OsOPR1. Biosci. Biotechnol. Biochem. 2007, 71, 3110–3115. [Google Scholar] [CrossRef] [PubMed]
Cao, H.; Zhang, X.; Li, F.; Han, Z.; Ding, B. A point mutation in the IAA14 promoter enhances callus formation and regeneration. Physiol. Mol. Biol. Plants 2024, 30, 1253–1263. [Google Scholar] [CrossRef] [PubMed]
Smita, S.; Katiyar, A.; Chinnusamy, V.; Pandey, D.M.; Bansal, K.C. Transcriptional Regulatory Network Analysis of MYB Transcription Factor Family Genes in Rice. Front. Plant Sci. 2015, 6, 1157. [Google Scholar] [CrossRef] [PubMed]
Zhang, Y.; Luo, M.; Cheng, L.; Lin, Y.; Chen, Q.; Sun, B.; Gu, X.; Wang, Y.; Li, M.; Luo, Y.; et al. Identification of the Cytosolic Glucose-6-Phosphate Dehydrogenase Gene from Strawberry Involved in Cold Stress Response. Int. J. Mol. Sci. 2020, 21, 7322. [Google Scholar] [CrossRef]
Dolferus, R.; Klok, E.J.; Ismond, K.; Delessert, C.; Wilson, S.; Good, A.; Peacock, J.; Dennis, L. Molecular basis of the anaerobic response in plants. IUBMB Life 2001, 51, 79–82. [Google Scholar] [CrossRef]
Siebertz, B.; Logemann, J.; Willmitzer, L.; Schell, J. Cis-analysis of the wound-inducible promoter wun1 in transgenic tobacco plants and histochemical localization of its expression. Plant Cell 1989, 1, 961–968. [Google Scholar] [CrossRef]
Dhatterwal, P.; Basu, S.; Mehrotra, S.; Mehrotra, R. Genome wide analysis of W-box element in Arabidopsis thaliana reveals TGAC motif with genes down regulated by heat and salinity. Sci. Rep. 2019, 9, 1681. [Google Scholar] [CrossRef]
Wang, K.; Wang, Z.; Li, F.; Ye, W.; Wang, J.; Song, G.; Yue, Z.; Cong, L.; Shang, H.; Zhu, S.; et al. The draft genome of a diploid cotton Gossypium raimondii. Nat. Genet. 2012, 44, 1098–1103. [Google Scholar] [CrossRef]
Zhu, T.; Liang, C.; Meng, Z.; Sun, G.; Meng, Z.; Guo, S.; Zhang, R. CottonFGD: An integrated functional genomics database for cotton. BMC Plant Biol. 2017, 17, 101. [Google Scholar] [CrossRef]
Letunic, I.; Bork, P. 20 years of the SMART protein domain annotation resource. Nucleic Acids Res. 2018, 46, D493–D496. [Google Scholar] [CrossRef]
Jones, P.; Binns, D.; Chang, H.-Y.; Fraser, M.; Li, W.; McAnulla, C.; McWilliam, H.; Maslen, J.; Mitchell, A.; Nuka, G.; et al. InterProScan 5: Genome-scale protein function classification. Bioinformatics 2014, 30, 1236–1240. [Google Scholar] [CrossRef] [PubMed]
Hallgren, J.; Tsirigos, K.D.; Pedersen, M.D.; Almagro Armenteros, J.J.; Marcatili, P.; Nielsen, H.; Krogh, A.; Winther, O. DeepTMHMM predicts alpha and beta transmembrane proteins using deep neural networks. bioRxiv 2022. [Google Scholar] [CrossRef]
Crooks, G.E.; Hon, G.; Chandonia, J.M.; Brenner, S.E. WebLogo: A sequence logo generator. Genome Res. 2004, 14, 1188–1190. [Google Scholar] [CrossRef] [PubMed]
Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef]
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]
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]
Kumar, S.; Stecher, G.; Suleski, M.; Hedges, S.B. TimeTree: A Resource for Timelines, Timetrees, and Divergence Times. Mol. Biol. Evol. 2017, 34, 1812–1819. [Google Scholar] [CrossRef]
Darby, C.A.; Stolzer, M.; Ropp, P.J.; Barker, D.; Durand, D. Xenolog classification. Bioinformatics 2016, 33, 640–649. [Google Scholar] [CrossRef]
Hu, Y.; Chen, J.; Fang, L.; Zhang, Z.; Ma, W.; Niu, Y.; Ju, L.; Deng, J.; Zhao, T.; Lian, J.; et al. Gossypium barbadense and Gossypium hirsutum genomes provide insights into the origin and evolution of allotetraploid cotton. Nat. Genet. 2019, 51, 739–748. [Google Scholar] [CrossRef]
Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]

Figure 1. Phylogenetic relationship of 240 DMP genes from 24 species, including Arabidopsis and cotton, showing a classification of DMPs in various groups. (A) MEGA 11 software was used to construct phylogenetic trees using the neighbor-joining (NJ) method with 1000 bootstrap replications. Each color shows a different group (I–V). Species are denoted with different colors and symbols in the phylogenetic tree. (B) Classification and arrangement of DMP genes in various groups based on a phylogenetic tree. Plant species classes are highlighted in various colors on the left. Five groups on the right side are made based on the phylogenetic tree containing DMP genes.

Figure 2. Comparison of the gene structure, conserved motifs and domains in DMP genes between Arabidopsis, cotton and 22 other species. (A) The NJ phylogenetic tree was constructed in MEGA 11 software using DMP sequences from 24 species. (B) The conserved protein motif composition of all DMPs and motif 1–10 is shown in various colored boxes. (C) Representation of the domain of Arabidopsis and other species. Yellow boxes represent the DMP domain in the sequences. (D) Exon and Intron structure of DMP genes of all Selected species. The red box represents exon or coding region, and the grey line represents intron with the scale at the bottom for measurement length of sequences.

Figure 3. Promoter region analysis of DMP sequences from two monocot and two dicot species along with phylogenetic tree of the DMPs for each species separately. (A,B) = Dicots (Arabidopsis thaliana and Gossypium raimondi). (C,D) = Monocots (Sorghum bicolor and Zea mays).

Figure 4. Evolutionary events in the DMP gene family in plants. Numbers of gene losses are shown in red color after “+” and duplications in blue color after “−”in the phylogenetic tree on the left side.

Figure 5. Expression profiles of DMP genes in various organ tissue from two monocot and two dicot species. Raw RNA-seq data were procured from relevant databases. (A) Arabidopsis thaliana from TAIR database (TAIR, v.11, http://www.arabidopsis.org). (B) Avena sativa and (C) Glycine max from Phytozome (https://phytozome.jgi.doe.gov/pz/portal.html/, accessed on 5 March 2024). (D) Gossypium hirsutum from COTTONGEN (http://www.cottongen.org/, accessed on 7 March 2024) and the Cotton Functional Genomics Database (CottonFGD, https://cottonfgd.org/). DMP gene expression levels are represented with various colors, with blue color indicating the least expression and red indicates the highest expression. Eudicot species (A,C,D); monocot species (B).

Figure 6. Expression levels of DMP sequences in various tissues of cotton. GhActin was used as a reference. Error bars denote standard deviation measured from three experiments.

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Ahmad, Z.; Tian, D.; Li, Y.; Aminu, I.M.; Tabusam, J.; Zhang, Y.; Zhu, S. Characterization, Evolution, Expression and Functional Divergence of the DMP Gene Family in Plants. Int. J. Mol. Sci. 2024, 25, 10435. https://doi.org/10.3390/ijms251910435

AMA Style

Ahmad Z, Tian D, Li Y, Aminu IM, Tabusam J, Zhang Y, Zhu S. Characterization, Evolution, Expression and Functional Divergence of the DMP Gene Family in Plants. International Journal of Molecular Sciences. 2024; 25(19):10435. https://doi.org/10.3390/ijms251910435

Chicago/Turabian Style

Ahmad, Zeeshan, Dingyan Tian, Yan Li, Isah Mansur Aminu, Javaria Tabusam, Yongshan Zhang, and Shouhong Zhu. 2024. "Characterization, Evolution, Expression and Functional Divergence of the DMP Gene Family in Plants" International Journal of Molecular Sciences 25, no. 19: 10435. https://doi.org/10.3390/ijms251910435

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

Article Menu

Characterization, Evolution, Expression and Functional Divergence of the DMP Gene Family in Plants

Abstract

1. Introduction

2. Results

2.1. Identification of the DMP Proteins in Various Species

2.2. Phylogenetic Analysis of the DMP Gene Family

2.3. Gene Structure, Motifs and Domain Analysis of DMPs

2.4. Promoter Region Analysis in Dicots and Monocots

2.5. Duplication and Loss of DMP Genes in Plants during Evolution

2.6. Tissue Expression Levels of DMPs in Monocots and Eudicots

2.7. Tissues Expression of DMPs in Cotton

3. Discussion

4. Materials and Methods

4.1. Identification and Retrieval of DMP Gene Family

4.2. Amino-Acid Sequence Alignment and Phylogenetic Analysis

4.3. Conserved Motifs, Domain and Gene Structure Analysis

4.4. Promoter Analysis of Monocot and Eudicot Species

4.5. Gene Loss and Duplication Event Analysis

4.6. Expression Profiles of DMPs in Monocots and Eudicots

4.7. RNA Extraction and qRT-PCR Analysis

5. Conclusions

Supplementary Materials

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI