- Research article
- Open Access
Phenotypic and transcriptomic characterization of a wheat tall mutant carrying an induced mutation in the C-terminal PFYRE motif of RHT-B1b
BMC Plant Biologyvolume 18, Article number: 253 (2018)
As central regulators of the gibberellic acid (GA) signaling pathway in plants, DELLA proteins function as growth repressors and affect diverse biological processes. The wheat RHT-B1b and RHT-D1b semi-dwarfing alleles, which encode GA-insensitive DELLA proteins, have been widely adopted in modern wheat varieties to improve lodging tolerance and harvest index. However, the molecular mechanisms by which DELLA modulates these responses in wheat remain largely unknown.
We identified a tall tetraploid wheat mutant line carrying an induced missense mutation (E529K) in the PFYRE motif of RHT-B1b that partially suppressed the semi-dwarf phenotype. The height-increasing effect of RHT-B1bE529K relative to RHT-B1b (19 cm or 21% increase) was significantly smaller than the effect of RHT-B1a (33 cm or 34% increase) relative to RHT-B1b in the same field experiment. The RHT-B1bE529K mutation was also associated with length increases in coleoptiles, seedling shoots, and stem internodes relative to the RHT-B1b allele. We detected no significant differences in germination rate, seedling root length, tiller number, flag leaf size, spike length, or yield components. Using RNA-seq, we compared gene expression profiles of plants encoding RHT-B1b and RHT-B1bE529K in coleoptile, first leaf, and elongating peduncles. We detected limited overlap among tissues of the genes differentially regulated by the two genotypes, and more genes upregulated (77%) than downregulated (23%) in RHT-B1bE529K relative to RHT-B1b. These results suggest that the wheat DELLA protein affects the transcriptome in a tissue-specific manner and that the mutation mainly eliminates or reduces repression functions of the RHT-B1 protein. Our study identified distinct sets of potential DELLA direct or indirect target genes involved in cell wall and carbohydrate metabolisms, cell cycle/division, and hormone pathways.
We identified the hypomorphic RHT-B1bE529K allele that confers an intermediate plant height and coleoptile elongation. This allele can be useful in rain-fed wheat breeding programs where the strong reduction in height and biomass associated with RHT-B1b has detrimental effects. Transcriptomic characterization of different tissues from the plants encoding RHT-B1bE529K and RHT-B1b provided valuable information for identifying DELLA downstream GA response genes in wheat.
Gibberellins (GAs) are a family of plant hormones that regulate a diverse range of developmental processes including germination, leaf expansion, stem elongation, phase transition, and flowering . GA biosynthesis and deactivation are highly regulated by a number of metabolic enzymes encoded by multi-gene families, and the downstream GA signaling pathway is mainly regulated by DELLA proteins [2, 3]. Lacking a canonical DNA binding site, DELLA represses growth in the absence of GA through physical interaction with its target proteins (reviewed in ). Mechanisms include the sequestration of transcription factors and transcriptional regulators and the transactivation of growth repressor genes [5,6,7,8,9,10]. When bioactive GAs are perceived by the receptor protein GID1 (GIBBERELLIN INSENSITIVE DWARF1), conformational change occurs in the N-terminal region of GID1 and enhances the formation of the GA–GID1–DELLA complex [11, 12]. This enables the SCF (SKP1, CULLIN, F-box) E3 ubiquitin ligase complex to recognize and polyubiquitinate DELLA, which leads to its degradation via the 26S proteasome pathway . DELLA targets are then released, promoting downstream GA-mediated growth responses.
DELLA proteins belong to the GRAS (GAI, RGA, and SCR) protein family and are composed of two major domains – the N-terminal regulatory domain (DELLA, LExLE, and TVHYNP motifs) and the C-terminal functional domain (LHRI, VHIID, LHRII, PFYRE, and SAW motifs). The N-terminal, or DELLA, domain, is involved in GID1 binding and confers functional specificity that distinguishes DELLA from other members of the GRAS subfamily [14, 15]. The C-terminal, or GRAS, domain is highly conserved among the GRAS subfamilies. The first leucine heptad repeat (LHRI) motif of the GRAS domain is essential for DELLA’s interaction with its target proteins (e.g. PIF4 , JAZ1 , and AtIDD3 ). The LHRI–VHIID–LHRII structure is highly conserved among GRAS subfamily proteins and has been shown to play similar roles . The adjacent PFYRE motif is less conserved and its specific role in protein-protein interaction is currently not clear [4, 15]. Finally, the distal SAW motif is crucial for the DELLA’s interaction with BZR1  and GAF1 .
In hexaploid wheat (Triticum aestivum; 2n = 42, AABBDD), the DELLA proteins are encoded by three RHT1 (REDUCED HEIGHT 1) homeologs – RHT-A1, RHT-B1, and RHT-D1 on chromosome 4AL, 4BS, and 4DS, respectively. During the ‘Green Revolution’ in the 1960s, the semi-dwarfing Rht-B1b and Rht-D1b alleles were introduced into breeding programs to develop high-yielding wheat varieties . The short and strong stem of these modern varieties greatly improved productivity by preventing lodging under high fertilizer input and enhancing the efficiency of photosynthate partitioning into grains [18,19,20]. Both Rht-B1b and Rht-D1b alleles have a point mutation inducing a premature stop codon in the N-terminal DELLA domain  that prevents GID1 binding. It has been suggested that translation reinitiates after these stop codons, to produce N-terminally truncated DELLA proteins with the active GRAS domain retaining its ability to repress growth. These truncated DELLAs cannot detect the GID1-mediated GA signals, thus conferring GA-insensitive semi-dwarfism .
Rht-B1b and Rht-D1b are the major semi-dwarfing alleles and are present in over 70% of wheat cultivars worldwide [22, 23]. Although their effects on yield are positive in optimal environments, negative pleiotropic effects such as low seedling vigor [24, 25] and susceptibility to Fusarium head blight [26, 27] have been reported. In order to improve allelic diversity, find alternative RHT1 alleles, and understand better the function of the different DELLA subdomains, several natural and induced variants in the RHT1 genes have been studied previously. Additional GA-insensitive dwarfing alleles found in natural germplasm include Rht-B1c (in-frame insertion of 30 amino acids between the DELLA and LExLE motifs [28, 29]), Rht-B1e and Rht-B1p (nonsense mutations located three and four amino acids prior to the Rht-B1b mutation, respectively [30, 31]), and Rht-D1c (increased copy number of Rht-D1b ). Also, natural polymorphic sites including insertions, deletions and point mutations have been identified in the coding and untranslated regions of all three RHT1 homeologs by sequencing  or ecotype targeting induced local lesions in genomes (EcoTILLING) . In addition to natural variation, suppressor screens using chemically mutagenized dwarf wheat identified derivative Rht-B1c alleles in the LHRI, VHIID, PFYRE, and SAW motifs that confer overgrowth phenotypes [34, 35].
In this study, we identified a mutant line exhibiting an increased height phenotype from a tetraploid wheat (Triticum turgidum; 2n = 28, AABB) TILLING population mutagenized with ethyl methanesulfonate (EMS). This line carries a missense mutation in the PFYRE motif of the RHT-B1b protein suppressing the dwarfing effect of Rht-B1b. We describe its effect on major agronomic traits and analyze transcriptomic changes induced in different tissues using RNA-seq.
Mutations in the C-terminal GRAS domain of RHT-B1b partially suppress the semi-dwarfing effect of RHT-B1b
A tall mutant line (T4-934) was identified in this study from an EMS-mutagenized population of the durum wheat variety Kronos developed before . Because the parental variety of this population (Kronos) carries the semi-dwarfing Rht-B1b allele, we hypothesized that the increased height of T4-934 may result from a second-site mutation in Rht-B1b suppressing the effects of this allele. We sequenced the C-terminal region of RHT-B1 and found that T4-934 carries a G-to-A mutation at position 1,585 from the start codon (G1585A), changing a negatively charged glutamate (E, the last E in the PFYRE motif designation) to a positively charged lysine (K) (hereafter RHT-B1bE529K, Fig. 1). The SIFT (Sorting Intolerant From Tolerant ) score of the RHT-B1bE529K mutation is zero, indicating that it changes a conserved amino acid residue and has a high probability of producing deleterious effects (Additional file 1: Table S1). We used the previously published exome sequencing of the complete Kronos TILLING population  to confirm that, except for E529K, there were no additional mutations in the RHT-B1 coding sequence in T4-934.
In the F2 and BC1F2 populations from the cross between T4-934 and the non-mutagenized Kronos (recurrent parent), plants homozygous for RHT-B1bE529K were significantly taller than plants homozygous for RHT-B1b (P < 0.0001) when grown under field conditions, with an estimated additive effect of 10.0 cm (F2) and 9.4 cm (BC1F2, Table 1). Estimated degrees of dominance  for the F2 (D = 0.18) and BC1F2 (D = 0.30) indicated partial dominance of the height-increasing allele (Table 1).
In the same field experiment, the BC1F2 population Kronos*2/T4-934 was grown together with the BC1F2 population Kronos*2/Gredho, which segregated for the wild-type RHT-B1a and the semi-dwarfing RHT-B1b alleles. In the Kronos*2/T4-934 population, homozygous RHT-B1bE529K plants were 19 cm taller (21%) than homozygous RHT-B1b plants (Fig. 2a), whereas in the Kronos*2/Gredho population homozygous RHT-B1a plants were 33 cm taller (34%) than homozygous RHT-B1b plants (Fig. 2b). A two-way factorial ANOVA including populations and RHT-B1 alleles as factors showed a highly significant interaction (P =0.0024), confirming that the height difference between RHT-B1a and RHT-B1b was significantly larger than that between RHT-B1bE529K and RHT-B1b. This result suggests that the E529K mutation is responsible for a partial suppression of the dwarfing effect of RHT-B1b.
The screening of the Kronos TILLING database  yielded six additional lines carrying nonsynonymous mutations in the distal part of the C-terminal GRAS domain of RHT-B1b including the SAW and part of PFYRE motifs (Fig. 1). Interestingly, one of these lines, T4-1405, carried the same E529K mutation as T4-934. Comparison of exome sequencing data in several independent genes using the TILLING database (https://dubcovskylab.ucdavis.edu/wheat_blast) showed that E529K was the only common mutation between the two lines, confirming their independent origin. We grew M4 populations (n = 16–18) segregating for each non-synonymous mutation under greenhouse conditions to determine their effect on plant height. Induced mutations identified in T4-1405 (E529K) and T4-3545 (W605*) showed highly significant increases (P = 0.0007 and <0.0001, respectively) in plant height compared to wild-type Kronos (Fig. 2c, d; Additional file 1: Table S1). Plants homozygous for these two mutant alleles were significantly taller (95.6 cm and 97.5 cm for T4-1405 and T4-3545, respectively) than wild-type Kronos (70.0 cm), and of similar height to the homozygous RHT-B1bE529K plants (94.0 cm). No significant effects on plant height were detected for the other four mutations (Fig. 1; Additional file 1: Table S1).
Taken together, these results demonstrate that the RHT-B1bE529K mutation in the PFYRE motif is sufficient to suppress significantly (but not completely) the semi-dwarf phenotype of lines carrying the RHT-B1b allele and that the mutant RHT-B1bE529K allele is partially dominant for increased plant height.
Effect of RHT-B1bE529K on other agronomic traits
In field experiments, all above-ground internodes were significantly longer (P < 0.0001) in homozygous RHT-B1bE529K plants than in homozygous RHT-B1b plants in both the F2 and BC1F2 populations (Table 1; Fig. 3), demonstrating that RHT-B1bE529K enhances stem growth throughout the elongation period. For all internodes, we found positive values for the degree of dominance, indicating partial dominance of the allele for increased length (Table 1; Fig. 3). RHT-B1bE529K had no significant effect on internode number, tiller number per plant, flag leaf length and width, days to heading, spike length, spikelet and grain number per spike, and grain weight (Table 1).
We found no significant difference in germination rate between homozygous RHT-B1bE529K and RHT-B1b plants in BC1F3 sister lines (Fig. 4a). However, in two-week-old seedlings, both coleoptiles and shoots were significantly longer in homozygous RHT-B1bE529K plants than RHT-B1b plants (Fig. 4b and c). Compared to homozygous RHT-B1b sister lines, coleoptiles and shoots of two-week-old homozygous RHT-B1bE529K plants were 17% (0.5 cm) and 22% (3.6 cm) longer (P < 0.0001), respectively. No significant difference was observed in seedling root length (Fig. 4b).
Taken together, these data show that RHT-B1bE529K enhances coleoptile, seedling shoot, and stem elongation without affecting the other agronomic traits evaluated in this study.
Effect of RHT-B1bE529K on GA sensitivity
In order to determine the effect of RHT-B1bE529K on GA sensitivity, we measured the response of homozygous RHT-B1bE529K BC1F3 plants to different GA3 concentrations (Fig. 5). As controls, we included GA-insensitive wild-type Kronos carrying RHT-B1b, and two GA-sensitive lines: an EMS-induced mutant T4-3822 previously revealed to have a 1.9 Mb deletion encompassing the whole RHT-B1 gene , and the tall tetraploid variety Gredho carrying the RHT-B1a allele.
As expected, the GA-insensitive Kronos showed no significant changes in coleoptile length under any tested GA3 concentration, while the GA-sensitive controls T4-3822 (RHT-B1-null) and Gredho (RHT-B1a) showed significant increases in coleoptile length under GA3 concentrations of ≥0.1 μM and ≥1.0 μM, respectively (Fig. 5). Homozygous RHT-B1bE529K lines showed significant increase in coleoptile length only under the highest GA3 concentration (10 μM). The level of GA sensitivity was determined using the GA response index (GRI) as described in “Materials and Methods” (Fig. 5). GRI of homozygous RHT-B1bE529K lines was slightly higher than that of wild-type Kronos, but was lower than T4-3822 and Gredho. Our results indicate that the RHT-B1bE529K mutation in the C-terminal functional domain increases seedling sensitivity to exogenous GA in spite of the presence of the Q64* stop codon mutation in the N-terminal regulatory domain of RHT-B1b (Fig. 1).
Transcriptome profile comparing different tissues from homozygous RHT-B1b and RHT-B1bE529K plants
To determine the effects of RHT-B1bE529K on the transcriptome, we performed an RNA-seq experiment using plant tissues that exhibited significant phenotypic differences between homozygous RHT-B1b and RHT-B1bE529K plants. In addition, these tissues are actively undergoing cell expansion, and are thus more likely to be active sites of GA signaling. We collected tissues from coleoptile, first leaf, and peduncles at early- and mid-stages (Z49 and Z52) of elongation from BC1F3 sister lines homozygous for the RHT-B1b or RHT-B1bE529K alleles. We constructed and sequenced 32 RNA-seq libraries (two genotypes × four tissues × four biological replicates), generating on average 35.9 million 100 bp single-end reads per library (Additional file 1: Table S2). For subsequent read count and differentially expressed gene (DEG) analyses, we used only reads that were uniquely mapped to a single annotated gene (IWGSC RefSeq v1.0; Additional file 1: Table S2, 70.3% of raw reads, average 25.3 million reads per library). Analysis of the reads mapping to the genomic location of the RHT-B1bE529K mutation (at 30,862,966 bp on chromosome 4B; IWGSC RefSeq v1.0) confirmed the expected genotype of each RNA-seq library.
In both RHT-B1b and RHT-B1bE529K genotypes, coleoptile tissue expressed the highest number of genes (70,912 and 69,688, respectively) followed by Z52 peduncle (68,127 and 68,208), Z49 peduncle (64,553 and 64,373), and first leaf (63,260 and 64,003). A heat-map clustering analysis using the expression data of all 32 samples showed clear separation among the four different tissues, which indicated consistent differences in the transcriptome profiles among the samples from different tissues (Fig. 6a). As expected, Z49 and Z52 peduncle samples clustered more closely together than other samples. Within each tissue, however, there was no clear separation by genotype. A similar pattern was observed in a principal component analysis (Fig. 6b), where the four tissues were separated into distinct clusters while plants encoding RHT-B1b and RHT-B1bE529K clustered closely together within each tissue. These results indicated that large transcriptomic differences exist between different tissues, while a relatively smaller number of genes are differentially expressed between plants carrying the different RHT-B1 alleles.
Differentially expressed genes between plants encoding RHT-B1b and RHT-B1bE529K
We compared homozygous RHT-B1b and RHT-B1bE529K plants to identify DEGs between genotypes in each of the four tissues. We conducted two different statistical analyses using the DESeq2 and EdgeR packages and considered a gene as differentially expressed only when they were significant (FDR-adj P < 0.05) in both tests (Fig. 7a; Additional file 1: Table S3 and S4; Additional file 2: Figure S1). A small number of DEGs (n = 229) were identified between the homozygous RHT-B1b and RHT-B1bE529K plants for the four tissues. As a comparison, nearly 120 times more DEGs (n = 27,362) were detected between Z49 and Z52 peduncles (6 days difference) using the same statistical criteria. No significant difference was detected in the expression of RHT-B1 between homozygous RHT-B1b and RHT-B1bE529K plants (Additional file 2: Figure S2), indicating that the RHT-B1bE529K mutation does not affect RHT-B1 transcript levels.
Among the 229 DEGs identified between homozygous RHT-B1b and RHT-B1bE529K plants in four different tissues (Fig. 7a), the majority (76.9%) were upregulated in the plants homozygous for the RHT-B1bE529K allele. Only one gene (TraesCS3A01G407200 encoding a putative NAD(P)H-quinone oxidoreductase subunit N) was differentially expressed in two different tissues, suggesting that distinct sets of genes are regulated by DELLA in a tissue-specific manner. The full list of 229 DEGs is provided in Additional file 1: Table S3. Table 2 presents 81 selected DEGs for which we found more detailed functional information in the literature, which are described below by tissue.
The seven upregulated DEGs included two genes encoding putative fatty acyl-coenzyme A reductases (FARs) that may play a role in synthesizing cell wall components [41, 42] (Table 2). A gene encoding plant defense related Hessian fly responsive-2 (Hfr-2)-like protein  was also upregulated. The two downregulated DEGs were a gene encoding a putative lateral organ boundaries (LOB) domain protein, a plant-specific transcription factor that regulates diverse developmental events , and a gene encoding a putative methionine aminopeptidase required for normal plant development .
Among 23 DEGs identified in the first leaf, 20 were upregulated and three were downregulated in the homozygous RHT-B1bE529K plants (Fig. 7a). Similar to the coleoptile, upregulated DEGs in the mutant plants included putative defense related genes encoding NBS-LRR proteins , a disease resistant protein RMP1 , and heavy metal transport/detoxification proteins  (Table 2). Many genes encoding kinase proteins were also upregulated, including receptor kinase genes potentially controlling plant development and defense mechanisms . Downregulated DEGs in the mutant plants included a gene encoding a UDP-glycosyltransferase potentially involved in glycosylation of plant hormones [50, 51], and a gene encoding an IAA-amino acid hydrolase likely involved in auxin regulation .
Z49 peduncle (early-stage elongation)
For the 119 upregulated and 44 downregulated DEGs identified in Z49 peduncle of the homozygous RHT-B1bE529K plants (Fig. 7a), we first conducted a GO enrichment analysis within the “Biological Process” terms to classify the DEGs (Fig. 7b). Enriched GO terms among upregulated DEGs were mainly associated with transmembrane transport and catabolic processes (Fig. 7b). DEGs associated with transmembrane transport included genes encoding a phosphate transporter, a sulfate transporter, two amino acid permeases, and a putative Na+/H+ antiporter (Table 2). DEGs associated with catabolic processes included genes encoding three putative RING/U-box proteins, a carboxypeptidase, two chitinase-like proteins, and a carotenoid cleavage dioxygenase. Five additional genes encoding putative arabinogalactan proteins (AGPs), which play important roles in cell expansion and stem elongation [53,54,55], were also upregulated. Two genes involved in the auxin pathway (i.e., an auxin influx transporter and an auxin-responsive protein) and two genes encoding LOB domain-containing proteins were upregulated. Downregulated DEGs in the homozygous RHT-B1bE529K plants were largely enriched in GO terms associated with cell division/cycle such as DNA replication, DNA unwinding, DNA repair, and chromosome organization (Fig. 7b). DEGs related to these GO terms included genes encoding putative helicases, DNA polymerases, DNA mismatch repair proteins, histone modifiers, and cyclin proteins (Table 2). Among the downregulated genes was one encoding a GRF1-interacting factor (GIF)-like protein, which is known to promote cell proliferation .
Z52 peduncle (mid-stage elongation)
Among 35 DEGs identified in Z52 peduncle, 30 were upregulated and five were downregulated in the homozygous RHT-B1bE529K plants (Fig. 7a). Many upregulated DEGs were involved in carbohydrate metabolism, i.e., genes encoding three carbonic anhydrases (CAs) , two fructose-bisphosphate aldolases (FBAs) , a sucrose-phosphate synthase (SPS) , a glyceraldehyde-3-phosphage dehydrogenase (GAPDH) , a pyrophosphate:fructose 6-phosphate 1-phosphotransferase (PFP) , and a dihydrolipoyl dehydrogenase  (Table 2). Other upregulated DEGs in the homozygous mutant plants included two genes encoding SPIRAL1 proteins that potentially control anisotropic cell expansion . Several genes associated with defense mechanisms were also upregulated, including genes encoding an aspartyl protease [64, 65], a chitinase , and an NBS-LRR protein . Downregulated DEGs included a gene encoding phospholipase 1, which catalyzes the initial step in jasmonic acid (JA) biosynthesis , and genes encoding defense-related protein such as a putative 2-oxoglutarate Fe (II) oxygenase  and a wound-induced protease inhibitor.
RHT-B1bE529K as an intermediate dwarfing resource
In addition to its strong effect on plant height, the introgression of the RHT-B1b allele into wheat germplasm has been associated with reductions in coleoptile length and seedling vigor, increased grain number per spike, and decreased grain weight relative to RHT-B1a [18, 19, 24, 69]. Our study shows that the RHT-B1bE529K mutation partially suppressed the dwarfing effect of RHT-B1b. The comparison of two BC1F2 backcross populations developed using the same Kronos recurrent parent and grown in the same field experiment, showed that homozygous RHT-B1bE529K plants were significantly taller than the homozygous RHT-B1b plants (21%), but that the effect was smaller than the height increase associated with RHT-B1a relative to RHT-B1b plants (34%) in the other population (Fig. 2a, b).
In addition, the effect of RHT-B1bE529K was also smaller than the 49% height increase observed in Kronos mutant line T4-3822 (81.4 cm to 120.9 cm), which carries a deletion encompassing the complete Rht-B1b gene (RHT-B1-null) . Although these two Kronos mutant lines were not grown in the same field experiments, the length of their coleoptiles at different GA concentrations were compared here in the same assay (Fig. 5). Coleoptiles were significantly shorter in homozygous RHT-B1bE529K than in homozygous RHT-B1-null plants, both in the absence and presence of GA3. The mutation in RHT-B1b prevents its degradation by GA, resulting in a constitutively active repressor. The RHT-B1bE529K mutation likely reduces but does not completely suppress, the repression ability of RHT-B1b resulting in a residual dwarfing ability.
In addition to its effect on final plant height, RHT-B1bE529K also increased the lengths of coleoptiles (17%) and seedling shoots (22%) relative to sister lines encoding the RHT-B1b allele (Table 1; Fig. 4). Therefore, RHT-B1bE529K might be useful as an alternative dwarfing source in rainfed environments, where taller semi-dwarf plants frequently perform better than shorter plants [70,71,72]. We have initiated the transfer of RHT-B1bE529K to hexaploid common wheat to test this hypothesis.
Although variation in DELLA can affect seed germination , the GA-insensitive gai-1 mutant in Arabidopsis showed no differences in germination . This result is consistent with previous reports on the effect of RHT-B1b suppressing stem elongation without affecting germination rate [73, 74], and with the results presented here. In addition, the RHT-B1bE529K mutation did not significantly affect seedling root length, tiller number, heading time, flag leaf size, spike length, spikelet number per spike, grain number per spike, and grain weight. These results are also consistent with previous reports on the effect of RHT-B1b suppressing stem elongation without affecting seedling root growth , flag leaf size , spike length , and spikelet number per spike . Contrasting effects of the dwarfing alleles on grain weight have been observed in different genotypes , so the lack of differences in grain weight detected in this study are not surprising. By contrast, consistent increases in grain number per spike have been associated with RHT1 dwarfing alleles in multiple genotypes [18, 19]. These increases were likely driven by increased fertility because no differences in spikelet number per spike were detected in these studies. In our F2 and BC1F2 populations, the semi-dwarf homozygous RHT-B1b plants also showed more grains per spike than the tall plants homozygous for RHT-B1bE529K (Table 1), but the differences were not significant (Table 1). We do not know if this lack of significant differences was a result of the partial suppression effect of RHT-B1bE529K plants on plant height, or because the PFYRE motif plays a smaller role in the DELLA effect on grain number and fertility.
RHT1-mediated GA signaling in wheat
We found three lines carrying mutations in RHT-B1b that increased height; two independent lines carrying the RHT-B1bE529K point mutation and a mutant line carrying a premature stop codon mutation (W605*) that truncated the last 17 amino acids of RHT-B1b (Figs. 1 and 2). A previous study using RHT-B1c (a stronger dwarfing allele than RHT-B1b) identified five independent mutations in a similar region of the C-terminal GRAS domain that suppressed the dwarfing effect of this allele, including one line carrying the W605* mutation identified in this study. However, these mutants were 2% to 24% shorter than the homozygous RHT-B1a plants, suggesting a partial suppression effect . Taken together, these results illustrate the importance of different motifs within the C-terminal GRAS domain in modulating GA-mediated growth responses. The characterization of the mutants described in this study can help define which amino acids within these motifs are involved in different repression functions of RHT-B1.
The GA insensitivity of GA insensitive DELLA domain mutants in Arabidopsis, wheat, barley and maize are all explained by the inability of the mutant DELLA proteins to interact with the GA receptor GID1 . Furthermore, mutations in the DELLA domain in RHT-B1b abolish the interaction with wheat GID1 [29, 79] and GID2 . Since the RHT-B1bE529K mutation reported in this study is in the C-terminal region of the protein and does not restore the truncated DELLA domain, it is unlikely that it will restore the interaction between RHT-B1b and GID1 and its sensitivity to GA. A simpler hypothesis is that the E529K mutation weakened some of the repression functions of RHT-B1b allowing expression of some of its target genes. This hypothesis is consistent with the significantly higher number of upregulated (77%) than downregulated (23%) genes in the RHT-B1bE529K mutant relative to the wild-type.
DELLA proteins mediate diverse GA responses mainly through protein-protein interactions, and different subdomains of DELLA interact with distinct subsets of target proteins [4, 81, 82]. It would be interesting to test if the EMS mutations identified here in the different RHT-B1b domains (Additional file 1: Table S5) and in other studies in RHT-B1c [34, 35, 74] affect different subsets of pleiotropic effect. Since no DELLA interacting proteins have been reported so far to specifically interact with the PFYRE domain , the mechanism by which the E529K mutation in this motif affects only a subset of the RHT-B1 pleiotropic effects remains unknown.
The RHT-B1b suppressor alleles identified in our study also provide an indirect confirmation of the hypothesis that the semi-dwarf phenotype of the Rht-B1b allele arises from translational re-initiation . The reduced GA signaling of Rht-B1b is thought to be the result of the re-initiation of translation at one of the subsequent methionine amino acids (i.e., M67, M69, or M71) present after the premature Q64* stop codon. The resulting protein lacking the N-terminal DELLA domain required for the formation of GA–GID1–DELLA complex is hypothesized to be responsible for the GA-insensitive phenotype conferred by RHT-B1b [11, 12, 29]. The suppression of the dwarfing effect of RHT-B1b by our induced mutations in the C-terminal domain provides indirect support for this hypothesis, since this result would not have been possible if there were no translation of RHT-B1b.
Transcriptomic changes induced by RHT-B1bE529K
RNA-seq experiments comparing sister lines encoding RHT-B1b and RHT-B1bE529K showed that only one out of 229 DEGs between the two genotypes (FDR-adj P < 0.05 in both DEseq2 and EdgeR) overlapped in two or more tissues (Additional file 1: Table S3). A similar result was observed when using a less stringent criterion, considering DEGs that were significant in at least one of the two statistical tests (FDR-adj P < 0.05 in DEseq2 or EdgeR; Additional file 1: Table S4). Only 4.4 % of the 4,963 DEGs using this criterion overlapped in two or more tissues (Additional file 2: Figure S3). Our results parallel previous reports in Arabidopsis showing that DELLAs regulate distinct sets of genes in different organs and developmental stages [81, 82].
DELLA proteins do not possess canonical DNA binding domains and control the expression of downstream GA response genes mainly through physical interactions with transcription factors and transcriptional regulators [4, 83]. Therefore, it is logical to speculate that the transcriptomic changes induced by RHT-B1bE529K are probably associated to disrupted protein-protein interaction. Weaker interaction between DELLA and these proteins would result in increased expression of the downstream GA response genes, thus stimulating growth. It would be interesting to investigate the effect of the E529K mutation on the strength of yeast-two hybrid interactions with known DELLA partners (reviewed in ).
GA promotes growth by enhancing the expression of genes involved in cell expansion and cell proliferation, and DELLA negatively regulates these processes in organ- and developmental stage-dependent manners [82, 84,85,86]. In agreement with previously published results, we found that two genes encoding FARs potentially involved in cell wall lipid synthesis [41, 42] were upregulated in the coleoptiles of homozygous RHT-B1bE529K plants. Similarly, five genes encoding cell wall protein AGPs involved in cell expansion and stem elongation [53,54,55] were upregulated at Z49 in peduncles of plants encoding the RHT-B1bE529K allele. By contrast, a number of genes promoting cell proliferation were downregulated in the same plants. These results suggest that the increased peduncle length observed in plants encoding RHT-B1bE529K is likely associated with increased cell expansion rather than cell proliferation.
Carbohydrate metabolism modulated by GA signaling is also an important determinant in plant growth [87, 88]. This was reflected in the upregulation of many genes involved in photosynthesis and respiration in the peduncles of homozygous RHT-B1bE529K plants at Z52. These include genes encoding CA required for photosynthetic CO2 fixation , FBA involved in Calvin cycle and glycolysis [58, 89], SPS involved in catalyzing sucrose biosynthesis , GAPDH and PFP involved in glycolysis [60, 61], and dihydrolipoyl dehydrogenase involved in the Krebs cycle  (Table 2). Interestingly, the modulation of the expression of some of these genes is known to induce changes in plant growth. For example, transgenic tobacco plants overexpressing FBA exhibited increased stem length and biomass , while FBA knockdown in tomato decreased stem diameter, biomass, and seed size . Also, overexpression of SPS in rice  and tobacco  was associated with a significant increase in plant height.
Finally, DELLAs mediate signals from multiple hormone pathways by interacting with multiple hormone signaling molecules . For example, DELLAs modulate auxin signaling by interacting with auxin response factors that regulate expression of auxin response genes . DELLAs also affect strigolactone (SL) signaling by interacting with the SL receptor DWARF14 . In this study, we observed increased transcript levels of auxin and SL pathway-related genes (e.g. ILR1, auxin influx transporter, and carotenoid cleavage dioxygenase; Table 2) in the plants encoding RHT-B1bE529K. These results reflect DELLA’s role in integrating signals from different hormone pathways. GA and JA signaling jointly regulate growth and defense balance responding to developmental and environmental cues, and this is largely mediated by the DELLA’s interaction with JAZ proteins, which are negative regulators of JA signaling . Many defense-related genes differentially expressed in plants encoding RHT-B1bE529K in our study (Table 2) may reflect the roles of DELLA in the integration of JA and GA signaling.
The absence of known GA targets among the wheat DEGs between RHT-B1b and RHT-B1E529K is not surprising considering that the E529K mutation showed only a mild response to GA3. The coleoptile experiments at different GA3 concentration showed significant coleoptile length increase in RHT-B1E529K only under the highest GA3 concentration, and the effects were smaller than in the plants carrying the wild-type RHT-B1a or the RHT-B1-null alleles (Fig. 5). A meta-transcriptome analysis using 12 Arabidopsis transcriptome data sets found limited DEG overlap among different data sets, likely reflecting the tissue specificity of GA responses . In spite of the limited overlap, this meta-transcriptome analysis observed some shared DEGs involved in GA metabolism/signaling, cell cycle, cell wall, and auxin metabolism / transport / signaling. Several of the wheat DEGs detected in this study between RHT-B1 and RHT-B1bE529K are involved in similar processes (cell cycle, cell wall, auxin pathway) suggesting some conservation in the general processes affected by DELLAs in different species. The study of the functions of the DEGs (both shared and wheat-specific) identified in this study would be facilitated by the availability of new reverse genetics tools such as the wheat sequenced TILLING populations  and the implementation of CRISPR/Cas9 in wheat [96, 97].
In summary, this study identified an induced wheat mutant encoding the RHT-B1bE529K hypomorphic allele, which partially suppressed the dwarfing effect of RHT-B1b. The intermediate plant height and coleoptile length resulting from this allele may be of interest to wheat breeding programs for water-limited regions. In addition, the significant effects of the E529K mutation provides indirect evidence that the RHT-B1b translation is re-initiated after the premature stop codon in the N-terminal DELLA domain. Our RNA-seq experiment showed that most of the changes in gene transcript levels between plants encoding RHT-B1bE529K or RHT-B1b were tissue specific, and identified distinct sets of potential DELLA down-steam target genes involved in cell wall and carbohydrate metabolisms, cell cycle/division, and hormone pathways.
Plant material and phenotype evaluation
We identified a tall mutant line ‘T4-934’ during the generation advancement of an EMS-induced tetraploid wheat (var. Kronos) TILLING population  grown in the field at the University of California, Davis, CA (38° 32’ N, 121° 46’ W) in 2013. Candidate sequencing of the C-terminal region of RHT-B1 (see “RHT-B1 sequencing and marker development” section below) showed that line T4-934 carries an induced mutation that results in a change from glutamate (E) to lysine (K) at position 529 (RHT-B1bE529K) in the PFYRE motif.
Using marker-assisted selection, T4-934 was backcrossed with the non-mutagenized variety ‘Kronos’ carrying Rht-B1b, to generate F2 (n = 84) and BC1F2 (n = 119) populations segregating for the alleles encoding RHT-B1b and RHT-B1bE529K. As a control, we developed another BC1F2 (n = 158) population from the cross Kronos*2/Gredho (PI 532239, RHT-B1a) segregating for the GA-sensitive RHT-B1a and insensitive RHT-B1b alleles. The three populations were sown in one-meter rows (five plants per row) in the field at the University of California, Davis, in November 2014 and were harvested in June 2015. Days to heading was measured as the number of days from sowing to the full emergence of the first spike. Fifteen days after all plants finished heading, the number of above-ground internodes, internode length, flag leaf length and width were measured as the average values from the three tallest tillers in each plant. Tiller number per plant was also measured. Plant height was measured as length from the ground to the top of the spike of the tallest tiller, excluding awns. Upon maturity, spikes from the three tallest tillers were harvested to measure spike length, spikelet number per spike, and grain number per spike. The averages from three spikes were used to represent an individual plant. Grain weight from the three spikes was used to calculate 1,000 grain weight.
The effect of RHT-B1bE529K on germination rate, coleoptile length, seedling shoot and root lengths were evaluated using BC1F3 sister lines homozygous for the RHT-B1bE529K or RHT-B1b alleles. For germination rate, 50 seeds of each BC1F3 sister line (five homozygous RHT-B1b and five homozygous RHT-B1bE529K lines) were sown on petri dishes with distilled water, kept at 4°C for 24 h. Germinating seeds were moved to a growth chamber with 16 h of light (intensity 230 μM m−2 s−1) at 22°C and 8 h night condition at 18°C. The number of seeds with at least 2 mm radicle emergence were counted every 8 h. The same BC1F3 sister lines were sown 2.5 cm below the top of germination paper (26 cm × 13 cm) moistened with distilled water, rolled and kept upright with the bottom 4 cm of the paper soaked in distilled water at 4°C for 24 h, and moved to the growth chamber with the same conditions as described above. Coleoptile, shoot and root lengths were measured 2 weeks after sowing. Coleoptile length was determined as the distance between the embryo and the tip of coleoptile. Shoot and root lengths were determined as the distances between the embryo and the tips of the first leaf and the longest root, respectively.
The GA sensitivity assay was conducted using non-mutagenized Kronos (RHT-B1b), a tetraploid variety ‘Gredho’ (RHT-B1a), an M4 Kronos mutant line ‘T4-3822’ homozygous for a deletion encompassing the entire RHT-B1 gene , and BC1F3 sister lines of ‘Kronos*2/T4-934’ homozygous for the RHT-B1bE529K mutation. The seeds were sown on germination paper as described above, kept at 4°C for 48 h, and moved to GA3 (Sigma-Aldrich, St. Louis, MO, USA) solutions with different concentrations (0, 0.1, 1.0, and 10 μM) under room temperature. After 10 days, coleoptile and shoot length was measured. The GA response index (GRI) was determined as the percent length increase under the highest GA concentration (10 μM) relative to the control (0 μM). The experiment was conducted as a randomized complete block design with four blocks, one replication (eight subsamples) per block and treatment combination.
By searching the sequenced tetraploid wheat TILLING database (; https://dubcovskylab.ucdavis.edu/wheat_blast), six additional nonsynonymous mutations were identified in the distal region of the C-terminal RHT-B1b including the SAW and part of PFYRE motifs. M4 lines (n = 16 – 18) segregating for each of the six mutations were grown in cone-shaped pots with 6.9 cm diameter and 25.4 cm depth in the greenhouse. Plant height was measured upon maturity as described above.
One-way ANOVAs were conducted on each agronomic trait evaluated in the field for the F2 and BC1F2 populations, using the RHT-B1 genotype as a single factor. Mean comparisons were conducted using Tukey’s multiple comparison tests. Similarly, one-way ANOVAs were conducted on plant height for the M4 populations segregating for each additional RHT-B1 mutation. T-tests were conducted to analyze germination and seedling traits from the experiment conducted in the growth chamber. For the GA sensitivity assay, an ANOVA was conducted using GA concentration as a fixed factor and blocks as a random factor, followed by Dunnett’s multiple comparison tests against the control without GA. For all ANOVAs, data violating ANOVA assumptions (normality of residuals by Shapiro-Wilk test, and homogeneity of variances by Levene’s test) were transformed using power transformations. After the statistical analyses, means for graphs and tables were de-transformed to the original scale. All statistical analyses were carried out using SAS 9.4 Software (SAS Institute, Cary, NC).
RHT-B1 sequencing and marker development
A pair of genome-specific primers (forward: 5’-GACTCCTCCTGCAGCACCTA-3’; reverse: 5’-AACCCGGCGTTGCCGAGG-3’) were designed to sequence the C-terminal region of RHT-B1 (433 – 1,697 bp region from the start codon). PCR was conducted in 20 μl total reaction volume with 50–100 ng template DNA, 0.25 μM of forward and reverse primers, 0.1 mM of each dNTP, 0.8 μl Taq polymerase, 1.5 mM MgCl2, 5% DMSO, and 10x PCR buffer. PCR was comprised of the initial denaturation at 94°C (5 min), 40 cycles of 94°C (20 sec), 60°C (30 sec), and 72°C (1 min), followed by the final extension at 72°C (7 min). Before sequencing, 5 μl PCR product was purified as described previously . The purified PCR product was sequenced using BigDye Terminator v3.1 Sequencing Standard Kit (Applied Biosystems, Foster City, CA) and an ABI-3730 DNA Sequencer (Applied Biosystems, Foster City, CA).
A CAPS (cleaved amplified polymorphic sequences) marker was designed to screen for the mutant RHT-B1bE529K allele. The PCR product, amplified using the same primers and conditions described above, was digested with the restriction enzyme StyI. The wild-type allele produced an undigested 1,265 bp fragment and the RHT-B1bE529K allele produced two fragments of 1,153 bp and 112 bp, which were visualized by agarose gel electrophoresis (Additional file 2: Figure S4).
To genotype the six additional RHT-B1 mutations, an additional reverse primer was designed in the 3’ UTR (5’-CTTCTTCTTCTTCAAGAGCG-3’) and used along with the same forward primer described above, to amplify the 1,434 bp distal region of Rht-B1b and 115 bp of the 3’-UTR. PCR and sequencing were conducted as described above.
BC1F3 seeds derived from a single heterozygous BC1F2 plant (RHT-B1b/RHT-B1bE529K) were used for the RNA-seq experiments. The seeds were sown on germination paper as described above, kept at 4°C for 48 h, and moved to a growth chamber with 16 h light (230 μM m−2 s−1 intensity) at 22°C and 8 h dark at 18°C. Ten days after sowing (DAS), seedlings were transplanted into plastic pots and grown in the same growth chamber conditions.
We sampled tissue at four different developmental points as defined in the Zadoks’ scale : Z07 coleoptile (7 DAS), Z11 first leaf (10 DAS), Z49 peduncle when first awn appears (40 DAS), and Z52 peduncle when first spikelet appears (46 DAS). Using the CAPS marker described above, homozygous RHT-B1b and homozygous RHT-B1bE529K sister lines were selected and used for the RNA-seq experiments.
Four biological replicates were used per genotype at each time point. Samples were ground to a fine powder in liquid nitrogen, and total RNA was extracted using the SpectrumTM Plant Total RNA Kit (Sigma-Aldrich, St. Louis, MO). RNA-seq libraries were constructed using the TruSeq RNA Sample Preparation Kit v2 (Illumina, San Diego, CA) following the manufacturer’s instructions, and each library was indexed using NEXTflex-96™ DNA Barcode (Bioo Scientific, Austin, TX) for multiplexed sequencing reactions. The libraries were sequenced on the HiSeq 3000 platform (Illumina, San Diego, CA, USA) using the single-end 100 bp (SE 100) module at the UC Davis Genome Center.
Sequencing data were analyzed as previously described . Briefly, raw reads were trimmed using ‘Sickle’ (version 1. 33, https://github.com/najoshi/sickle with default parameters except -l 25 -q 25) and ‘Scythe’ (version 0. 991, https://github.com/vsbuffalo/scythe with default parameters). The processed reads were mapped to the IWGSC RefSeq v1.0 Chinese Spring reference genome (https://urgi.versailles.inra.fr/download/iwgsc/IWGSC_RefSeq_Assemblies/v1.0/) using ‘GSNAPl’ (version 2016-11-07, with default parameters except -m 5 -n 1 -N 1 -A sam, ). We excluded the D genome sequences from the reference before mapping the processed reads of the tetraploid lines (A and B genomes). Raw counts per transcript were calculated using ‘htseq-count’ (version 0.6.1, with default parameters except -m union -a 30 -t gene, ) with the GFF3 (General Feature Format version 3) file from the IWGSC RefSeq v1.0 annotation (199,680 features including high and low confidence genes and lncRNAs) and the SAM (Sequence Alignment/Mapping) files from each sample. Only transcripts with three or more reads in at least two biological replicates were retained using an R package ‘noleaven’ (https://github.com/topherconley/noleaven).
Raw counts were normalized using ‘DESeq2’ (version 1.18.1, ), and differentially expressed genes (DEGs) were identified between homozygous RHT-B1b and RHT-B1bE529K plants at each tissue/time point using ‘DESeq2’ and ‘EdgeR’ (version 3.20.6, ). Significant DEGs were selected using an adjusted false discovery rate threshold of P < 0.05 in both methods (FDR-adj P; ). Functional annotation and gene ontology (GO) enrichment analysis for DEGs were conducted as previously described . DEGs shared among different tissues were identified and visualized with UpSet plots  using an R package ‘UpSetR’ (version 1.4.0, ).
differentially expressed gene
FATTY ACYL-COENZYME A REDUCTASE
false discovery rate
GIBBERELLIN INSENSITIVE DWARF1
GAI, RGA, and SCR
GA response index
leucine heptad repeat
lateral organ boundary
PYROPHOSPHATE:FRUCTOSE 6-PHOSPHATE 1-PHOSPHOTRANSFERASE
targeting induced local lesions in genomes
Hauvermale AL, Ariizumi T, Steber CM. Gibberellin signaling: a theme and variations on DELLA repression. Plant Physiol. 2012;160(1):83–92.
Yamaguchi S. Gibberellin metabolism and its regulation. Annu Rev Plant Biol. 2008;59:225–51.
Hedden P, Thomas SG. Gibberellin biosynthesis and its regulation. Biochem J. 2012;444:11–25.
van de Velde K, Ruelens P, Geuten K, Rohde A, van der Straeten D. Exploiting DELLA signaling in cereals. Trends Plant Sci. 2017;22(10):880–93.
de Lucas M, Davière JM, Rodríguez-Falcón M, Pontin M, Iglesias-Pedraz JM, Lorrain S, Fankhauser C, Blázquez MA, Titarenko E, Prat S. A molecular framework for light and gibberellin control of cell elongation. Nature. 2008;451(7177):480–4.
Feng SH, Martinez C, Gusmaroli G, Wang Y, Zhou JL, Wang F, Chen LY, Yu L, Iglesias-Pedraz JM, Kircher S, Schäfer E, Fu XD, Fan LM, Deng XW. Coordinated regulation of Arabidopsis thaliana development by light and gibberellins. Nature. 2008;451(7177):475–9.
Bai MY, Shang JX, Oh E, Fan M, Bai Y, Zentella R, Sun TP, Wang ZY. Brassinosteroid, gibberellin and phytochrome impinge on a common transcription module in Arabidopsis. Nat Cell Biol. 2012;14(8):810–7.
Hou XL, Lee LYC, Xia KF, Yen YY, Yu H. DELLAs modulate jasmonate signaling via competitive binding to JAZs. Dev Cell. 2010;19(6):884–94.
Nakamura H, Xue YL, Miyakawa T, Hou F, Qin HM, Fukui K, Shi X, Ito E, Ito S, Park SH, Miyauchi Y, Asano A, Totsuka N, Ueda T, Tanokura M, Asami T. Molecular mechanism of strigolactone perception by DWARF14. Nat Commun. 2013;4:2613.
Fukazawa J, Teramura H, Murakoshi S, Nasuno K, Nishida N, Ito T, Yoshida M, Kamiya Y, Yamaguchi S, Takahashi Y. DELLAs function as coactivators of GAI-ASSOCIATED FACTOR1 in regulation of gibberellin homeostasis and signaling in Arabidopsis. Plant Cell. 2014;26(7):2920–38.
Ueguchi-Tanaka M, Nakajima M, Motoyuki A, Matsuoka M. Gibberellin receptor and its role in gibberellin signaling in plants. Annu Rev Plant Biol. 2007;58:183–98.
Murase K, Hirano Y, Sun TP, Hakoshima T. Gibberellin-induced DELLA recognition by the gibberellin receptor GID1. Nature. 2008;456(7221):459–63.
Sasaki A, Itoh H, Gomi K, Ueguchi-Tanaka M, Ishiyama K, Kobayashi M, Jeong DH, An G, Kitano H, Ashikari M, Matsuoka M. Accumulation of phosphorylated repressor for gibberellin signaling in an F-box mutant. Science. 2003;299(5614):1896–8.
Sun XL, Jones WT, Harvey D, Edwards PJB, Pascal SM, Kirk C, Considine T, Sheerin DJ, Rakonjac J, Oldfield CJ, Xue B, Dunker AK, Uversky VN. N-terminal domains of DELLA proteins are intrinsically unstructured in the absence of interaction with GID1/gibberellic acid receptors. J Biol Chem. 2010;285(15):11557–71.
Sun XL, Jones WT, Rikkerink EHA. GRAS proteins: the versatile roles of intrinsically disordered proteins in plant signalling. Biochem J. 2012;442:1–12.
Yoshida H, Hirano K, Sato T, Mitsuda N, Nomoto M, Maeo K, Koketsu E, Mitani R, Kawamura M, Ishiguro S, Tada Y, Ohme-Takagi M, Matsuoka M, Ueguchi-Tanaka M. DELLA protein functions as a transcriptional activator through the DNA binding of the INDETERMINATE DOMAIN family proteins. Proc Natl Acad Sci U S A. 2014;111(21):7861–6.
Hedden P. The genes of the Green Revolution. Trends Genet. 2003;19(1):5–9.
Youssefian S, Kirby EJM, Gale MD. Pleiotropic effects of the GA-insensitive Rht dwarfing genes in wheat. 2. Effects on leaf, stem, ear and floret growth. Field Crops Res. 1992;28(3):191–210.
Flintham JE, Börner A, Worland AJ, Gale MD. Optimizing wheat grain yield: Effects of Rht (gibberellin-insensitive) dwarfing genes. J Agr Sci. 1997;128:11–25.
Evenson RE, Gollin D. Assessing the impact of the Green Revolution, 1960 to 2000. Science. 2003;300(5620):758–62.
Peng J, Richards DE, Hartley NM, Murphy GP, Devos KM, Flintham JE, Beales J, Fish LJ, Worland AJ, Pelica F, Sudhakar D, Christou P, Snape JW, Gale MD, Harberd NP. ‘Green revolution’ genes encode mutant gibberellin response modulators. Nature. 1999;400(6741):256–61.
Evans LT. Feeding the ten billion. Plant and population growth. Cambridge: Cambridge University Press; 1998.
Guedira M, Brown-Guedira G, van Sanford D, Sneller C, Souza E, Marshall D. Distribution of Rht genes in modern and historic winter wheat cultivars from the Eastern and Central USA. Crop Sci. 2010;50(5):1811–22.
Rebetzke GJ, Richards RA, Fettell NA, Long M, Condon AG, Forrester RI, Botwright TL. Genotypic increases in coleoptile length improves stand establishment, vigour and grain yield of deep-sown wheat. Field Crops Res. 2007;100(1):10–23.
Rebetzke GJ, Appels R, Morrison AD, Richards RA, McDonald G, Ellis MH, Spielmeyer W, Bonnett DG. Quantitative trait loci on chromosome 4B for coleoptile length and early vigour in wheat (Triticum aestivum L.). Aust J Agr Res. 2001;52(11-12):1221–34.
Srinivasachary GN, Steed A, Simmonds J, Leverington-Waite M, Wang Y, Snape J, Nicholson P. Susceptibility to Fusarium head blight is associated with the Rht-D1b semi-dwarfing allele in wheat. Theor Appl Genet. 2008;116(8):1145–53.
Srinivasachary GN, Steed A, Hollins TW, Bayles R, Jennings P, Nicholson P. Semi-dwarfing Rht-B1 and Rht-D1 loci of wheat differ significantly in their influence on resistance to Fusarium head blight. Theor Appl Genet. 2009;118(4):695–702.
Wen W, Deng QY, Jia HY, Wei LZ, Wei JB, Wan HS, Yang LM, Cao WJ, Ma ZQ. Sequence variations of the partially dominant DELLA gene Rht-B1c in wheat and their functional impacts. J Exp Bot. 2013;64(11):3299–312.
Pearce S, Saville R, Vaughan SP, Chandler PM, Wilhelm EP, Sparks CA, Al-Kaff N, Korolev A, Boulton MI, Phillips AL, Hedden P, Nicholson P, Thomas SG. Molecular characterization of Rht-1 dwarfing genes in hexaploid wheat. Plant Physiol. 2011;157(4):1820–31.
Divashuk MG, Vasilyev AV, Bespalova LA, Karlov GI. Identity of the Rht-11 and Rht-B1e reduced plant height genes. Russ J Genet. 2012;48(7):761–3.
Bazhenov MS, Divashuk MG, Amagai Y, Watanabe N, Karlov GI. Isolation of the dwarfing Rht-B1p (Rht17) gene from wheat and the development of an allele-specific PCR marker. Mol Breed. 2015;35(11):213.
Wilhelm EP, Mackay IJ, Saville RJ, Korolev AV, Balfourier F, Greenland AJ, Boulton MI, Powell W. Haplotype dictionary for the Rht-1 loci in wheat. Theor Appl Genet. 2013;126(7):1733–47.
Li A, Yang W, Lou X, Liu D, Sun J, Guo X, Wang J, Li Y, Zhan K, Ling HQ, Zhang A. Novel natural allelic variations at the Rht-1 loci in wheat. J Integr Plant Biol. 2013;55(11):1026–37.
Chandler PM, Harding CA. ‘Overgrowth’ mutants in barley and wheat: new alleles and phenotypes of the ‘Green Revolution’ DELLA gene. J Exp Bot. 2013;64(6):1603–13.
Derkx AP, Harding CA, Miraghazadeh A, Chandler PM. Overgrowth (Della) mutants of wheat: development, growth and yield of intragenic suppressors of the Rht-B1c dwarfing gene. Funct Plant Biol. 2017;44(5):525–37.
Uauy C, Paraiso F, Colasuonno P, Tran RK, Tsai H, Berardi S, Comai L, Dubcovsky J. A modified TILLING approach to detect induced mutations in tetraploid and hexaploid wheat. BMC Plant Biol. 2009;9:115.
Kumar P, Henikoff S, Ng PC. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc. 2009;4(7):1073–82.
Krasileva KV, Vasquez-Gross HA, Howell T, Bailey P, Paraiso F, Clissold L, Simmonds J, Ramirez-Gonzalez RH, Wang XD, Borrill P, Fosker C, Ayling S, Phillips AL, Uauy C, Dubcovsky J. Uncovering hidden variation in polyploid wheat. Proc Natl Acad Sci U S A. 2017;114(6):E913–21.
Falconer D. Introduction to quantitative genetics. London: Oliver and Boyd; 1960.
Mo YJ, Howell T, Vasquez-Gross H, de Haro LA, Dubcovsky J, Pearce S. Mapping causal mutations by exome sequencing in a wheat TILLING population: a tall mutant case study. Mol Genet Genomics. 2018;293(2):463–77.
Domergue F, Vishwanath SJ, Joubés J, Ono J, Lee JA, Bourdon M, Alhattab R, Lowe C, Pascal S, Lessire R, Rowland O. Three Arabidopsis fatty acyl-coenzyme A reductases, FAR1, FAR4, and FAR5, generate primary fatty alcohols associated with suberin deposition. Plant Physiol. 2010;153(4):1539–54.
Rowland O, Domergue F. Plant fatty acyl reductases: Enzymes generating fatty alcohols for protective layers with potential for industrial applications. Plant Sci. 2012;193:28–38.
Puthoff DP, Sardesai N, Subramanyam S, Nemacheck JA, Williams CE. Hfr-2, a wheat cytolytic toxin-like gene, is up-regulated by virulent Hessian fly larval feeding. Mol Plant Pathol. 2005;6(4):411–23.
Xu C, Luo F, Hochholdinger F. LOB domain proteins: beyond lateral organ boundaries. Trends Plant Sci. 2016;21(2):159–67.
Ross S, Giglione C, Pierre M, Espagne C, Meinnel T. Functional and developmental impact of cytosolic protein N-terminal methionine excision in Arabidopsis. Plant Physiol. 2005;137(2):623–37.
McHale L, Tan X, Koehl P, Michelmore RW. Plant NBS-LRR proteins: adaptable guards. Genome Biol. 2006;7(4):212.
Boyes DC, Nam J, Dangl JL. The Arabidopsis thaliana RPM1 disease resistance gene product is a peripheral plasma membrane protein that is degraded coincident with the hypersensitive response. Proc Natl Acad Sci U S A. 1998;95(26):15849–54.
Fukuoka S, Saka N, Koga H, Ono K, Shimizu T, Ebana K, Hayashi N, Takahashi A, Hirochika H, Okuno K, Yano M. Loss of function of a proline-containing protein confers durable disease resistance in rice. Science. 2009;325(5943):998–1001.
Shiu SH, Bleecker AB. Plant receptor-like kinase gene family: diversity, function, and signaling. Sci STKE. 2001;2001(113):re22.
Li Y, Baldauf S, Lim EK, Bowles DJ. Phylogenetic analysis of the UDP-glycosyltransferase multigene family of Arabidopsis thaliana. J Biol Chem. 2001;276(6):4338–43.
Ross J, Li Y, Lim E, Bowles DJ. Higher plant glycosyltransferases. Genome Biol. 2001;2(2):reviews3004.1.
LeClere S, Tellez R, Rampey RA, Matsuda SP, Bartel B. Characterization of a family of IAA-amino acid conjugate hydrolases from Arabidopsis. J Biol Chem. 2002;277(23):20446–52.
Willats WG, Knox JP. A role for arabinogalactan-proteins in plant cell expansion: evidence from studies on the interaction of β-glucosyl Yariv reagent with seedlings of Arabidopsis thaliana. Plant J. 1996;9(6):919–25.
Park MH, Suzuki Y, Chono M, Knox JP, Yamaguchi I. CsAGP1, a gibberellin-responsive gene from cucumber hypocotyls, encodes a classical arabinogalactan protein and is involved in stem elongation. Plant Physiol. 2003;131(3):1450–9.
Lee KJ, Sakata Y, Mau SL, Pettolino F, Bacic A, Quatrano RS, Knight CD, Knox JP. Arabinogalactan proteins are required for apical cell extension in the moss Physcomitrella patens. Plant Cell. 2005;17(11):3051–65.
Kim JH, Tsukaya H. Regulation of plant growth and development by the GROWTH-REGULATING FACTOR and GRF-INTERACTING FACTOR duo. J Exp Bot. 2015;66(20):6093–107.
Badger MR, Price GD. The role of carbonic anhydrase in photosynthesis. Annu Rev Plant Phys. 1994;45:369–92.
Marsh JJ, Lebherz HG. Fructose-bisphosphate aldolases: an evolutionary history. Trends Biochem Sci. 1992;17(3):110–3.
Huber SC, Huber JL. Role and regulation of sucrose-phosphate synthase in higher plants. Annu Rev Plant Phys. 1996;47:431–44.
Tristan C, Shahani N, Sedlak TW, Sawa A. The diverse functions of GAPDH: Views from different subcellular compartments. Cell Signal. 2011;23(2):317–23.
Nielsen TH, Rung JH, Villadsen D. Fructose-2,6-bisphosphate: a traffic signal in plant metabolism. Trends Plant Sci. 2004;9(11):556–63.
Tovar-Méndez A, Miernyk JA, Randall DD. Regulation of pyruvate dehydrogenase complex activity in plant cells. Eur J Biochem. 2003;270(6):1043–9.
Nakajima K, Furutani I, Tachimoto H, Matsubara H, Hashimoto T. SPIRAL1 encodes a plant-specific microtubule-localized protein required for directional control of rapidly expanding Arabidopsis cells. Plant Cell. 2004;16(5):1178–90.
Palma JM, Sandalio LM, Corpas FJ, Romero-Puertas MC, McCarthy I, del Río LA. Plant proteases, protein degradation, and oxidative stress: role of peroxisomes. Plant Physiol Biochem. 2002;40(6-8):521–30.
Li YR, Kabbage M, Liu WD, Dickman MB. Aspartyl protease-mediated cleavage of BAG6 is necessary for autophagy and fungal resistance in plants. Plant Cell. 2016;28(1):233–47.
Kasprzewska A. Plant chitinases - regulation and function. Cell Mol Biol Lett. 2003;8(3):809–24.
Ishiguro S, Kawai-Oda A, Ueda J, Nishida I, Okada K. The DEFECTIVE IN ANTHER DEHISCENCE1 gene encodes a novel phospholipase A1 catalyzing the initial step of jasmonic acid biosynthesis, which synchronizes pollen maturation, anther dehiscence, and flower opening in Arabidopsis. Plant Cell. 2001;13(10):2191–209.
van Damme M, Huibers RP, Elberse J, van den Ackerveken G. Arabidopsis DMR6 encodes a putative 2OG-Fe (II) oxygenase that is defense-associated but required for susceptibility to downy mildew. Plant J. 2008;54(5):785–93.
Schillinger WF, Donaldson E, Allan RE, Jones SS. Winter wheat seedling emergence from deep sowing depths. Agron J. 1998;90(5):582–6.
Richards RA. The effect of dwarfing genes in spring wheat in dry environments. I. Agronomic characteristics. Aust J Agr Res. 1992;43(3):517–27.
Butler JD, Byrne PF, Mohammadi V, Chapman PL, Haley SD. Agronomic performance of Rht alleles in a spring wheat population across a range of moisture levels. Crop Sci. 2005;45(3):939.
Zhang J, Gizaw SA, Bossolini E, Hegarty J, Howell T, Carter AH, Akhunov E, Dubcovsky J. Identification and validation of QTL for grain yield and plant water status under contrasting water treatments in fall-sown spring wheats. Theor Appl Genet. in press.
Gooding MJ, Uppal RK, Addisu M, Harris KD, Uauy C, Simmonds JR, Murdoch AJ. Reduced height alleles (Rht) and Hagberg falling number of wheat. J Cereal Sci. 2012;55(3):305–11.
van de Velde K, Chandler PM, van der Straeten D, Rohde A. Differential coupling of gibberellin responses by Rht-B1c suppressor alleles and Rht-B1b in wheat highlights a unique role for the DELLA N-terminus in dormancy. J Exp Bot. 2017;68(3):443–55.
Wojciechowski T, Gooding MJ, Ramsay L, Gregory PJ. The effects of dwarfing genes on seedling root growth of wheat. J Exp Bot. 2009;60(9):2565–73.
McCaig TN, Morgan JA. Root and shoot dry matter partitioning in near-isogenic sheat lines differing in height. Can J Plant Sci. 1993;73(3):679–89.
Youssefian S, Kirby EJM, Gale MD. Pleiotropic effects of the GA-insensitive Rht dwarfing genes in wheat. 1. Effects on development of the ear, stem and leaves. Field Crops Res. 1992;28(3):179–90.
Willige BC, Ghosh S, Nill C, Zourelidou M, Dohmann EMN, Maier A, Schwechheimer C. The DELLA domain of GA INSENSITIVE mediates the interaction with the GA INSENSITIVE DWARF1A gibberellin receptor of Arabidopsis. Plant Cell. 2007;19(4):1209–20.
Wu J, Kong X, Wan J, Liu X, Zhang X, Guo X, Zhou R, Zhao G, Jing R, Fu X, Jia J. Dominant and pleiotropic effects of a GAI gene in wheat results from a lack of interaction between DELLA and GID1. Plant Physiol. 2011;157(4):2120-30.
Lou XY, Li X, Li AX, Pu MY, Shoaib M, Liu DC, Sun JZ, Zhang AM, Yang WL. Molecular characterization of three GIBBERELLIN-INSENSITIVE DWARF2 homologous genes in common wheat. PLoS One. 2016;11(6):e0157642.
Locascio A, Blázquez MA, Alabadí D. Genomic analysis of DELLA protein activity. Plant Cell Physiol. 2013;54(8):1229–37.
Claeys H, De Bodt S, Inzé D. Gibberellins and DELLAs: central nodes in growth regulatory networks. Trends Plant Sci. 2014;19(4):231–9.
Sun TP. The molecular mechanism and evolution of the GA–GID1–DELLA signaling module in plants. Curr Biol. 2011;21(9):R338–45.
Ubeda-Tomás S, Swarup R, Coates J, Swarup K, Laplaze L, Beemster GTS, Hedden P, Bhalerao R, Bennett MJ. Root growth in Arabidopsis requires gibberellin/DELLA signalling in the endodermis. Nat Cell Biol. 2008;10(5):625–8.
Achard P, Gusti A, Cheminant S, Alioua M, Dhondt S, Coppens F, Beemster GTS, Genschik P. Gibberellin signaling controls cell proliferation rate in Arabidopsis. Curr Biol. 2009;19(14):1188–93.
Gonzalez N, de Bodt S, Sulpice R, Jikumaru Y, Chae E, Dhondt S, van Daele T, de Milde L, Weigel D, Kamiya Y, Stitt M, Beemster GTS, Inze D. Increased leaf size: different means to an end. Plant Physiol. 2010;153(3):1261–79.
Iqbal N, Nazar R, Khan MIR, Masood A, Khan NA. Role of gibberellins in regulation of source–sink relations under optimal and limiting environmental conditions. Curr Sci. 2011;100(7):998–1007.
Ribeiro DM, Araújo WL, Fernie AR, Schippers JHM, Mueller-Roeber B. Action of gibberellins on growth and metabolism of Arabidopsis plants associated with high concentration of carbon dioxide. Plant Physiol. 2012;160(4):1781–94.
Flechner A, Gross W, Martin WF, Schnarrenberger C. Chloroplast class I and class II aldolases are bifunctional for fructose-1,6-biphosphate and sedoheptulose-1,7-biphosphate cleavage in the Calvin cycle. Febs Lett. 1999;447(2-3):200–2.
Uematsu K, Suzuki N, Iwamae T, Inui M, Yukawa H. Increased fructose 1,6-bisphosphate aldolase in plastids enhances growth and photosynthesis of tobacco plants. J Exp Bot. 2012;63(8):3001–9.
Cai B, Li Q, Liu F, Bi H, Ai X. Decreasing fructose 1,6-bisphosphate aldolase activity reduces plant growth and tolerance to chilling stress in tomato seedlings. Physiol Plant. 2017. https://doi.org/10.1111/ppl.12682.
Ishimaru K, Ono K, Kashiwagi T. Identification of a new gene controlling plant height in rice using the candidate-gene strategy. Planta. 2004;218(3):388–95.
Park AY, Canam T, Kang KY, Ellis DD, Mansfield SD. Over-expression of an arabidopsis family A sucrose phosphate synthase (SPS) gene alters plant growth and fibre development. Transgenic Res. 2008;17(2):181–92.
Davière JM, Achard P. A pivotal role of DELLAs in regulating multiple hormone signals. Mol Plant. 2016;9(1):10–20.
Oh E, Zhu JY, Bai MY, Arenhart RA, Sun Y, Wang ZY. Cell elongation is regulated through a central circuit of interacting transcription factors in the Arabidopsis hypocotyl. Elife. 2014;3:e03031.
Shan QW, Wang YP, Li J, Gao CX. Genome editing in rice and wheat using the CRISPR/Cas system. Nat Protoc. 2014;9(10):2395–410.
Zhang Y, Liang Z, Zong Y, Wang YP, Liu JX, Chen KL, Qiu JL, Gao CX. Efficient and transgene-free genome editing in wheat through transient expression of CRISPR/Cas9 DNA or RNA. Nat Commun. 2016;7:12617.
Rohland N, Reich D. Cost-effective, high-throughput DNA sequencing libraries for multiplexed target capture. Genome Res. 2012;22(5):939–46.
Zadoks JC, Chang TT, Konzak CF. A decimal code for the growth stages of cereals. Weed Res. 1974;14:415–21.
Pearce S, Kippes N, Chen A, Debernardi JM, Dubcovsky J. RNA-seq studies using wheat PHYTOCHROME B and PHYTOCHROME C mutants reveal shared and specific functions in the regulation of flowering and shade-avoidance pathways. BMC Plant Biol. 2016;16(1):141.
Wu TD, Nacu S. Fast and SNP-tolerant detection of complex variants and splicing in short reads. Bioinformatics. 2010;26(7):873–81.
Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc B Met. 1995;57(1):289–300.
Lex A, Gehlenborg N, Strobelt H, Vuillemot R, Pfister H. UpSet: visualization of intersecting sets. Ieee T Vis Comput Gr. 2014;20(12):1983–92.
Conway JR, Lex A, Gehlenborg N. UpSetR: an R package for the visualization of intersecting sets and their properties. Bioinformatics. 2017;33(18):2938–40.
We thank the International Wheat Genome Sequencing Consortium for pre-publication access to the RefSeq v1.0 Chinese Spring reference sequences.
This project was supported by the Agriculture and Food Research Initiative Competitive Grant 2017-67007-25939 (WheatCAP) from the USDA National Institute of Food and Agriculture and the Howard Hughes Medical Institute. Y. Mo is a Howard Hughes Medical Institute’s International Student Research fellow and a Monsanto Beachell-Borlaug International scholar. The funding bodies played no role in the design, execution, interpretation of the data, or writing of the manuscript.
Availability of data and materials
RNA-seq reads generated in this study are available in NCBI’s Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra, Bioproject PRJNA435752). Seeds of the RHT-B1bE529K mutant have been deposited in the USDA-ARS National Plant Germplasm System (NPGS) as accession PI 687144 (https://npgsweb.ars-grin.gov/gringlobal/accessiondetail.aspx?id=1961292).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Effects of additional induced mutations in the C-terminal of RHT-B1 on plant height. Table S2. RNA-seq data summary of 32 libraries. Table S3. Differentially expressed genes between homozygous RHT-B1b and RHT-B1bE529K plants (FDR-adj P < 0.05 in both DESeq2 and EdgeR). Table S4. Differentially expressed genes between homozygous RHT-B1b and RHT-B1bE529K plants (FDR-adj P < 0.05 in DESeq2 or EdgeR). Table S5. EMS-induced nonsynonymous mutations in conserved RHT-B1 motifs. (XLSX 551 kb)
Figure S1. Number of differentially regulated genes (DEGs) between plants encoding RHT-B1b and RHT-B1bE529K identified by DESeq2 or EdgeR in (a) coleoptile, (b) first leaf, (c) Z49 peduncle, and (d) Z52 peduncle. Upregulated and downregulated DEGs are in red and blue, respectively. Figure S2. Expression of RHT-B1 (TraesCS4B01G043100) in plants encoding RHT-B1b and RHT-B1bE529K in coleoptile, first leaf, Z49 and Z52 peduncle. Normalized counts were obtained from the RNA-seq data. NS: not significant according to the t-test. Error bars indicate ±1 standard error. Figure S3. UpSet plots showing the numbers of differentially expressed genes (DEGs) between homozygous RHT-B1b and RHT-B1bE529K plants that are unique or overlapping among four different tissues. DEGs significant (FDR-adj P < 0.05) in either DESeq2 or EdgeR were counted. (a) Number of upregulated DEGs. (b) Number of downregulated DEGs. Figure S4. An example genotyping result of the CAPS marker developed to detect the RHT-B1bE529K mutation. W: homozygous RHT-B1b, M: homozygous RHT-B1bE529K, H: heterozygote (DOCX 312 kb)