Whether Gametophytes Are Reduced or Unreduced in Angiosperms Might Be Determined Metabolically

In angiosperms, meiotic failure coupled with the formation of genetically unreduced gametophytes in ovules (apomeiosis) constitute major components of gametophytic apomixis. These aberrant developmental events are generally thought to be caused by mutation. However, efforts to locate the responsible mutations have failed. Herein, we tested a fundamentally different hypothesis: apomeiosis is a polyphenism of meiosis, with meiosis and apomeiosis being maintained by different states of metabolic homeostasis. Microarray analyses of ovules and pistils were used to differentiate meiotic from apomeiotic processes in Boechera (Brassicaceae). Genes associated with translation, cell division, epigenetic silencing, flowering, and meiosis characterized sexual Boechera (meiotic). In contrast, genes associated with stress responses, abscisic acid signaling, reactive oxygen species production, and stress attenuation mechanisms characterized apomictic Boechera (apomeiotic). We next tested whether these metabolic differences regulate reproductive mode. Apomeiosis switched to meiosis when premeiotic ovules of apomicts were cultured on media that increased oxidative stress. These treatments included drought, starvation, and H2O2 applications. In contrast, meiosis switched to apomeiosis when premeiotic pistils of sexual plants were cultured on media that relieved oxidative stress. These treatments included antioxidants, glucose, abscisic acid, fluridone, and 5-azacytidine. High-frequency apomeiosis was initiated in all sexual species tested: Brassicaceae, Boechera stricta, Boechera exilis, and Arabidopsis thaliana; Fabaceae, Vigna unguiculata; Asteraceae, Antennaria dioica. Unreduced gametophytes formed from ameiotic female and male sporocytes, first division restitution dyads, and nucellar cells. These results are consistent with modes of reproduction and types of apomixis, in natural apomicts, being regulated metabolically.

Fourteen DEGs identified between the YG stages of B. microphylla and B. × formosa were selected for quantitative reverse transcription PCR (qRT-PCR). RNA was converted to cDNA using Superscript III First-Strand Synthesis kits (Invitrogen). Primers were designed from A. thaliana data using Primer3 (bioinfo.ut.ee/primer3/), quality checked by PrimerSelect (dnastar.com), and searched (blast.ncbi.nlm. nih.gov/Blast.cgi) for like sequences. Sequences predicted to form internal loops or dimers or to be homologous to other genes were discarded. Primers were tested using Boechera gDNA. qRT-PCR reactions (25 µL) contained QuantiTect SYBR Green (Qiagen), 200-250 nM forward and reverse primers and 3-5 µL of dilute cDNA. Reaction cycles were: 50 • C, 2 min; 95 • C, 15 min; 40 cycles of 94 • C for 15 s, 58 • C for 30 s, and 72 • C for 30 s. Melt curve analyses were performed from 65 • C to 95 • C in 5 sec 0.5 • C increments. A DNA Engine (Opticon 2, Continuous Fluorescence Detection System) with MicroAmp Optical 96-well plates (Bio-Rad, Hercules, CA, USA) was used. Three reps with two technical reps each were conducted. Expression values and standard deviations were calculated using the comparative C T method using separate sets of housekeeping replicates for each gene. Delta Ct values were normalized (delta-delta Ct method) and subjected to Student's t-tests (p ≤ 0.05).

Drought and Heat-Stress Experiments
Seedlings (ca. 60 days post-planting) were randomly assigned to three groups, well-watered, drought-stressed, and heat plus drought-stressed (8-12 plants per taxon per group) and vernalized (60 days, 4 • C). After vernalization, plants targeted for the well-watered and droughted treatments were grown in the greenhouse, and plants targeted for heat plus drought were grown in a growth chamber (32 • C, 500 µmol m −2 s −1 photosynthetic photon flux, 16/8 h day/night photoperiod). To minimize surface evaporation, perlite was added to the surface of pots, and field capacity weights were obtained (weighed 2 h after drenching). Enough water was added each day to the well-watered plants to bring them to field capacity (based on pot weight). Water was added to droughted plants each day to equal 50% of the transpiration rate of well-watered plants. Upon flowering, pistils were collected for cytology and expression profiling. Ovule parameters were compared by chi-square tests [25].

