Peng-Fei Qi

The γ-gliadin multigene family in common wheat (Triticum aestivum) and its closely related species

BackgroundThe unique properties of wheat flour primarily depend on gluten, which is the most important source of protein for human being. γ-Gliadins have been considered to be the most ancient of the wheat gluten family. The complex family structure of γ-gliadins complicates the determination of their function. Moreover, γ-gliadins contain several sets of celiac disease epitopes. However, no systematic research has been conducted yet.ResultsA total of 170 γ-gliadin genes were isolated from common wheat and its closely related species, among which 138 sequences are putatively functional. The ORF lengths of these sequences range from 678 to 1089 bp, and the repetitive region is mainly responsible for the size heterogeneity of γ-gliadins. The repeat motif P(Q/L/S/T/I/V/R/A)F(S/Y/V/Q/I/C/L)P(R/L/S/T/H/C/Y)Q1–2(P(S/L/T/A/F/H)QQ)1–2is repeated from 7 to 22 times. Sequence polymorphism and linkage disequilibrium analyses show that γ-gliadins are highly diverse. Phylogenic analyses indicate that there is no obvious discrimination between Sitopsis and Ae. tauschii at the Gli-1 loci, compared with diploid wheat. According to the number and placement of cysteine residues, we defined nine cysteine patterns and 17 subgroups. Alternatively, we classified γ-gliadins into two types based on the length of repetitive domain. Amino acid composition analyses indicate that there is a wide range of essential amino acids in γ-gliadins, and those γ-gliadins from subgroup SG-10 and SG-12 and γ-gliadins with a short repetitive domain are more nutritional. A screening of toxic epitopes shows that γ-gliadins with a pattern of C9 and γ-gliadins with a short repetitive domain almost lack any epitopes.Conclusionγ-Gliadin sequences in wheat and closely related Aegilops species are diverse. Each group/subgroup contributes differently to nutritional quality and epitope content. It is suggested that the genes with a short repetitive domain are more nutritional and valuable. Therefore, it is possible to breed wheat varieties, the γ-gliadins of which are less, even non-toxic and more nutritional.

Biochemical and molecular characterization of gliadins

Gliadins account for about 40–50% of the total proteins in wheat seeds and play an important role in the nutritional and processing quality of flour. Usually, gliadins can be divided into α-(α/β), γ-, and ω-groups, whereas the low-molecular-weight (LMW) gliadins are novel seed storage proteins. The low-molecular-weight glutenin subunits (LMW-GSs) are also designated as gliadins in a few publications. The genes encoding gliadins are mainly located on the short arms of group 6 and group 1 chromosomes, and not evenly distributed. Repetitive sequences cover most of the uncoding regions, which attributed greatly to the evolution of wheat genome. The primary structure of each gliadin is divided into several domains, and the long repetitive domains consist of peptide motifs. Conserved cysteine residues mainly form intramolecular disulfide bonds. The rare potential intermolecular disulfide bonds and the long repetitive domains play an important role in the quality of wheat flour. There is a general idea that gliadin genes, even prolamin genes, have a common origin and subsequent divergence leads to gene polymorphism. The γ-gliadins are considered to be the most ancient of the wheat prolamin family. Several elements in the 5′-flanking (e.g., CAAT and TATA box) and the 3′-flanking sequences have been detected, which has been shown to be necessary for the proper expression of gliadins.

Molecular characterization of α-gliadin genes from wild emmer wheat (Triticum dicoccoides)

According to the two distal and conserved regions of known α-gliadin genes, gene-specific primers for α-gliadin were designed, and the coding regions of four gliadin genes (i.e. GliTd-1, GliTd-2, GliTd-3 and GliTd-4) with the length of about 800 bp were isolated from the genomic DNA of wild emmer wheat (Triticum dicoccoides). No introns were observed. Sequence comparison indicated that these genes should be classified as α-gliadins. GliTd-3 (GenBank accession No.DQ140351) and GliTd-4 (DQ140352) were potentially functional, whereas GliTd-1 (DQ140349) and GliTd-2 (DQ140350) were both pseudogenes by the definition of in-frame stop codons and frameshifts. Six conserved cysteine residues were observed. Sequence analysis suggested that the motif units of repetitive domain for the four newly detected genes were different from the known genes, and the QQQP sequence before the position 60 was more toxic to coeliac patients. Codons for proline were strongly biased. Codons (CAG and CAA) for glutamine were clustered into the specific regions, and the high percentage of pseudogenes resulted from the mutation of CAG → TAG.

Genome-Wide Quantitative Trait Locus Mapping Identifies Multiple Major Loci for Brittle Rachis and Threshability in Tibetan Semi-Wild Wheat (Triticum aestivum ssp. tibetanum Shao)

