Next Article in Journal
Post-Bloom CPPU Application Is Effective at Improving Fruit Set and Suppressing Coloration but Ineffective at Increasing Fruit Size in Litchi
Next Article in Special Issue
Pectin, Lignin and Disease Resistance in Brassica napus L.: An Update
Previous Article in Journal
Gibberellic Acid Concentrations and Storage of Caryocar brasiliense (Caryocaraceae) Seeds Propagated in Tubes
Previous Article in Special Issue
Systematic Characterization of Brassica napus HIR Gene Family Reveals a Positive Role of BnHIR2.7 in Sclerotinia sclerotiorum Resistance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Structure Variations and Expression Characteristics of DMP8 and DMP9 Genes in Brassicaceae

Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Zhongguancun Nan Da Jie No.12, Haidian District, Beijing 100081, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Horticulturae 2022, 8(11), 1095; https://doi.org/10.3390/horticulturae8111095
Submission received: 20 October 2022 / Revised: 14 November 2022 / Accepted: 17 November 2022 / Published: 21 November 2022
(This article belongs to the Special Issue Advances in Brassica Crops Genomics and Breeding, 2nd Edition)

Abstract

:
Doubled haploid (DH) technology based on in vivo haploid induction (HI), which is used to obtain true-breeding lines within a single generation, is a technique that significantly increases modern crop-breeding efficiency. Recently, dicot Arabidopsis thaliana lines containing mutations in DMP8/9 were used as haploid inducer lines, but the use of this new HI mechanism is limited in Brassicaceae species, which include many important vegetable, oil, and fodder crops. Here, we investigated the phylogenetic distribution of the DMP8 and DMP9 homologous genes from 26 sequenced Brassicaceae species. We found that DMP8 only exists in the tribe Arabideae, while multiple copies of the DMP9 gene are presenting in all the investigated Brassicaceae species. The syntenic DMP9 genes were divided into two groups derived from the S genomic block and R genomic block, respectively. We further investigated the duplication, structure variations, and expression of the DMP9 genes in Brassica species that had undergone an extra whole-genome triplication. Our results revealed that DMP9 was lost in the most fractionated (MF2) subgenome, and the retained DMP9s in the least fractionated (LF) subgenome and medium fractionated (MF1) subgenome showed diversified expression patterns, indicating their functional diversification. Our results will be useful for obtaining the target DMP genes for the establishing of HI lines in Brassicaceae crops.

1. Introduction

Haploid induction (HI) technology is used to obtain homozygous lines in a single generation. Compared with the traditional breeding approach, which requires multiple generations of selfing or backcrossing, the approach employing doubled haploids increases crop-breeding efficiency [1]. As important cultivated crops, haploids of Brassica species could be obtained through in vitro anther or microspore cultures; however, the procedure is time-consuming and complex, and the induction is genotype-dependent [2,3].
In vivo HI technology, which has been widely used in monocot breeding, has revolutionized traditional crop breeding, especially in maize (Zea mays). The loss of function of ZmPLA1/MTL/NLD (patatin-like phospholipase A1/MATRILINEAL/NOT LIKE DAD) in maize, which encodes a sperm cell-expressed phospholipase [4,5,6], can trigger maternal HI, and the HI rate (HIR) can reach 2–3%. This new approach has been extended to other monocots, including rice [7] and wheat [8,9], but has not been extensively applied to dicots owing to the lack of identifiable ZmPLA1/MTL/NLD orthologues [10]. In contrast, another maternal haploid inducer gene identified in maize, Zea mays DOMAIN OF UNKNOWN FUNCTION 679 membrane protein (ZmDMP), belonging to an unknown plant-specific protein family, has a conserved amino acid sequence in both monocots and dicots. The HIR in ZmDMP mutants is 0.1–0.3%, which can be increased 5–6-fold when zmdmp is combined with mtl/zmpla1/nld alleles. In Arabidopsis, the DMP protein family has ten members (DMP1 to 10), and all DMP proteins have four transmembrane domains. However, a phylogenetic analysis showed that AtDMP8 and AtDMP9 had the closest relation with ZmDMP. The loss of function of AtDMP8 and AtDMP9 triggers in vivo maternal HI in Arabidopsis thaliana [10].
A previous study found that AtDMP8 and AtDMP9, expressed explicitly in both generative and sperm cells, are involved in a gamete interaction that leads to correct double fertilization in flowering plants [11]. To date, haploid inducer systems with DMP genes have been reported for several other dicot plants, including tomato (Solanum lycopersicum), tobacco (Nicotiana tabacum), Medicago truncatula, Brassica napus, and Brassica oleracea [12,13,14,15].
The Brassicaceae family, containing 340 genera and 4636 species, is a very important economic plant family [16]. The Brassicaceae family, which includes many important vegetable, oil, and fodder crops [17], provides nutrients and phytochemicals for human health. The Brassicaceae family has undergone three major ancient whole-genome duplication (WGD) events: a paleohexaploidy (γ) event and two paleotetraploidy (β and α) events, which are referred to as the At-γ, At-β, and At-α of the core Brassicaceae [3,18]. A key agricultural genus of the Brassicaceae family, the Brassica genus, had undergone an extra whole-genome triplication (WGT) event. These duplication events, followed by rapid DNA sequence divergence and fractionation of the ancestral gene sequence, reduced the gene similarity and increased the adaptive ability to survive and diversify in the changing environment [18,19]. Furthermore, WGDs tend to cause expression divergence within and between species, contributing to morphological diversity and expanding regulatory networks [20,21,22].
Because of the complex genome of Brassicaceae, it is essential to analyze the gene evolution and expression characteristics, especially in those crops for which it is difficult to obtain transgenic plants, such as Brassica rapa. Herein, we investigated the evolutionary history of the DMP8 and DMP9 genes in 26 sequenced Brassicaceae species by constructing phylogenetic trees and examining the syntenic relationship of the DMP9 genes. Specifically, we tested whether the duplicated DMP9 genes represented expression divergence within or between Brassica species. It is essential to study their evolution to accurately find the functional target of DMP9 genes, and such a study also paves the way for applying this HI system in B. rapa.

2. Materials and Methods

2.1. Plant Materials

The germinated seeds of B. rapa “Chiifu”, B. rapa “B9008”, and B. oleracea “JZS” were sown in pots in a greenhouse without climate control in Beijing. Plants were grown during the spring of 2021 under normal growth conditions. Pollen was collected from these plants for DNA isolation, RNA isolation, and RNA-seq.

2.2. Identification and Phylogenetic Analysis of DMP9 Homologous Genes

