Dan-Dan Liu1*†
Da-Ru Wang2†
Xuan-Yu Yang1
Chang-Hui Zhao1Shao-Hua Li1Guang-Li Sha3Rui-Fen Zhang3Hong-Juan Ge3Xian-Song Tong4
Chun-Xiang You2*- 1College of Agriculture, Yunnan University, Kunming, Yunnan, China
- 2National Key Laboratory of Crop Biology, National Research Center for Apple Engineering and Technology, College of Horticulture Science and Engineering, Shandong Agricultural University, Tai-An, Shandong, China
- 3Qingdao Academy of Agricultural Sciences, Qingdao, Shandong, China
- 4Fu-ning Popularizing Agricultural Techniques Center, Fu-ning, Yunnan, China
Apomixis is the asexual reproduction through seeds that leads to the production of genetically uniform progeny. It has become an important tool in plant breeding because it facilitates the retention of genotypes with desirable traits and allows seeds to be obtained directly from mother plants. Apomixis is rare in most economically important crops, but it occurs in some Malus species. Here, the apomictic characteristics of Malus were examined using four apomictic and two sexually reproducing Malus plants. Results from transcriptome analysis showed that plant hormone signal transduction was the main factor affecting apomictic reproductive development. Four of the apomictic Malus plants examined were triploid, and pollen was either absent or present in very low densities in the stamen. Variation in the presence of pollen was associated with variation in the apomictic percentage; specifically, pollen was absent in the stamens of tea crabapple plants with the highest apomictic percentage. Furthermore, pollen mother cells failed to progress normally into meiosis and pollen mitosis, a trait mostly observed in apomictic Malus plants. The expression levels of meiosis-related genes were upregulated in apomictic plants. Our findings indicate that our simple method of detecting pollen abortion could be used to identify apple plants that are capable of apomictic reproduction.
1 Introduction
Apomictic seeds are the products of asexual reproduction and give rise to maternal clones; meiosis and double fertilization are not required for the production of apomictic seeds (Hand and Koltunow, 2014). Apomixis results in the production of genetically identical offspring. It has become an important breeding tool in agriculture because it permits the generation of clones of mother plants as well as the retention of crop genotypes with desirable features. Studies of apomixis have been conducted for over 150 years. Female and male heteromorphic plants of Achornea ilicifolia were shown to produce seeds without pollination by Smith in 1841 (Smith, 1841). Apomixis has been documented to occur in over 400 species in 61 families, including Poaceae, Rosaceae, and Asteraceae. Apomictic reproduction is not widespread in crops, but it is common in horticultural crops, such as citrus, crabapple, walnut, mango, and Chinese chive (Cosmulescu et al., 2012; Yamashita et al., 2012; Liu et al., 2014; Kuhn et al., 2017; Wu et al., 2021).
In apomictic plants, female gametes are formed without meiosis (apomeiosis) and the embryo develops without fertilization; the formation of functional endosperm is also mediated by several developmental adaptations in apomictic plants (Koltunow and Grossniklaus, 2003; Hao and Liu, 2011). Apomixis can be classified as obligate or facultative depending on whether plants can reproduce sexually. Apomixis can also occur via diplospory, apospory, and adventitious embryo reproduction depending on the developmental origin of apomictically derived embryos. Diplospory and apospory are both forms of gametophytic apomixis, wherein unreduced cells give rise to a megagametophyte via mitosis, and adventitious embryo reproduction is a form of sporophytic apomixis, wherein unreduced cells give rise directly to an embryo (Tucker and Koltunow, 2009). Apomixis, especially facultative apomixis, is common in apple trees in the family Rosaceae, and natural pollination often leads to the production of a small number of sexual hybrids (Liu et al., 2014).
Several apomixis-related genes involved in the regulation of embryo sac genes have been identified in Arabidopsis thaliana, petunia, maize, and other plants in recent years. For example, SWITCH1/DYAD promotes female gametophyte formation by causing cells to bypass meiosis (Mercier et al., 2001; Mercier et al., 2003). SOMATIC EMBRYOGENESIS RECEPTOR KINASE 1 (SERK1), LEAF COTYLEDON 1 (LEC1), LEC2, and BABYBOOM (BBM) have been shown to induce somatic embryogenesis in plants (Boutilier et al., 2002; Ranganath, 2004; Santos and Aragão, 2009). In Pliosella, which is a subgenus of Hieracium, apomixis and the autonomous development of seeds are regulated by two independent dominant loci, LOSS OF APOMEIOSIS (LOA) and LOSS OF PARTHENOGENESIS (LOP), respectively (Catanach et al., 2006; Okada et al., 2007; Koltunow et al., 2011). CitRWP, which contains an RWP-RK domain, regulates polyembryonic characters and various morphological aspects of the plant body during citrus embryogenesis (Wang et al., 2017; Wang et al., 2022). In Zanthoxylum bungeanum, the expression of AGL11 is correlated with nucellar embryo development, and the ectopic expression of this gene can result in abnormal flower development and apomixis phenotypes in Arabidopsis (Fei et al., 2021). Various studies have attempted to genetically engineer plants with gametophytic apomixis. The MiMe (Mitosis instead of Meiosis) system was first created in Arabidopsis via the simultaneous mutation of three key meiotic genes (SPO11-1, REC8, and OSD1). MiMe, coupled with the expression of BBM1 (a parthenogenesis gene) in the egg cell, can also induce the production of clonal progeny in hybrid rice and mediate the retention of genome-wide parental heterozygosity (Khanday et al., 2019). Nevertheless, there are still several shortcomings associated with synthetic apomixis.
Apomictic reproduction has been documented in a total of 10 Malus species (Shu, 1999; Li, 2001). Six of these species are native to China, and these species comprise a group containing 23 species (e.g., Malus hupehensis, Malus sikkimensis, Malus rockii, Malus sieboldii, Malus platycarpa, Malus sargentii, Malus xiaojinensis, Malus lancifolia, and Malus corolla). Most angiospermous show facultative apomixis because they can set seeds both sexually and apomictically depending on their genotype and environmental conditions (Carman et al., 2011; Liu et al., 2014). Apospory is widespread among polyploid species in the genus Malus (e.g., tea crabapple) (Schmidt, 1977). In apospory, the embryo sac develops directly from ovule cells and the embryo develops from the unreduced egg cells (Liu et al., 2014). Tea crabapple (M. hupehensis var. pingyiensis) is facultatively apomictic, and both parthenogenesis and apomeiosis have been documented in this species. The apomictic capacity of this species is controlled by separate genetic loci, and these loci affect both apomeiosis and parthenogenesis (Liu et al., 2014).
Much research has focused on developing approaches to facilitate the production of apomictic seeds and produce seedlings with uniform characteristics from apomictic crabapple rootstocks. Although much progress has been made in the genetic engineering of Malus species, there are still many challenges that need to be overcome for the genetic improvement of apomictic Malus species. A deep understanding of the mechanism of apomixis is essential for ensuring that genetic modifications have their desired effect. Tea crabapple, which is referred to as PYTC in this study, is a triploid plant that is facultatively apomictic; it is often used as a rootstock in production. Previous analysis of the development and genetic characteristics of apomixis in PYTC has shown that the bypassing of pollination and fertilization is mainly achieved through a delay in embryo sac development (female late-in-date), and this mediates the production of asexual seeds (Liu et al., 2014). Here, we evaluated the apomictic capacity of PYTC and then characterized several genes involved in facultative apomixis by conducting a comparative transcriptomic analysis of the flowers of a hybrid population. The reference transcriptome generated in our study will aid future studies aimed at elucidating the regulatory networks and pathways involved in apomixis in Malus. We also evaluated the apomictic capacity in Malus; most of the material examined was triploid, and pollen development was abnormal. Overall, our study provides a simple method that could be used to identify apomictic apple rootstocks.
2 Materials and methods
2.1 Plant materials
Apomictic crabapple species examined in our study included the triploid tea crabapple (M. hupehensis Redh. var. pingyiensis, PYTC, 2n=3X=51, facultatively apomictic), triploid M. sargentii (SJHT, 2n=3X=51, facultatively apomictic), M. sikkimensis (SYHT, 2n=3X=51, facultatively apomictic), and M. sieboldii (XJHT, 2n=3X=51, facultatively apomictic); both Malus domestica “Jonagold” (QNJ, 2n=3X=51, triploid but can reproduce sexually) and Malus baccata (SDZ, 2n=3X=34, diploid and can reproduce sexually) were used as controls. All plants were grown in the experimental orchard of Shandong Agricultural University and Yunnan University.
2.2 Evaluation of apomixis and estimation of ploidy
The apomictic capacity of the crabapple species examined was evaluated following the methods of Liu et al. (2014). To determine the parthenogenesis percentage, pistils were decapitated in the spring before the flowers opened. The parthenogenesis percentage was calculated as the number of parthenocarpic fruits that could produce seeds in the absence of fertilization relative to the number of decapitated flowers. To determine the apomeiosis percentage, asexual seeds were harvested from parthenocarpic fruits. The apomeiosis percentage was calculated as the number of triploid seeds relative to the total number of seeds tested.
Ploidy was determined using flow cytometry assays, which were conducted using the Partec CyStain UV Precise T reagent kit (PARTEC, Cod. 05-5003) per the manufacturer’s instructions. First, young leaves or embryos from mature seeds were cut into smaller pieces with a razor blade in nuclei extraction buffer. Staining buffer was then added with the one-fifth volume of nuclei extraction buffer following the manufacturer’s instructions. Next, debris was removed, and nuclei were collected by filtering the nuclear suspensions two times through a 20-µm nylon mesh. An arc lamp-based flow cytometer (PARTEC PA, Germany) was used to measure the fluorescence intensity of the nuclei. The fluorescence intensity of 30,000 nuclei per cytogram was measured. At least three different plants for each line were measured, and each sample was analyzed in triplicate.
2.3 Morphological observations and pollen viability experiment
Morphological observations of the anthers were made by taking photographs with a Nikon SMZ18 microscope, and a scalpel was used to dissect the anthers after the photographs were taken.
The anthers were peeled, cut open, and placed on a slide, and forceps were used to extract the pollen grains. Two to three drops of I2–KI solution were added; after mixing thoroughly, the slide was covered with a coverslip for staining. One of the potassium iodide was then dissolved in a small amount of distilled water. After the potassium iodide and distilled water were fully dissolved, 0.5 g of iodine was added. The mixture was then shaken to dissolve all the components; it was then diluted to 150 ml and stored in a brown glass bottle.
2.4 Histological analysis
Anthers were fixed in formaldehyde–alcohol–acetic acid fixative (90 ml of 70% alcohol + 5 ml of formalin + 5 ml of glacial acetic acid). Paraffin sectioning was conducted following the method described in a previous study to characterize the histology of the fixed samples (Yang et al., 2013).
2.5 RNA sequencing and data analysis
The Trizol kit (Invitrogen, USA) was used to extract total RNA from all the plant materials, and the RNA sequencing (RNA-seq) library was constructed using the RNA obtained. After sequencing, the raw data were filtered to obtain high-quality clean reads. Differentially expressed genes (DEGs) were identified using the DEGseq R package with the following criteria: false discovery rate (FDR) < 0.001 and |log2(fold change (FC))| ≥ 1. The GOSeq R software package (p-value ≤ 0.05) was used to conduct a Gene Ontology (GO) analysis of the DEGs.
2.6 qRT-PCR assays
An ABI 7500 Real-time PCR system was used to conduct qRT-PCR analyses of the extracted RNA following the procedures described in a previous study (Wang et al., 2022). The 2−ΔΔCt method was used to determine expression levels. Three biological replicates were performed for each sample. Supplementary Table 1 shows the primers that were used in qRT-PCR analyses.
2.7 Statistical analysis
Data processing system software was used to analyze the data (Wang et al., 2022). The significance of differences between groups was evaluated using a one-way analysis of variance (ANOVA), followed by the Tukey–Kramer test; the threshold for statistical significance was p-value < 0.05.
3 Results
3.1 The apomictic capacity of Malus plants varied
The apomictic percentages of Malus plants were evaluated from 2019 to 2021. The lowest apomictic percentages were observed for QNJ (1.3%) and SDZ (2.8%). The apomictic percentages of PYTC, SJHT, SYHT, and XJHT were 94.7%, 86.1%, 41.5%, and 52.7%, respectively (Figure 1A). The apomictic percentage of PYTC was the highest among all apple species. Additionally, ploidy analyses revealed that SDZ was diploid, and PYTC, SJHT, SYHT, XJHT, and QNJ were triploid (Figure 1B), suggesting that all detected apomictic Malus plants were polyploid. Nonetheless, the sexual line QNJ was also found to be polyploid, suggesting that polyploidy may not be the only driving factor leading to apomixis.
Figure 1 Apomictic percentage and ploidy analysis of different Malus plants. (A) Apomictic percentage analysis. (B) Ploidy analysis of apomictic Malus plants.
3.2 DEGs in Malus with high and low apomictic percentages
Plants in the genus Malus are cross-pollinated, and the progeny show high heterozygosity. Hybrid progeny were obtained via artificial pollination, and these hybrid progeny populations were used to clarify the mechanism underlying apomixis in Malus. Apomictic reproduction is performed facultatively in tea crabapple (PYTC) in its natural environment, and PYTC exhibits the highest apomictic percentage among all Malus plants. Sexually reproducing hybrids were obtained via artificial pollination of PYTC pistils with pollen from another wild apple species (Malus pumila Mill. var. “Maypole”). Floral development comprises six stages (Liu et al., 2014). The ovules of hybrids derived from tea crabapple and Maypole were collected at stage 2, which is the pre-meiotic development stage, from two populations differing in parthenogenesis percentages: P1 (21.3 ± 1.6%; low parthenogenesis percentage) and P2 (76.7 ± 5.3%; high parthenogenesis percentage) (Liu et al., 2014).
To identify apomixis-related genes in Malus, RNA-seq analysis was conducted on P1 and P2. Overall, a total of 24,126,591 reads were obtained in P1, and a total of 22,785,112 reads were obtained in P2. In P1, the GC content was 47.11%, the total number of clean reads obtained was 22,785,112, and the total number of bases was 7,237,977,400; in P2, the GC content was 46.89%, and the total number of bases was 6,835,533,600 (Supplementary Table 2). A total of 3,154 DEGs were identified using the following criteria: FDR < 0.001 and |log2(FC) ≥ 1|. The expression of 1,459 and 1,695 DEGs was downregulated and upregulated, respectively (Figure 2A). A volcano plot was constructed to visualize variation in the number of DEGs among comparison groups. The differential expression of genes was higher the closer they were located at the ends of the volcano plot (Figure 2B). Cluster analysis was conducted on the DEGs (Figure 2C).
Figure 2 Preliminary analysis of transcriptomic data from P1 and P2. (A) Statistical analysis of downregulated and upregulated DEGs in the P1 vs. P2 comparison group. (B) Volcano plots of DEGs in the P1 vs. P2 comparison group. (C) Heatmap of the expression levels of DEGs in the P1 vs. P2 comparison group.
3.3 Functional classification of DEGs
GO analysis was conducted to characterize the functions of DEGs in P1 and P2 (Figure 3A). DEGs were assigned GO terms in three categories: cellular component (CC), molecular function (MF), and biological process (BP). In the CC category, DEGs were mainly enriched in the following GO terms: “cell part,” “cell,” and “organelle.” In the MF category, DEGs were mainly enriched in the following GO terms: “catalytic activity” and “binding.” In the BP category, most DEGs were enriched in the following GO terms: “metabolic process,” “cellular process,” and “single-organism process.” Most of these processes are involved in pollen development.
Figure 3 GO and KEGG analysis of DEGs in the P1 vs. P2 comparison group. (A) GO analysis of DEGs in the P1 vs. P2 comparison group. (B) KEGG analysis of DEGs in the P1 vs. P2 comparison group.
Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was conducted on the DEGs to further clarify their functions. The main pathways were divided into the following categories: organismal systems, genetic information processing, environmental information processing, metabolism, and cellular processes (Figure 3B). In genetic information processing, the most enriched pathways were “ribosome,” “DNA replication,” and “ribosome biogenesis in eukaryotes.” The most enriched pathway in environmental information processing was “plant hormone signal transduction.” The most enriched pathway in metabolism was “purine metabolism.” These findings suggest that hormone signal transduction played a role in regulating apomixis.
3.4 Pollen development was abnormal in apomictic Malus plants
All detected apomictic Malus plants were polyploid. We studied the pollen development process and measured pollen viability in apomictic Malus plants to determine whether pollen development was affected by polyploidy. No pronounced differences in anther morphology were observed between apomictic and sexually reproducing Malus plants. Pollen was absent in the anthers of PYTC and SJHT, which were the two Malus species with the highest apomictic percentage (Figure 4A). In addition, the abundance of pollen in the apomictic Malus species SYHT and XJHT was less than that in the sexually reproducing species QNJ and SDZ. The activity of pollen grains in SYHT, XJHT, QNJ, and SDZ was evaluated via I2–KI staining (Figure 4B). The amount of pollen was significantly lower in each field of SYHT and XJHT than in QNJ and SDZ. Most of the pollen grains in QNJ and SDZ were stained blue or blue-black, suggesting that the degree of starch accumulation was high in QNJ and SDZ. By contrast, most of the pollen grains in SYHT and XJHT were stained yellow or yellow-brown, and only a few were stained blue or blue-black. This suggests that starch accumulated in only a few of the pollen grains and thus that most of the pollen grains were not viable pollen.
Figure 4 Morphological observations of anthers and evaluation of the activity of pollen grains in apple plants. (A) Morphological structure of the anthers of PYTC, SJHT, SYHT, XJHT, QNJ, and SDZ and the presence of pollen after dissection. Bar, 500 μm. (B) Pollen viability assays for SYHT, XJHT, QNJ, and SDZ. Bar, 200 μm.
3.5 The expression patterns of genes involved in pollen development are altered in apomictic Malus plants
We identified genes involved in pollen development using the RNA-seq data to clarify the relationship between pollen development and apomixis (Figure 5A). We identified a total of 31 significant DEGs in the P1 vs. P2 comparison group. Next, qRT-PCR analyses were conducted to measure the expression levels of six genes (MD16G1014300, MD07G1233000, MD01G1236300, MD13G1013400, MD10G1337500, and MD14G1229400) in QNJ, PYTC, and SJHT (Figure 5B). The homologous gene of MD07G1233000 (ETG1) in Arabidopsis plays a key role in the establishment of sister chromatid cohesion in mitosis and meiosis. MD01G1236300 and MD13G1013400 (ATA1) may regulate the development of the tapetum, which provides nutrients for the development of pollen. MD16G1014300 regulates the development of flowers, MD10G1337500 encodes a heavy metal transport/detoxification superfamily protein, and MD14G1229400 plays a key role in the development of the cytoplasm of pollen. The expression of all six genes was lower in QNJ than in PYTC and SJHT, and differences in the expression of this gene among Malus species were consistent with differences in the apomictic percentage. These findings suggest that there is a relationship between apomixis and pollen development.
Figure 5 Analysis of DEGs involved in pollen development in the P1 vs. P2 comparison group. (A) Cluster heatmap of genes involved in pollen development in the P1 vs. P2 comparison group. (B) The expression patterns of six of the DEGs in (A) were analyzed by qRT-PCR. Groups with different lowercase letters are significantly different according to ANOVA, followed by the Tukey–Kramer test (p < 0.05). Data are the mean ± standard deviation of three independent replicates.
3.6 Meiotic abnormalities result in the premature termination of pollen development in apomictic Malus plants
Paraffin sectioning was used to identify the stage at which pollen development ceases and determine why pollen development was abnormal in apomictic Malus plants. Pollen development in PYTC and SJHT ceased at the meiosis stage of pollen mother cells (Figures 6F, I), and defects in tapetal development were the ultimate cause of pollen abortion. No mature pollen grains were detected in the anthers at the final maturation stage (Figures 6H, L). Pollen development was normal in the sexually reproducing QNJ (Figures 6A–D).
Figure 6 Paraffin sectioning and quantification of the expression of genes involved in meiosis at different stages of pollen development in the anthers of QNJ, PYTC, and SJHT. (A–D) QNJ anthers at different developmental stages. (A) Microspore mother cell stage; (B) early uninucleate stage; (C) late uninucleate stage; and (D) mature stage. (E–H) PYTC anthers at different developmental stages. (E, F) Microspore mother cell stage; (G) callose degradation; (H) and mature stage. (I–L) SJHT anthers at different developmental stages. (I) Microspore mother cell stage; (J, K) callose degradation; (K, L) mature stage. (M) The expression levels of meiosis-related genes were determined by qRT-PCR analyses. Groups with different lowercase letters are significantly different according to ANOVA, followed by the Tukey–Kramer test (p < 0.05). Data are the mean ± standard deviation of three independent replicates.
Pollen development did not proceed normally in PYTC and SJHT because of the failed meiosis of microspore mother cells. To identify candidate genes that contribute to the observed meiotic abnormalities in PYTC and SJHT, we measured the expression levels of some of the meiosis-related genes (Supplementary Figure 1). Genes MD10G1311800 and MD05G1336300 may be homologs of Arabidopsis ATK1, a gene that plays a key role in the assembly of the spindle during male meiosis and may code for a kinesin, while MD15G1344200 may correspond to Arabidopsis CYCA1;1, a gene that may play a role in the mitotic cell cycle phase transition, and MD15G1140600 negatively regulates cell biosynthesis. We also found that the expression of these four genes was higher in PYTC and SJHT than in QNJ (Figure 6M). The differential expression of these genes might contribute to the abnormal meiosis of PYTC and SJHT.
4 Discussion
Generally, plants undergo sexual reproduction to form seeds. However, some plants have evolved alternative modes of reproduction that allow them to bypass the meiotic and double fertilization processes of sexual reproduction before seed formation. This process of asexual seed production is called apomixis (Bicknell and Koltunow, 2004). Several studies have attempted to generate apomictic plants via genetic approaches, but no breakthroughs in the generation of apomictic plants have yet been made (Leblanc et al., 2009; Rodriguez-Leal and Vielle-Calzada, 2012). One of the major challenges is that apomixis is a highly complex trait, and apomeiosis, parthenogenesis, and developmental adaptation are required for the formation of the functional endosperm. Several genes are involved in regulating each of these processes (Koltunow and Grossniklaus, 2003). Differences in the expression of plant hormone signaling genes have been observed between sexually reproducing and apomictic Arabidopsis and Paspalum notatum (Figueiredo and Köhler, 2018; Ortiz et al., 2019). In this study, transcriptome analysis of the hybrid progeny of PYTC revealed that hormone signal transduction was the main factor affecting apomictic reproductive development, and additional studies are needed to clarify the relative importance of apomixis in Malus plants.
Nearly all asexually reproducing plants (including both herbaceous and woody plants) identified to date are polyploid (Sax, 1959). All the Malus plants in our study were polyploid, and the apomictic percentage of these plants varied, indicating that there might be a relationship between apomixis and polyploidy (Figure 1). PYTC, SJHT, SYHT, and XJHT were all triploid and had high apomictic percentages; SDZ, which does not reproduce via apomixis, was diploid. Although QNJ was triploid, it is not capable of apomictic reproduction, suggesting that genes related to apomixis have been lost in this apple species.
Pollen was absent in the anthers of PYTC and SJHT (Figure 4A). The amount of pollen in the anthers of SYHT and XJHT was low (Figure 4B). Analysis of the expression of genes involved in pollen development revealed that the expression of these genes in PYTC and SJHT significantly differed from their expression patterns in QNJ (Figure 5B). Pollen development is a complex process, and abnormalities in meiosis, the tapetal layer, and pollen sac formation can result in abortive or abnormal development (Jeong et al., 2014; Wu et al., 2015; Somaratne et al., 2017; Zhao et al., 2018; Yang et al., 2019). In apomictic Malus plants, pollen development ceased at the meiosis stage, which suggests that apomeiosis failed in pollen development. Thus, pollen arrest could be used to identify apomictic plants; however, polyploidy does not necessarily result in pollen abortion. In Melastomataceae species, apomictic plants with no viable pollen or with pollen with low viability do not receive visits from pollinators and likely produce strictly apomictic fruits, whereas apomictic and sexually reproducing plants with high pollen viability might produce fruits and seeds through both sexual and apomictic reproduction (Maia et al., 2016). Pollen viability levels are important determinants of flower visitation regardless of whether plants reproduce sexually or apomictically. In Malus plants, pollen viability seems linked to a high occurrence of apomixis.
Polyploidy and hybridization might induce shifts in mechanisms controlling reproduction, and the repatterning of gene expression might result in the evolution of gametophytic apomixis (Carman, 2001; Grimanelli et al., 2001; Carman, 2007). Pollen development in plants is a complex process that requires meiosis, the division of haploid nuclei, haploid nuclear division, cellular differentiation of male gametophytes, and other processes (Goldberg et al., 1993; Hafidh et al., 2016; Liu and Wang, 2021). In diploid plants, pollen is formed via the meiosis and mitosis of pollen mother cells. We conducted a transcriptome analysis to clarify the relationship between apomixis and pollen development in a population with a low apomictic percentage (P1) and a population with a high apomictic percentage (P2). KEGG analysis revealed DEGs that were significantly enriched in DNA replication, ribosome, plant hormone signal transduction, and purine metabolism (Figure 3B). This suggests that P2 plants show abnormalities in meiosis and hormone metabolism during pollen development and metabolism. Many studies have shown that there is a close relationship between hormone metabolism and pollen development. Auxin levels in pollen grains affect the normal development of stamens and anthers in Arabidopsis (Salinas-Grenet et al., 2018). The programmed death of microspore mother cells during meiosis is partly induced by abscisic acid and indole-3-acetic acid (Kovaleva et al., 2018).
To further clarify the causes of pollen abortion, we studied the pollen development of QNJ, PYTC, and SJHT using paraffin sectioning (Figures 6A–L). We found that meiosis failed in pollen mother cells in PYTC and SJHT. Meiosis is a critically important stage in the process of pollen development, and meiotic abnormalities can result in pollen abortion. Several genes are involved in the regulation of meiosis. In Arabidopsis, approximately 50 genes have been shown to directly regulate meiosis, and many have been shown to indirectly regulate meiosis (Sanchez-Moran and Armstrong, 2014). We found that the transcript level of the meiosis-related genes MD05G1336300 and MD15G1344200 in Malus plants was upregulated in PYTC and SJHT, leading to abnormal division of pollen mother cells, suggesting that meiotic abnormalities were the cause of abnormal pollen development. The regulation of meiosis in Malus species might be more complex given that they are perennial woody plants.
5 Conclusion
In sum, Malus plants that reproduce apomictically may be identified by determining whether the morphology of pollen is abnormal and its viability is low. Thus, the methods used in our study may facilitate the identification of apple plants that are capable of apomictic reproduction. Further characterization of genes involved in pollen mother cell division such as MD05G1336300 (ATK1) and MD15G1344200 (CYCA1;1) may improve our understanding of meiosis in apples for breeding purposes.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Author contributions
Data curation, D-RW. Formal analysis, X-YY, C-HZ and S-HL. Funding acquisition, D-DL and C-XY. Investigation, D-DL and C-XY. Resources, G-LS, R-FZ and H-JG. Visualization, D-RW and X-ST. Writing—original draft, D-RW and D-DL. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Key R&D Program of China (2018YFD1000106), the National Natural Science Foundation of China (32260739 and 31660566), Yunnan Province Applied Basic Research Programs of Science and Technology Department (202001BB050047), and the Yunnan Xingdian Talent Plan (YNWR-QNBJ-2019-058).
Acknowledgments
We thank the professor of Yu-Jin Hao for his assistance with the project. We thank TopEdit (www.topeditsci.org) for linguistic assistance during the preparation of this manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1065032/full#supplementary-material
References
Bicknell, R. A., Koltunow, A. M. (2004). Understanding apomixis: recent advances and remaining conundrums. Plant Cell 16 Suppl, S228–S245. doi: 10.1105/tpc.017921
Boutilier, K., Offringa, R., Sharma, V. K., Kieft, H., Ouellet, T., Zhang, L., et al. (2002). Ectopic expression of BABY BOOM triggers a conversion from vegetative to embryonic growth. Plant Cell 14 (8), 1737–1749. doi: 10.1105/tpc.001941
Carman, J. G. (2001). The gene effect: Genome collisions and apomixis. Flowering Apomixis: From Mech. to Genet. Eng 95–110.
Carman, J. G., Jamison, M., Elliott, E., Dwivedi, K. K., Naumova, T. N. (2011). Apospory appears to accelerate onset of meiosis and sexual embryo sac formation in sorghum ovules. BMC Plant Biol. 11, 9. doi: 10.1186/1471-2229-11-9
Catanach, A. S., Erasmuson, S. K., Podivinsky, E., Jordan, B. R., Bicknell, R. (2006). Deletion mapping of genetic regions associated with apomixis in hieracium. Proc. Natl. Acad. Sci. U.S.A. 103 (49), 18650–18655. doi: 10.1073/pnas.0605588103
Cosmulescu, S. N., Botu, M., Achim, G. (2012). Determination of apomictic fruit set ratio in several Romanian walnut (Juglans regia l.) cultivars. Notulae Botanicae Horti Agrobotanici Cluj-Napoca 40 (1), 229–233. doi: 10.15835/nbha4016407
Fei, X., Shi, Q., Qi, Y., Wang, S., Lei, Y., Hu, H., et al. (2021). ZbAGL11, a class d MADS-box transcription factor of zanthoxylum bungeanum, is involved in sporophytic apomixis. Hortic. Res. 8 (1), 23. doi: 10.1038/s41438-020-00459-x
Figueiredo, D. D., Köhler, C. (2018). Auxin: a molecular trigger of seed development. Genes Dev. 32 (7-8), 479–490. doi: 10.1101/gad.312546.118
Goldberg, R. B., Beals, T. P., Sanders, P. M. (1993). Anther development: basic principles and practical applications. Plant Cell 5 (10), 1217–1229. doi: 10.1105/tpc.5.10.1217
Grimanelli, D., Leblanc, O., Perotti, E., Grossniklaus, U. (2001). Developmental genetics of gametophytic apomixis. Trends Genet. 17 (10), 597–604. doi: 10.1016/S0168-9525(01)02454-4
Hafidh, S., Fíla, J., Honys, D. (2016). Male Gametophyte development and function in angiosperms: a general concept. Plant Reprod. 29 (1-2), 31–51. doi: 10.1007/s00497-015-0272-4
Hand, M. L., Koltunow, A. M. (2014). The genetic control of apomixis: asexual seed formation. Genetics 197 (2), 441–450. doi: 10.1534/genetics.114.163105
Hao, Y.-J., Liu, D.-D. (2011). Apomixis in plants–10000 scientific problems 381–383. (China: Science press).
Jeong, H. J., Kang, J. H., Zhao, M., Kwon, J. K., Choi, H. S., Bae, J. H., et al. (2014). Tomato Male sterile 1035 is essential for pollen development and meiosis in anthers. J. Exp. Bot. 65 (22), 6693–6709. doi: 10.1093/jxb/eru389
Khanday, I., Skinner, D., Yang, B., Mercier, R., Sundaresan, V. (2019). A male-expressed rice embryogenic trigger redirected for asexual propagation through seeds. Nature 565 (7737), 91–95. doi: 10.1038/s41586-018-0785-8
Koltunow, A. M., Grossniklaus, U. (2003). Apomixis: a developmental perspective. Annu. Rev. Plant Biol. 54, 547–574. doi: 10.1146/annurev.arplant.54.110901.160842
Koltunow, A. M., Johnson, S. D., Rodrigues, J. C., Okada, T., Hu, Y., Tsuchiya, T., et al. (2011). Sexual reproduction is the default mode in apomictic hieracium subgenus pilosella, in which two dominant loci function to enable apomixis. Plant J. 66 (5), 890–902. doi: 10.1111/j.1365-313X.2011.04556.x
Kovaleva, L. V., Voronkov, A. S., Zakharova, E. V., Andreev, I. M. (2018). ABA and IAA control microsporogenesis in petunia hybrida l. Protoplasma 255 (3), 751–759. doi: 10.1007/s00709-017-1185-x
Kuhn, D. N., Bally, I. S. E., Dillon, N. L., Innes, D., Groh, A. M., Rahaman, J., et al. (2017). Genetic map of mango: A tool for mango breeding. Front. Plant Sci. 8, 577. doi: 10.3389/fpls.2017.00577
Leblanc, O., Grimanelli, D., Hernandez-Rodriguez, M., Galindo, P. A., Soriano-Martinez, A. M., Perotti, E. (2009). Seed development and inheritance studies in apomictic maize-tripsacum hybrids reveal barriers for the transfer of apomixis into sexual crops. Int. J. Dev. Biol. 53 (4), 585–596. doi: 10.1387/ijdb.082813ol
Liu, D.-D., Fang, M.-J., Dong, Q.-L., Hu, D.-G., Zhou, L.-J., Sha, G.-L., et al. (2014). Unreduced embryo sacs escape fertilization via a ‘female-late-on-date’ strategy to produce clonal seeds in apomictic crabapples. Scientia Hortic. 167, 76–83. doi: 10.1016/j.scienta.2013.12.035
Liu, L., Wang, T. (2021). Male Gametophyte development in flowering plants: A story of quarantine and sacrifice. J. Plant Physiol. 258-259, 153365. doi: 10.1016/j.jplph.2021.153365
Maia, F. R., Varassin, I. G., Goldenberg, R. (2016). Apomixis does not affect visitation to flowers of Melastomataceae, but pollen sterility does. Plant Biol. (Stuttgart Germany) 18 (1), 132–138. doi: 10.1111/plb.12364
Mercier, R., Armstrong, S. J., Horlow, C., Jackson, N. P., Makaroff, C. A., Vezon, D., et al. (2003). The meiotic protein SWI1 is required for axial element formation and recombination initiation in arabidopsis. Development 130 (14), 3309–3318. doi: 10.1242/dev.00550
Mercier, R., Vezon, D., Bullier, E., Motamayor, J. C., Sellier, A., Lefèvre, F., et al. (2001). SWITCH1 (SWI1): a novel protein required for the establishment of sister chromatid cohesion and for bivalent formation at meiosis. Genes Dev. 15 (14), 1859–1871. doi: 10.1101/gad.203201
Okada, T., Catanach, A. S., Johnson, S. D., Bicknell, R. A., Koltunow, A. M. (2007). An hieracium mutant, loss of apomeiosis 1 (loa1) is defective in the initiation of apomixis. Sexual Plant Reprod. 20 (4), 199–211. doi: 10.1007/s00497-007-0057-5
Ortiz, J. P. A., Leblanc, O., Rohr, C., Grisolia, M., Siena, L. A., Podio, M., et al. (2019). Small RNA-seq reveals novel regulatory components for apomixis in paspalum notatum. BMC Genomics 20 (1), 487. doi: 10.1186/s12864-019-5881-0
Ranganath, R. M. (2004). Harnessing the developmental potential of nucellar cells: barriers and opportunities. Trends Biotechnol. 22 (10), 504–510. doi: 10.1016/j.tibtech.2004.08.010
Rodriguez-Leal, D., Vielle-Calzada, J. P. (2012). Regulation of apomixis: learning from sexual experience. Curr. Opin. Plant Biol. 15 (5), 549–555. doi: 10.1016/j.pbi.2012.09.005
Salinas-Grenet, H., Herrera-Vásquez, A., Parra, S., Cortez, A., Gutiérrez, L., Pollmann, S., et al. (2018). Modulation of auxin levels in pollen grains affects stamen development and anther dehiscence in arabidopsis. Int. J. Mol. Sci. 19 (9), 2480. doi: 10.3390/ijms19092480
Sanchez-Moran, E., Armstrong, S. J. (2014). Meiotic chromosome synapsis and recombination in arabidopsis thaliana: new ways of integrating cytological and molecular approaches. Chromosome Res. 22 (2), 179–190. doi: 10.1007/s10577-014-9426-8
Santos, M. O., Aragão, F. J. (2009). Role of SERK genes in plant environmental response. Plant Signal Behav. 4 (12), 1111–1113. doi: 10.4161/psb.4.12.9900
Sax, K. (1959). The cytogenetics of facultative apomixis in malus species. J. Arnold Arboretum 40 (3), 289–297. doi: 10.5962/p.186033
Schmidt, H. (1977). Contributions on the breeding of apomictic apple stocks. Z Pflanzenzucht 78, 3–12.
Smith, J. (1841). Notice of a plant which produces seeds without any apparent action of pollen. Trans. Linnaean Soc. London (meeting June 18 1839) 1841, 18, (4), 509–512. doi: 10.1111/j.1095-8339.1838.tb00200.x
Somaratne, Y., Tian, Y., Zhang, H., Wang, M., Huo, Y., Cao, F., et al. (2017). ABNORMAL POLLEN VACUOLATION1 (APV1) is required for male fertility by contributing to anther cuticle and pollen exine formation in maize. Plant J. 90 (1), 96–110. doi: 10.1111/tpj.13476
Tucker, M. R., Koltunow, A. M. G. (2009). Sexual and asexual (apomictic) seed development in flowering plants: molecular, morphological and evolutionary relationships. Funct. Plant Biol. 36 (6), 490–504. doi: 10.1071/FP09078
Wang, N., Song, X., Ye, J., Zhang, S., Cao, Z., Zhu, C., et al. (2022). Structural variation and parallel evolution of apomixis in citrus during domestication and diversification. Natl. Sci. Rev 9 (10). doi: 10.1093/nsr/nwac114
Wang, D., Wang, X., Zhang, C., Yang, K., Wang, X., Cui, J., et al. (2022). Genome-wide identification, expression, and functional analysis of MdMSI genes in apples (Malus domestica borkh.). Front. Genet. 13, 846321. doi: 10.3389/fgene.2022.846321
Wang, X., Xu, Y., Zhang, S., Cao, L., Huang, Y., Cheng, J., et al. (2017). Genomic analyses of primitive, wild and cultivated citrus provide insights into asexual reproduction. Nat. Genet. 49 (5), 765–772. doi: 10.1038/ng.3839
Wang, D. R., Yang, K., Wang, X., Lin, X. L., Rui, L., Liu, H. F., et al. (2022). Overexpression of MdZAT5, an C2H2-type zinc finger protein, regulates anthocyanin accumulation and salt stress response in apple calli and arabidopsis. Int. J. Mol. Sci. 23 (3), 1897. doi: 10.3390/ijms23031897
Wu, J., Shahid, M. Q., Chen, L., Chen, Z., Wang, L., Liu, X., et al. (2015). Polyploidy enhances F1 pollen sterility loci interactions that increase meiosis abnormalities and pollen sterility in autotetraploid rice. Plant Physiol. 169 (4), 2700–2717. doi: 10.1104/pp.15.00791
Wu, G. A., Sugimoto, C., Kinjo, H., Azama, C., Mitsube, F., Talon, M., et al. (2021). Diversification of mandarin citrus by hybrid speciation and apomixis. Nat. Commun. 12 (1), 4377. doi: 10.1038/s41467-021-24653-0
Yamashita, K., Nakazawa, Y., Namai, K., Amagai, M., Tsukazaki, H., Wako, T., et al. (2012). Modes of inheritance of two apomixis components, diplospory and parthenogenesis, in Chinese chive (Allium ramosum) revealed by analysis of the segregating population generated by back-crossing between amphimictic and apomictic diploids. Breed Sci. 62 (2), 160–169. doi: 10.1270/jsbbs.62.160
Yang, Z., Liu, L., Sun, L., Yu, P., Zhang, P., Abbas, A., et al. (2019). OsMS1 functions as a transcriptional activator to regulate programmed tapetum development and pollen exine formation in rice. Plant Mol. Biol. 99 (1-2), 175–191. doi: 10.1007/s11103-018-0811-0
Yang, J., Tian, L., Sun, M. X., Huang, X. Y., Zhu, J., Guan, Y. F., et al. (2013). AUXIN RESPONSE FACTOR17 is essential for pollen wall pattern formation in arabidopsis. Plant Physiol. 162 (2), 720–731. doi: 10.1104/pp.113.214940
Keywords: apple, apomixis, transcriptome analysis, pollen abortion, apomeiosis
Citation: Liu D-D, Wang D-R, Yang X-Y, Zhao C-H, Li S-H, Sha G-L, Zhang R-F, Ge H-J, Tong X-S and You C-X (2023) Apomictic Malus plants exhibit abnormal pollen development. Front. Plant Sci. 14:1065032. doi: 10.3389/fpls.2023.1065032
Received: 09 October 2022; Accepted: 30 January 2023;
Published: 20 February 2023.
Edited by:
Ryo Fujimoto, Kobe University, JapanReviewed by:
Pablo Bolaños-Villegas, University of Costa Rica, Costa RicaChinedu Charles Nwafor, University of Nebraska-Lincoln, United States
Copyright © 2023 Liu, Wang, Yang, Zhao, Li, Sha, Zhang, Ge, Tong and You. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dan-Dan Liu, liudandan@ynu.edu.cn; Chun-Xiang You, youchunxiang@sdau.edu.cn
†These authors have contributed equally to this work