Next Article in Journal
Oxidative Stress in Type 2 Diabetes: Impacts from Pathogenesis to Lifestyle Modifications
Next Article in Special Issue
Physiological and Transcriptomic Analyses Reveal Commonalities and Specificities in Wheat in Response to Aluminum and Manganese
Previous Article in Journal
Ex Vivo Chromosomal Radiosensitivity Testing in Patients with Pathological Germline Variants in Breast Cancer High-Susceptibility Genes BReast CAncer 1 and BReast CAncer 2
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
CIMB
Volume 45
Issue 8
10.3390/cimb45080419
cimb-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Wheat Water-Soluble Carbohydrate Remobilisation under Water Deficit by 1-FEH w3

by
Nusrat Khan
1,2,
Jingjuan Zhang
1,
Shahidul Islam
1,2,
Rudi Appels
3 and
Bernard Dell
1,*
1
Centre for Crop and Food Innovation, Food Futures Institute, Murdoch University, 90 South Street, Murdoch, WA 6163, Australia
2
Department of Plant Sciences, North Dakota State University, Fargo, ND 58102, USA
3
Faculty of Science, University of Melbourne, Parkville, VIC 3010, Australia
*
Author to whom correspondence should be addressed.
Curr. Issues Mol. Biol. 2023, 45(8), 6634-6650; https://doi.org/10.3390/cimb45080419
Submission received: 18 July 2023 / Revised: 3 August 2023 / Accepted: 8 August 2023 / Published: 11 August 2023
(This article belongs to the Special Issue Advanced Research in Wheat Genome and Breeding)
Browse Figures
Versions Notes

Abstract

:
Fructan 1-exohydrolase (1-FEH) is one of the major enzymes in water-soluble carbohydrate (WSC) remobilisation for grains in wheat. We investigated the functional role of 1-FEH w1, w2, and w3 isoforms in WSC remobilisation under post-anthesis water deficit using mutation lines derived from the Australian wheat variety Chara. F1 seeds, developed by backcrossing the 1-FEH w1, w2, and w3 mutation lines with Chara, were genotyped using the Infinium 90K SNP iSelect platform to characterise the mutated region. Putative deletions were identified in FEH mutation lines encompassing the FEH genomic regions. Mapping analysis demonstrated that mutations affected significantly longer regions than the target FEH gene regions. Functional roles of the non-target genes were carried out utilising bioinformatics and confirmed that the non-target genes were unlikely to confound the effects considered to be due to the influence of 1-FEH gene functions. Glasshouse experiments revealed that the 1-FEH w3 mutation line had a slower degradation and remobilisation of fructans than the 1-FEH w2 and w1 mutation lines and Chara, which reduced grain filling and grain yield. Thus, 1-FEH w3 plays a vital role in reducing yield loss under drought. This insight into the distinct role of the 1-FEH isoforms provides new gene targets for water-deficit-tolerant wheat breeding.

1. Introduction

Water deficit during grain filling reduces photosynthesis and increases the dependency of wheat crops on water-soluble carbohydrates (WSC) stored in the stem for grain filling. Wheat fructans, a major component of WSC, must be depolymerised (hydrolysed) before remobilisation to the grain, and this is regulated by a class of enzymes called fructan exohydrolases (FEHs) [1]. So far, two major types of FEHs have been characterised and cloned: 1-FEHs that preferentially degrade inulin-type fructans having β-(2,1) linkage [2,3] and 6-FEHs that degrade levan-type fructans with β-(2,6) linkage [2,4,5].
1-FEH catalyses fructan breakdown and is considered a β-(2,1)-trimmer during active fructan biosynthesis [2]. Elevated 1-FEH activities lead to the breakdown of fructan, thus releasing more sucrose and fructose when the demand for sucrose is high for grain filling [2,6,7,8,9]. Three isoforms of the 1-FEH gene have been reported, namely 1-FEH w1 (1-FEH-6A), 1-FEH w2 (1-FEH-6D), and 1-FEH w3 (1-FEH-6B). Previous findings from our laboratory suggested that the dominant expression of 1-FEH w3 might facilitate higher fructan remobilisation under terminal water deficit [10]. However, the individual functions of the three isoforms are still unknown due to the complex nature of the plant’s response to water-deficit stress.
Knockout mutation lines of each of the three isoforms (1-FEH w1, w2, and w3) were developed by heavy ion bombardment at CSIRO, Canberra, Australia, in the background of Chara [11], which is a moderately water-deficit-tolerant bread wheat cultivar. Mutating particular gene/s of interest is a common genetic approach to studying biological phenomena and defining gene function [12]. Knocking out genes is a straightforward method of producing gene mutation lines, and this approach has been used in many studies [13]. While mutation lines are essential for forward genetic studies of gene function, recent advancement in sequencing the wheat genome [14] also enables their application in reverse genetics. However, in polyploidy plant species like wheat, which contain multiple copies of the same gene (homologous gene) encoded by each of the ancestral genomes, mutational gene inactivation/deletion is very important to determine whether different copies of the same or related genes are functionally active or redundant [11].
The use of mutational gene inactivation/deletion in polyploid plant species is complicated [15], and for wheat, the challenge is that it possesses three ‘subgenomes’, designated ‘A’, ‘B’, and ‘D’, which are approximately 17 Gb in total size and 40 times the size of the rice genome. In addition to its large size, the wheat genome contains a high percentage of repetitive segments (80–90%) [16,17]. Generally, most wheat genes have three homoeologous copies on the A, B, and D subgenomes, leading to functional redundancy [15,18]. This gene redundancy causes ambiguities in gene knockout strategies, and, as a result, wheat is tolerant of mutations. Although such tolerance to mutation is generally desirable to keep the overall phenotype functional through the generations, genomic redundancy is challenging for producing loss-of-function mutants for a gene of interest [19]. Moreover, heavy ion bombardment causes different mutation types, including base substitutions, deletions, and chromosomal alteration [20,21]. Therefore, it is essential to investigate how gene inactivation works, for example, whether inactivation is by DNA sequence alterations or complete deletion and whether the mutation has caused any polymorphism in the coding sequence or an indirect modification of gene function. Most importantly, any closely linked genes, which might influence WSC remobilisation, need to be identified.
Whole-genome sequencing is a straightforward and powerful tool to directly identify mutagen-induced nucleotide changes linked to the causal mutation and recognise specific genomic locations of the mutated region [22]. High-density SNP arrays have been optimised and successfully used for genetic studies of several economically important crops [23,24,25,26,27,28]. In wheat, the construction of high-resolution genetic maps based on the microarray hybridisation-based technique called Diversity Arrays Technology (DArT) was developed. It enabled high-throughput genotyping without relying on sequence information [29]. Advancements in next-generation sequencing (NGS) technology have made it possible to screen for SNPs which have been used to develop high-throughput SNP-typing platforms such as BeadExpress [30], KASPar [31], Infinium [32], and designer SNP arrays [33]. A 9K SNP wheat chip detected genomic regions targeted by breeding and improvement selection [32,34]. Wang et al. [35] reported a 90K SNP iSelect array that comprises approximately 90,000 gene-associated SNPs with the potential to provide dense coverage of the wheat genome. This technology offers a resource for diversity studies and high-resolution dissection of complex traits in wheat.
The present study used a 90K iSelect array to survey SNPs across the mutation lines and the parent cultivar Chara. Data from the 90K analysis were compared with the whole wheat genome assembly to investigate the 1-FEH gene mutations, map the mutation region in the chromosome, and annotate the mutational region. Accordingly, the response of each mutation line to WSC remobilisation and how that is influenced by terminal water stress was characterised. The knowledge generated from this research regarding the specific roles of the 1-FEH gene isoforms provides the gene target for wheat breeding, particularly post-anthesis water-deficit conditions.

2. Materials and Methods