Genes and genome datasets for 26 Brassicaceae species, including Brassica rapa (v3.0), Brassica oleracea, Brassica nigra, Brassica juncea, Brassica napus, Brassica carinata, Boechera retrofracta, Boechera stricta, Crucihimalaya himaliaca, Capsella grandiflora, Capsella rubella, Camelina sativa, Arabidopsis lyrata, Arabidopsis halleri, Arabidopsis thaliana, Descurainia sophia, Cardamine hirsuta, Leavenworthia alabamica, Raphanus sativus, Sisymbrium irio, Isatis indigotica, Arabis alpina, Schrenkiella parvula, Thlaspi arvense, Thellungiella halophila, and Aethionema arabicum, were downloaded from the BRAD (http://brassicadb.cn, accessed on 10 November 2022). The 18 B. rapa genomes included in the pan-genome were obtained from Cai et al. [23]. Based on the results of all 26 Brassicaceae species and a B. rapa pan-genome multiple sequence alignment, a neighbor-joining tree was constructed using a bootstrap test with 1000 replicates. The phylogenetic tree was visualized in MEGA 7.0 [24]. For identifying syntenic genes, we used the function of “Syntenic Gene @ Subgenomes” of the BRAD database (http://brassicadb.cn/syntenic-gene/, accessed on 10 November 2022), which identifies syntenic genes in Brassicaceae species based on SynOrths [25,26].

2.3. Transcriptome Sequencing

To investigate the expression patterns of DMP9 homologous genes, we collected flower buds at differential development stages from B. rapa “Chiifu”, B. rapa “B9008”, and B. oleracea “JZS” plants, each having three biological replicates. Pollen was collected from flower buds of different lengths and observed under a microscope (Carl Zeiss Primo Star 415500-0057-000).
Pollen development was divided into six stages: diploid nucleus, haploid nucleus, microspore, uninuclear pollen, bicellular pollen, and tricellular pollen (mature pollen), according to observation by a microscope. In total, 54 samples were collected and used for RNA-seq via the Illumina platform.

2.4. Conserved Domain and Gene Structure Analysis

The gene structures of DMP9 copies in the B. rapa pan-genome were generated using TBtools software [27]. The parameters were set to the default values [27].

2.5. RNA-Seq Data Analysis

The total RNA was extracted with a Trizol reagent from the samples of flower buds at the six pollen developmental stages collected from B. rapa “B9008”, B. rapa “Chiifu”, and B. oleracea “JZS”. The RNA was used for RNA-seq via the Illumina platform. All raw reads were filtered by fastp [28] using the parameter “-z 4 -q 20 -u 30 -n 5”. All clean reads were mapped onto the corresponding genome using Hisat2 (version 2.2.0) [29]. We calculated the TPM (transcripts per million) value of each gene using StringTie (version 2.1.3b) [30]. Data from different experimental sets were analyzed for statistical significance using one-way analysis of variance (ANOVA) with GraphPad Prism 7.0 software. A p-value < 0.05 was considered significant.

3. Results

3.1. AtDMP8 and AtDMP9 Homologous Genes in 26 Brassicaceae Genomes

The full-length DNA sequences of AtDMP8 and AtDMP9 were used for a BLASTP search (http://brassicadb.cn/#/BLAST/, accessed on 10 November 2022) to identify homologous genes with the identity ≥ 85% in Brassicaceae species. The results showed that a few homologous segments, showing a high identity with the DMP9 gene, were not annotated in the genomes (Supplementary Table S1). Thereafter, we performed a syntenic analysis for the DMP8 and DMP9 homologous genes using SynOrths [25]. The homologous genes were divided into two parts: syntenic genes and non-syntenic genes. Syntenic genes from different species are orthologs located on syntenic fragments. Therefore, they often share similar functions and originate from a common ancestor [25,31]. Here, we analyzed the evolution of syntenic homologous genes of DMP8 and DMP9 in the 26 species of the Brassicaceae family.
A previous study showed that the DMP8 in the Brassica genus was lost [15]. Our results showed that the DMP8 genes located in the A block, showing a good syntenic relationship, only existed in A. thaliana and A. lyrata (Supplementary Table S1, Figure 1a). The DMP8 duplication is a unique duplication event that occurs in the genus Arabidopsis. Moreover, we identified 50 AtDMP9-like genes whose identities were greater than 85% in the 26 sequenced Brassicaceae genomes (Supplementary Table S1). This result is consistent with the expectation of multiple DMP9 copies in the Brassicaceae species, derived from multiple whole-genome duplications of the Brassicaceae ancestor and the triploidization event during the genome evolution of Brassiceae. The number of homologous DMP9 genes varied in different Brassicaceae species, with two DMP9 homologs in 18 species, one DMP9 homolog in five species, and four DMP9 homologs in three species, accounting for 69.23, 19.23, and 11.54% of all the Brassicaceae species, respectively.
Taking the maize ZmDMP gene as an outgroup, we constructed a phylogenetic tree of the DMP8 and DMP9 homologous genes of the 26 sequenced Brassicaceae species based on their DNA sequences (Figure 1a). The phylogenetic tree was visualized in MEGA 7.0 [24]. The phylogeny results showed that the homologous genes of DMP8 and DMP9 clustered into two groups: Group 1 and Group 2. Almost all sequenced Brassicaceae species in Group 1 (excluding R. sativus) retained DMP9 genes, while only A. thaliana and A. lyrata retained DMP8 genes. In contrast, in Group 2, only the genus Brassica and C. himaliaca possessed DMP9 genes (Figure 1a). Overall, the phylogenetic analysis revealed that all sequenced Brassicaceae species retained DMP9 genes.
Based on the 24 basic genomic blocks (GBs, A-X) of Brassicaceae genomes, we performed the syntenic analysis of DMP9 homologous genes using SynOrths [25,32] (Figure 1b). Our study revealed that the DMP9 homologous genes in Brassicaceae were clustered into two groups: Group 1 (green branch) and Group 2 (red branch). In addition, the genes in Group 1 were mainly derived from the S block, and those in Group 2 were mainly derived from the R block. Alternatively, the numbers of the DMP9 homologous genes varied across the genomic blocks of 26 Brassicaceae species. Most Brassicaceae species possessed DMP9 genes in the S block (excluding A. halleri, R. sativus, and A. alpina).
In the R block, the DMP9 homologous genes only presented in the genus Brassica, R. sativus, and C. himaliaca. Alternatively, 13 DMP9 homologous genes had no syntenic relationship located in Group 1, and one of B. carinata was located in Group 2 (Figure 1). There was only one exception for DMP9 in B. carinata, where tandem duplicated copies located in the R block rather than homoeologous copies were observed. Among the 26 Brassicaceae species, only R. sativus retained one DMP9 paralogous copy after the WGD and WGT events.

3.2. The Evolution of DMP9 Homologous Genes between Species with an Extra WGT

As the genus Brassica had undergone WGT, and the genome possesses three subgenomes—LF (least fractionated), MF1 (medium fractionated 1), and MF2 (most fractionated 2)—it was expected that there would be three DMP9 copies in the genomes of Brassica species [33]. To verify the impact of WGT on the DMP9 evolution in Brassicaceae, we analyzed the loss of DMP9 copies in each of the three subgenomes in the genus Brassica (B. rapa, B. oleracea, B. nigra, B. juncea, B. napus, and B. carinata), and R. sativus. The results showed that the DMP9 copy in the MF2 subgenome had been lost in all the analyzed Brassica species. All the diploid Brassica genus species, B. rapa (AA), B. oleracea (CC), and B. nigra (BB), retained two DMP9 copies, while in R. sativus, it remained a single duplicate. Among the three allotetraploid species, both B. juncea (AABB) and B. napus (AACC) contained four DMP9 copies, and each of the two diploid ancestors contributed two copies located in the LF and MF1 subgenomes; however, B. carinata (BBCC) kept four copies only in the B subgenome, including one copy located in the LF subgenome, two copies in a tandem duplication located in the MF1 subgenome, and the last one being a nonsyntenic gene.
Interestingly, we found that the DMP9 homologous genes of species in the genus Brassica located in the S block belonged to the LF subgenome, and the copy located in the R block belonged to the MF1 subgenome. Specifically, the LF and MF1 subgenomes retained only one S block and one R block in the diploid Brassica species (B. rapa, B. oleracea, and B. nigra), respectively (Figure 1b). This result demonstrated that the WGT event accelerated the loss of DMP9 copies.
Moreover, during the formation of the DMP9 loci in B. carinata, the DMP9 gene was amplified by a tandem duplication (TD) event and a segmental duplication, which were smaller-scale duplications (including TD, segmental duplication, and gene transposition duplication) [34,35]. These distributions are shown in Figure 1.
Note that all these analyses on the evolutionary history of DMP9 homologous genes in the Brassicaceae species are based on one reference genome per species, and that the gene copy number and functionality of some gene copies may be different within a species.

3.3. DMP9 Homologous Genes in B. rapa

To investigate the impact of the DMP9 gene evolution within a species, we used B. rapa as a model. B. rapa shares a common ancient hexaploid ancestor with Brassica plants and has undergone an additional WGT, which has played an essential role in the speciation and morphotype diversification of Brassica plants [2].
B. rapa has a rich variety of morphotypes showing extreme traits, such as the leafy heads of Chinese cabbage, the enlarged roots or stem tubers of turnips, and the enlarged stems of caixin. In this study, we examined the B. rapa pan-genome derived from the genome assembly of the 18-accession evolution of DMP9 genes. The result showed that there were 41 DMP9 homologous genes in the B. rapa pan-genome. Among the 18 B. rapa genomes, 14 contained two DMP9 homologs, accounting for 83.33%. Three genomes contained three DMP9 homologous genes, and one genome (TBA) contained four DMP9 homologs (Table 1, Figure 2a). We noticed that the homologs of DMP9 were tandemly duplicated in the R block of two genomes (Figure 2b). The results showed that the DMP9 homologous genes had a relatively complex evolution route in B. rapa.
Taking the A. thaliana AtDMP9 gene as an outgroup, we constructed a phylogenetic tree of the DMP9 homologous genes identified from the 18 B. rapa genomes based on their DNA sequences (Figure 2a). A syntenic relationship analysis showed that the DMP9 homologous genes were retained in the LF (S block) and MF1 (R block) subgenomes, but were lost in the MF2 subgenome (Figure 2b). The homologous genes were more conserved in the R block of the MF1 subgenome compared to the S block in the LF subgenome. Although most of the B. rapa genomes retained two DMP9 homologous genes, four (CCA, BRO, OIB, and TBA) retained more DMP9 copies during their evolution. However, it seemed that the tandem duplication of DMP9 was not subspecies-specific, as the two Chinese cabbage genomes were different at this locus. CCA and BRO retained a tandem duplicated gene, while CCB did not. TBA retained four copies, but only two of them were annotated. To explore the features and potential functions of DMP9 copies, we performed a gene structural analysis based on the gene sequence similarity of the encoded proteins (Figure 2c). The results showed that most DMP9 copies indicated high conservation. The tandem duplication in B. rapa CCA and BRO shared a high similarity, but there was a 5745 bp insertion in the tandemly duplicated copy of CCA. In short, we proposed that the DMP9 gene not only diversified among the Brassicaceae species, but also within species, such as in the B. rapa pan-genome.

3.4. Transcriptome Analysis of DMP9 Homologous Genes in B. rapa and B. oleracea

Expression divergence between duplicated genes has long been of interest to geneticists and evolutionary biologists because it is considered as the first step in functional divergence between duplicate genes, which increases the chance of the retention of duplicated genes in a genome [36,37]. In a previous study, it has been found that DMP8 was lost and DMP9 retained two homologous genes in B. oleracea. The BolDMP9.1 gene was nonfunctional, owing to a 1-bp deletion in the exon, while BolDMP9.2 mutants could induce a maternal haploid with an HIR of about 2.35% [15]. B. oleracea and B. rapa showed a close relationship with A. thaliana. Based on the evolution analysis and previous studies of the DMP9 homologous genes in Brassicaceae, we used B. oleracea “JZS”, B. rapa “Chiifu” (B. rapa reference genome Chiifu), and B. rapa “B9008” (CCA genome) to further study the expression patterns of DMP9 genes in B. rapa (Table 2 and Figure 3a).
To further compare the differences of DMP9 copies between the two B. rapa species, we analyzed the corresponding relationship and amino acid sequences with A. thaliana. The results showed that although the DMP9 copies contained the conserved domain, in the third and fourth transmembrane (TM) spans, the proteins encoded by the tandem duplicated copies CCADMP9.1 and CCADMP9.2 showed a greater structural divergence than others (Figure 3b).
To explore the expression patterns of DMP9 homologous genes, we analyzed RNA-seq data generated from different pollen developmental stages, including the diploid nucleus, haploid nucleus, microspore, uninuclear pollen, bicellular pollen, and tricellular pollen (mature pollen) (Figure 3c). An expression analysis of the DMP9 homologous genes revealed that the DMP9 homologous genes of B. rapa located in the LF (BraDMP9.1 and CCADMP9.1) and MF2 (BraDMP9.2 and CCADMP9.3) subgenomes were highly expressed in mature pollen, but were expressed at much lower levels in immature pollen of different stages, as in A. thaliana (Figure 3c) [15]. In addition, the expression level of DMP9 copies in the MF1 subgenome was higher than that of the genes in the LF subgenome. In contrast, the tandem duplicated DMP9 homologous gene (CCADMP9.2) in B. rapa “B9008” was most highly expressed at the haploid nucleus stage, indicating that the DMP9 gene diversified in expression patterns. In B. oleracea, the two DMP9 copies (BolDMP9.1 and BolDMP9.2) were highly expressed in mature pollen, which was consistent with that observed in A. thaliana and B. rapa, and the MF1 copy (BolDMP9.2) had a higher expression level than the LF copy (BolDMP9.1).

4. Discussion

Doubled haploid (DH) technology based on in vivo HI has led to a new approach to crop breeding. Significant breakthroughs in DH production have been made in the last few years by identifying maize genes that induce maternal haploid embryos in vivo. In A. thaliana, the loss function of the homologous genes of maize, ZmDMP, AtDMP8, and AtDMP9, could trigger a maternal haploid [10]. Significant breakthroughs in DH production have been made in the last few years by identifying maize genes that induce maternal haploid embryos in vivo.
Recently, more studies have reported that a loss of function of the DMP8 and DMP9 homologous genes could trigger HI in Medicago truncatula, tomato, Brassica napus, Brassica oleracea, and Nicotiana tabacum [12,13,14,38]. However, the use of this new HI mechanism is still limited in Brassicaceae species.
The Brassicaceae is a medium-sized family that includes many important economical crops. These crops have undergone WGDs that were followed by gene losses during diploidization. The Brassica genus underwent an extra WGT event, resulting in a complex relationship among the duplicated DMP genes. Therefore, the diversification of the DMP8 and DMP9 genes is central to the HI function in Brassicaceae species. In this study, we took advantage of 26 sequenced Brassicaceae genomes to investigate the evolution and diversification of the DMP8 and DMP9 homologous genes. We proposed that the DMP8 gene was retained in the model plant A. thaliana and its congener A. lyrata, but that it cannot be found in other Brassicaceae species. Although it is likely that DMP8 was a unique gene duplication in the common ancestor of A. thaliana and A. lyrata, we did not exclude the possibility that it may have been lost in other Brassicaceae species. Therefore, the DMP9 homologous genes should be the major focus when we plan to produce HI in Brassica crops. The importance of DMP9 has been demonstrated by the successful generation of HI in B. oleracea [15].
Thereafter, we identified DMP9 homologous genes in 26 sequenced Brassicaceae species through a BLASTP search. The results showed that DMP9 copies were diverse. Among the 51 DMP9 copies identified from the Brassicaceae genomes, four copies were not annotated on the genome. The synteny analysis showed that the DMP9 homologous genes were divided into syntenic genes and non-syntenic genes. The syntenic DMP9 genes are located in two genomic blocks: the S block and the R block. Most of the DMP9 copies were located in the S block, while only the genus Brassica, R. sativus, and C. himaliaca retained one copy in the R block. We also found that the directions of DMP9 copies were different in the genomes. In addition to the gene loss and gene fragment duplication of DMP9 copies, we also observed the tandemly duplicated DMP9 in the B. carinata genome. In summary, the diversification of DMP9 copies makes it more difficult to obtain the haploid inducer lines in Brassicaceae. Based on the currently sequenced genome data, we analyzed and obtained valuable information. We believe that after the assembly of these genomes has been improved, we can determine the locations of those DMP genes that are not allocated.
At present, haploid inducer systems with dmp9 genes have been reported for B. napus and B. oleracea in the Brassica genus [14,15,38]. For example, in B. napus, one genome retained three DMP9 copies with a gene loss on the LF subgenome [14], while another genome retained four DMP9 copies [37]. This result is consistent with what we found in the B. rapa pan-genome. The copy number variations detected in the pan-genome were all due to DMP9 gene duplication on the LF subgenome, while the MF subgenome seemed more stable than the LF subgenome. In B. oleracea, Zhao et al. [15] revealed that two DMP9 copies were retained. However, the DMP9 copy on the LF subgenome lost its function due to a 1-bp deletion in the exon. These results indicated that DMP9 genes varied in the copy number, and it was necessary to determine the appropriate copy as the target gene in the construction of HI in Brassica crops.
Regarding the loss of function of the DMP9 duplicates in B. napus and B. oleracea, we examined the DMP9 duplicate sequences and expression patterns in the B. rapa pan-genome. The results showed that the DMP9 copies on the LF did not contain the deletion of that of B. oleracea in 18 B. rapa genomes. However, the B. rapa genomes retained tandem duplicates. An expression analysis revealed that the tandem duplicates showed expressional divergence, indicating the functional divergence between duplicate genes. Local gene duplication generates tandem duplicated genes and is ubiquitous during genome evolution [39]. In Brassica species, although the haploid inducer lines have been obtained in B. napus and B. oleracea, whether all the varieties or subspecies can obtain the haploid inducer needs further analysis.
In addition to the role in haploid induction, it is reported that AtDMP8 and AtDMP9 regulate HAP2/GCS1 trafficking for egg–sperm fusion, and that this function of sperm cell activation is conserved in seed plants [40]. The DMP9 protein is involved in gamete interactions that lead to correct double fertilization in flowering plants [11,41]. The increased number of DMP9 genes in Brassica and their divergence in gene expression indicate their functional differentiation. Therefore, our research on the evolutionary relationship of DMP8 and DMP9 genes is of great value in elucidating the gene function of DMP8 and DMP9 duplicates.
In this study, the DMP9 evolution analysis has brought a lot of valuable information regarding the HI line in Brassicaceae species, especially in crops such as B. rapa, for which it is difficult to obtain transgenic plants. The target genes can be found through an evolutionary analysis, which can greatly improve the efficiency of obtaining haploid inducer lines, especially for crops with multi-copy genes. Although the genetic transformation of B. rapa has not been well developed, we can use other biological methods to further evaluate the gene function of DMP9 homologous genes, such as in situ reverse transcription polymerase chain reaction (in situ RT-PCR) or insertion–deletion (InDel) markers utilizing high-resolution melting (HRM) curve analysis, which have been exploited and used to detect the expression position and the population genetic analysis [42,43]. On the whole, the results from our present study lay the foundation for establishing the B. rapa HI line.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/horticulturae8111095/s1, Table S1: List of DMP8 and DMP9 homologous genes in the 26 Brassicaceae species [44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64].

Author Contributions

T.Z. and X.C. analyzed and interpreted the data. T.Z. drafted and revised the manuscript. T.Z. grew plants, collected tissues, and extracted RNA. J.L., L.Z., J.W. and X.W. improved the manuscript. J.W. and X.W. conceived the research, supervised the experiment and data analysis, and modified the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Beijing Joint Research Program for Germplasm Innovation and New Variety Breeding under Grant No. G20220628003.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The RNA-seq data can be downloaded from genome database of BRAD (Link: http://39.100.233.196:82/download_genome/datasets/DMPpaper/, accessed on 20 November 2022).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jacquier, N.; Gilles, L.M.; Pyott, D.E.; Martinant, J.P.; Widiez, T. Puzzling out plant reproduction by haploid induction for innovations in plant breeding. Nat. Plants 2020, 6, 610–619. [Google Scholar] [CrossRef] [PubMed]
  2. Cheng, F.; Wu, J.; Wang, X.W. Genome triplication drove the diversification of Brassica plants. Hortic Res. 2014, 1, 14024. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Warwick, S.I.; Francis, A.; Al-Shehbaz, I.A. Brassicaceae: Species checklist and database on CD-Rom. Plant Syst. Evol. 2006, 259, 249–258. [Google Scholar] [CrossRef]
  4. Gilles, L.M.; Khaled, A.; Laffaire, J.B.; Chaignon, S.; Gendrot, G.; Laplaige, J.; Berges, H.; Beydon, G.; Bayle, V.; Barret, P.; et al. Loss of pollen-specific phospholipase NOT LIKE DAD triggers gynogenesis in maize. EMBO J. 2017, 36, 707–717. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Kelliher, T.; Starr, D.; Richbourg, L.; Chintamanani, S.; Delzer, B.; Nuccio, M.L.; Green, J.; Chen, Z.Y.; McCuiston, J.; Wang, W.L.; et al. MATRILINEAL, a sperm-specific phospholipase, triggers maize haploid induction. Nature 2017, 542, 105. [Google Scholar] [CrossRef] [PubMed]
  6. Liu, C.X.; Li, X.; Meng, D.X.; Zhong, Y.; Chen, C.; Dong, X.; Xu, X.W.; Chen, B.J.; Li, W.; Li, L.; et al. A 4-bp Insertion at ZmPLA1 Encoding a Putative Phospholipase A Generates Haploid Induction in Maize. Mol. Plant 2017, 10, 520–522. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Yao, L.; Zhang, Y.; Liu, C.X.; Liu, Y.B.; Wang, Y.L.; Liang, D.W.; Liu, J.T.; Sahoo, G.; Kelliher, T. OsMATL mutation induces haploid seed formation in indica rice. Nat. Plants 2018, 4, 530–533. [Google Scholar] [CrossRef]
  8. Liu, C.X.; Zhong, Y.; Qi, X.L.; Chen, M.; Liu, Z.K.; Chen, C.; Tian, X.L.; Li, J.L.; Jiao, Y.Y.; Wang, D.; et al. Extension of the in vivo haploid induction system from diploid maize to hexaploid wheat. Plant Biotechnol. J. 2020, 18, 316–318. [Google Scholar] [CrossRef] [Green Version]
  9. Liu, H.Y.; Wang, K.; Jia, Z.M.; Gong, Q.; Lin, Z.S.; Du, L.P.; Pei, X.W.; Ye, X.G. Efficient induction of haploid plants in wheat by editing of TaMTL using an optimized Agrobacterium-mediated CRISPR system. J. Exp. Bot. 2020, 71, 1337–1349. [Google Scholar] [CrossRef]
  10. Zhong, Y.; Chen, B.J.; Li, M.R.; Wang, D.; Jiao, Y.Y.; Qi, X.L.; Wang, M.; Liu, Z.K.; Chen, C.; Wang, Y.W.; et al. A DMP-triggered in vivo maternal haploid induction system in the dicotyledonous Arabidopsis. Nat. Plants 2020, 6, 466–472. [Google Scholar] [CrossRef]
  11. Takahashi, T.; Mori, T.; Ueda, K.; Yamada, P.; 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] [Green Version]
  12. Wang, N.; Xia, X.Z.; Jiang, T.; Li, L.L.; Zhang, P.C.; Niu, L.F.; Cheng, H.M.; Wang, K.J.; Lin, H. In planta haploid induction by genome editing of DMP in the model legume Medicago truncatula. Plant Biotechnol. J. 2022, 20, 22–24. [Google Scholar] [CrossRef] [PubMed]
  13. Zhong, Y.; Chen, B.J.; Wang, D.; Zhu, X.J.; Li, M.R.; Zhang, J.Z.; Chen, M.; Wang, M.; Riksen, T.; Liu, J.C.; et al. In vivo maternal haploid induction in tomato. Plant Biotechnol. J. 2022, 20, 250–252. [Google Scholar] [CrossRef] [PubMed]
  14. 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]
  15. Zhao, X.; Yuan, K.; Liu, Y.; Zhang, N.; Yang, L.; Zhang, Y.; Wang, Y.; Ji, J.; Fang, Z.; Han, F.; et al. In vivo maternal haploid induction based on genome editing of DMP in Brassica oleracea. Plant Biotechnol. J. 2022, 20, 22–24. [Google Scholar] [CrossRef]
  16. Wu, J.; Liang, J.; Lin, R.; Cai, X.; Zhang, L.; Guo, X.; Wang, T.; Chen, H.; Wang, X. Investigation of Brassica and its relative genomes in the postgenomics era. Hortic. Res. 2022, 9, uhac182. [Google Scholar] [CrossRef]
  17. Zhu, B.; Liang, Z.; Zang, Y.; Zhu, Z.; Yang, J. Diversity of glucosinolates among common brassicaceae vegetables in China. Hortic. Plant J. 2022. [Google Scholar] [CrossRef]
  18. Franzke, A.; Lysak, M.A.; Al-Shehbaz, I.A.; Koch, M.A.; Mummenhoff, K. Cabbage family affairs: The evolutionary history of Brassicaceae. Trends Plant Sci. 2011, 16, 108–116. [Google Scholar] [CrossRef]
  19. Tang, H.B.; Wang, X.Y.; Bowers, J.E.; Ming, R.; Alam, M.; Paterson, A.H. Unraveling ancient hexaploidy through multiply-aligned angiosperm gene maps. Genome Res. 2008, 18, 1944–1954. [Google Scholar] [CrossRef] [Green Version]
  20. Carroll, S.B. Endless forms: The evolution of gene regulation and morphological diversity. Cell 2000, 101, 577–580. [Google Scholar] [CrossRef]
  21. Gu, Z.L.; Rifkin, S.A.; White, K.P.; Li, W.H. Duplicate genes increase gene expression diversity within and between species. Nat. Genet. 2004, 36, 577–579. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Lynch, M.; Conery, J.S. The evolutionary fate and consequences of duplicate genes. Science 2000, 290, 1151–1155. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Cai, X.; Chang, L.C.; Zhang, T.T.; Chen, H.X.; Zhang, L.; Lin, R.M.; Liang, J.L.; Wu, J.; Freeling, M.; Wang, X.W. Impacts of allopolyploidization and structural variation on intraspecific diversification in Brassica rapa. Genome Biol. 2021, 22, 166–189. [Google Scholar] [CrossRef] [PubMed]
  24. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  25. Cheng, F.; Wu, J.; Fang, L.; Wang, X.W. Syntenic gene analysis between Brassica rapa and other Brassicaceae species. Front. Plant Sci. 2012, 3, 198. [Google Scholar] [CrossRef] [Green Version]
  26. Chen, H.; Wang, T.; He, X.; Cai, X.; Lin, R.; Liang, J.; Wu, J.; King, G.; Wang, X. BRAD V3.0: An ungraded Brassicaceae database. Nucleic Acids Res. 2021, 50, D1432–D1441. [Google Scholar] [CrossRef]
  27. Chen, C.J.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.H.; Xia, R. TBtools: An Integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  28. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef]
  29. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [Green Version]
  30. Kovaka, S.; Zimin, A.V.; Pertea, G.M.; Razaghi, R.; Salzberg, S.L.; Pertea, M. Transcriptome assembly from long-read RNA-seq alignments with StringTie2. Genome Biol. 2019, 20, 1–13. [Google Scholar] [CrossRef]
  31. Lyons, E.; Pedersen, B.; Kane, J.; Alam, M.; Ming, R.; Tang, H.B.; Wang, X.Y.; Bowers, J.; Paterson, A.; Lisch, D.; et al. Finding and comparing syntenic regions among Arabidopsis and the outgroups papaya, poplar, and grape: CoGe with rosids. Plant Physiol. 2008, 148, 1772–1781. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Schranz, M.E.; Lysak, M.A.; Mitchell-Olds, T. The ABC’s of comparative genomics in the Brassicaceae: Building blocks of crucifer genomes. Trends Plant Sci. 2006, 11, 535–542. [Google Scholar] [CrossRef] [PubMed]
  33. Wang, X.W.; Wang, H.Z.; Wang, J.; Sun, R.F.; Wu, J.; Liu, S.Y.; Bai, Y.Q.; Mun, J.H.; Bancroft, I.; Cheng, F.; et al. The genome of the mesopolyploid crop species Brassica rapa. Nat. Genet. 2011, 43, 1035–1157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Edger, P.P.; Pires, J.C. Gene and genome duplications: The impact of dosage-sensitivity on the fate of nuclear genes. Chromosome Res. 2009, 17, 699–717. [Google Scholar] [CrossRef] [Green Version]
  35. Fang, L.; Cheng, F.; Wu, J.; Wang, X. The impact of genome triplication on tandem gene evolution in Brassica rapa. Front. Plant Sci. 2012, 3, 261. [Google Scholar] [CrossRef] [Green Version]
  36. Engel, W.; Hof, J.O.; Wolf, U. Gene duplication by polyploid evolution: The isoenzyme of the sorbitol dehydrogenase in herring- and salmon-like fishes (Isospondyli). Humangenetik 1970, 9, 157–163. [Google Scholar] [CrossRef]
  37. Ferris, S.D.; Whitt, G.S. Evolution of the differential regulation of duplicate genes after polyploidization. J. Mol. Evol. 1979, 12, 267–317. [Google Scholar] [CrossRef]
  38. Li, Y.; Li, D.; Xiao, Q.; Wang, H.; Wen, J.; Tu, J.; Shen, J.; Fu, T.; Yi, B. An in planta haploid induction system in Brassica napus. J. Integr. Plant Biol. 2022, 64, 1140–1144. [Google Scholar] [CrossRef]
  39. Kane, J.; Freeling, M.; Lyons, E. The Evolution of a High Copy Gene Array in Arabidopsis. J. Mol. Evol. 2010, 70, 531–544. [Google Scholar] [CrossRef] [Green Version]
  40. Wang, W.; Xiong, H.; Zhao, P.; Sun, M. DMP8 and 9 regulate HAP2/GCS1 trafficking for the timely acquisition of sperm fusion competence. Proc. Natl. Acad. Sci. USA 2022, 45, 119. [Google Scholar] [CrossRef]
  41. Philipp, C.; Maria, L.; Stefanie, S. Gamete fusion is facilitated by two sperm cell-expressed DUF679 membrane proteins. Nat. Plant 2019, 5, 253–257. [Google Scholar] [CrossRef]
  42. Song, Y.; Guo, X.; Wu, Y.; Liang, J.; Lin, R.; Yan, Z.; Wang, X. An optimized protocol for detecting guard cell specific gene expression by in situ RT-PCR in Brassica rapa. Hortic. Plant J. 2022, 8, 311–318. [Google Scholar] [CrossRef]
  43. Chen, R.; Chang, L.; Cai, X.; Wu, J.; Lin, R.; Song, Y.; Wang, X. Development of InDel markers for Brassica rapa based on a Highp-resolution melting curve. Hortic. Plant J. 2021, 7, 31–37. [Google Scholar] [CrossRef]
  44. Kaul, S.; Koo, H.L.; Jenkins, J.; Rizzo, M.; Rooney, T.; Tallon, L.J.; Feldblyum, T.; Nierman, W.; Benito, M.I.; Lin, X.Y.; et al. Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature 2000, 408, 796–815. [Google Scholar] [CrossRef] [Green Version]
  45. Hu, T.T.; Pattyn, P.; Bakker, E.G.; Cao, J.; Cheng, J.F.; Clark, R.M.; Fahlgren, N.; Fawcett, J.A.; Grimwood, J.; Gundlach, H.; et al. The Arabidopsis lyrata genome sequence and the basis of rapid genome size change. Nat. Genet. 2011, 43, 476–481. [Google Scholar] [CrossRef]
  46. Kliver, S.; Rayko, M.; Komissarov, A.; Bakin, E.; Zhernakova, D.; Prasad, K.; Rushworth, C.; Baskar, R.; Smetanin, D.; Schmutz, J.; et al. Assembly of the Boechera retrofracta genome and evolutionary analysis of apomixis-associated genes. Genes 2018, 9, 185. [Google Scholar] [CrossRef] [Green Version]
  47. Li, J.; Bi, C.; Tu, J.; Lu, Z. The complete mitochondrial genome sequence of Boechera stricta. Mitochondrial DNA B Resour. 2018, 3, 896–897. [Google Scholar] [CrossRef] [Green Version]
  48. Zhang, T.; Qiao, Q.; Novikova, P.Y.; Wang, Q.; Yue, J.; Guan, Y.; Ming, S.; Liu, T.; De, J.; Liu, Y.; et al. Genome of Crucihimalaya himalaica, a close relative of Arabidopsis, shows ecological adaptation to high altitude. Proc. Natl. Acad. Sci. USA 2019, 116, 7137–7146. [Google Scholar] [CrossRef] [Green Version]
  49. Slotte, T.; Hazzouri, K.M.; Agren, J.A.; Koenig, D.; Maumus, F.; Guo, Y.L.; Steige, K.; Platts, A.E.; Escobar, J.S.; Newman, L.K.; et al. The Capsella rubella genome and the genomic consequences of rapid mating system evolution. Nat. Genet. 2013, 45, 831–835. [Google Scholar] [CrossRef] [Green Version]
  50. Kagale, S.; Koh, C.; Nixon, J.; Bollina, V.; Clarke, W.E.; Tuteja, R.; Spillane, C.; Robinson, S.J.; Links, M.G.; Clarke, C.; et al. The emerging biofuel crop Camelina sativa retains a highly undifferentiated hexaploid genome structure. Nat. Commun. 2014, 5, 3706. [Google Scholar] [CrossRef]
  51. Akama, S.; Shimizu-Inatsugi, R.; Shimizu, K.K.; Sese, J. Genome-wide quantification of homeolog expression ratio revealed nonstochastic gene regulation in synthetic allopolyploid Arabidopsis. Nucleic Acids Res. 2014, 42, e46. [Google Scholar] [CrossRef] [PubMed]
  52. Gan, X.; Hay, A.; Kwantes, M.; Haberer, G.; Hallab, A.; Ioio, R.D.; Hofhuis, H.; Pieper, B.; Cartolano, M.; Neumann, U.; et al. The Cardamine hirsuta genome offers insight into the evolution of morphological diversity. Nat. Plants 2016, 2, 16167. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Haudry, A.; Platts, A.E.; Vello, E.; Hoen, D.R.; Leclercq, M.; Williamson, R.J.; Forczek, E.; Joly-Lopez, Z.; Steffen, J.G.; Hazzouri, K.M.; et al. An atlas of over 90,000 conserved noncoding sequences provides insight into crucifer regulatory regions. Nat. Genet. 2013, 45, 891–898. [Google Scholar] [CrossRef] [PubMed]
  54. Zhang, L.; Cai, X.; Wu, J.; Liu, M.; Grob, S.; Cheng, F.; Liang, J.; Cai, C.; Liu, Z.; Liu, B.; et al. Erratum: Author Correction: Improved Brassica rapa reference genome by single-molecule sequencing and chromosome conformation capture technologies. Hortic. Res. 2019, 6, 124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Liu, S.; Liu, Y.; Yang, X.; Tong, C.; Edwards, D.; Parkin, I.A.; Zhao, M.; Ma, J.; Yu, J.; Huang, S.; et al. The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. Nat. Commun. 2014, 5, 3930. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Paritosh, K.; Pradhan, A.K.; Pental, D. A highly contiguous genome assembly of Brassica nigra (BB) and revised nomenclature for the pseudochromosomes. BMC Genom. 2020, 21, 887. [Google Scholar] [CrossRef]
  57. Yang, J.; Liu, D.; Wang, X.; Ji, C.; Cheng, F.; Liu, B.; Hu, Z.; Chen, S.; Pental, D.; Ju, Y.; et al. The genome sequence of allopolyploid Brassica juncea and analysis of differential homoeolog gene expression influencing selection. Nat. Genet. 2016, 48, 1225–1232. [Google Scholar] [CrossRef]
  58. Rousseau-Gueutin, M.; Belser, C.; Da Silva, C.; Richard, G.; Istace, B.; Cruaud, C.; Falentin, C.; Boideau, F.; Boutte, J.; Delourme, R.; et al. Long-read assembly of the Brassica napus reference genome Darmor-bzh. Gigascience 2020, 9, giaa137. [Google Scholar] [CrossRef]
  59. Song, X.; Wei, Y.; Xiao, D.; Gong, K.; Sun, P.; Ren, Y.; Yuan, J.; Wu, T.; Yang, Q.; Li, X.; et al. Brassica carinata genome characterization clarifies U’s triangle model of evolution and polyploidy in Brassica. Plant Physiol. 2021, 186, 388–406. [Google Scholar] [CrossRef]
  60. Kitashiba, H.; Li, F.; Hirakawa, H.; Kawanabe, T.; Zou, Z.; Hasegawa, Y.; Tonosaki, K.; Shirasawa, S.; Fukushima, A.; Yokoi, S.; et al. Draft sequences of the radish (Raphanus sativus L.) genome. DNA Res. 2014, 21, 481–490. [Google Scholar] [CrossRef]
  61. Kang, M.; Wu, H.; Yang, Q.; Huang, L.; Hu, Q.; Ma, T.; Li, Z.; Liu, J. A chromosome-scale genome assembly of Isatis indigotica, an important medicinal plant used in traditional Chinese medicine: An Isatis genome. Hortic. Res. 2020, 7, 18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Willing, E.M.; Rawat, V.; Mandakova, T.; Maumus, F.; James, G.V.; Nordstrom, K.J.; Becker, C.; Warthmann, N.; Chica, C.; Szarzynska, B.; et al. Genome expansion of Arabis alpina linked with retrotransposition and reduced symmetric DNA methylation. Nat. Plants 2015, 1, 14023. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Dorn, K.M.; Fankhauser, J.D.; Wyse, D.L.; Marks, M.D. A draft genome of field pennycress (Thlaspi arvense) provides tools for the domestication of a new winter biofuel crop. DNA Res. 2015, 22, 121–131. [Google Scholar] [CrossRef] [PubMed]
  64. Yang, R.; Jarvis, D.E.; Chen, H.; Beilstein, M.A.; Grimwood, J.; Jenkins, J.; Shu, S.; Prochnik, S.; Xin, M.; Ma, C.; et al. The Reference genome of the halophytic plant Eutrema salsugineum. Front. Plant Sci. 2013, 4, 46. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Synteny and phylogenetic evolution analysis of DMP9 homologous genes in 26 Brassicaceae species. (a) An unrooted neighbor-joining phylogenetic tree of 54 DMP9 homologous genes from 26 Brassicaceae species and maize. DMP9 homologous genes in Brassicaceae species clustered in two groups: Group 1 (green branch) and Group 2 (red branch). The blue background represents the gene of A block. The green background represents the gene of S block. The red background represents the gene of R block. The gray background represents the gene that has no syntenic relationship. The ZmDMP gene of maize was used as the outgroup. Numbers on branches are bootstrap support percentages based on 1000 replicates. The phylogenetic tree was constructed with MEGA7.0 using the maximum likelihood algorithm with default parameters [24]. (b) Schematic representation of synteny of DMP9 homologous genes among 26 Brassicaceae species. The background colors represent different genomic blocks; green is the S block, red is the R block, and gray is the non-syntenic genes. The white box with a full line represents a syntenic gene. The white box with a dotted line represents no gene.