Expression Profiling
We recently documented reproduction in 44 Boechera taxa. These ranged from 100% sexual in species to nearly 100% apomictic, by TDG formation, in hybrids. Between these extremes were facultative apomicts that expressed combinations of meiosis and TDG and/or Hieracium-type gametophyte (HAG) formation [16]. Here, we report transcriptome analyses for four of these taxa: sexual B. stricta, strongly diplosporous B. × formosa, facultatively aposporous, and diplosporous B. microphylla (Cache County, UT), and facultatively diplosporous B. lignifera. Pistil lengths, array normalization results, array vs. qRT-PCR comparisons, and primers used for qRT-PCR are shown in Figures S1-S3 and Table S4, respectively.
Comparisons among the MMC (megasporocyte), AM (active meiocyte), and YG (young gametophyte) stages within diplosporous B. lignifera and aposporous B. microphylla produced few DEGs (Table S5). This near absence of DEGs may have occurred because pistil length intervals for the respective samples differed only minimally ( Figure S1). However, few DEGs were detected at these stages even when comparisons were made across taxa. Hence, we combined the MMC, AM, and YG stages, within taxa, to make a single B. lignifera vs. B. microphylla comparison, which yielded 283 DEGs. At the MG stage (mature gametophyte), 609 DEGs were detected (Figure 1a, Table S5).
Genes 2020, 11, x FOR PEER REVIEW 6 of 37 respective samples differed only minimally ( Figure S1). However, few DEGs were detected at these stages even when comparisons were made across taxa. Hence, we combined the MMC, AM, and YG stages, within taxa, to make a single B. lignifera vs. B. microphylla comparison, which yielded 283 DEGs. At the MG stage (mature gametophyte), 609 DEGs were detected (Figure 1a, Table S5). Across all comparisons, 3723 distinct DEGs were detected ( Figure 1a; Table S5), which are linked to 3980 distinct GO terms (Figure 1b; Table S6). Overrepresented GOs (OGOs) occurs when more DEGs are linked to a GO term than expected by chance (pantherdb.org/). Among our OGO comparisons (Figure 1c), 930 distinct OGOs were detected (Table S7). EGOs occur when more genes of a GO term are up or down-regulated (significant or not) than expected by chance (pantherdb.org/). Since all expression data are used, EGO analyses can detect shifts in gene expression that may be missed by OGO analyses [31,32]. Among our EGO comparisons (Figure 1d), 1773 distinct EGOs were detected (Table S8). Collectively, 4822 distinct GO terms were identified. To more readily visualize large metabolism-related shifts in gene expression, we partitioned the 4822 GO terms into 13 userdefined subgroups (Tables S1, S2). Subgroup frequencies per taxon were then compared among taxa by chi-square analyses (Figures 2 and 3; Table S3; see Materials and Methods). Across all comparisons, 3723 distinct DEGs were detected (Figure 1a; Table S5), which are linked to 3980 distinct GO terms (Figure 1b; Table S6). Overrepresented GOs (OGOs) occurs when more DEGs are linked to a GO term than expected by chance (pantherdb.org/). Among our OGO comparisons (Figure 1c), 930 distinct OGOs were detected (Table S7). EGOs occur when more genes of a GO term are up or down-regulated (significant or not) than expected by chance (pantherdb.org/). Since all expression data are used, EGO analyses can detect shifts in gene expression that may be missed by OGO analyses [31,32]. Among our EGO comparisons (Figure 1d), 1773 distinct EGOs were detected (Table S8). Collectively, 4822 distinct GO terms were identified. To more readily visualize large metabolism-related shifts in gene expression, we partitioned the 4822 GO terms into 13 user-defined subgroups (Tables S1 and S2). Subgroup frequencies per taxon were then compared among taxa by chi-square analyses (Figures 2 and 3; Table S3; see Materials and Methods).  . Normal and exploded subgroups tended to be more frequently observed in less and more meiotic taxa, respectively. Starred subgroups differed significantly in frequency (p ≤ 0.05): green stars, normal; red stars, exploded; values next to pies are total numbers of GOs observed.
Only 5-20% of ovules of the TDG forming B. lignifera were facultatively meiotic [16]. Hence, we tested whether fewer meiosis related GO terms (all GOs of DEGs) occurred in B. lignifera compared to the TDG and HAG forming B. microphylla, where ca. 40% of ovules were meiotic [16]. For the combined MMC to YG comparison, meiosis and gametophyte related GOs occurred nearly exclusively in the more meiotic (and aposporous) B. microphylla (Table S6, lines 1-35 vs. 1098-1151, see also lines 2164-92 vs. 4112-62). These meiosis and gametophyte related GOs belong to the 'cell division' subgroup, which collectively represented a significantly larger percentage of GO terms in the aposporous B. microphylla compared to the diplosporous B. lignifera (Figure 2a, red star).
We used EGO analyses to compare B. lignifera to B. microphylla at each stage individually. The MMC stage produced the most EGOs ( Figure 1d). Here, signaling, stress, bioenergetics, metabolism, catabolism, and transporter EGOs occurred more frequently in ovules of the TDG forming B. lignifera (least meiotic) (Figure 2b). These included carbohydrate, fatty acid, and protein metabolism processes, responses to organic and inorganic stimuli, activation of homeostasis processes, regulation of ROS and other stresses, and ABA signaling (Table S8, lines 11-79, 84-91, 95-144). In contrast, transcription, translation, development, and cell division EGOs were more frequent in ovules of the more meiotic HAG forming B. microphylla (Figure 2b). These included 13 meiosis associated terms and numerous ribosome-associated terms (Table S8,   . Normal and exploded subgroups tended to be more frequently observed in less and more meiotic taxa, respectively. Starred subgroups differed significantly in frequency (p ≤ 0.05): green stars, normal; red stars, exploded; values next to pies are total numbers of GOs observed.
Only 5-20% of ovules of the TDG forming B. lignifera were facultatively meiotic [16]. Hence, we tested whether fewer meiosis related GO terms (all GOs of DEGs) occurred in B. lignifera compared to the TDG and HAG forming B. microphylla, where ca. 40% of ovules were meiotic [16]. For the combined MMC to YG comparison, meiosis and gametophyte related GOs occurred nearly exclusively in the more meiotic (and aposporous) B. microphylla (Table S6, lines 1-35 vs. 1098-1151, see also lines 2164-92 vs. 4112-62). These meiosis and gametophyte related GOs belong to the 'cell division' subgroup, which collectively represented a significantly larger percentage of GO terms in the aposporous B. microphylla compared to the diplosporous B. lignifera (Figure 2a, red star).
We used EGO analyses to compare B. lignifera to B. microphylla at each stage individually. The MMC stage produced the most EGOs ( Figure 1d). Here, signaling, stress, bioenergetics, metabolism, catabolism, and transporter EGOs occurred more frequently in ovules of the TDG forming B. lignifera (least meiotic) (Figure 2b). These included carbohydrate, fatty acid, and protein metabolism processes, responses to organic and inorganic stimuli, activation of homeostasis processes, regulation of ROS and other stresses, and ABA signaling (Table S8, lines 11-79, 84-91, 95-144). In contrast, transcription, translation, development, and cell division EGOs were more frequent in ovules of the more meiotic HAG forming B. microphylla (Figure 2b). These included 13 meiosis associated terms and numerous ribosome-associated terms (Table S8, lines  . Fewer differences were observed at the AM, YG, and MG stages (Table S8, lines 290-471; Figure 2a,b). Collectively, these findings revealed major metabolism-related differences. Gene expression in the TDG forming B. lignifera favored signaling, stress, and catabolism. In the meiotic HAG forming B. microphylla, gene expression favored transcription, translation, development, and cell division.
Meiosis associated EGO subgroups at the MMC stage (exploded pie sections, Figure 2b) comprised ca. 10% and 80% of the diplosporous B. lignifera and the more meiotic but aposporous B. microphylla, respectively. Interestingly, nearly the same percentages were observed when diplosporous B. lignifera and aposporous B. microphylla were compared individually to the more diplosporous B. × formosa (Figure 2f,j). These results reveal a continuum where meiosis associated gene expression was highest in B. microphylla and lowest in B. × formosa with B. lignifera in between (compare Figure 2b,f,j). One of our objectives was to determine if these distinct shifts in gene expression are responsible for reproductive mode and apomixis type or are simply coincidental (discussed below). Analyses involving EGOs, OGOs, and all GOs of DEGs at the YG or AM/YG stages generally mirrored those at the AM stage for comparisons between the nearly obligate apomictic B. × formosa and the facultatively apomictic B. lignifera and B. microphylla (Figure 2f-i, j-m).
For this study, we collected Boechera species from native habitats in the Great Basin and Rocky Mountain Cordillera of North America [16]. However, embryological analyses revealed apomixis for taxa we thought, at the time of collection, were sexual. We eventually found sexual species, but resources for ovule dissections were no longer available. Instead, we created expression profiles from whole pistils of sexual B. stricta, each pistil of which contained ca. 160-180 ovules. We then compared ovule profiles of the apomicts to the same pistil profile from sexual B. stricta ( Figure 3). We first considered differential gene expression between the nearly obligate TDG forming B. × formosa versus the completely sexual B. stricta. This comparison produced the most DEGs, OGOs, all GOs of DEGs and EGOs ( Figure 1). Importantly, it produced highly similar frequency patterns for meiosis-associated vs. apomeiosis-associated EGO subgroups (Figure 3a) as had been observed for the more vs. the less meiotic apomicts (Figure 2b,f,j). Similar GO subgroup patterns were observed for OGOs and all GOs of DEGs (Figure 3a-c) and for the facultative apomicts compared to sexual B. stricta (Figure 3d-g).
Carbon starvation downregulates translation by inactivating TARGET OF RAPAMYCIN (TOR) [33]. DEGs in the red and green boxes (Figure 4a) are translation machinery genes ( Table S9). Most of these were upregulated in B. stricta (red box), which suggests that diplosporous B. × formosa was experiencing carbon starvation. Further evidence for this involved upregulated levels of ARABIDOPSIS NAC DOMAIN CONTAINING PROTEIN 2 (ATAF1), a transcription factor (TF) that upregulates TREHALASE1 (TRE1). TRE1 promotes the formation of a sugar starvation metabolome by catabolizing the signaling molecule trehalose-6-phosphate, thus further reducing TOR signaling by increasing SUCROSE-NON-FERMENTING-RELATED PROTEIN KINASE-1 (SnRK1) activity [34][35][36]. Carbon starvation also upregulates the vacuolar glucose exporter EARLY RESPONSIVE TO DEHYDRATION-LIKE 6 gene (ERDL6) [37] and the DUAL-SPECIFICITY PROTEIN PHOSPHATASE 4 gene (SEX4) [38], the products of which replenish monosaccharide levels to maintain primary metabolism. SEX4 might also inactivate SnRK1, thus reinitiating TOR activity in the presence of restored energy levels [39]. These three-carbon starvation genes were upregulated in diplosporous B. × formosa (Figure 4a). In contrast, SNF1-RELATED PROTEIN KINASE REGULATORY β SUBUNIT 1 (AKINβ1), which encodes an SnRK1 subunit that directs SnRK1 to the nucleus for gene regulation but also silences SnRK1 upon N-myristoylation [40], was downregulated in B. × formosa (Figure 4a), possibly as a feedback mechanism responding to TOR inactivation.  (Table S9). Carbon starvation downregulates translation by inactivating TARGET OF RAPAMYCIN (TOR) [33]. DEGs in the red and green boxes (Figure 4a) are translation machinery genes ( Table S9). Most of these were upregulated in B. stricta (red box), which suggests that diplosporous B. x formosa was experiencing carbon starvation. Further evidence for this involved upregulated levels of ARABIDOPSIS NAC DOMAIN CONTAINING PROTEIN 2 (ATAF1), a transcription factor (TF) that upregulates TREHALASE1 (TRE1). TRE1 promotes the formation of a sugar starvation metabolome by catabolizing the signaling molecule trehalose-6-phosphate, thus further reducing TOR signaling by increasing SUCROSE-NON-FERMENTING-RELATED PROTEIN KINASE-1 (SnRK1) activity [34][35][36]. Carbon starvation also upregulates the vacuolar glucose exporter EARLY RESPONSIVE TO DEHYDRATION-LIKE 6 gene (ERDL6) [37] and the DUAL-SPECIFICITY PROTEIN PHOSPHATASE 4 gene (SEX4) [38], the products of which replenish monosaccharide levels to maintain primary metabolism. SEX4 might also inactivate SnRK1, thus reinitiating TOR activity in the presence of restored energy levels [39]. These three-carbon starvation genes were upregulated in diplosporous B. x formosa (Figure 4a). In contrast, SNF1-RELATED PROTEIN KINASE REGULATORY β SUBUNIT 1 (AKINβ1), which encodes an SnRK1 subunit that directs SnRK1 to the nucleus for gene regulation but also silences SnRK1 upon N-myristoylation [40], was downregulated in B. x formosa (Figure 4a), possibly as a feedback mechanism responding to TOR inactivation. The early flowering gene EMBRYONIC FLOWER 2 (EMF2) encodes a polycomb repressive complex member that suppresses flowering [41]. It was downregulated in B. x formosa (Figure 4b), suggesting that this apomict may exhibit an early flowering phenotype, which is common for apomicts in nature [4]. Additionally, genes responsible for flower formation were mostly downregulated in B. x formosa (Figure 4b, red box). Taken together, these results suggest that the earliness of flowering was favored in apomictic B. x formosa but that several epigenetically regulated The early flowering gene EMBRYONIC FLOWER 2 (EMF2) encodes a polycomb repressive complex member that suppresses flowering [41]. It was downregulated in B. × formosa (Figure 4b), suggesting that this apomict may exhibit an early flowering phenotype, which is common for apomicts in nature [4]. Additionally, genes responsible for flower formation were mostly downregulated in B. × formosa (Figure 4b, red box). Taken together, these results suggest that the earliness of flowering was favored in apomictic B. × formosa but that several epigenetically regulated programs of reproduction were only weakly expressed in this taxon.
In apomictic strains of Saccharomyces cerevisiae (brewer's yeast), heat stress switches apomeiotic divisions, which normally produce dyads of unreduced spores, to meiotic divisions, which produce tetrads of reduced spores. However, when mitochondrial transcription in apomictic yeast is interrupted by erythromycin, the heat-stress switch from apomeiosis to meiosis no longer occurs, i.e., the yeast remains apomictic. This finding, plus research in other organisms, suggests that meiosis depends on stress-induced nucleomitochondrial communication [3,42,43]. In the present study, genes for two nuclear TFs, NAC DOMAIN CONTAINING PROTEIN 17 (ANACO17) and 53 (ANACO53), which participate in drought-induced mitochondrial retrograde signaling, were upregulated in sexual B. stricta (downregulated in the TDG-forming 1st division restitution B. × formosa, Figure 4c). This is consistent with a mitochondrial role in the initiation and completion of meiosis I, a role that is downregulated during 1st division restitution apomeioses. ANAC017 normally induces senescence, suppresses auxin and chloroplast function, reprograms mitochondrial metabolism for lower energy conditions, and possibly activates SnRK1 [44]. It also upregulates ANACO53, which upregulates ROS production genes and initiates unfolded protein responses (UPRs) that attenuate proteotoxic conditions in the ER [45]. MMS ZWEI HOMOLOGUE 2 (MMZ2), which was strongly upregulated in B. stricta, is upregulated by DNA damage and UPRs, and it normally participates in post replicative DNA repair. Since these genes were upregulated in sexual B. stricta, their expression is consistent with the view that meiosis evolved as a DNA repair mechanism [15,[46][47][48] at about the same time that early eukaryotes were domesticating their oxygen respiring endocytically obtained organelles [46]. The downregulation of meiosis genes and mitochondrial stress response genes in apomictic B. × formosa (Figure 4c,d) is further evidence that reproductive processes, from flower formation through meiosis, were only weakly supported at the molecular level in apomictic B. × formosa.
DEGs of two additional signaling related categories, gene silencing and BR signaling, also provide clues concerning sex apomixis switching ( Figure 5). Chromatin remodeling by methylation and demethylation occurs throughout plant development. In sexual cell lineages, RNA-directed DNA methylation (RdDM) is required to accomplish meiosis-specific gene splicing [49]. RdDM may also occur in tissues that produce sex cells because it silences potentially harmful transposons [49]. In the present study, silencing processes by RdDM and by polycomb repressive complex activity were upregulated in the ovules of sexual B. stricta vs. apomictic B. × formosa (Figure 5a-c). This again is consistent with other upregulated reproductive processes in B. stricta as described above.
ALFIN-LIKE 7 (AL7) is a nuclear-localized H3K4me3 binding protein that directs polycomb repressive complex (PcG) 1s (PCR1s) to transcriptionally active chromatin where the H3K4me3 expression marks are converted to H2Aub1 repression marks. These are then stabilized by PCR2 H3K27me3 silencing marks [50]. ULTRAPETALA1 (ULT1) is a trithorax factor that functions as an anti-repressor to counteract PcG silencing at 1000s of loci, including many genes where specific up or downregulation is required for normal reproductive development [51,52]. The upregulation in B. stricta of AL7 and ULT1 genes (Figure 5a-c) is further evidence that floral development at the molecular level was strongly supported in sexual B. stricta but weakly supported in apomictic B. × formosa. Overall, it appears that B. stricta ovules embarked on epigenetic paths of sexual development that apomeiotic ovules were reluctant to follow. provide clues concerning sex apomixis switching ( Figure 5). Chromatin remodeling by methylation and demethylation occurs throughout plant development. In sexual cell lineages, RNA-directed DNA methylation (RdDM) is required to accomplish meiosis-specific gene splicing [49]. RdDM may also occur in tissues that produce sex cells because it silences potentially harmful transposons [49]. In the present study, silencing processes by RdDM and by polycomb repressive complex activity were upregulated in the ovules of sexual B. stricta vs. apomictic B. x formosa (Figure 5a-c). This again is consistent with other upregulated reproductive processes in B. stricta as described above.  BR interacts with other processes to regulate growth and development as well as abiotic and biotic stresses [53,54]. The major hub of BR signaling is the kinase BRASSINOSTEROID-INSENSITIVE 2 (BIN2), which phosphorylation inactivates the major BR TFs BRASSINAZOLE-RESISTANT 1 (BZR1) and BRI1-EMS-SUPPRESSOR 1 (BES1). When BIN2 is inhibited, BZR1, and BES1 are activated, and this activation upregulates cascades of gene expression involving multiple downstream TFs and ca. 2000 downstream genes [55]. BR is a growth-promoting hormone. Hence, it is generally upregulated by TOR and downregulated by ABA. Consistent with this pattern, the two BR regulated signaling genes identified in our B. × formosa vs. B. stricta comparison (Table S9) were upregulated in B. stricta ( Figure 5d): CYCLIN D3;1 (CYCD3) is a cyclin that initiates cell division [56], and BR-ENHANCED EXPRESSION 1 (BEE1) combines with a blue light photoreceptor, CRYPTOCHROME 2 (CRY2), to initiate photoperiod induced flowering [57]. Several additional BR associated genes were upregulated in B. stricta. These included two BR synthesis genes ( Figure 5d), which, like TOR (Figure 4a), are upregulated under favorable conditions: HALF FILLED (CESTA) is a TF that upregulates multiple BR biosynthesis genes [58], including CONSTITUTIVE PHOTOMORPHOGENIC DWARF (CPD), which is critical for BR biosynthesis [59]. Additionally upregulated in B. stricta were two BR receptor-associated proteins, including BRI1-ASSOCIATED RECEPTOR KINASE (BAK1) [54] and SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE 1 (SERK1) [60]. Both of these proteins enhance BR signal transduction across membranes by forming protective complexes with the cell surface BR receptor kinase BR INSENSITIVE 1 (BRI1). One of these, SERK1, also binds EXCESS MICROSPOROCYTES 1 (EMS1) during reproductive development in processes that convert somatic cells to reproductive cells (microsporocytes) in anthers [61]. The downregulation of these BR genes in B. × formosa is consistent with other downregulated reproduction processes in B. × formosa (Figures 4 and 5).

Molecular Pathways Upregulated in Apomeiotic Taxa
Several BR associated genes were upregulated in B. × formosa (Figure 5d). One of these, BRASSINOSTEROID-6-OXIDASE 2 (CYP85A2), is a circadian rhythm regulated gene that encodes an enzyme that catalyzes the last step in BR biosynthesis. When upregulated, CYP85A2 increases BR but decreases ABA signaling [62]. BRZ-INSENSITIVE-LONG HYPOCOTYLS 4 (BIL4) was also upregulated in B. × formosa. It encodes a BRI1-interacting transmembrane protein that is involved in cell elongation. BIL4 suppresses BRI1 degradation and physically maintains BRI1 close to the plasma membrane [63]. Two protein phosphatase 2C (PP2C) genes, ABA INSENSITIVE 1 (ABI1) and PROTEIN PHOSPHATASE 2CA (PP2CA) were also upregulated in B. × formosa. These genes upregulate BZR1 and BES1 by phosphatase inactivating BIN2 (Figure 5d). They are also upregulated by ABA as part of a feedback inhibition mechanism where (i) stress induces ABA biosynthesis, (ii) ABA inactivates ABI1 and PP2CA, (iii) ABA signaling upregulates ABI1 and PP2CA biosynthesis, and (iv) additional accumulation of ABI1 and PP2CA suppresses ABA signaling ( Figure 6). These variable findings are not surprising given the abundant crosstalk that occurs among tissues and cells during hormonal and environmental signaling [54].
Genes 2020, 11, x FOR PEER REVIEW 13 of 37 Figure 6. ABA signaling associated differentially expressed genes (DEGs) observed between diplosporous Boechera x formosa and sexual Boechera stricta (see Tables S5 and S9). Superscripts identify taxa in which DEGs are upregulated; F, B. x formosa; S, B. stricta; red and green highlights, up and down-regulation in B. stricta, respectively; arrows indicate upregulation or activation; blunt end lines represent downregulation or suppression.
A very different state of cellular homeostasis exists during stress than during favorable conditions. During stress, ABA-inactivated PP2Cs permit the self-activation by autophosphorylation of SnRK2,3s as well as BIN2. Activated BIN2 phosphorylation inactivates the growth-enhancing BR TFs but phosphorylation activates the NAC TF RESPONSIVE TO DESICCATION 26 (RD26), which upregulates multiple stress response genes [53,64]. RD26 was upregulated 4.2 fold in B. x formosa vs. B. stricta ( Figure 5d, Table S9). Activated SnRK2,3s upregulate ABA specific stress response genes [65] (Figure 6). These data suggest that apomictic ovules experienced states of cellular homeostasis that were maintained by enhanced stress perception with feedback suppression of ABA signaling occurring through CYP85A2, ABI1, and PP2CA [53,62,64,65].
The molecular pathways presented to this point suggest that apomeiosis is caused by carbon starvation (Figure 4a, low-level TOR activity), ABA ( Figure 6), and upregulated drought response genes ( Figure 5d). However, these pathways do not explain how a stress-like state of metabolic homeostasis induces apomixis, especially when stressed apomicts in many cases revert to sex [15,66]. Recognizing this paradox, we hypothesized that the stress-like state of metabolic homeostasis observed in ovules of apomicts is a chronic condition that is maintained by imbalances in wild type gene expression that occur as a result of hybridization, polyploidization, or other chromosome aberrations [13,67,68]. We further hypothesized that this chronic stress response state enhances stress tolerance, and it is this enhanced stress tolerance that simultaneously suppresses sex while inducing apomixis. Accordingly, we searched our DEGs for candidates that might be responsible for shifting metabolic homeostasis toward a chronic dual state of stress-response but even stronger stress- Figure 6. ABA signaling associated differentially expressed genes (DEGs) observed between diplosporous Boechera × formosa and sexual Boechera stricta (see Tables S5 and S9). Superscripts identify taxa in which DEGs are upregulated; F, B. × formosa; S, B. stricta; red and green highlights, up and down-regulation in B. stricta, respectively; arrows indicate upregulation or activation; blunt end lines represent downregulation or suppression.
A very different state of cellular homeostasis exists during stress than during favorable conditions. During stress, ABA-inactivated PP2Cs permit the self-activation by autophosphorylation of SnRK2,3s as well as BIN2. Activated BIN2 phosphorylation inactivates the growth-enhancing BR TFs but phosphorylation activates the NAC TF RESPONSIVE TO DESICCATION 26 (RD26), which upregulates multiple stress response genes [53,64]. RD26 was upregulated 4.2 fold in B. × formosa vs. B. stricta (Figure 5d, Table S9). Activated SnRK2,3s upregulate ABA specific stress response genes [65] (Figure 6). These data suggest that apomictic ovules experienced states of cellular homeostasis that were maintained by enhanced stress perception with feedback suppression of ABA signaling occurring through CYP85A2, ABI1, and PP2CA [53,62,64,65].
The molecular pathways presented to this point suggest that apomeiosis is caused by carbon starvation (Figure 4a, low-level TOR activity), ABA (Figure 6), and upregulated drought response genes ( Figure 5d). However, these pathways do not explain how a stress-like state of metabolic homeostasis induces apomixis, especially when stressed apomicts in many cases revert to sex [15,66]. Recognizing this paradox, we hypothesized that the stress-like state of metabolic homeostasis observed in ovules of apomicts is a chronic condition that is maintained by imbalances in wild type gene expression that occur as a result of hybridization, polyploidization, or other chromosome aberrations [13,67,68]. We further hypothesized that this chronic stress response state enhances stress tolerance, and it is this enhanced stress tolerance that simultaneously suppresses sex while inducing apomixis. Accordingly, we searched our DEGs for candidates that might be responsible for shifting metabolic homeostasis toward a chronic dual state of stress-response but even stronger stress-attenuation (Figures 6 and 7). We first considered the possibility that the biosynthesis of ABA itself, which is upregulated by stress [69], might be chronically upregulated in the apomicts. ABA DEFICIENT 2 (ABA2), an ABA biosynthesis gene, was upregulated in aposporous B. microphylla vs. sexual B. stricta. However, this gene was also upregulated in B. microphylla vs. B. x formosa (Table S9), which is inconsistent with ABA synthesis being a major contributor. Furthermore, the ABA and salt stress-induced TF DIVARICATA 2 (DIV2), which negatively regulates ABA biosynthesis by downregulating ABA1,2 [70], was upregulated in B. x formosa vs. B. stricta (Table S9). Hence, a chronic stress-like phenotype being caused by upregulated ABA biosynthesis is unlikely.
Decreased catabolism of ABA is another possibility. Here, a drought-upregulated histone deacetylase gene, HISTONE DEACETYLASE 9 (HDA9), was upregulated 1.9 fold in B. x formosa vs. B. stricta ( Figure 6). HDA9 forms a complex with ABA INSENSITIVE 4 (ABI4) that downregulates the major ABA catabolism genes CYTOCHROME P450, FAMILY 707, SUBFAMILY A, POLYPEPTIDE 1,2 (CYP707A1,2) [71]. Moreover, the salt stress-induced PROTEIN WITH THE RING DOMAIN AND TMEMB_185A 1 (PPRT1) gene [72], a C3HC4 zinc-finger ubiquitin E3 ligase that upregulates the ABA catabolism genes CYP707A1,3 [73], was downregulated in B. x formosa vs. B. stricta. Hence, the up and downregulation of HDA9 and PPRT1, respectively, in B. x formosa likely contributed to its stronger ABA signaling. Hence, from our expression profiles, suppressing ABA catabolism is a reasonable candidate for upregulating ABA signaling and inducing a homeostatic state of stress attenuation.
Upregulating the biosynthesis of ABA receptor genes is another way in which ABA signaling in B. x formosa might have been enhanced. Normally, stress-induced ABA feedback inhibits ABA receptor gene biosynthesis, and this may have occurred for the ABA receptor genes PYRABACTIN RESISTANCE 1-LIKE 5,6 (PYL5,6) (downregulated 67% in B. x formosa vs. B. stricta; PYL5 was also downregulated 67% in B. lignifera and B. microphylla vs. B. stricta; Table S9). However, the downregulation of ABA receptor genes by feedback mechanisms is suppressed by the ABA-INDUCED TRANSCRIPTION REPRESSOR 2 gene (AITR2) [74,75], which was upregulated 2.1-fold in B. x formosa vs. B. stricta. Additionally, the C2-DOMAIN ABA-RELATED 9 gene (CAR9), an ABA receptor ( Figure 6), was upregulated 4.3-fold in B. x formosa vs. B. stricta. CAR9 positively regulates ABA signaling by interacting with ABA and various PYRABACTIN RESISTANCE (PYR) and PYL ABA receptors at the plasma membrane [76]. Furthermore, the ABA receptor gene GPCR-TYPE G  Tables S5 and S9). Superscripts identify taxa in which DEGs were upregulated; F, B. × formosa; S, B. stricta; red and green highlights indicate up and down-regulation in B. stricta, respectively; arrows indicate upregulation or activation; blunt end lines represent downregulation or suppression.
We first considered the possibility that the biosynthesis of ABA itself, which is upregulated by stress [69], might be chronically upregulated in the apomicts. ABA DEFICIENT 2 (ABA2), an ABA biosynthesis gene, was upregulated in aposporous B. microphylla vs. sexual B. stricta. However, this gene was also upregulated in B. microphylla vs. B. × formosa (Table S9), which is inconsistent with ABA synthesis being a major contributor. Furthermore, the ABA and salt stress-induced TF DIVARICATA 2 (DIV2), which negatively regulates ABA biosynthesis by downregulating ABA1,2 [70], was upregulated in B. × formosa vs. B. stricta (Table S9). Hence, a chronic stress-like phenotype being caused by upregulated ABA biosynthesis is unlikely.
Upregulating the biosynthesis of ABA receptor genes is another way in which ABA signaling in B. × formosa might have been enhanced. Normally, stress-induced ABA feedback inhibits ABA receptor gene biosynthesis, and this may have occurred for the ABA receptor genes PYRABACTIN RESISTANCE 1-LIKE 5,6 (PYL5,6) (downregulated 67% in B. × formosa vs. B. stricta; PYL5 was also downregulated 67% in B. lignifera and B. microphylla vs. B. stricta; Table S9). However, the downregulation of ABA receptor genes by feedback mechanisms is suppressed by the ABA-INDUCED TRANSCRIPTION REPRESSOR 2 gene (AITR2) [74,75], which was upregulated 2.1-fold in B. × formosa vs. B. stricta.
Downregulating proteins that target ABA receptors for degradation, could also enhance ABA signaling. Here, the DET1 AND DDB1 ASSOCIATED 1 gene (DDA1), which downregulates ABA signaling by targeting PYL4,8,9 for degradation [78], was downregulated in all three apomicts (Table S9). Hence, DDA1 is a reasonable candidate for shifting metabolic homeostasis toward ABA signaling. In general, the majority of ABA receptors were, for various reasons, upregulated in apomictic B. × formosa ( Figure 6).
Many of the ABA response DEGs ( Figure 6) are also upregulated by heat, salt, or osmotic stress ( Table S9), suggesting that hypersensitivity to these conditions may have also contributed to the chronic stress induction/tolerance homeostasis observed in B. × formosa. Many of these were late embryogenesis abundant (LEA) genes, the proteins of which protect cells and proteins from the deleterious effects of ROS, drought, and freezing [87].  [91]; and the protein-serine/threonine phosphatase gene AP2C1, upregulated 12.1 and 6.1-fold in B. × formosa vs. B. stricta and B. lignifera, respectively. AP2C1 modulates K + homeostasis by dephosphorylation inactivating CBL-INTERACTING PROTEIN KINASE 9 (CIPK9), which normally enhances K + ion uptake through K + channels, particularly under low K + conditions [92]. It also moderates responses to pathogens by dephosphorylation inactivating MAP KINASE 4,6 (MPK4,6) in the cytoplasm and nucleus. These MAPKs affect the expression of hundreds of genes by modifying the activity of 88 TFs from 21 TF families [93].
ABA signaling upregulates ROS synthesis, ROS signaling, and large suites of peroxisomal and other antioxidant network genes [91,94]. The first step in ROS signaling is an upregulation of ROS biosynthesis by peroxidases or RBOHs. Five ROS biosynthesis genes were highly upregulated in the TDG forming B. × formosa. PEROXIDASE CA (PRXCA) (Figure 7) encodes a cell wall peroxidase that produces H 2 O 2 in the apoplast [95]. It was upregulated in B. × formosa 11.4, 10.0, and 8.5-fold vs. the meiosis exhibiting B. stricta, B. microphylla, and B. lignifera, respectively (Table S9). NEET is a ROS homeostasis regulator that was upregulated in B. × formosa 7.8, 4.9, and 10.7-fold vs. B. stricta, B. microphylla, and B. lignifera. AT-NEET (NEET) encodes a protein that regulates Fe and ROS homeostasis by controlling the biogenesis of Fe-S clusters that function in electron transport. In the absence of NEET, excess Fe-S clusters form and react with oxygen to produce ROS [96]. PROLINE-RICH RECEPTOR-LIKE PROTEIN KINASE 4 (PERK4) enhances ROS accumulation by enabling ABA-induced expression of RBOH genes, including RBOHC [97]. PERK4 was upregulated in B. × formosa 6.7, 5.6, and 6.5-fold vs. B. stricta, B. microphylla, and B. lignifera. RBOHD produces superoxide in response to heat and wounding and its gene was upregulated in B. × formosa 2.9, 5.1, and 8.2-fold vs. B. stricta, B. microphylla, and B. lignifera. Both RBOHD and RBOHF are upregulated by ABA, and their expression, through ROS signaling, is critical to cell survival under ER stress [98]. RBOHF was upregulated in B. × formosa 2.9, 4.2, and 3.5-fold vs. B. stricta, B. microphylla, and B. lignifera.

Metabolically Induced Switching from Apomeiosis to Meiosis
We hypothesized above that apomicts express chronic stress responses coupled with chronic, overcompensating stress attenuation mechanisms that suppress meiosis and induce apomeiosis. If this is correct, meiosis should be inducible in apomicts by applying stresses that exceed the apomict's stress attenuation capabilities. In 1951, Bocher [117] reported shifts from 1st division restitutions in microsporocytes of apomictic Boechera, which produced unreduced pollen, to complete meioses in weak plants that flowered sporadically late in the season. He noted that the processes that produce unreduced male and female spores, 1st division restitution, occurred in anthers and ovules. He also suggested that stress had induced the observed shifts from apomeiosis to meiosis, late in the season, and he speculated that stress might affect ovules in the same way. We viewed Bocher's speculations as an opportunity to gain further insights into these processes of sex/apomixis switching in facultative apomicts. Accordingly, we exposed randomly selected groups of Boechera plants to three growth conditions: well-watered, droughted, and droughted plus heat. Taxa studied were sexual B. stricta, facultative TDG forming B. lignifera, triploid, and highly TDG forming B. cf. gunnisoniana, HAG forming B. microphylla from Millard County Utah (meiotic and aposporous, but not diplosporous), and facultative TDG forming B. exilis × retrofracta [16]. The combination of drought and heat increased meiosis frequencies (tetrad formation) from 15 to 75% in B. lignifera and from 2 to 77% in B. cf. gunnisoniana (Figure 8a). Since most dyads in these droughted and heat-stressed plants (23 to 25% of ovules, Figure 8a) were possibly produced by meiosis, actual frequencies of sexual meioses in these stressed TDG-forming apomicts may have been closer to 90%. Drought alone in TDG forming B. lignifera and B. cf. gunnisoniana was less effective in switching apomeiosis to meiosis (Figure 8a). Disturbed meioses that produced pentads and hexads of imbalanced megaspores also occurred among the stress-induced meioses of triploid B. cf. gunnisoniana (Figure 8c,d). Likewise, Bocher also observed microspore pentads and hexads in the anthers of late-flowering Boechera triploids [117].
Due to lethality, data were not obtained for the double stressed (drought and heat) B. stricta, B. exilis × retrofracta, and B. microphylla. Drought alone had minor effects on the reproductive mode in sexual B. stricta and facultative B. exilis × retrofracta. The sexual dyad stage of meiosis is brief and difficult to catch cytologically. Accordingly, it is possible that drought accelerated the dyad stage (e.g., more rapid onset of M II ) such that fewer dyads were observed in ovules of B. stricta (Figure 8a). Drought also shifted frequencies of gametophyte formation from aposporous to sexual in the Millard County B. microphylla (Figure 8b). Figure 8e-g provides examples of sexual tetrads, diplosporous dyads, and early onset stages of aposporous gametophyte formation as observed in this experiment.
Our greenhouse experiment exposed plants to chronic drought or drought plus heat for the 10-40 days required for flowering to occur following vernalization. Accordingly, we asked whether long-term chronic stress is required to switch apomeiosis to meiosis or whether this switch might occur rapidly, i.e., following acute stresses applied to pistils within days or even hours of apomeiosis. To answer this question, we developed a simple MS-based [26] pistil culture protocol that supported ovule development in cultured pistils from the late pre-MMC stage (1-I) through the enlarging MMC (2-II), mature MMC (2-III), meiosis (2-IV), late meiosis to functional megaspore (2-V), and early vacuolate multinucleate embryo sac (3-III) stages. Ovules of pistils cultured on basal medium (without pharmacological additives) progressed through these developmental stages within 48-72 h of culture initiation. However, development beyond the early gametophyte stage was irregular, in controls and treatments, which prevented us from evaluating in vitro the effects of metabolic modifications on later stages of gametophyte development and seed formation.
The TDG forming triploid B. cf. gunnisoniana responded strongly to chronic drought (Figure 8a). Hence, we chose it for acute drought experiments that involved culturing 1.5-1.9 mm long pistils, each containing 40-50 ovules at the early to late MMC stage, on low water potential media (Figure 9a). Even on control media (−0.50 MPa), frequencies of meiotic tetrads were elevated above the 2-3% background level. Apparently, excising and culturing these young desiccation sensitive pistils was itself stressful enough to switch apomeiosis to meiosis in some ovules, presumably by overwhelming stress attenuation processes in ovules of the cultured pistils. Meiotic tetrad frequencies were highest (~50%) in the −0.56 MPa treatment, and the highest tetrad frequencies occurred in pistils that contained late MMC staged ovules at the time of culture initiation. Hence, the duration of drought exposure may be less important than drought exposure occurring at the late MMC stage when meiosis or apomeiosis fates are determined. Ovules with less mature MMCs at the time of culture initiation possibly had more time in culture to restore metabolic homeostasis. Importantly, this experiment revealed that extended periods of stress are not required to reprogram genomes for meiosis. Instead, this reprogramming occurred within hours of stress induction. stressed TDG-forming apomicts may have been closer to 90%. Drought alone in TDG forming B. lignifera and B. cf. gunnisoniana was less effective in switching apomeiosis to meiosis (Figure 8a). Disturbed meioses that produced pentads and hexads of imbalanced megaspores also occurred among the stress-induced meioses of triploid B. cf. gunnisoniana (Figure 8c,d). Likewise, Bocher also observed microspore pentads and hexads in the anthers of late-flowering Boechera triploids [117]. Due to lethality, data were not obtained for the double stressed (drought and heat) B. stricta, B. exilis x retrofracta, and B. microphylla. Drought alone had minor effects on the reproductive mode in sexual B. stricta and facultative B. exilis x retrofracta. The sexual dyad stage of meiosis is brief and difficult to catch cytologically. Accordingly, it is possible that drought accelerated the dyad stage (e.g., more rapid onset of MII) such that fewer dyads were observed in ovules of B. stricta (Figure 8a). Drought also shifted frequencies of gametophyte formation from aposporous to sexual in the Millard County B. microphylla (Figure 8b). Figure 8e-g provides examples of sexual tetrads, diplosporous dyads, and early onset stages of aposporous gametophyte formation as observed in this experiment.
Our greenhouse experiment exposed plants to chronic drought or drought plus heat for the 10-40 days required for flowering to occur following vernalization. Accordingly, we asked whether longterm chronic stress is required to switch apomeiosis to meiosis or whether this switch might occur rapidly, i.e., following acute stresses applied to pistils within days or even hours of apomeiosis. To answer this question, we developed a simple MS-based [26] pistil culture protocol that supported ovule development in cultured pistils from the late pre-MMC stage (1-I) through the enlarging MMC (2-II), mature MMC (2-III), meiosis (2-IV), late meiosis to functional megaspore (2-V), and early vacuolate multinucleate embryo sac (3-III) stages. Ovules of pistils cultured on basal medium (without pharmacological additives) progressed through these developmental stages within 48-72 h Comparisons of the more apomeiotic vs. the more meiotic taxa revealed reduced molecular support for translation in ovules of apomicts (Figure 2b,f,j and Figure 4a). Moreover, these ovules were more metabolically active, were responding to ABA, and were more actively attenuating stress (Figures 6 and 7). Hence, it appears that our chronic and acute drought stresses overwhelmed the stress attenuation capabilities of these apomicts and switched apomeiosis to meiosis. We next asked whether severe carbon starvation, which generates oxidative stress and terminates translation [118], might also increase apomeiosis to meiosis conversions. To address this question, we cultured pistils of B. cf. gunnisoniana (containing early to late-stage MMCs) on basal medium either with or without sucrose. Again, younger pistils on the control medium (with sucrose) produced fewer tetrads than more mature pistils (Figure 9b, compare 2-III to 2-IV results). Presumably, this occurred because cells of the younger 2-III staged pistils had more time to adjust their stress shocked metabolism to pre-excision levels prior to reaching the point of stress tolerance that defines facultativeness (apomeiosis/meiosis switching) in apomicts. In the no sucrose treatment, TDG (1st division restitution dyad) frequencies decreased by ca. 50%. This was matched by corresponding increases in meiotic tetrads, but it was also matched by HAG formation onset (Figure 9g). Frequencies of both phenomena approached 50% in ovules of the more mature pistils (Figure 9b). Hence, under severe carbon starvation, apomeiosis switched to meiosis, but the meioses in these ovules were often accompanied by HAG formation (apospory) with little or no net decrease in apomixis. We next tested whether oxidative stress (exogenously applying H 2 O 2 ) would induce apomeiosis to meiosis conversions as well as diplospory to apospory conversions, i.e., switching nucellar cell programmed cell death (the norm) to HAG formation. For this experiment, we used pistils (with ovules at the MMC stage) from three apomicts, the facultative TDG forming B. lignifera (1.0-1.4 mm long pistils), the nearly obligate TDG forming B. cf. gunnisoniana (1.5-1.9 mm long pistils), and the facultative TDG forming B. retrofracta × stricta (1.6-2.1 mm long pistils). Pistils were excised and soaked for 5 min in H 2 O 2 solutions prior to culture (Figure 9c). These H 2 O 2 pretreatments induced up to 60% apomeiosis to meiosis conversions. However, all megaspores from many of the induced meiotic tetrads degenerated, and in many cases, this degeneration was accompanied by HAG formation (Figure 9h-i). These conversions were again most prominent in ovules of more mature pistils, suggesting that younger ovules had more time (ca. 12-24 h) to detoxify the exogenously added H 2 O 2 before committing themselves to very different developmental paths.
The TDG forming triploid B. cf. gunnisoniana responded strongly to chronic drought (Figure 8a). Hence, we chose it for acute drought experiments that involved culturing 1.5-1.9 mm long pistils, each containing 40-50 ovules at the early to late MMC stage, on low water potential media ( Figure  9a). Even on control media (−0.50 MPa), frequencies of meiotic tetrads were elevated above the 2-3% background level. Apparently, excising and culturing these young desiccation sensitive pistils was itself stressful enough to switch apomeiosis to meiosis in some ovules, presumably by overwhelming stress attenuation processes in ovules of the cultured pistils. Meiotic tetrad frequencies were highest (~50%) in the −0.56 MPa treatment, and the highest tetrad frequencies occurred in pistils that contained late MMC staged ovules at the time of culture initiation. Hence, the duration of drought exposure may be less important than drought exposure occurring at the late MMC stage when meiosis or apomeiosis fates are determined. Ovules with less mature MMCs at the time of culture initiation possibly had more time in culture to restore metabolic homeostasis. Importantly, this experiment revealed that extended periods of stress are not required to reprogram genomes for meiosis. Instead, this reprogramming occurred within hours of stress induction. Figure 9. Effects of pharmacological treatments and pistil maturity on conversion from Taraxacumtype gametophyte formation (TDG) to (i) meiosis, (ii) Hieracium-type gametophyte formation (HAG), and iii) Antennaria-type gametophyte formation (ADG): (a) drought (media water potential); (b) energy starvation; (c) reactive oxygen species stress (H2O2); (d) abscisic acid (ABA); (e) brassinazole (BRZ, brassinosteroid synthesis inhibitor); (f) 5-azacytidine (DNA methylation inhibitor); (g) sucrose starved Boechera cf. gunnisoniana ovule with a Taraxacum-type diplosporous dyad (black arrows) and three 1 nucleate vacuolate HAGs (red arrows; v, vacuoles); (h-i) H2O2 induced conversions from TDG formation to meiosis, but with all four megaspores degenerating (between black arrows), and from The experiments reported so far illustrate the importance of plants maintaining highly stable states of metabolic homeostasis that restrict reproduction to a single sexual or apomictic type, and they suggest that apomixis occurs when oxidative stress is overly attenuated (Figure 7). Continuing with this theme, we asked whether upregulating stress attenuation by increasing ABA signaling (Figures 6 and 7) might prevent meiotic prophase and shift TDG formation to ADG formation in B. cf. gunnisoniana (Figure 9d). Such treatments induced ADG formation in 25-35% of ovules (gametophyte formation directly from MMCs; Figure 9j-k). These findings place oxidative stress and oxidative stress attenuation at the center of reproductive decisions in Boechera, with strong and moderate stress attenuation causing ADG and TDG formation and oxidative stress causing nucellar cell destabilization followed by HAG formation.
Our expression profiling studies suggested that ABA signaling in the TDG forming B. × formosa reduces BR signaling at BIN2 (Figure 5d). To further evaluate this, we cultured TDG forming B. cf. gunnisoniana pistils, with ovules at the MMC stage, on media containing the BR biosynthesis inhibitor brassinazole (BRZ). This more thorough reduction in BR signaling increased frequencies of meiosis and apospory (Figure 9e), possibly by further suppressing mechanisms of oxidative stress attenuation.
Meiosis is a cell-lineage-specific process that requires a meiosis-specific transcriptome that includes meiosis-specific de novo RNA-directed DNA methylation (RdDM) [49]. Here, we asked whether disrupting RdDM using 5-azaC, a methyltransferase inhibitor that reduces DNA methylation [119], might prevent meiotic prophase in TDG forming B. cf. gunnisoniana but still permit gametophyte formation directly from MMCs (ADG formation). Though not as efficient as ABA, 5-azaC induced ADG formation at frequencies that were correlated with 5-azaC concentrations (Figure 9f). That both ABA and 5-azaC induce ADG formation (Figure 9d,f,j-k) is consistent with ABA-induced oxidative stress attenuation silencing meiosis-specific RdDM in germline cells (Figures 4-7).

Metabolically Induced Switching from Meiosis to Apomeiosis
We demonstrated above that the apomicts studied herein chronically express stress attenuation mechanisms, but when these mechanisms are overwhelmed, apomeiosis reverts to meiosis. Here, we show that the reverse also occurs, i.e., when stress attenuation mechanisms are enhanced in sexual plants, meiosis reverts to apomeiosis. For these experiments, we chose four sexual diploid species: A. thaliana, B. stricta, and B. exilis of the Brassicaceae family and V. unguiculata of the Fabaceae family. Pistils containing MMC staged ovules were used. Pistil lengths were 0.5-1.2, 1.3-1.9, 1.5-2.3, and 1.5-1.9 mm, respectively. Dioecious A. dioica (Asteraceae) was used to study microsporogenesis.
Treatments applied to excised A. thaliana pistils included a 60 min pretreatment emersion. In the controls (Figure 10a), percentages of tetrads to dyads were as expected. However, 17% of control ovules (presoaked for 60 min) contained vacuolate HAGs. This response may have been caused by anoxia induced ROS (e.g., Figure 9c). To test this, we pre-soaked sets of pistils for 7 min in water or 100 mmol L −1 H 2 O 2 and then cultured them on basal medium or medium amended with 1.0 µmol L −1 epiBL. Water alone did not induce HAG formation, but such formation did occur in 12-15% of ovules pretreated with H 2 O 2 (Figure 10b). Diplospory was absent in these treatments. Hence, exogenous H 2 O 2 may have reinforced meiosis while destabilizing nucellar cells and either inducing or allowing aposporous development. The TF BZR1, which is activated by endogenous H 2 O 2 [120], may have also been involved. These experiments indicate that nucellar cells in ovules of diplosporous and sexual plants respond similarly to H 2 O 2 treatments in terms of inducing apospory.
To determine whether apomeiosis is inducible by reducing oxidative stress, we excised early MMC staged A. thaliana pistils, soaked them in the antioxidant (S)-2-aminobutane-1,4-dithiol hydrochloride (DTBA) for 60 min, and cultured them on DTBA containing medium for 4 days. In 15 of 41 scored ovules, a first division meiotic restitution occurred, which produced a dyad of megaspores (Figure 10a). The micropylar dyad member in some of these was degenerating and the chalazal dyad member was becoming vacuolate and undergoing endomitosis to form a multinucleate gametophyte (Figure 10c,d).
An additional 7 of 41 scored ovules completed meiosis, but all four megaspores in these ovules were degenerating and being replaced by HAGs (Figure 10e-f). controls (Figure 10a), percentages of tetrads to dyads were as expected. However, 17% of control ovules (presoaked for 60 min) contained vacuolate HAGs. This response may have been caused by anoxia induced ROS (e.g., Figure 9c). To test this, we pre-soaked sets of pistils for 7 min in water or 100 mmol L −1 H2O2 and then cultured them on basal medium or medium amended with 1.0 μmol L −1 epiBL. Water alone did not induce HAG formation, but such formation did occur in 12-15% of ovules pretreated with H2O2 (Figure 10b). Diplospory was absent in these treatments. Hence, exogenous H2O2 may have reinforced meiosis while destabilizing nucellar cells and either inducing or allowing aposporous development. The TF BZR1, which is activated by endogenous H2O2 [120], may have also been involved. These experiments indicate that nucellar cells in ovules of diplosporous and sexual plants respond similarly to H2O2 treatments in terms of inducing apospory.  We also asked whether BR alone would induce apomeiosis. In 24 of 50 ovules, BR induced a 1st division restitution and subsequent TDG formation. Low-level HAG formation (8%) also occurred. Adding an H 2 O 2 pretreatment to BR-treated pistils prevented TDG formation possibly by reinforced meiosis (Figure 10b). Nevertheless, the combination of BR and H 2 O 2 supported post-meiotic abortion coupled with HAG formation. We next tested whether 5-azaC would terminate meiosis and induce ADG formation in A. thaliana as it had done in B. cf. gunnisoniana (Figure 9d,j,k). Here, we soaked pre-MMC staged A. thaliana pistils in 5-azaC for 60 min and cultured them on control or 5-azaC medium. In 5-azaC treated pistils, 14 of 106 appropriately staged ovules were forming ADGs (Figure 10a,g,h). Of all treated pistils (Figure 10a), 64% (180 of 282) initiated either TDG, ADG, or HAG formation, i.e., all three major types of apomeiosis were induced in A. thaliana by metabolic modifications.
To determine if apomeiosis is inducible in other angiosperms, we exposed pistils of sexual B. stricta and sexual V. unguiculata to treatments similar to those used for A. thaliana (Figure 11a,b). The A. thaliana and B. stricta treatments were performed in duplicate, but with glucose and sucrose used individually as energy sources. The sugar main effect was not significant in A. thaliana, but TDG formation (Figure 11c-d) increased significantly in B. stricta when glucose was used (38 of 152 ovules vs. 26 of 164 ovules). ADG formation was only rarely observed in B. stricta but was induced by all treatments in V. unguiculata (Figure 11b,j-k). The V. unguiculata study included 74 pistils, but only 34 ovules were at the MMC stage at culture initiation. Of these 34, 33 had been exposed to apomeiosis-inducing treatments, of which 24 produced ADGs (73%). One MMC staged ovule had been exposed to the sucrose control, and it was meiotic. formation (Figure 11c-d) increased significantly in B. stricta when glucose was used (38 of 152 ovules vs. 26 of 164 ovules). ADG formation was only rarely observed in B. stricta but was induced by all treatments in V. unguiculata (Figure 11b,j-k). The V. unguiculata study included 74 pistils, but only 34 ovules were at the MMC stage at culture initiation. Of these 34, 33 had been exposed to apomeiosisinducing treatments, of which 24 produced ADGs (73%). One MMC staged ovule had been exposed to the sucrose control, and it was meiotic. ABA signaling includes oxidative stress induction by respiratory burst oxidase homologs and superoxide dismutases and subsequent oxidative stress attenuation by antioxidants. The evolution of this compensatory stress tolerance process facilitated the evolution of land plants. Here, we wondered if terminating ABA synthesis might induce apomeiosis by downregulating PP2C production, hence upregulating SnRK2,3 activities and the subsequent production of stress attenuation mechanisms (Figures 6 and 7). To evaluate this, we exposed pistils of sexual B. exilis to 5.0 µmol L −1 fluridone treatments for 24-48 h before meiosis initiation (early to the mid-MMC stage). Fluridone, which inhibits ABA synthesis [121,122], induced 17%, 22%, and 2% of ovules to undergo ADG, HAG, and TDG formation, respectively (Figure 12a-g). In some ovules, ADG and HAG formation occurred simultaneously (Figure 12d-f). Evolutionarily, the nucellus is sporangial, like fern sporangia [5]. By altering ABA signaling, nucellar cells may have been released from programmed cell death, their normal fate, to pursue an ancestral fate, that of gametophyte formation.
We also tested whether chemicals that induce apomeiosis in megasporocytes might also induce apomeiosis in microsporocytes. Here we placed cut stems of male dioecious A. dioica (Asteraceae) in liquid MS media amended with fluridone. As was observed for MMCs, the ameiotic PMCs developed precociously into unreduced male gametophytes (Figure 12h-j). Note that the four-fold meiotically replicated and expanding tetrads of microspores crowded the anther locule in the controls (Figure 12i). In contrast, the meiotically skipped, precociously forming unreduced pollen grains of the fluridone treated plants were widely spaced throughout the locule (Figure 12j). 5.0 μmol L −1 fluridone treatments for 24-48 h before meiosis initiation (early to the mid-MMC stage). Fluridone, which inhibits ABA synthesis [121,122], induced 17%, 22%, and 2% of ovules to undergo ADG, HAG, and TDG formation, respectively (Figure 12a-g). In some ovules, ADG and HAG formation occurred simultaneously (Figure 12d-f). Evolutionarily, the nucellus is sporangial, like fern sporangia [5]. By altering ABA signaling, nucellar cells may have been released from programmed cell death, their normal fate, to pursue an ancestral fate, that of gametophyte formation. (j) multinucleate and vacuolate ameiotic 1-3 nucleate pollen grains that formed directly from microsporocytes in A. dioica anthers from stems that were cultured in fluridone amended medium. Numbers of informative ovules or anther meiocytes analyzed appear next to treatment bars in (a) and (h). Pooled chi-squares were significant (p < 0.001). red arrows, degenerating haploid megaspores; n, nucleus; v, vacuole; cytology bars, 20 μm.
We also tested whether chemicals that induce apomeiosis in megasporocytes might also induce apomeiosis in microsporocytes. Here we placed cut stems of male dioecious A. dioica (Asteraceae) in (i) normal microspore tetrads in anthers of control A. dioica plants; (j) multinucleate and vacuolate ameiotic 1-3 nucleate pollen grains that formed directly from microsporocytes in A. dioica anthers from stems that were cultured in fluridone amended medium. Numbers of informative ovules or anther meiocytes analyzed appear next to treatment bars in (a) and (h). Pooled chi-squares were significant (p < 0.001). red arrows, degenerating haploid megaspores; n, nucleus; v, vacuole; cytology bars, 20 µm.

Sucrose Non-Fermenting-Related Protein Kinases May Be Central to Meiosis Apomeiosis Switching
The evolution of land plants required new water retaining and transporting morphologies, new stress detection and signaling pathways, and new processes for catabolizing toxins induced by terrestrial stresses [123][124][125]. Here, the evolution of ABA and SnRK2,3 stress response pathways was critical for increasing ROS attenuation capabilities [124]. ROS form constantly as byproducts of respiration, photosynthesis, and the activities of peroxidases and peroxisomes. However, in land plants, ROS are also produced to propagate waves of ROS signals from stressed cells to more distant cells [126]. Two ABA-induced pathways accomplish this: (i) propagation of ROS signaling through tissues via RBOHs and SODs ( Figure 7) and (ii) upregulation of stress attenuation mechanisms through SnRK2,3 ABI5 signaling ( Figure 6). During such signaling, H 2 O 2 activates or inactivates redox-regulated protein phosphatases and other enzymes [127,128]. It is critical that excess H 2 O 2 is catabolized by ROS attenuation mechanisms [94,129,130]. Metabolic processes that establish such redox homeostasis function by balancing ROS generation, from multiple sources, with ROS catabolism. We discovered, and report herein, that apomeiosis occurs when this redox balance is chronically tilted toward H 2 O 2 catabolism, i.e., when SnRK2,3s and RD26 are activated ( Figure 13). We further show that meiosis apomeiosis interconversions are inducible by applying chemicals that acutely disturb these taxa specific states of redox homeostasis. catabolized by ROS attenuation mechanisms [94,129,130]. Metabolic processes that establish such redox homeostasis function by balancing ROS generation, from multiple sources, with ROS catabolism. We discovered, and report herein, that apomeiosis occurs when this redox balance is chronically tilted toward H2O2 catabolism, i.e., when SnRK2,3s and RD26 are activated ( Figure 13). We further show that meiosis apomeiosis interconversions are inducible by applying chemicals that acutely disturb these taxa specific states of redox homeostasis. In addition to the plant-specific SnRK2,3s (Figure 6), SnRK1 was also more active in the apomicts (Figure 4a). In eukaryotes, SnRK1 regulates growth no-growth decisions. When phosphorylated, SnRK1 dissociates the TORC1 complex [33,36,124,[131][132][133][134][135], and this suppresses translation and BR regulated cell proliferation [33,120,124,[136][137][138][139]. Translation downregulation in the more apomeiotic taxa (14 of 16 comparisons during early ovule development, Figure 2a,e-i and Figure 3a-g) indicates that TOR was less active in the apomicts due to increased SnRK1 activity. That stress or low energy conditions in the apomicts activated SnRK1, which suppresses TOR, is further shown by upregulated genes in the apomicts for fatty acid, amino acid, carbohydrate, and other macromolecule catabolism (Figure 2a,f,g,l, and Figure 3a-d,g).
In the apomicts, profiling revealed upregulated genes associated with signaling, stress responses, metabolism, photosynthesis, and bioenergetics (Figure 1a,e-j,l, and Figure 3a-g). Here, the B. × formosa population may have contained triploids, which is the common polyploid level among Boechera apomicts [17]. Polyploids often show distinct patterns of gene expression [140] as well as increased stress attenuation capabilities [15,141]. Increased metabolic activity seems counterintuitive given the low energy induced suppression of TOR in this taxon. However, the evolution of SnRK2,3s enabled land plants to avoid toxic levels of ROS, due to enhanced ROS attenuation processes, while maintaining growth and other metabolic processes [124,125]. Under favorable conditions, all three categories of SnRKs (SnRK1,2,3s) were largely phosphatase inactivated in the more meiotic plants by ABI1,2 and the HIGHLY ABA-INDUCED PP2C GENE 2 (AIP1). In contrast, these SnRKs appeared to remain partially active in the apomicts, which likely explains their greater ROS attenuation capacity. Switching from apomeiosis to meiosis (facultativeness) in response to stress (Figures 8 and 9) suggests that ROS accumulation, due to stress treatments, exceeded ROS catabolism ( Figure 13).

Metabolically Inducing Meiosis Apomeiosis Interconversions
Stress attenuation in plants involves the activities of ascorbic acid, ascorbate peroxidases, glutathiones, thioredoxins, catalases, superoxide dismutases (SODs), tocopherols, flavonoids, carotenoids, and ion homeostasis regulators. These are upregulated by oxidative stress, ABA-activated SnRKs, and ABA insensitive SnRKs (activated by salt, osmotic, or other stress) [124,128]. Genes encoding many of these were upregulated in ovules of the more apomeiotic taxa (Tables S2 and S5). In this respect, Shah et al. [142] provided evidence that chronic stress tolerances in apomictic Boechera might persist throughout the life cycle. They further suggested that this increased tolerance contributes to geographical parthenogenesis (superior stress tolerance in apomicts). We provide evidence herein that the chronic states of stress tolerance in apomicts suppress meiosis and induce apomeiosis. They might also suppress syngamy and induce parthenogenesis. This notion recently gained traction by the demonstration that parthenogenesis in rice is inducible by altering the location of expression of a wild type APETALA2/ERF gene [143], a member of a TF family that regulates development and stress by integrating redox, hormone, and environmental signals [144].
Using pharmacological approaches, we tested whether the differences in molecular pathways observed between apomeiotic and meiotic taxa (Figures 2-7) are responsible for inducing meiosis or apomeiosis or are simply coincidental. To do this, we identified chemicals that perturb specific elements of metabolic homeostasis. Here, apomeiosis switched to meiosis when the chosen treatments increased oxidative stress. These included drought, energy starvation, topical H 2 O 2 applications, and BR synthesis inhibition (Figure 9a-c,e; Figure 13a-d,f). In contrast, meiosis switched to apomeiosis when the chosen treatments relieved oxidative stress. These included glucose, BR, antioxidants, ABA, and the ABA synthesis inhibitor fluridone (Figures 10-12 and Figure 13e,g-n). Additionally, 5-azaC terminated meiosis in sexual plants (Figures 10 and 11), terminated meiotic prophase in TDG forming plants (Figure 9f), and permitted or induced ADG formation, directly from ameiotic MMCs, in sexual and TDG forming plants. This suggests that while RdDM may be required for meiosis [49,145], it is not required for the early stages of gametophyte formation.
When used as the sole carbon source, glucose blocked meiosis, and induced TDG formation in B. stricta and ADG formation in V. unguiculata (Figure 11b-d,j,k). In S. cerevisiae, glucose, which inhibits Snf1 kinase (SnRK1 in plants) and possibly activates TOR [146], also blocks meiosis [147]. Thus, in plants and yeast, meiosis is correlated with Snf1/SnRK1 inactivity and possibly TOR activity (Figure 4a). In plants, meiosis was further correlated with SnRK2,3 inactivity (reduced ABA signaling). In contrast, the more apomeiotic taxa expressed elevated levels of ABA signaling and correspondingly higher levels of antioxidant production (Figures 6 and 7). Since glucose in plants also upregulates ABA synthesis [148], glucose upregulated apomeiosis in B. stricta and V. unguiculata ( Figure 11) is consistent with upregulated ABA signaling followed by increased stress attenuation ( Figure 13). Surprisingly, carefully timed H 2 O 2 treatments not only converted apomeiosis to meiosis in TDG forming Boechera (apomixis to sex) but also induced HAG formation (sex to apomixis) in many of the same ovules (Figure 9c). H 2 O 2 also induced HAG formation in sexual A. thaliana (Figure 10b). We suggested that H 2 O 2 induces HAG formation directly. However, treating pistils with the antioxidant DTBA also induced HAG formation (Figure 10a,e,f). Hence, it is also possible that the exogenous H 2 O 2 treatments triggered an amplification of ROS attenuation mechanisms by upregulating ABA signaling (Figures 6 and 7). Upregulated ROS catabolism may have then freed nucellar cells from redox induced programmed cell death, thus allowing them to pursue gametophyte formation.
In addition to inducing ADG and HAG formation in sexual B. exilis (Figure 12a-g), fluridone also induced gametophyte formation directly from PMCs (Figure 12h-j). Since ABA upregulates PP2C formation, the suppression of ABA synthesis may have indirectly increased antioxidant production downstream of SnRK2,3s ( Figure 6). PMCs, MMCs, and nucellar cells share an evolutionary origin that includes homosporous meioses, e.g., as occur in sporangia of many ferns [12,149,150]. Hence, it is not surprising that metabolic disturbances that induce ADG formation from MMCs and HAG formation from nucellar cells (Figure 12a-g) can also induce gametophyte formation from PMCs (Figure 12h-j).
The effects of fluridone persisted throughout the 3-4 days duration of the experiment. Accordingly, most nucellar cells and some integument cells initiated HAG formation (Figure 12d-g).

Apomixis Genes: Horses of a Different Color
From the early 1900s to the 1940s, common views concerning the inheritance of apomixis stemmed from Strasburger's work, which suggested that certain sexual plants possess propensities for apomixis that are released by hybridization and polyploidization [4,151]. Starting in the late 1930s, genetic studies suggested that apomixis genes exist, which were later found to be isolated in nonrecombining genomic regions [10]. Aside from these more common views, a short-lived hypothesis (in the 1920s) was introduced by Ernst [152] who concluded that apomixis was caused by hybridization. Importantly, Ernst emphasized that apomixis propensities in parent plants are not required. In the 1990 s, Carman [13] comprehensively reviewed the genomics of known apomicts and reported highly significant correlations between apomixis and hybridization, polyploidization, and a relative absence of recent paleopolyploidization (low chromosome base numbers). Carman also concluded that hybridization was singularly important to apomixis evolution and that apomixis specific genes are not necessary. Ernst thought that phylogenetic distance between parental species determined apomixis [152]. In contrast, Carman cited physiological and biochemical examples of asynchronously expressed duplicate genes in hybrids and hypothesized that apomixis is caused by asynchronous expression of germline development genes in the ovules of certain hybrids.
The asynchrony hypothesis suggests that ecotype-driven divergences in wild type genes occurred during speciation along ecological gradients and that these new alleles modified ovule development timing. Following secondary contact hybridization, germline asynchronies were resolved in certain hybrids by eliminating meiosis and syngamy, thus producing apomixis [13]. The present paper updates this hypothesis. Here, asynchronies in germline development are replaced by dosage-dependent disparities in metabolic homeostasis ( Figure 13). In terms of genetic mapping, this revision simply shifts the types of genes expected to be found in apomixis associated genomic regions from meiosis/syngamy genes to metabolic homeostasis genes. Interestingly, apomixis is much more common in land plants than in many other eukaryotes. The addition of a major ROS attenuation mechanism (ABA/SnRK2,3 signaling) possibly made apomixis-conducive states of redox homeostasis more readily attainable through hybridization and polyploidization.
It is important to address the question as to why "occasional apomixis" in sexual plants-sensu Asker and Jerling [4]-is rare. The simplest explanation is that the chronic states of stress attenuation responsible for apomixis ( Figure 13, left side) are not environmentally inducible in sexual plants. Instead, such states of homeostasis are achieved only rarely through genomic aberrations associated with hybridization, recombination, or polyploidization. In contrast, states of stress attenuation in natural apomicts, induced by genomic aberrations, are variable ( Figure 13, right side) such that some plants are weakly apomictic while others remain nearly obligately apomictic even when stressed.
Several interesting examples of homeostatic states of chronic stress attenuation, which possibly enforce apomixis, exist in apomicts of other kingdoms. Apomeiosis in bdelloid and monogonont rotifers (Metazoa) is mitotic-like, and the 2n eggs develop into embryos parthenogenetically. Bdelloid rotifers are obligately apomictic, and they exhibit constitutive anti-oxidant protection provided by extensively duplicated anti-oxidant genes [153,154]. Generally, in cyclical animal apomicts, it is the onset of overwhelming stress, from various factors, that cause apomictic females to apomictically produce clonal but sexually functional males and females (genders induced environmentally). These breed to produce resting eggs that hatch into apomictic females. This pattern is observed in monogonont rotifers, water fleas (Daphnia), aphids (Arthropoda, Metazoa) [155,156], diatoms (Chromista) [157], and green alga (Plantae) from the genera Volvox (Plantae). As in bdelloid rotifers, a major class of oxidative stress detoxifying enzymes, glutathione S-transferases, are upregulated during the apomictic female phase in the pea aphid Acyrthosiphon pisum [156]. Major oxidative stresses are also required to induce the switch from apomictic to sexual reproduction in Daphnia [155] and Volvox [158][159][160], and to induce apomixis to sex conversions in facultatively apomictic slime molds (Protozoa) [161][162][163] and brewer's yeast (Fungi) [3]. In our experiments, attenuating oxidative stress in sexual sporocytes of Boechera, Arabidopsis, Vigna, and Antennaria caused meiosis to switch to apomeiosis. In contrast, when we overwhelmed ovules of facultatively to nearly obligately apomictic Boechera with oxidative stress, apomeiosis switched to meiosis. If our interpretations are correct, states of metabolic homeostasis that induce apomeiosis in extant eukaryotes are generally rare, are characterized by strong anti-oxidant properties, and are usually caused by genomic aberrations that grossly disturb the expression of wild type genes.

Conclusions
Switching from apomeiosis to meiosis in apomictic plants and from meiosis to apomeiosis in sexual plants was readily achieved through appropriate treatments. However, the treatment effects dissipated rapidly, with taxon-specific levels of cellular homeostasis apparently returning quickly to pretreatment levels. For example, the effectiveness of treatments applied 24-48 h before apomeiosis were diminished by up to 90% compared to treatments applied only a few hours prior to apomeiosis ( Figure 9). Accordingly, achieving levels of cellular homeostasis that assure complete penetrance of apomeiosis (and parthenogenesis) will likely require genetic engineering modifications that acutely alter cellular levels of metabolic homeostasis throughout ovule and early seed development. Candidate genes and pathways for achieving such modifications are described herein.
Our studies shed light on two important barriers that have frustrated the development of apomixis as a technology for stabilizing hybrid vigor in crops: (i) a correct understanding of what apomixis is, and (ii) a correct understanding of what biological properties trigger meiosis/apomeiosis interconversions. Concerning the first, we have demonstrated that apomeiosis does not require mutations, meiotic or otherwise, but is instead a polyphenism of meiosis, a conserved polyphenism that may be as ancient as meiosis itself [5]. Concerning the second, we have demonstrated that the polyphenic shifts from apomeiosis to meiosis or vice versa are metabolically regulated. Here, MMCs of apomictic TDG forming taxa reverted to meiosis at high frequencies when exposed to meiosis promoting metabolic conditions. Likewise, MMCs, nucellar cells, and/or microsporocytes from five sexual taxa-representing three angiospermous families-reverted to apomeiosis at high frequencies when exposed to apomeiosis promoting metabolic conditions.
The experiments reported herein can be easily modified to further elucidate metabolic and molecular details of meiosis apomeiosis switching using florets from a single genotype, or even from a single plant. Through such studies, it may be possible to identify wild type genes that faithfully induce apomixis in crops when their levels and locations of expression are transgenically manipulated.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4425/11/12/1449/s1, Figure S1: Mean (±SD) lengths of pistils by ovule development stage, Figure S2: Microarrays that passed quality control analyses, Figure S3: qRT-PCR verification of microarray expression profiling, Table S1: Rules used to assign gene ontology (GO) terms to groups and subgroups, Table S2: Gene ontology (GO) groups and subgroups for all GO terms of differentially expressed genes, Table S3: Calculations for chi-square tests of independence and for the determination of standard and adjusted residuals, Table S4: Primers used for qRT-PCR verification of microarray data, Table S5: Differentially expressed genes for 17 microarray comparisons, Table S6: All gene ontology categories (GOs) of all differentially expressed genes (DEGs) by comparison, Table S7: GO terms by comparison using DEGs that were overrepresented or underrepresented, Table S8: Enriched gene ontology terms (EGOs) of comparisons, Table S9: Genes by comparison identified by keyword searches of DEG descriptions.