2.1. Plant Material

A set of mutation lines of the three isoforms of 1-FEH (1-FEH w1, w2, and w3), previously developed in the Australian wheat cultivar Chara, was bulked up and used in this experiment (Table 1). Four germplasms were included in the study: one line each of w1 mutant, w2 mutant, and w3 mutant was selected together with the parent Chara.

2.2. Backcrossing for 90K SNP Analysis

The mutated lines were backcrossed with the parental cultivar Chara in a Murdoch University glasshouse to check the mutation’s stability and reduce background genetic noise due to mutations elsewhere in the genome. The parent Chara was used as the female parent. The F1 seeds were collected on maturation at 65 days after pollination and were used to develop the iSelect 90K SNP array.

2.3. Genotypic DNA Extraction

DNA was extracted using an Illustra DNA Extraction Kit Phytopure (illustra Nucleon Phytopure Genomic DNA Extraction Kits, GE Healthcare, Amersham Place, Buckinghamshire, UK) following the standard protocol supplied by the manufacturer (Illustra Nucleon Phytoprue for small sample, 0.1 g RPN8510). One microlitre of extracted DNA was run on 1% agarose gel and standard lambda DNA (Axyzen Life Sciences, Arizona, AZ, USA) to check the quality. The DNA concentration was quantified using NanoDrop (ND 1000 Spectrophotometer, Thermo Fisher Scientific, Waltham, MA, USA). The DNA samples were then diluted to 20 ng µL−1 in 50 µL EB buffer (10 mM Tris-HCl pH 8.5) that stabilised the DNA.

2.4. The 90K SNP Chip Analysis

The Infinium 90K SNP iSelect platform was used for genotyping using BeadStation and iScan according to the manufacturer’s protocol from Illumina. Deletions in the HIB (Heavy Ion Bombardment) mutants were identified from output files using Illumina’s GenomeStudio v2011.1 (Illumina Inc., San Diego, CA, USA) and custom MS Excel macros and/or perl scrip, where shifts away from the euploid SNP cluster position indicate a deleted locus. Standard cluster files provided with Infinium products identify the expected intensity levels of genotype classes for each SNP. Raw intensity data for the three FEH HIBs and three FEH HIB × Chara backcrosses (as well as pre-existing data for Chara and other unrelated CSIRO HIBs) were loaded into Genomestudio, and the NormTheta (x-axis coordinate) and NormR (y-axis coordinate) values from the SNP cluster plots were exported. Shifts in cluster position were calculated between the euploid Chara and each mutant sample. The differences in Norm Theta and in NormR were calculated following the formula: Theta [Chara] − Theta [samples].
Shifts in cluster position were also calculated between the average HIB position (96 unrelated HIB mutants genotyped of Chara previously used in a different project and each of the FEH mutant samples). The 96 HIBs were used as an extra tool for helping to identify putative deletions due to the presence of Chara biotypes (i.e., different versions of ‘Chara’). Differences in Norma Theta and in Norma R were calculated using the formula: Theta [Av. HIB] − Theta [sample]. The 96 HIBs were used as an additional tool for helping identify deletions due to the presence of Chara biotypes (i.e., different versions of “Chara”). SNP loci putatively tracking deleted segments were identified in each FEH HIB mutant. SNPs with a minimum cluster shift between Chara and the HIB mutant of 0.1 along the Norm Theta axis were considered to be markers for deletion.
In order to help clarify the position of the FEH genes in chromosomes 6A, 6B, and 6D, sequences similar to the FEH gene sequence (FJ184991) were aligned in Gydle against the Chinese Spring genome assembly IWGSC. Regions containing 1-FEH-like sequences were identified from the alignments. SNPs were retained if they showed evidence for deletion and/or if they were mapped to 6A, 6B, or 6D. Markers considered as flanking deleted segments were identified for each HIB line by looking at the shift in cluster position. All SNPs were considered, even those with cluster shifts <0.2 Norm Theta and/or <0.5 Norm R. These markers were visually assessed in GenomeStudio to assess polymorphisms relative to Chara. To validate that the deletion detected was in fact located on the subgenome of interest (e.g., 6A) rather than a homoeologoue (e.g., 6B or 6D), those single-site SNPs located near either end of the deletion were visually assessed in GenomeStudio. A ‘null’ allele, where there is no fluorescence due to deletion of the hybridisation site, was considered to be indicative of the deletion occurring on the subgenome.

2.5. The 90K SNP Data Analysis