Tibetan semi-wild wheat (Triticum aestivum ssp. tibetanum Shao) is a semi-wild hexaploid wheat resource that is only naturally distributed in the Qinghai-Tibet Plateau. Brittle rachis and hard threshing are two important characters of Tibetan semi-wild wheat. A whole-genome linkage map of T. aestivum ssp. tibetanum was constructed using a recombinant inbred line population (Q1028×ZM9023) with 186 lines, 564 diversity array technology markers, and 117 simple sequence repeat markers. Phenotypic data on brittle rachis and threshability, as two quantitative traits, were evaluated on the basis of the number of average spike rachis fragments per spike and percent threshability in 2012 and 2013, respectively. Quantitative trait locus (QTL) mapping performed using inclusive composite interval mapping analysis clearly identified four QTLs for brittle rachis and three QTLs for threshability. However, three loci on 2DS, 2DL, and 5AL showed pleiotropism for brittle rachis and threshability; they respectively explained 5.3%, 18.6%, and 18.6% of phenotypic variation for brittle rachis and 17.4%, 13.2%, and 35.2% of phenotypic variation for threshability. A locus on 3DS showed an independent effect on brittle rachis, which explained 38.7% of the phenotypic variation. The loci on 2DS and 3DS probably represented the effect of Tg and Br1, respectively. The locus on 5AL was in very close proximity to the Q gene, but was different from the predicted q in Tibetan semi-wild wheat. To our knowledge, the locus on 2DL has never been reported in common wheat but was prominent in T. aestivum ssp. tibetanum accession Q1028. It remarkably interacted with the locus on 5AL to affect brittle rachis. Several major loci for brittle rachis and threshability were identified in Tibetan semi-wild wheat, improving the understanding of these two characters and suggesting the occurrence of special evolution in Tibetan semi-wild wheat.

The molecular diversity of α-gliadin genes in the tribe Triticeae

Many of the unique properties of wheat flour are derived from seed storage proteins such as the α-gliadins. In this study these α-gliadin genes from diploid Triticeae species were systemically characterized, and divided into 3 classes according to the distinct organization of their protein domains. Our analyses indicated that these α-gliadins varied in the number of cysteine residues they contained. Most of the α-gliadin genes were grouped according to their genomic origins within the phylogenetic tree. As expected, sequence alignments suggested that the repetitive domain and the two polyglutamine regions were responsible for length variations of α-gliadins as were the insertion/deletion of structural domains within the three different classes (I, II, and III) of α-gliadins. A screening of celiac disease toxic epitopes indicated that the α-gliadins of the class II, derived from the Ns genome, contain no epitope, and that some other genomes contain much fewer epitopes than the A, S(B) and D genomes of wheat. Our results suggest that the observed genetic differences in α-gliadins of Triticeae might indicate their use as a fertile ground for the breeding of less CD-toxic wheat varieties.

Characterization of high-molecular-weight glutenin subunits from Eremopyrum bonaepartis and identification of a novel variant with unusual high molecular weight and altered cysteine residues

We characterized two high-molecular-weight glutenin subunit (HMW-GS) variants from Eremopyrum bonaepartis, determined their complete open reading frames, and further expressed them in a bacterial system. The variants have many novel structural features compared with typical subunits encoded by Glu-1 loci: 1Fx3.7 and 1Fy1.5 exhibit hybrid properties of x- and y-type subunits. In addition, unusual molecular mass and altered number and distribution of cysteine residues were unique features of HMW-GSs encoded by Glu-F1 from E. bonaepartis. The mature 1Fx3.7 subunit has a full length of 1,223 amino acid residues, making it the largest subunit found thus far, while 1Fy1.5 is just 496 residues. In addition, the mutated PGQQ repeat motif was found in the repetitive region of 1Fx3.7. Although it has a similar molecular mass to that previously reported for 1Dx2.2, 1Dx2.2* and 1Sshx2.9 subunits, 1Fx3.7 appears to have had a different evolutionary history. The N-terminal and repetitive regions have a total of four additional cysteine residues, giving 1Fx3.7 a total of eight cysteines, while 1Fy1.5 has only six cysteines because the GHCPTSPQQ nonapeptide at the end of the repetitive region is deleted. With its extra cysteine residues and the longest repetitive region, features that are relevant to good wheat quality, the 1Fx3.7 subunit gene could be an excellent candidate for applications in wheat quality improvement.

Molecular Characterization of the pina Gene in Einkorn Wheat