Figure 1. Synteny and phylogenetic evolution analysis of DMP9 homologous genes in 26 Brassicaceae species. (a) An unrooted neighbor-joining phylogenetic tree of 54 DMP9 homologous genes from 26 Brassicaceae species and maize. DMP9 homologous genes in Brassicaceae species clustered in two groups: Group 1 (green branch) and Group 2 (red branch). The blue background represents the gene of A block. The green background represents the gene of S block. The red background represents the gene of R block. The gray background represents the gene that has no syntenic relationship. The ZmDMP gene of maize was used as the outgroup. Numbers on branches are bootstrap support percentages based on 1000 replicates. The phylogenetic tree was constructed with MEGA7.0 using the maximum likelihood algorithm with default parameters [24]. (b) Schematic representation of synteny of DMP9 homologous genes among 26 Brassicaceae species. The background colors represent different genomic blocks; green is the S block, red is the R block, and gray is the non-syntenic genes. The white box with a full line represents a syntenic gene. The white box with a dotted line represents no gene.
Horticulturae 08 01095 g001aHorticulturae 08 01095 g001b
Figure 2. Evolution of DMP9 homologous genes in B. rapa pan-genome. (a) Phylogeny relationships of DMP9 genes in B. rapa pan-genome. An unrooted neighbor-joining phylogenetic tree of 41 DMPs from B. rapa pan-genome. DMPs belonged to two branches. The green background represents the S block. The red background represents the R block. The gray background represents the non-syntenic gene of the A block. The AtDMP9 gene was used as the outgroup. (b) The syntenic relationship of DMP9 genes in B. rapa pan-genome. The green background represents the S block. The red background represents the R block. The white box with a full line represents the syntenic gene. The white box with the dotted line represents no gene. (c) Gene structures of syntenic DMP9 genes in the B. rapa pan-genome. The green boxes represent CDS, and the yellow boxes indicate an untranslated region.
Figure 2. Evolution of DMP9 homologous genes in B. rapa pan-genome. (a) Phylogeny relationships of DMP9 genes in B. rapa pan-genome. An unrooted neighbor-joining phylogenetic tree of 41 DMPs from B. rapa pan-genome. DMPs belonged to two branches. The green background represents the S block. The red background represents the R block. The gray background represents the non-syntenic gene of the A block. The AtDMP9 gene was used as the outgroup. (b) The syntenic relationship of DMP9 genes in B. rapa pan-genome. The green background represents the S block. The red background represents the R block. The white box with a full line represents the syntenic gene. The white box with the dotted line represents no gene. (c) Gene structures of syntenic DMP9 genes in the B. rapa pan-genome. The green boxes represent CDS, and the yellow boxes indicate an untranslated region.
Horticulturae 08 01095 g002
Figure 3. (a) The syntenic relationship of A. thaliana and B. rapa “Chiifu” and “B9008”. (b) Alignment of amino acid sequences encoded by DMP9 genes of B. oleracea and two B. rapa species. The green background represents the DMP9 homologous genes of the LF subgenome. The red background represents the DMP9 homologous genes of the MF1 subgenome. (c) Relative expression profiles of DMP9 homologous genes in B. rapa “Chiifu”, B. rapa “B9008”, and B. oleracea “JZS” at the different pollen developmental stages. The error bars represent the standard deviations of three independent biological repeats. The asterisks represent significant differences via t-tests (* p ≤ 0.05; ** p ≤ 0.01; **** p ≤ 0.0001). Pollen at different developmental stages was observed by DAPI staining.
Figure 3. (a) The syntenic relationship of A. thaliana and B. rapa “Chiifu” and “B9008”. (b) Alignment of amino acid sequences encoded by DMP9 genes of B. oleracea and two B. rapa species. The green background represents the DMP9 homologous genes of the LF subgenome. The red background represents the DMP9 homologous genes of the MF1 subgenome. (c) Relative expression profiles of DMP9 homologous genes in B. rapa “Chiifu”, B. rapa “B9008”, and B. oleracea “JZS” at the different pollen developmental stages. The error bars represent the standard deviations of three independent biological repeats. The asterisks represent significant differences via t-tests (* p ≤ 0.05; ** p ≤ 0.01; **** p ≤ 0.0001). Pollen at different developmental stages was observed by DAPI staining.
Horticulturae 08 01095 g003
Table 1. List of DMP9 homologous genes of 18 B. rapa varieties and subspecies.
Table 1. List of DMP9 homologous genes of 18 B. rapa varieties and subspecies.
SampleNameDescriptionE-ValueIdentityGenomic BlockAccession
CCAssp. pekinensisChinese cabbage080.30%SCCADMP9.1
068.30%SCCADMP9.2
085.8%RCCADMP9.3
CCBssp. pekinensisChinese cabbage091.00%SCCBDMP9.1
085.60%RCCBDMP9.2
Chiifu #ssp. pekinensisChinese cabbage091.00%SBraDMP9.1
085.80%RBraDMP9.2
CXAssp. parachinensisCaixin090.30%SCXADMP9.1
085.60%RCXADMP9.2
CXBssp. parachinensisCaixin090.30%SCXBDMP9.1
085.80%RCXBDMP9.2
OIAssp. oleiferaOil seeds091.00%SOIADMP9.1
085.80%ROIADMP9.2
OIBssp. oleiferaRapid cycling090.30%SOIBDMP9.1
090.30%SOIBDMP9.2
085.60%ROIBDMP9.3
OICssp. oleiferaRapid cycling090.30%SOICDMP9.1
085.60%ROICDMP9.2
Z1ssp. oleiferaSarson type 090.30%SZ1DMP9.1
085.60%RZ1DMP9.2
TUAssp. rapaTurnip090.70%STUADMP9.1
085.80%RTUADMP9.2
TUEssp. rapaTurnip091.00%STUEDMP9.1
085.60%RTUEDMP9.2
PCAssp. chinensisPak choi090.30%SPCADMP9.1
085.80%RPCADMP9.2
PCBssp. chinensisPak choi090.30%SPCBDMP9.1
085.70%RPCBDMP9.2
WTCssp. narinosaWutacai091.00%SWTCDMP9.1
085.60%RWTCDMP9.2
TCAssp. chinensis var.tai-tsaiTaicai091.00%STCADMP9.1
085.80%RTCADMP9.2
BROBroccolietoBroccolleto068.30%SBroDMP9.1
090.70%SBroDMP9.2
085.80%RBroDMP9.3
MIZssp. nipposinicaMizuna091.00%SMIZDMP9.1
085.70%RMIZDMP9.2
TBA (WLD)unknownFrom Tibet, China090.70%STBADMP9.1
085.80%RTBADMP9.2
085.60%unknownTBADMP9.3 *
085.60%unknownTBADMP9.4 *
Notes: # The accession was used as B. rapa reference genome. * indicates non-syntenic genes.
Table 2. DMP9 homologous gene information of B. oleracea “JZS”, B. rapa “Chiifu”, and B. rapa “B9008”.
Table 2. DMP9 homologous gene information of B. oleracea “JZS”, B. rapa “Chiifu”, and B. rapa “B9008”.
Gene NameChrGene IDStartEndGene BlockSpeciesVarieties/
Subspecies
AtDMP9A05At5g396501587519915876116SA. thalianaA. thaliana
BraDMP9.1A04BraA04g012030.3C95241839524908SB. rapaChiifu
BraDMP9.2A03BraA03g003970.3C16916981692426RB. rapaChiifu
CCADMP9.1 *A04A04p12560.1_BraCCA89994899000040SB. rapaB9008
CCADMP9.2 *A04A04p12570.1_BraCCA90005099001147SB. rapaB9008
CCADMP9.3A03A03p04040.1_BraCCA17469051747628RB. rapaB9008
BolDMP9.1C04BolC04g044930.2J4636183846362575SB. oleraceaJZS
BolDMP9.2C03BolC03g004320.2J21184352119158RB. oleraceaJZS
Note: * indicates tandem duplication.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, T.; Liang, J.; Cai, X.; Zhang, L.; Wu, J.; Wang, X. Analysis of Structure Variations and Expression Characteristics of DMP8 and DMP9 Genes in Brassicaceae. Horticulturae 2022, 8, 1095. https://doi.org/10.3390/horticulturae8111095

AMA Style

Zhang T, Liang J, Cai X, Zhang L, Wu J, Wang X. Analysis of Structure Variations and Expression Characteristics of DMP8 and DMP9 Genes in Brassicaceae. Horticulturae. 2022; 8(11):1095. https://doi.org/10.3390/horticulturae8111095

Chicago/Turabian Style

Zhang, Tingting, Jianli Liang, Xu Cai, Lei Zhang, Jian Wu, and Xiaowu Wang. 2022. "Analysis of Structure Variations and Expression Characteristics of DMP8 and DMP9 Genes in Brassicaceae" Horticulturae 8, no. 11: 1095. https://doi.org/10.3390/horticulturae8111095

APA Style

Zhang, T., Liang, J., Cai, X., Zhang, L., Wu, J., & Wang, X. (2022). Analysis of Structure Variations and Expression Characteristics of DMP8 and DMP9 Genes in Brassicaceae. Horticulturae, 8(11), 1095. https://doi.org/10.3390/horticulturae8111095

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