For mapping and analysing the SNP data, UGENE—Integrated Bioinformatics Tools (http://ugene.net/06/04/2018, accessed on 15 July 2023), UniprotKB (http://www.uniprot.org/, accessed on 15 July 2023), and CLC Main Workbench—QIAGEN Bioinformatics were used. CDS sequences and exons of 1-FEHw1, w3, and w2 were used for blast analysis against the IWGSC chromosomes 6A, 6B, and 6D. The public RNAseq database exVIP (http://www.wheat-expression.com/, accessed on 9 October 2022) was used to study the gene expression of these gene copies. Sequence annotations, including GO terms and Pfam domain information along with structural predictions, were performed using the BioMart in EnsmblPlants (http://plants.ensembl.org/biomart/martview/, accessed on 15 July 2023).

2.6. Glasshouse Experimental Setup

The experiment was set up as a Randomised Complete Block Design (RCBD) with 3 factors (Line × Treatment × Sampling stage) and 3 replications (Table 1), resulting in a total of 324 pots. Pots were initially randomised following the DiGGer design of R statistical software. Pots were rotated every alternate day to ensure equal solar radiation to all plants. The potting mix was the same as that used to grow wheat in glasshouses at Murdoch University [10]. The average temperature thresholds for night and day were 15 °C and 25 °C, respectively, and for humidity, they were 80% and 50%, respectively. The day length was not controlled. Hence, the photoperiod changed with the season. The day lengths ranged from 12 h in June to 14 h in November. Plants were watered using an irrigation mat for 5 min each, 3 times daily before anthesis.

2.7. Imposing Water Deficit

The main stem was identified at the mid-tillering stage, and the flowering date was recorded for every plant by tagging them on the first flowering date. When 50% of the plants were flowering, half were exposed to water deficit until grain maturity by turning off the water supplied via irrigation mats, and the other half were kept well-watered. In this stage, pots were arranged on four benches (two for water-deficit plants and two for well-watered plants). Weights of the unwatered pots were taken each second day to determine the extent of water loss, and the well-watered plants were weighed once a week to confirm the well-watered condition. Water was turned off to the well-watered pots when plants became yellow on maturity.

2.8. Sampling

The first sample was collected at the AR5 stage [36], which was followed by the second sample collected when the ear fully emerged, and the following six samplings for the preliminary and the main experiments were carried out at 4-day intervals after initiation of the water-deficit treatment (Table 1). The final sample was collected at maturity. Sampling for WSC analysis was undertaken between 11.00 and 17.00 h to study the shoot carbohydrate reserve accumulation and remobilisation pattern during grain filling.
At each sampling stage, three organs (the stem, leaf sheath, and blade) for each sample were collected from the main tiller of three plants of each pot as one sample and immediately put on dry ice. The remaining two plants were left for maturity and seed collection. Three biological replicates of each line were collected at each sampling stage. Samples were stored at −20 °C.
Harvest data were collected at maturity from the main stem and the tillers of 9 plants from 3 pots per line per treatment. The following parameters were recorded: kernel number (KN) per main spike, KN per plant, and spike number per plant. Grain was oven dried for 10 h at 60 °C, and grain weight per main spike, grain weight per plant, and thousand grain weight (TGW) of the main spike were recorded.

2.9. Carbohydrate Analysis

Samples for carbohydrate analysis were freeze dried (−50 °C) for 72 h and then oven-dried (70 °C) for 24 h. Frozen samples were chopped manually and ground using a tissue-lyser into less than 5 mm pieces in a cold room and manually divided into two parts. One part was stored at −80 °C for the RNA study, and the other was stored at −20 °C for WSC analysis. Carbohydrate analysis was performed on 170 mg of each ground sample using boiling deionised water for carbohydrate extraction followed by colorimetry using the anthrone reagent [10,37].

2.10. Calculation of WSC Remobilisation

The WSC concentration was highest at 12 DAA in all mutation lines and Chara. The WSC started declining from 12 DAA, and the rate of decline was elevated up to 22 DAA. Therefore, WSC remobilisation was calculated from 12 to 22 DAA.
WSC remobilisation rate (%dw × day−1) = (Stem WSC concentration (%dw) at 12 DAA − Stem WSC concentration (%dw) at 22 DAA)/10

2.11. Data Analysis

Analysis of variance was performed using IBM SPSS Statistics V21.0. Post hoc comparisons were conducted using Tukey’s Multiple Range Test at p < 0.05. Graphs were generated using Sigmaplot 13.0.

3. Results

3.1. Detecting the Mutation Region

From iSelect illumine 90K SNP array analysis, putative deletions were identified in the FEH-Heavy Ion Bombardment (HIB) mutations encompassing the FEH genomic regions (Table 2). The 6A mutants were larger deletions than the 6B and 6D ones. For each HIB line, markers at each end of the deletion were identified to estimate the size of the deletion.
SNPs with a minimum cluster shift between Chara and the HIB mutant of 0.1 along the Norm theta axis or 0.5 along the NormR axis were considered to be putative markers for a deletion. Deletions were confirmed by the Chara X mutant backcross, which were located at a ‘HET’ position, i.e., halfway between the mutant and the euploid Chara. The Infinium assay probes were bound to multiple sites in the genome due to the presence of homologues (e.g., 6A, 6B, and 6D) and paralogues, and polymorphism on any one of these sites was detected based on cluster position. The deletion of a segment of DNA in an HIB mutant affected the number of loci to which the probe could hybridise. This resulted in a shift in the NormTheta and/or NormR axes relative to the euploid sample, i.e., Chara (Supplementary Figure S1 as an example). When the probe had only one hybridisation site in the genome (e.g., hybridises to 6A but not 6B or 6D), the deletion of a hybridisation site resulted in a ‘null’ allele, i.e., a lack of fluorescent signal (Supplementary Figure S2 as an example). Those single-site SNPs, located near either end of the deletion, were visually assessed in Genome Studio. A ’null’ allele was considered to be indicative of the deletion occurring on the subgenome. Biotypes were observed between the Chara used to make HIBs and Chara used to backcross with the HIBs (Supplementary Figure S3 as an example).

3.2. Mapping of Mutation Regions in Chromosomes

High-density polymorphism clusters, considered to be mutated regions, were found on chromosomes 6A, 6B, and 6D, respectively (Table 2). By detecting the approximate genomic location of the mutation, the affected genes were subsequently identified and mapped. The deletion size was 28,246,961 bp (from 12,836,428 to 41,083,389 bp) on chromosome 6A (Supplementary Figure S4); 4,876,122 bp (from 54,419,374 to 59,295,496 bp) on chromosome 6B (Supplementary Figure S5) and 1,600,754 bp (from 28,682,554 to 30,283,308 bp) (Supplementary Figure S6) on chromosome 6D.
The entire list of genes in these intervals is provided in Supplementary Tables S1–S3.

3.3. Gene Annotation at the Mutated Region and Functional Analysis

The genes for both the targets and non-targets within the mutation regions of chromosomes 6A, 6B, and 6D were annotated with IWGSC Refseq v2.1 assembly using CLC Genomic Workbench. For the mutation region of 1-FEH w1, 479 high-confidence (HC) genes were identified on chromosome 6A, including the 1-FEH w1 (Supplementary Table S1, Supplementary Figure S4), 90 HC genes were identified in the 1-FEH w3 mutation region on chromosome 6B including the 1-FEH w3 gene (Supplementary Table S2, Supplementary Figure S5) and 144 HC genes were identified in chromosome 6D’s mutated region along with 1-FEH w2 (Supplementary Table S3, Supplementary Figure S6). Functional annotation of the identified genes was carried out using the BioMart tool available on EnsemblPlants and by manually blasting against the IWGSC gene functional annotation V.1. A total of six, two, and five genes have been identified on the mutated regions of chromosome 6A, 6B and 6D, with similar function to the 1-FEH w1, 1-FEH w3 and 1-FEH w2, respectively. Further functional confirmation of uncharacterised non-target genes became essential to validate that the phenotypic change in the mutated plant is due to the deletion of the target genes, namely 1-FEHw1, w2, w3. Thus, the function prediction of uncharacterised proteins was carried out following computational biology approaches. However, most of the genes were identified as encoding uncharacterised proteins based on sequence blasts in the available databases. To further predict the possible interference of the uncharacterised suite of genes with the function of 1-FEH, gene expressional behaviours were compared from the available datasets at ExpVIP (http://www.wheat-expression.com, accessed on 9 October 2022). Results demonstrated that all of the genes had either very little or no expression compared to the corresponding 1- FEH isoform at the reproductive stage (Supplementary Figures S7–S9).
Detailed interpretation of the functional properties of non-target genes indicated that none of them were directly related to 1-FEH w1, w2, and w3 gene functions (Table 3). So, it can be concluded that the phenotypic differences in the mutated lines in relation to WSC mobilisation can be considered to be due to the deletion of 1-FEH isoforms.

3.4. Observation of the Water-Deficit Regime in Mutated Lines

By 14 days after the water-deficit treatment was imposed, the soil moisture content had fallen to 22%, and plants showed signs of mild stress (first visible at 10 days) from leaf wilting. Plants reached the permanent wilting point (13% soil moisture content) after 20 days of water being withheld. Water deficit promoted early senescence compared to the well-watered plants (Figure 1).

3.5. Water Deficit Promoted Early Remobilisation of WSC

In well-watered plants, stem WSC levels peaked at 27 days after anthesis (DAA). There were no significant differences among parent and mutation lines (w1, w2, and w3) (Figure 2). In contrast, under water deficit, WSC peaked at 12 DAA, which is 15 days earlier than in the well-watered plants. Accordingly, the highest stem WSC level appeared significantly lower in water-deficit plants compared to the well-watered ones. In well-watered plants, the highest level of WSC was in Chara, with 35%, which was followed by w3, w2, and w1 mutation lines with 33.9, 29.9, and 34.5%, respectively. On the other hand, in water-deficit plants, the highest level of stem WSC under water deficit was 31% in Chara followed by 29.9, 28.5, and 26.3% in w2, w1, and w3 mutation lines.

3.6. Timing of Decline in WSC Differed among Mutation Lines

There was a clear variation in the pattern of decline in WSC concentration in plants subjected to water deficit. The WSC remobilisation rate was slower in the 1-FEH w3 mutant line compared to Chara and the other two mutation lines. This was evident from the fact that the 1-FEH w3 mutation line maintained (at 22 DAA) a significantly higher (11.1%) level of stem WSC than Chara (5.7%, p < 0.07), 1-FEH w2 (6.7%) and 1-FEH w1 (9.5%) mutant lines (Figure 2). However, there were no significant mean differences between Chara and the mutation lines in well-watered plants. Accordingly, under water-deficit conditions, between 12 DAA (peak value for WSC) and 22 DAA (permanent wilting point), the remobilisation rate of WSC was significantly (p < 0.05) lower in the w3 mutation line (15.2%) compared to Chara (25.2%) and was slightly slower (not significantly) than mutation lines w2 (21.8%) and w1 (20.5%) (Figure 3).

3.7. WSC Concentration and Grain Development

There was a negative correlation between stem WSC concentration and developing grain weight under water deficit (between 12 and 32 DAA), but it was not the case in well-watered plants (Figure 2). In water-deficit plants, the WSC concentration began to decrease from 12 DAA, and the developing grain weight started to increase, indicating WSC was being remobilised to the grain. Even though at 12 DAA, there was no significant difference in the main spike seed weight, by 32 DAA, the seed weight was significantly (p < 0.05) lower in the 1-FEH w3 mutation line (0.56 g/ear) compared to Chara (0.87 g/ear) (Figure 2). However, there were no significant differences in 1-FEH w2 and 1-FEH w1 mutation lines compared with Chara for these developmental stages.

3.8. Grain Weight at Maturity Was Influenced by 1-FEH Isoforms under Water-Deficit Condition

Water deficit reduced (p < 0.05) the grain weight of the 1-FEH w3 mutant line compared to Chara and the other two mutation lines (Figure 4). Under water deficit, grain wt/main spike was the lowest in the w3 mutation line (0.41 g ± 0.03) and the highest in Chara (0.71 g ± 0.03 SE). The other two mutation lines 1-FEH w2 and 1-FEH w1 produced 0.63 g ± 0.05 and 0.58 g ± 0.08 grain/main spike, respectively. In contrast, in the well-watered treatment, the 1-FEH w3 mutant line and Chara had similar seed wt/main spike (1.43 g ± 0.09 and 1.5 g ± 0.04, respectively) followed by w2 (1.02 g ± 0.02) and w1 (0.96 g ± 0.16). The TGW (thousand grain weight), a major determinant of the yield, was significantly lower (p < 0.05) in the w3 mutation line (17.47 g) compared to Chara (22.53 g) under water deficit, but the w2 and w1 mutation lines did not show any significant difference with Chara. Under the well-watered condition, there were no significant differences among the mutation lines and Chara.

3.9. Correlation Analysis of Grain Weight with WSC Level and Remobilisation Rate

To understand potential factors influencing grain weight under water deficit, two correlation analyses were carried out: (1) main stem seed weight vs. the highest level of WSC; and (2) grain weight per main stem vs. WSC remobilisation rate. The correlation co-efficient of grain weight per main stem with sugar remobilisation rate was 0.78 (r = 0.78) (p < 0.01), while it was 0.66 (r = 0.66) (p < 0.05) (Supplementary Figure S10) in the case of the highest WSC level. Thus, grain weight per main stem was more strongly correlated with sugar remobilisation than the WSC level.

4. Discussion

4.1. Characterisation of 1-FEH Gene Mutations

Characterisation of the 1-FEH mutation lines was carried out using high-density SNP genotyping array and bioinformatics approaches that provided detailed information regarding the physical modification of mutated segments on the chromosomes and the nature of the mutation. This information was crucial to explain the functional loss of the 1-FEH genes in the mutation lines. It is important to note that the backcrossing step was necessary for the comparative analysis of mutants obtained from all mutagenesis screens, irrespective of the type of mutant identification strategy used, or type of mutagen used, or the organism used. Deletion mutations typically do not revert. To identify mutation regions, we focused on a high-density variation on a single chromosome for each mutant which corresponds to regions of high mutagen-induced damage that were not removed during backcrossing and therefore most likely genetically linked to the causal mutation [38]. We therefore focused our attention on these physical regions to identify candidate mutations. This approach offered the opportunity to compile and compare the parent line, mutation lines, and backcrossed lines and to identify the mutation region.

4.2. Chromosomal Mapping of Mutation Regions and Candidate Genes

To validate the phenotypic expression of gene function loss due to mutation, we characterised the structural changes that happened in the genes and adjacent regions of the chromosomes. The high-density SNPs clusters identified on chromosomes 6A, 6B, and 6D, respectively, clearly demonstrated the physical location of the mutated regions of the 1-FEH mutation lines. The advancement of wheat whole genome sequencing [14] has furthermore enhanced the capacity for a detailed characterisation the chromosomal mutated regions by comparing the mutated line with the parent line. The sizes of the mutated/deleted regions were variable among the mutation lines. The largest deletion region was identified in the 1-FEH w1 mutation line, which was 28,246,961 bp on chromosome 6A. Deletion sizes of mutation lines 1-FEH w3 and 1-FEH w2 were 4,876,122 bp on chromosome 6B and 1,600,754 bp on chromosome 6D, respectively. Generally, mutation by heavy ion irradiation mutates a larger region compared to other ion radiation or chemical mutagenesis, and it can be highly variable. Knocking out 1-FEH isoforms by using CRISPR-Cas9 genome editing technology or functional interference by using RNAi technology could provide the opportunity of precise and accurate investigation. However, a very high level of sequence similarity between 1-FEH isoforms made it challenging to use these technologies. Moreover, several other members of the FEH gene family have high sequence similarities with 1-FEH isoforms. These increase the possibility of off-target mutation/knockout using CRISPR-Cas9 or RNAi.

4.3. Annotation of the Deleted Region

Mapping analysis demonstrated that the mutation affected significantly longer regions than the target gene regions of 1-FEH w1, w3, and w2 on chromosomes 6A, 6B, and 6D, respectively. Thus, annotating the whole mutated region was crucial to understanding the effect of deletion on the measured phenotypes. This analysis provided us with the functional information of target and non-target mutated genes. The investigation was also carried out to identify whether any non-target genes located within the mutated region influenced the target gene’s function.
A total of 379, 89, and 144 non-target high-confidence genes were identified in the mutated regions of chromosomes 6A, 6B, and 6D, respectively. These findings indicated that the mutation targeting the functional loss of 1-FEH w1, w2, and w3 affected a considerable number of genes at the adjacent chromosomal location.

4.4. Function Prediction for Non-Target Genes in Deleted Regions

Following the assembly of the mutation regions, further investigation on the functional annotation of non-target genes became essential to validate whether the phenotypic change in the mutated plant is due to the deletion of the target genes 1-FEH-w1, -w2, -w3. This could be investigated by using either the bioinformatics approach or the molecular biology approach. With the release of wheat, whole genome sequences and robust gene annotation data, bioinformatics has become more reliable than ever. Moreover, the development of several powerful bioinformatic tools and databases has made this platform reliable to researchers as has been evident in the large number of published articles [14]. On the other hand, employing a molecular biology approach for the functional confirmation or expressional behaviour of large numbers of non-target genes is not feasible. Thus, this study used robust bioinformatic tools to confirm the functional behaviour of the non-target genes.
The IWGSC refseqv1.0 Functional Annotation v1 was the prime database used for the investigation of gene function. Furthermore, function prediction of non-targeted genes was carried out following computational biology approaches. The usual approach uses similarity-based searches within annotation databases. Several tools and databases were used to predict information about gene function using orthologue genes with high sequence similarity and proteins with significant structural homology.
The isoforms of 1-FEH are hydrolase enzymes that catalyse the hydrolysis of chemical bonds. However, different hydrolase enzymes work on different kinds of bonds, meaning they break down different molecules. 1-FEHs break down the β-(2-1) linkage of fructan, which is a storage carbohydrate in cereals. From the annotation, we identified 1-FEH w1, w2 and w3 in the mutated regions of chromosome 6A, 6D and 6B, respectively, whose main functions are breaking down the sugar bond. Within the deletion regions of chromosomes 6A, 6B, and 6D, besides 1-FEH w1, w3, and w2, another 5, 1, and 4 glycoside hydrolase genes were identified. However, their expressional analysis across publicly available databases clearly demonstrated that their expression level was negligible compared to the isoforms of 1-FEH (Supplementary Figures S4–S6). Thus, we considered it unlikely that these genes had any influence on the sugar transport system.

4.5. Influence of 1-FEH Gene Mutations on WSC Remobilisation

This study investigated the performance of 1-FEH isoforms mutation lines in sugar remobilisation under water-deficit and well-watered conditions. Mutation lines demonstrated functional variation of the isoforms in WSC remobilisation and grain filling. WSC content in the stem was measured at different stages during the grain filling, since it is the direct way to monitor differences in the accumulation and mobilisation of stem WSC.

4.6. The 1-FEH Gene Mutations Showed Functional Variation Only under Water-Deficit Condition Not under Well-Watered Condition

The results of this study clearly demonstrated that the functional role of different isoforms of the 1-FEH gene was influenced by water availability. Under well-watered conditions, there was no difference in the WSC concentration, remobilisation rate, and grain weights between the mutation lines compared to the parent Chara under well-watered conditions. However, a significant variation was observed in the WSC remobilisation rate and grain weights under water deficit. Thus, it is evident that the 1-FEH enzyme plays a significant role in stem WSC regulation in water-deficit plants, leading to an effect on grain weight. In contrast, the equal performance of mutation lines and the parent in well-watered plants indicates that the 1-FEH enzyme did not play a significant functional role in stem WSC accumulation or remobilisation under well-watered conditions. A previous study [39] indicated differences in these isoforms of 1-FEH at the transcription level in water-deficit wheat. This study used the 1-FEH mutation lines for the first time, which confirmed this observation at the phenotype level.

4.7. Under Water-Deficit Conditions, WSC Remobilisation Rate and Grain Filling Were Slowed Down by the Deletion of 1-FEH w3

Under water deficit, WSC remobilisation was slower in the 1-FEH w3 mutation line compared to 1-FEH w1 and 1-FEH w2 mutation lines and the parent Chara. In water-deficit plants, the stem WSC peaked at 12 DAA, indicating the starting point of WSC remobilisation. Generally, the stage when the concentration of WSC reaches the highest level is considered as the start of carbohydrate remobilisation, since the process leads to a decline in the WSC concentration [40]. On the other hand, plants reached the permanent wilting point at 22 DAA, which would have reduced metabolic activity due to the scarcity of water. Thus, the period between 12 and 22 DAA in this trial was the most crucial for WSC remobilisation under water deficit. Usually, seed development occurs over 6 weeks after anthesis, and during the first two weeks, only cell division takes place in the endosperm and relatively little dry weight is gained. But from two weeks after anthesis, starch and protein begin to accumulate rapidly in the kernel and the dry weight increases almost in a linear manner. This is the crucial time for accumulating most of the final kernel weight [41]. So, the remobilisation of WSC during the mid-grain filling stage is crucial to maintain grain weight [42].
Noticeably, compared with the parent Chara and the other two mutant lines, at 22 DAA, the stem WSC concentration remained the highest in the 1-FEH w3 mutation line (10% DW) despite the fact that it had the lowest stem WSC content (26.25%) at the peak at 12 DAA. This clearly indicates that the average remobilisation rate of WSC was the slowest in the 1-FEH w3 mutation line. Indeed, WSC remobilisation was 39.8% slower in the 1-FEH w3 mutation line than in Chara. In contrast, mutation lines 1-FEH w1 (18.8% slower) and 1-FEH w2 (13.6% slower) maintained the statistically non-significant differences in remobilisation as the parent Chara. So, it is clearly evident that 1-FEH w3 is playing an important role in the remobilisation of WSC under terminal water-deficit conditions.

4.8. Components of Grain Yield Were Reduced in the 1-FEH w3 Mutation Line but Not in 1-FEH w2 or 1-FEH w1 Mutation Lines

The results indicated that the slow remobilisation rate of WSC in the 1-FEH w3 mutation line directly affected seed development and ultimately resulted in a significant decrease in the final yield. Accordingly, the higher remobilisation in the other two mutation lines (1-FEH w1 and 1-FEH w2) and the parent Chara, which carries the 1-FEH w3 wild type, helped to reduce the grain weight loss from water-deficit stress. This study clearly shows that wheat without 1-FEH w3 would be disadvantaged when grown under post-anthesis water deficit.
It is evident that the three isoforms of 1-FEH do not play equal roles in stem WSC remobilisation under water deficit. The absence of w1 and w2 did not have any significant effect on WSC remobilisation and seed weight development compared to the parent Chara. However, the absence of 1-FEH w3 disrupts sugar transportation to the developing grain in water-deficit plants, resulting in lower seed weight. Chara and the mutation lines 1-FEH w1 and 1-FEH w2 had the 1-FEH w3 isoform, which helped them to minimise the yield loss due to water deficit, which agreed with the previous studies [43,44]. It is worth mentioning that the three isoforms of 1-FEH could have an additive effect in the resulting phenotype. A further investigation by utilising single, double and triple mutant/s of 1-FEH isoforms would be able to provide insight on the additive effect of 1-FEH isoforms.
This study reported the individual functions of the three isoforms of the 1-FEH gene by using gene mutation through heavy ion irradiations. An extensive amount of information from 90K SNP-based genotyping, a wide range of genomic and bioinformatic analyses, and a glasshouse-based drought experiment supported the functional difference of 1-FEH isoforms. Future studies by overexpressing the individual isoforms in wheat or in the model plant species and employing extensive molecular biology approaches will provide a further level of functional confirmation of the isoforms in regard to the WSC remobilisation under terminal drought.

5. Conclusions

We used a gene mutation approach to investigate the individual function of the three isoforms of 1-FEH in wheat. The 90K SNP-based characterisation of gene mutation demonstrated that the deleted regions on all three chromosomes were larger than the targeted genomic regions. However, gene functional prediction analysis confirmed that none of the non-target mutated genes had a significant influence on sugar transportation or remobilisation. An investigation of WSC remobilisation rates, grain-filling rate, and grain weight demonstrated that 1-FEH w3 plays the most crucial role in stem WSC remobilisation under post-anthesis water deficit compared to the other two isoforms.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cimb45080419/s1, Supplementary Tables S1–S3: List of HC genes in deleted region in 6A, 6B, and 6D. Supplementary Figures S1–S3: Examples of SNP detecting a deletion in all homeologues (Screenshot from GenomeStudio). Supplementary Figures S4–S6: Deletion regions on chromosomes 6A, 6B, and 6D. Supplementary Figures S7–S9: Expression comparison of 1-FEH w1, 1-FEH w3 and 1-FEH w2 with the functionally close genes in the mutated region on 6A, 6B and 6D. Supplementary Figure S10: Correlation between seed weight/main stem and WSC remobilisation and WSC concentration.

Author Contributions

Conceptualisation, J.Z., R.A. and B.D.; Methodology, N.K. and J.Z.; Software, S.I.; Validation, J.Z., R.A. and B.D.; Formal Analysis, N.K. and S.I.; Investigation, N.K.; Resources, J.Z. and B.D.; Data Curation, S.I.; Writing—Original Draft Preparation, N.K.; Writing—Review and Editing, J.Z., S.I., R.A. and B.D.; Visualisation, N.K. and S.I.; Supervision, J.Z., R.A. and B.D.; Project Administration, N.K.; Funding Acquisition, N.K. and B.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Grain Research Development Corporation (GRDC), Australia. Grant No.: GRS 10796.

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

Data will be available on request.

Acknowledgments

Angela Juhasz, Senior Research fellow, at Murdoch University, Australia helped with data analysis. State Agricultural Biotechnology Centre (SABC), Murdoch University, Australia, for providing research facilities.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yáñez, A.; Tapia, G.; Guerra, F.; Del Pozo, A. Stem carbohydrate dynamics and expression of genes involved in fructan accumulation and remobili-zation during grain growth in wheat (Triticum aestivum L.) genotypes with contrasting tolerance to water stress. PLoS ONE 2017, 12, e0177667. [Google Scholar] [CrossRef] [Green Version]
  2. Ende, W.V.D.; Clerens, S.; Vergauwen, R.; Van Riet, L.; Van Laere, A.; Yoshida, M.; Kawakami, A. Fructan 1-Exohydrolases. β-(2,1)-Trimmers during Graminan Biosynthesis in Stems of Wheat? Purification, Characterization, Mass Mapping, and Cloning of Two Fructan 1-Exohydrolase Isoforms. Plant Physiol. 2003, 131, 621–631. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Ende, W.V.D.; Michiels, A.; Van Wonterghem, D.; Clerens, S.P.; De Roover, J.; Van Laere, A.J. Defoliation Induces Fructan 1-Exohydrolase II in Witloof Chicory Roots. Cloning and Purification of Two Isoforms, Fructan 1-Exohydrolase IIa and Fructan 1-Exohydrolase IIb. Mass Fingerprint of the Fructan 1-Exohydrolase II Enzymes. Plant Physiol. 2001, 126, 1186–1195. [Google Scholar] [CrossRef] [Green Version]
  4. Doyle, A.; Fischer, R. Dry matter accumulation and water use relationships in wheat crops. Aust. J. Agric. Res. 1979, 30, 815–829. [Google Scholar] [CrossRef]
  5. Van Riet, L.; Nagaraj, V.; Ende, W.V.D.; Clerens, S.; Wiemken, A.; Van Laere, A. Purification, cloning and functional characterization of a fructan 6-exohydrolase from wheat (Triticum aestivum L.). J. Exp. Bot. 2006, 57, 213–223. [Google Scholar] [CrossRef] [PubMed]
  6. Kawakami, A.; Yoshida, M. Fructan: Fructan 1-fructosyltransferase, a key enzyme for biosynthesis of graminan oli-gomers in hardened wheat. Planta 2005, 223, 90–104. [Google Scholar] [CrossRef]
  7. Valluru, R.; Van den Ende, W. Plant fructans in stress environments: Emerging concepts and future prospects. J. Exp. Bot. 2008, 59, 2905–2916. [Google Scholar] [CrossRef] [Green Version]
  8. Van Riet, L.; Altenbach, D.; Vergauwen, R.; Clerens, S.; Kawakami, A.; Yoshida, M.; Ende, W.V.D.; Wiemken, A.; Van Laere, A. Purification, cloning and functional differences of a third fructan 1-exohydrolase (1-FEHw3) from wheat (Triticum aestivum). Physiol. Plant. 2008, 133, 242–253. [Google Scholar] [CrossRef] [PubMed]
  9. Yang, J.; Zhang, J.; Wang, Z.; Zhu, Q.; Liu, L. Activities of fructan- and sucrose-metabolizing enzymes in wheat stems subjected to water stress during grain filling. Planta 2004, 220, 331–343. [Google Scholar] [CrossRef] [PubMed]
  10. Zhang, J.; Dell, B.; Conocono, E.; Waters, I.; Setter, T.; Appels, R. Water deficits in wheat: Fructan exohydrolase (1-FEH) mRNA expression and relationship to soluble carbohydrate concentrations in two varieties. New Phytol. 2009, 181, 843–850. [Google Scholar] [CrossRef] [Green Version]
  11. Fitzgerald, T.L.; Kazan, K.; Li, Z.; Morell, M.K.; Manners, J.M. A high-throughput method for the detection of homoeologous gene deletions in hexaploid wheat. BMC Plant Biol. 2010, 10, 264. [Google Scholar] [CrossRef] [Green Version]
  12. Lloyd, J.; Meinke, D. A Comprehensive Dataset of Genes with a Loss-of-Function Mutant Phenotype in Arabidopsis. Plant Physiol. 2012, 158, 1115–1129. [Google Scholar] [CrossRef] [Green Version]
  13. Knoll, J.E.; Ramos, M.L.; Zeng, Y.; Holbrook, C.C.; Chow, M.; Chen, S.; Maleki, S.; Bhattacharya, A.; Ozias-Akins, P. TILLING for allergen reduction and improvement of quality traits in peanut (Arachis hypogaea L.). BMC Plant Biol. 2011, 11, 81. [Google Scholar] [CrossRef] [Green Version]
  14. The International Wheat Genome Sequencing Consortium (IWGSC); Appels, R.; Eversole, K.; Feuillet, C.; Keller, B.; Rogers, J.; Stein, N.; Pozniak, C.J.; Stein, N.; Choulet, F.; et al. Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science 2018, 361, eaar7191. [Google Scholar] [CrossRef] [Green Version]
  15. Fitzgerald, T.L.; Kazan, K.; Manners, J.M. The Application of Reverse Genetics to Polyploid Plant Species. Crit. Rev. Plant Sci. 2012, 31, 181–200. [Google Scholar] [CrossRef]
  16. Šafář, J.; Šimková, H.; Kubaláková, M.; Číhalíková, J.; Suchánková, P.; Bartoš, J.; Doležel, J. Development of Chromosome-Specific BAC Resources for Genomics of Bread Wheat. Cytogenet. Genome Res. 2010, 129, 211–223. [Google Scholar] [CrossRef] [PubMed]
  17. Wanjugi, H.; Coleman-Derr, D.; Huo, N.; Kianian, S.F.; Luo, M.-C.; Wu, J.; Anderson, O.; Gu, Y.Q. Rapid development of PCR-based genome-specific repetitive DNA junction markers in wheat. Genome 2009, 52, 576–587. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Lai, K.; Duran, C.; Berkman, P.J.; Lorenc, M.T.; Stiller, J.; Manoli, S.; Hayden, M.J.; Forrest, K.L.; Fleury, D.; Baumann, U.; et al. Single nucleotide polymorphism discovery from wheat next-generation sequence data. Plant Biotechnol. J. 2012, 10, 73–749. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Wang, T.L.; Uauy, C.; Robson, F.; Till, B. TILLING in extremis. Plant Biotechnol. J. 2012, 10, 761–772. [Google Scholar] [CrossRef] [PubMed]
  20. Cecchini, E.; Mulligan, B.J.; Covey, S.N.; Milner, J.J. Characterization of gamma irradiation-induced deletion mutations at a selectable locus in Arabidopsis. Mutat. Res. Mol. Mech. Mutagen. 1998, 401, 199–206. [Google Scholar] [CrossRef] [PubMed]
  21. Morita, R.; Kusaba, M.; Iida, S.; Yamaguchi, H.; Nishio, T.; Nishimura, M. Molecular characterization of mutations induced by gamma irradiation in rice. Genes Genet. Syst. 2009, 84, 361–370. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Nordström, K.J.V.; Albani, M.C.; James, G.V.; Gutjahr, C.; Hartwig, B.; Turck, F.; Paszkowski, U.; Coupland, G.; Schneeberger, K. Mutation identification by direct comparison of whole-genome sequencing data from mutant and wild-type individuals using k-mers. Nat. Biotechnol. 2013, 31, 325–330. [Google Scholar] [CrossRef] [Green Version]
  23. Akond, M.; Liu, S.; Schoener, L.; Anderson, J.A.; Kantartzi, S.K.; Meksem, K.; Song, O.; Wang, D.; Wen, Z.; Lightfoot, D.A.; et al. A SNP-based genetic linkage map of soybean using the SoySNP6K Illumina Infinium BeadChip geno-typing array. Plant Genet. Genom. Biotechnol. 2013, 1, 80–89. [Google Scholar] [CrossRef]
  24. Cook, J.P.; McMullen, M.D.; Holland, J.B.; Tian, F.; Bradbury, P.; Ross-Ibarra, J.; Buckler, E.S.; Flint-Garcia, S.A. Genetic Architecture of Maize Kernel Composition in the Nested Association Mapping and Inbred Association Panels. Plant Physiol. 2012, 158, 824–834. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Ganal, M.W.; Durstewitz, G.; Polley, A.; Bérard, A.; Buckler, E.S.; Charcosset, A.; Clarke, J.D.; Graner, E.-M.; Hansen, M.; Joets, J.; et al. A Large Maize (Zea mays L.) SNP Genotyping Array: Development and Germplasm Genotyping, and Genetic Mapping to Compare with the B73 Reference Genome. PLoS ONE 2011, 6, e28334. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Sim, S.-C.; Durstewitz, G.; Plieske, J.; Wieseke, R.; Ganal, M.W.; Van Deynze, A.; Hamilton, J.P.; Buell, C.R.; Causse, M.; Wijeratne, S.; et al. Development of a Large SNP Genotyping Array and Generation of High-Density Genetic Maps in Tomato. PLoS ONE 2012, 7, e40563. [Google Scholar] [CrossRef] [PubMed]
  27. Wiedmann, R.T.; Smith, T.P.; Nonneman, D.J. SNP discovery in swine by reduced representation and high throughput pyrosequencing. BMC Genet. 2008, 9, 81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Zhao, K.; Tung, C.-W.; Eizenga, G.C.; Wright, M.H.; Ali, M.L.; Price, A.H.; Norton, G.J.; Islam, M.R.; Reynolds, A.; Mezey, J.; et al. Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa. Nat. Commun. 2011, 2, 467. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Akbari, M.; Wenzl, P.; Caig, V.; Carling, J.; Xia, L.; Yang, S.; Uszynski, G.; Mohler, V.; Lehmensiek, A.; Kuchel, H.; et al. Diversity arrays technology (DArT) for high-throughput profiling of the hexaploid wheat genome. Theor. Appl. Genet. 2006, 113, 1409–1420. [Google Scholar] [CrossRef]
  30. Trebbi, D.; Maccaferri, M.; de Heer, P.; Sørensen, A.; Giuliani, S.; Salvi, S.; Sanguineti, M.C.; Massi, A.; van der Vossen, E.A.G.; Tuberosa, R. High-throughput SNP discovery and genotyping in durum wheat (Triticum durum Desf.). Theor. Appl. Genet. 2011, 123, 555–569. [Google Scholar] [CrossRef]
  31. Allen, A.M.; Barker, G.L.; Berry, S.T.; Coghill, J.A.; Gwilliam, R.; Kirby, S.; Robinson, P.; Brenchley, R.C.; D’Amore, R.; McKenzie, N.; et al. Transcript-specific, single-nucleotide polymorphism discovery and linkage analysis in hexaploid bread wheat (Triticum aestivum L.). Plant Biotechnol. J. 2011, 9, 1086–1099. [Google Scholar] [CrossRef] [Green Version]
  32. Cavanagh, C.R.; Chao, S.; Wang, S.; Huang, B.E.; Stephen, S.; Kiani, S.; Forrest, K.; Saintenac, C.; Brown-Guedira, G.L.; Akhunova, A.; et al. Genome-wide comparative diversity uncovers multiple targets of selection for improvement in hexaploid wheat landraces and cultivars. Proc. Natl. Acad. Sci. USA 2013, 110, 8057–8062. [Google Scholar] [CrossRef]
  33. Keeble-Gagnère, G.; Pasam, R.; Forrest, K.L.; Wong, D.; Robinson, H.; Godoy, J.; Rattey, A.; Moody, D.; Mullan, D.; Walmsley, T.; et al. Novel Design of Imputation-Enabled SNP Arrays for Breeding and Research Applications Supporting Multi-Species Hybridization. Front. Plant Sci. 2021, 12, 756877. [Google Scholar] [CrossRef]
  34. Akhunov, E.; Nicolet, C.; Dvorak, J. Single nucleotide polymorphism genotyping in polyploid wheat with the Illumina GoldenGate assay. Theor. Appl. Genet. 2009, 119, 507–517. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Wang, S.; Wong, D.; Forrest, K.; Allen, A.; Chao, S.; Huang, B.E.; Maccaferri, M.; Salvi, S.; Milner, S.G.; Cattivelli, L.; et al. Characterization of polyploid wheat genomic diversity using a high-density 90 000 single nucleotide polymorphism array. Plant Biotechnol. J. 2014, 12, 787–796. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Islam, S.; Wang, P.; Vincent, D.; Khan, J.M.; Juhasz, A.; Diepeveen, D.; Lipscombe, R.J.; Ma, W.; Appels, R. Proteomic profiling of developing wheat heads under water-stress. Funct. Integr. Genom. 2020, 20, 695–710. [Google Scholar] [CrossRef] [PubMed]
  37. Fales, F. The assimilation and degradation of carbohydrates by yeast cells. J. Biol. Chem. 1951, 193, 113–124. [Google Scholar] [CrossRef] [PubMed]
  38. Zuryn, S.; Le Gras, S.; Jamet, K.; Jarriault, S. A Strategy for Direct Mapping and Identification of Mutations by Whole-Genome Sequencing. Genetics 2010, 186, 427–430. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Yang, J.; Zhang, J.; Huang, Z.; Zhu, Q.; Wang, L. Remobilization of Carbon Reserves Is Improved by Controlled Soil-Drying during Grain Filling of Wheat. Crop. Sci. 2000, 40, 1645–1655. [Google Scholar] [CrossRef]
  40. Cruz-Aguado, J.A.; Rodés, R.; Pérez, I.P.; Dorado, M. Morphological characteristics and yield components associated with accumulation and loss of dry mass in the internodes of wheat. Field Crop. Res. 2000, 66, 129–139. [Google Scholar] [CrossRef]
  41. Simmons, S.; Oelke, E.; Anderson, P. Growth and Development Guide for Spring Wheat; University of Minnesota Agricultural Extension Service: St. Paul, MN, USA, 1985. [Google Scholar]
  42. Zhang, J.; Dell, B.; Biddulph, B.; Khan, N.; Xu, Y.; Luo, H.; Appels, R. Vernalization gene combination to maximize grain yield in bread wheat (Triticum aestivum L.) in diverse environments. Euphytica 2014, 198, 439–454. [Google Scholar] [CrossRef]
  43. Khoshro, H.H.; Taleei, A.; Bihamta, M.R.; Shahbazi, M.; Abbasi, A.; Ramezanpour, S.S. Expression analysis of the genes involved in accumulation and remobilization of assimilates in wheat stem under terminal drought stress. Plant Growth Regul. 2014, 74, 165–176. [Google Scholar] [CrossRef]
  44. Zhang, J.; Xu, Y.; Chen, W.; Dell, B.; Vergauwen, R.; Biddulph, B.; Khan, N.; Luo, H.; Appels, R.; Van den Ende, W.; et al. A wheat 1-FEH w3 variant underlies enzyme activity for stem WSC remobilization to grain under drought. New Phytol. 2015, 205, 293–305. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Cimb 45 00419 g001
Figure 1. Effect of the water-deficit treatment on physical appearance of wheat plants over time. In each panel, the pot on the left-hand side is a well-watered reference, and the pot on right-hand side is the water-deficit treatment. The number as % under each panel is the soil moisture content of the water-deficit-treated pot. Appearance of wilted leaves and early leaf senescence were noted in the water-deficit-treated plants.
Figure 1. Effect of the water-deficit treatment on physical appearance of wheat plants over time. In each panel, the pot on the left-hand side is a well-watered reference, and the pot on right-hand side is the water-deficit treatment. The number as % under each panel is the soil moisture content of the water-deficit-treated pot. Appearance of wilted leaves and early leaf senescence were noted in the water-deficit-treated plants.
Cimb 45 00419 g001
Cimb 45 00419 g002
Figure 2. Change in stem water-soluble carbohydrate (WSC) concentration (a) and grain dry weight per ear with time after anthesis (b). Closed symbols and bars = well watered, and open symbols and bars = water deficit; from 10 days before anthesis to 32 days after anthesis, Bars with the same letter are not statistically different at p = 0.05 (n = 9), vertical bars indicate the ±SE of the mean of three replicates.
Figure 2. Change in stem water-soluble carbohydrate (WSC) concentration (a) and grain dry weight per ear with time after anthesis (b). Closed symbols and bars = well watered, and open symbols and bars = water deficit; from 10 days before anthesis to 32 days after anthesis, Bars with the same letter are not statistically different at p = 0.05 (n = 9), vertical bars indicate the ±SE of the mean of three replicates.
Cimb 45 00419 g002
Cimb 45 00419 g003
Figure 3. Decline in stem WSC (%DW) in three mutation lines and Chara under water deficit from 12 to 22 days after anthesis. Bars with the same letter are not statistically different at p = 0.05 (n = 9) vertical bars represent ±SE of the mean of three replicates.
Figure 3. Decline in stem WSC (%DW) in three mutation lines and Chara under water deficit from 12 to 22 days after anthesis. Bars with the same letter are not statistically different at p = 0.05 (n = 9) vertical bars represent ±SE of the mean of three replicates.
Cimb 45 00419 g003
Cimb 45 00419 g004
Figure 4. Main stem seed weight (g), seed no/spike and TGW at maturity in well-watered and water-deficit plants in all mutation lines and parent line Chara. Closed symbols = well-watered, and open symbols = water deficit; vertical bars represent ±SE of the mean of three replicates. Values with the same letter are not statistically different at p = 0.05 (n = 9) under water-deficit treatment.
Figure 4. Main stem seed weight (g), seed no/spike and TGW at maturity in well-watered and water-deficit plants in all mutation lines and parent line Chara. Closed symbols = well-watered, and open symbols = water deficit; vertical bars represent ±SE of the mean of three replicates. Values with the same letter are not statistically different at p = 0.05 (n = 9) under water-deficit treatment.
Cimb 45 00419 g004
Table 1. Experimental details.
Table 1. Experimental details.
Linesw1 mutant line -BW22-M3-22-w1
w3 mutant line -BW18b-M3-426-w3
w2 mutant line -BW21-M3-433-w2
Parent: Chara
TreatmentsTerminal water deficit (water deficit imposed from anthesis (50% plants at the flowering stage) until maturity)
Wellwatered
Sample collection growth stagesAR5 (Auricle distance 5 cm)
Fully emerged spike
7 days after anthesis
12 days after anthesis
17 days after anthesis
22 days after anthesis
27 days after anthesis
32 days after anthesis
At grain maturity
Table 2. Summary of deletions observed in each 1-FEH mutant.
Table 2. Summary of deletions observed in each 1-FEH mutant.
6A6A x C6B6B x C6D6D x C
ChrSize of the Deletion (bp)Region Containing FEH-like Sequence22C X22426Cx 426433Cx 433SNP IndexSNP NamePositionDescription of Observed PolyComment
6AS
28,246,961






 DH 76,092wsnp_CAP8_c6680_31368991,28,36,428---1,28,36,621PolySNP near top of del in 22
DH 3722BobWhite_c5092_4221,30,71,495-----1,30,71,594NullNull SNP near top of del in 22
FEH_6ADH
DH 71,378Tdurum_contig42906_7323,26,58,889----3,26,58,989NullExample of SNP located in FEH
DHdhdh55,811RAC875_c25556_12503,26,60,703---3,26,60,803PolyExample of SNP located in FEH
DH 1675BobWhite_c23193_1704,02,21,884---4,02,21,984Poly
DH 41,259Kukri_c14765_16554,10,83,842---4,10,83,933Poly
DH 3378BobWhite_c44549_834,10,83,289---4,10,83,389Poly
6BS
4,876,122


 DH 10,716BS00074183_515,44,19,374---5,44,19,474NullNull SNP near top of del in 426
FEH_6B DH
dhDHdh55,811RAC875_c25556_12505,72,96,208---5,72,86,308PolyExample of SNP located in FEH
DH..75,083tplb0055h14_4835,92,95,396---5,92,95,496PolySNP near bottom of del in 426
6DS
1,600,754


 dh DH18,485D_GBQ4KXB01B5NHZ_3362,86,82,554---2,86,82,785PolySNP near top of del in 433
FEH_6D
dhdhDH55,811RAC875_c25556_12503,02,82,758---3,02,82,858PolyExample of SNP located in FEH
sdhdhDH46,499Kukri_c5531_3583,02,83,208---3,02,83,308PolySNP near end of del in 433. Located in FEH
Note: D (Red): deletion detected on target chromosome; d: deletion detected on homoeologous chromosome. H (Green): deletion detected on target chromosome in heterozygous condition, h: deletion detected on homoeologous chromosome in heterozygous condition. Highlighted (yellow) area showing position for targeted gene. Polymorphisms (poly) assign genomic changes resulting from HIB to being either nulls or a change in sequence.
Table 3. Functional classification of the proteins positioned within the mutated segment.
Table 3. Functional classification of the proteins positioned within the mutated segment.
Major ClassificationChromosome 6AChromosome 6BChromosome 6D
HydrolaseTraesCS6A02G048200
TraesCS6A02G042300
TraesCS6A02G042400
TraesCS6A02G060700
(Fructan 1-exohydrolase w1 (Source: UniProtKB/Swiss Prot; Acc: Q84PN8))
TraesCS6A02G060800
TraesCS6A02G069000		TraesCS6B02G080800
TraesCS6B02G080700
Fructan 1-exohydrolase w3
(Source: UniProtKB/Swiss-Prot; Acc: B6DZC8)		TraesCS6D02G066800
TraesCS6D02G064100
TraesCS6D02G064200
TraesCS6D02G064300

TraesCS6D02G064400
Fructan 1-exohydrolase w2 (Source: UniProtKB/Swiss-Prot; Acc: Q84LA1)
LigaseTraesCS6A02G067500
TraesCS6A02G067600
TraesCS6A02G067700
TraesCS6A02G067800		 TraesCS6D02G065100
TraesCS6D02G065200
OxidoreductaseTraesCS6A02G034100
TraesCS6A02G041900
TraesCS6A02G055700		 TraesCS6D02G072600
TransferaseTraesCS6A02G059300
TraesCS6A02G059700
TraesCS6A02G066000
TraesCS6A02G026800
TraesCS6A02G037400
TraesCS6A02G039500
TraesCS6A02G051200
TraesCS6A02G051300
TraesCS6A02G059700
TraesCS6A02G066000		TraesCS6B02G079600
TraesCS6B02G080100
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Khan, N.; Zhang, J.; Islam, S.; Appels, R.; Dell, B. Wheat Water-Soluble Carbohydrate Remobilisation under Water Deficit by 1-FEH w3. Curr. Issues Mol. Biol. 2023, 45, 6634-6650. https://doi.org/10.3390/cimb45080419

AMA Style

Khan N, Zhang J, Islam S, Appels R, Dell B. Wheat Water-Soluble Carbohydrate Remobilisation under Water Deficit by 1-FEH w3. Current Issues in Molecular Biology. 2023; 45(8):6634-6650. https://doi.org/10.3390/cimb45080419

Chicago/Turabian Style

Khan, Nusrat, Jingjuan Zhang, Shahidul Islam, Rudi Appels, and Bernard Dell. 2023. "Wheat Water-Soluble Carbohydrate Remobilisation under Water Deficit by 1-FEH w3" Current Issues in Molecular Biology 45, no. 8: 6634-6650. https://doi.org/10.3390/cimb45080419

APA Style

Khan, N., Zhang, J., Islam, S., Appels, R., & Dell, B. (2023). Wheat Water-Soluble Carbohydrate Remobilisation under Water Deficit by 1-FEH w3. Current Issues in Molecular Biology, 45(8), 6634-6650. https://doi.org/10.3390/cimb45080419

Article Metrics

Back to TopTop