Fifty-six sequences encoding the pina protein were characterized from three species or subspecies of einkorn wheat. These sequences contained 1,595 nucleotides, including 1,270 conserved sites, 21 single nucleotide polymorphisms (SNPs), and 16 indels. The average frequency of SNPs and indels was one out of 76.1 and 99.9 bases, respectively. Five SNPs and no indels were found in the translated sequences. Fourteen haplotypes were defined, and the accessions in each haplotype ranged from 1 to 18. There were nine haplotypes in Triticum monococcum ssp. aegilopoides, eight in T. monococcum ssp. monococcum, and two in T. urartu. Phylogenetic analysis showed that pina genes from different species or subspecies could be clearly differentiated based on the open reading frame. Genes from T. urartu grouped together, whereas genes from T. monococcum ssp. aegilopoides and T. monococcum ssp. monococcum were shared by three and two clusters, respectively. Both the haplotype and phylogenetic analyses indicated that T. monococcum ssp. aegilopoides was more diverse. These results would contribute to the understanding of functional aspects and efficient utilization of pina genes.

Chitin synthase gene FgCHS8 affects virulence and fungal cell wall sensitivity to environmental stress in Fusarium graminearum

Fusarium graminearum is the major causal agent of Fusarium head blight (FHB) of wheat and barley and is considered to be one of the most devastating plant diseases worldwide. Chitin is a critical component of the fungal cell wall and is polymerized from UDP-N-acetyl-alpha-D-glucosamine by chitin synthase. We characterized FgCHS8, a new class of the chitin synthase gene in F. graminearum. Disruption of FgCHS8 resulted in reduced accumulation of chitin, decreased chitin synthase activity, and had no effect on conidia growth when compared with the wild-type isolate. ΔFgCHS8 had a growth rate comparable to that of the wild-type isolate in vitro. However, ΔFgCHS8 had reduced growth when grown on agar supplemented with either 0.025% SDS or 0.9 mM salicylic acid. ΔFgCHS8 produced significantly less deoxynivalenol and exhibited reduced pathogenicity in wheat spikes. Re-introduction of a functional FgCHS8 gene into the ΔFgCHS8 mutant strain restored the wild-type phenotypes. Fluorescence microscopy revealed that FgCHS8 protein was initially expressed in the septa zone, and then gradually distributed over the entire cellular membrane, indicating that FgCHS8 was required for cell wall development. Our results demonstrated that FgCHS8 is important for cell wall sensitivity to environmental stress factors and deoxynivalenol production in F. graminearum.

Development of chromosome 6D-specific markers for α-gliadin genes and their use in assessing dynamic changes at the Gli-2 loci

To develop chromosome 6D-specific point mutation (PM) markers for α-gliadin genes, 79 α-gliadin sequences cloned from Aegilops tauschii and another 40 α-gliadin genes with known chromosome locations were used in multi-sequence alignment and phylogenic analysis. Additional multiple alignment adjustments were performed manually to facilitate discovery of putative chromosome 6D-specific point mutations. A total of 85 PM primers were designed to detect 68 candidate chromosome 6D-specific point mutations. Experimental tests revealed 11 chromosome 6D-specific PM markers by using genomic DNA from homoeologous group 6 nullisomic–tetrasomic lines of Chinese Spring and putative diploid and tetraploid ancestors of hexaploid wheat as PCR templates. Detection of PM markers in one synthetic hexaploid wheat and its parental lines indicated that some α-gliadin genes were lost from Gli-2 loci during the formation of hexaploid wheat by amphidiploidization of the genomes of Triticum turgidum and Ae. tauschii. Detection of these PM markers in Ae. tauschii, T. aestivum and its four subspecies indicated that at least two genetically distinct sources of Ae. tauschii contributed germplasm to the D genome of T. aestivum.

Uncovering the Dispersion History, Adaptive Evolution and Selection of Wheat in China

Summary Wheat was introduced to China approximately 4500 years ago, where it adapted over a span of time to various environments in agro‐ecological growing zones. We investigated 717 Chinese and 14 Iranian/Turkish geographically diverse, locally adapted wheat landraces with 27 933 DArTseq (for 717 landraces) and 312 831 Wheat660K (for a subset of 285 landraces) markers. This study highlights the adaptive evolutionary history of wheat cultivation in China. Environmental stresses and independent selection efforts have resulted in considerable genome‐wide divergence at the population level in Chinese wheat landraces. In total, 148 regions of the wheat genome show signs of selection in at least one geographic area. Our data show adaptive events across geographic areas, from the xeric northwest to the mesic south, along and among homoeologous chromosomes, with fewer variations in the D genome than in the A and B genomes. Multiple variations in interdependent functional genes such as regulatory and metabolic genes controlling germination and flowering time were characterized, showing clear allelic frequency changes corresponding to the dispersion of wheat in China. Population structure and selection data reveal that Chinese wheat spread from the northwestern Caspian Sea region to South China, adapting during its agricultural trajectory to increasingly mesic and warm climatic areas.