Homologous genes of IMPDH in A. thaliana
A previously cloned IMPDH gene in A. thaliana (IMPDH1, At1g79470) (Collart et al. 1996) and its homologous gene (IMPDH2, At1g16350) (Witte and Herde 2020), share high similarity in their amino acid sequences (identity: 84%; similarity: 91%) (Fig. S1). A comparison of the amino acid sequences of AtIMPDH1 and AtIMPDH2 with those of human HsIMPDH1 and HsIMPDH2 showed that the binding sites for NAD+ and IMP and the cysteine residues essential for enzyme activity were all conserved in AtIMPDH1 and AtIMPDH2 (Fig. S1). However, the N-terminal Bateman domain, which controls enzyme activity, is less conserved between human and Arabidopsis IMPDHs, suggesting that the regulatory mechanism of these two may differ (Fig. S1). Next, we compared the expression levels of AtIMPDH1 and AtIMPDH2 (hereafter referred to as IMPDH1 and IMPDH2, respectively) based on publicly available RNA-seq data and found that IMPDH2 was more highly expressed than IMPDH1 in most tissues. The only exception was mature pollen, where IMPDH2 expression was barely detectable (Fig. S2). These data suggested that IMPDH2, rather than IMPDH1, has a crucial role in most organs.
Embryonic lethality in impdh1 impdh2 double mutants
To investigate the physiological roles of IMPDH1 and IMPDH2, we established Arabidopsis T-DNA insertion lines containing impdh1-1 (SAIL_5_B10C1) with an insertion in the third exon, as well as impdh2-1 (SALK_201269) and impdh2-2 (SALK_203056), both with an insertion in the second exon (Fig. 2a-b). The expression levels of IMPDH1 and IMPDH2 across these T-DNA insertion sites were not detected in any of the mutants (primer pairs C-E and a-c, respectively). These expression levels were generally lower than those observed in the controls before (primer pairs A-B and a-b) and after (primer pair d-e) the T-DNA insertion sites, except for those of transcripts between primers F-G in impdh1-1 (Fig. 2a, c). Next, we attempted to generate double mutants of both the genes. The impdh1-1 and impdh2 mutants showed no defects in seed development, but IMPDH1/impdh1-1 impdh2-1 and IMPDH1/impdh1-1 impdh2-2 (hereafter referred to as ‘impdh1(+/–) impdh2-1’ and ‘impdh1(+/–) impdh2-2’, respectively) developed aborted seeds, suggesting that impdh1 impdh2 double mutants are defective in gametogenesis and/or embryogenesis (Fig. 2d). We then tested the segregation ratios of impdh1(+/–) impdh2-1 and impdh1(+/–) impdh2-2, which confirmed the lethality of the impdh1 impdh2 double mutants (Table 1), with the ratio of IMPDH1/IMPDH1 to IMPDH1/impdh1-1 heterozygous plants clearly less than the expected ratio of 1:2 (Table 1). The segregation ratio of the reciprocal crossing between impdh1(+/–) impdh2-2 and impdh2-2 showed that although both the cases were negatively affected, the effect was more severe when the females were mutated in impdh1-1 impdh2-2 (Table 2). These results strongly suggested that both the genes are redundant and are important for development.
Table 1
Segregation of self progeny.
Parental genotypes | Genotypes of IMPDH1 (n) | | c2 test (p) |
| +/+ | +/- | -/- | Total (n) | (vs. 1:2) |
IMPDH1/impdh1-1 impdh2-1 | 285 | 136 | 0 | 421 | 1.37E-50 |
IMPDH1/impdh1-1 impdh2-2 | 226 | 96 | 0 | 322 | 1.01E-44 |
Table 2
Segregation of reciprocal crosses.
Parental genotypes | Genotypes of IMPDH1 (n) | | c2 test (p) |
(Female x Male) | +/+ | +/- | Total (n) | (vs. 1:1) |
IMPDH1/impdh1-1 impdh2-2 x impdh2-2 | 126 | 11 | 137 | 8.78E-23 |
impdh2-2 x IMPDH1/impdh1-1 impdh2-2 | 102 | 67 | 169 | 7.09E-03 |
Table 3. Frequencies of leaf shape phenotype. | |
| Leaf shape frequencies | Total (n) |
Genotype | Flat | (%) | Trumpet | (%) | Needle | (%) |
WT | 219 | 100.0 | 0 | 0.0 | 0 | 0.0 | 219 |
as2-1 | 240 | 100.0 | 0 | 0.0 | 0 | 0.0 | 240 |
rpl4d-3 | 245 | 100.0 | 0 | 0.0 | 0 | 0.0 | 245 |
as2-1 rpl4d-3 | 50 | 22.3 | 16 | 7.1 | 158 | 70.5 | 224 |
impdh2-1 | 256 | 100.0 | 0 | 0.0 | 0 | 0.0 | 256 |
i1(+/-) i2-1 | 252 | 100.0 | 0 | 0.0 | 0 | 0.0 | 252 |
as2-1 impdh2-1 | 227 | 100.0 | 0 | 0.0 | 0 | 0.0 | 227 |
as2-1 i1(+/-) i2-1a | 187 | 92.1 | 10 | 4.9 | 6 | 3.0 | 203 |
impdh2-2 | 221 | 100.0 | 0 | 0.0 | 0 | 0.0 | 221 |
i1(+/-) i2-2 | 210 | 100.0 | 0 | 0.0 | 0 | 0.0 | 210 |
as2-1 impdh2-2 | 223 | 100.0 | 0 | 0.0 | 0 | 0.0 | 223 |
as2-1 i1(+/-) i2-2b | 220 | 94.0 | 8 | 3.4 | 6 | 2.6 | 234 |
aas2-1 i1(+/–) i2-1: as2-1 IMPDH1/impdh1-1 impdh2-1 | | |
bas2-1 i1(+/–) i2-2: as2-1 IMPDH1/impdh1-1 impdh2-2 | | |
Developmental defects in impdh mutant plants
Because the function of IMPDH is essential, quantitative sufficiency of IMPDH function may be important for growth. Therefore, we examined the phenotypes (shoot fresh weight, chlorophyll content in cotyledons, and primary root length) of the impdh1-1 and impdh2 mutants in detail. Although impdh1-1 showed no apparent difference from the wild-type, impdh2-1 and impdh2-2 had shorter primary root lengths and lighter shoot fresh weight after 7 days of growth (Fig. 3a–c). These results are consistent with the higher expression of IMPDH2 than that of IMPDH1 in most tissues (Fig. S2). Additionally, the impdh2 mutants displayed transient light green cotyledons and lower chlorophyll content than that of the wild-type on the 4 days of growth; the chlorophyll content of the impdh2 mutants recovered to roughly the same level as that of the wild-type by the 7 days of growth (Fig. 3a, d). Next, we isolated impdh1(+/–) impdh2 mutants to test their phenotypes at reduced IMPDH levels. Upon observation of seedlings derived from impdh1(+/–) impdh2-1 and impdh1(+/–) impdh2-2, some plants had severe growth defects and a light leaf color, whereas others did not (Fig. S3a). Therefore, we checked the genotypes of 20 seedlings of each phenotype and found that all the growth-defective plants were impdh1(+/–) impdh2-1 or impdh1(+/–) impdh2-2, whereas all others were impdh2-1 or impdh2-2 (Fig. S3a–c). Subsequently, in the following experiments, we tried to distinguish impdh1(+/–) impdh2-1 and impdh1(+/–) impdh2-2 from impdh2-1 and impdh2-2 based on differences in the post-germinative growth. The growth of impdh1(+/–) impdh2-1 and impdh1(+/–) impdh2-2 showed even stronger growth defects in both primary root length and shoot fresh weight than those of impdh2-1 and impdh2-2, respectively (Fig. 3a–c). Even after 7 days of growth, although the leaf chlorophyll content of the impdh2 mutants had recovered to the level observed in the wild-type, those of impdh1(+/–) impdh2-1 and impdh1(+/–) impdh2-2 remained low (Fig. 3a, d). In addition, the frequency of impdh1(+/–) impdh2-1 and impdh1(+/–) impdh2-2 was high in seedlings with three or four cotyledons, which was not the case for the impdh1-1 or impdh2 mutants (Fig. S4a-b). These results indicated that sufficient amounts of IMPDH are important for the integrity of chlorophyll biosynthesis, growth, and cotyledon development.
Effect of mycophenolic acid (MPA) and guanosine treatments on impdh mutant plants
MPA, an inhibitor of IMPDH, causes ribosomal stress responses in animal cells, such as the suppression of rRNA synthesis, reduction of nucleolus size, and stabilization of p53 protein in specific cell types (Huang et al. 2008; Freedman et al. 2019; Kofuji et al. 2019). In plants, MPA treatment caused the translocation of several nucleolar-localized proteins to the nucleoplasm (Ahn et al. 2016; Choi et al. 2020). To estimate if plant IMPDH is also affected by MPA, we tested its effect on the impdh mutants. MPA treatment of wild-type and impdh1-1 at a concentration of 3 nM had no effect on primary root elongation, whereas MPA at a concentration of 30 nM inhibited primary root elongation by approximately half (Fig. 4A). In contrast, the growth of impdh2-1 and impdh2-2 was significantly inhibited compared to the control condition upon treatment with 3 nM MPA, which had no effect on the growth of the wild-type or impdh1-1; in addition, the growth was drastically inhibited upon treatment with 30 nM MPA (Fig. 4A). These results suggested that the relatively low abundance of IMPDH in impdh2-1 and impdh2-2 inhibited root growth at lower concentrations of MPA than that observed in the wild-type and impdh1-1.
There are two known pathways for GTP biosynthesis in animals and possibly in plants: de novo biosynthesis pathway involving IMPDH and a purine salvage pathway that does not involve IMPDH (Fig. 1, Kofuji and Sasaki 2020; Witte and Herde 2020). Accordingly, GTP can be synthesized from guanosine via the purine salvage pathway and does not require functional IMPDH. Therefore, we investigated whether the growth impairment caused by the impdh mutants used in this study due to insufficient GTP levels could be alleviated by guanosine treatment. The growth of impdh2-1 and impdh2-2 was significantly recovered upon treatment with 20 µM guanosine, whereas treatment with 20 µM guanosine had no effect on the growth of both the wild-type and impdh1-1 (Fig. 4b). Furthermore, in impdh1-1(+/–) impdh2-1 and impdh1-1(+/–) impdh2-2, treatment with 200 µM guanosine partially recovered the growth inhibition as well as the light color of the cotyledons (Fig. 4b). The results of these MPA and guanosine treatment experiments strongly suggested that the impdh1 and impdh2 mutants, especially the impdh1-1(+/–) impdh2 mutants, lack IMPDH enzyme activity and reduce GTP levels.
Involvement of ribosomal stress in impdh mutant plants
Since the inhibition of IMPDH function in animals causes ribosomal stress (Huang et al. 2008; Freedman et al. 2019; Kofuji et al. 2019), we investigated whether this stress also occurred in the plant impdh mutants. Because several Arabidopsis ribosome biogenesis-associated mutants were more resistant to translation inhibitors than the wild-type (Rosado et al. 2010; Hsu et al. 2014; Zhu et al. 2016; Maekawa et al. 2018), we also tested the growth of the impdh mutants on a medium containing translation inhibitor (chloramphenicol, spectinomycin, or streptomycin) and measured the primary root length. The primary root length of the impdh2 mutant in the control medium was lesser than that of the wild-type, which either became equivalent to or higher in the impdh2 mutants grown on the translation inhibitor-containing medium than that in the wild-type (Fig. 5). In other words, the impdh2-1 and impdh2-2 mutants were more resistant to the translation inhibitor than the wild-type and impdh1-1.
Since human IMPDH inhibition-induced ribosomal stress reduces the nucleolus size, which is an important compartment for ribosome biogenesis (Kofuji et al. 2019), we tested whether the impdh mutants had a reduced nucleolus size. Detection and comparison of the nucleolus size near the root tip showed that while impdh1-1 did not differ from the wild-type, impdh2-1 and impdh2-2 showed a reduction of ~ 40% in size, as compared with the median value (Fig. 6). Owing to the experimental setup, the impdh1(+/–) impdh2 mutants were too small to be used as test samples. These results indicated that a characteristic ribosomal stress response occurs in the impdh mutants, which depends on the level of IMPDH.
rRNA maturation in the impdh mutant plants
Because MPA treatment or suppression of IMPDH function also results in the inhibition of rRNA synthesis in animals (Huang et al. 2008; Kofuji et al. 2019), we attempted to detect rRNA profiles in impdh mutant plants grown for 7 days. First, quantitative PCR was used to quantify the amount of rRNA relative to that of cytosolic 18S rRNA. The amounts of plastid 16S and 23S rRNA in the impdh1(+/–) impdh2 mutants were clearly lower than those in the other samples, whereas the amounts of the other rRNAs were not different (Fig. 7a). The processing of cytosolic and plastid rRNA was detected using an RNA gel-blot using two probes (See Fig. 7b for rRNA processing steps and probe sites). There were no differences in cytosolic rRNA processing patterns between the samples (Fig. 7c). For plastid rRNA, the impdh2 mutants showed slightly more plastid 1.7 kb 16S rRNA levels than the wild-type, but no other apparent differences, and the impdh1-1 mutant did not differ from the wild-type. In contrast, the amount of mature 1.5 kb 16S rRNA and its ratio to its precursor 1.7 kb 16S rRNA (1.5 kb/1.7 kb) was clearly lower in the impdh1(+/–) impdh2 mutants than in the other samples. In addition, although it did not appear to be stuck in a particular processing step, the impdh1(+/–) impdh2 mutants generally had lower amounts of rRNA derived from rrn23 than the other samples (Fig. 7c). These results indicated that impdh mutations have an inhibitory effect on the synthesis and processing of plastid rRNA, rather than those of cytosolic rRNA.
Adaxial-abaxial leaf polarity alterations in as2-1 impdh1(+/–) impdh2 mutant plants
Mutations in plant ribosomal protein genes and ribosome biogenesis-related factors enhance the leaf abaxialization phenotype of the as2 mutant (Horiguchi et al. 2011; Machida et al. 2015). Genes that enhance the as2 phenotype are called ‘modifier genes’ (Machida et al. 2015). Therefore, we generated multiple mutant lines with as2-1 to examine whether impdh mutations had similar developmental effects. as2-1 rpl4d-3 showed a strong leaf abaxialization phenotype (Horiguchi et al. 2011) and was used as the positive control, which developed ~ 70% needle- and ~ 7% trumpet-like leaves (Fig. 8 and Table 3). Observations of several multiple mutant lines revealed that the as2-1 impdh1(+/–) impdh2 mutants developed trumpet- and needle-like leaves at a frequency of 3.4–4.9% and 2.6–3.0% respectively, which was not found in any other mutant line (Fig. 8 and Table 3). These results indicated that the impdh1(+/–) impdh2 mutations enhance the as2 mutation and that IMPDH serves as a 'modifier gene' of as2. Taken together, our results indicated that the impdh mutations causes both previously known and novel ribosomal stress responses, with a major impact on plastids.