Qian-Tao Jiang

Genome-wide identification and evaluation of novel internal control genes for Q-PCR based transcript normalization in wheat

To accurately quantify gene expression using quantitative PCR amplification, it is vital that one or more ideal internal control genes are used to normalize the samples to be compared. Ideally, the expression level of those internal control genes should vary as little as possible between tissues, developmental stages and environmental conditions. In this study, 32 candidate genes for internal control were obtained from the analysis of nine independent experiments which included 333 Affymetrix GeneChip Wheat Genome arrays. Expression levels of the selected genes were then evaluated by quantitative real-time PCR with cDNA samples from different tissues, stages of development and environmental conditions. Finally, fifteen novel internal control genes were selected and their respective expression profiles were compared using NormFinder, geNorm, Pearson correlation coefficients and the twofold-change method. The novel internal control genes from this study were compared with thirteen traditional ones for their expression stability. It was observed that seven of the novel internal control genes were better than the traditional ones in expression stability under all the tested cDNA samples. Among the traditional internal control genes, the elongation factor 1-alpha exhibited strong expression stability, whereas the 18S rRNA, Alpha-tubulin, Actin and GAPDH genes had very poor expression stability in the range of wheat samples tested. Therefore, the use of the novel internal control genes for normalization should improve the accuracy and validity of gene expression analysis.

Quantitative Trait Loci Associated with Micronutrient Concentrations in Two Recombinant Inbred Wheat Lines

Micronutrient malnutrition affects over three billion people worldwide, especially women and children in developing countries. Increasing the bioavailable concentrations of essential elements in the edible portions of crops is an effective resolution to address this issue. To determine the genetic factors controlling micronutrient concentration in wheat, the quantitative trait locus (QTL) analysis for iron, zinc, copper, manganese, and selenium concentrations in two recombinant inbred line populations was performed. In all, 39 QTLs for five micronutrient concentrations were identified in this study. Of these, 22 alleles from synthetic wheat SHW-L1 and seven alleles from the progeny line of the synthetic wheat Chuanmai 42 showed an increase in micronutrient concentrations. Five QTLs on chromosomes 2A, 3D, 4D, and 5B found in both the populations showed significant phenotypic variation for 2–3 micronutrient concentrations. Our results might help understand the genetic control of micronutrient concentration and allow the utilization of genetic resources of synthetic hexaploid wheat for improving micronutrient efficiency of cultivated wheat by using molecular marker-assisted selection.

The evolution pattern of rDNA ITS in Avena and phylogenetic relationship of the Avena species (Poaceae: Aveneae).

Ribosomal ITS sequences are commonly used for phylogenetic reconstruction because they are included in rDNA repeats, and these repeats often undergo rapid concerted evolution within and between arrays. Therefore, the rDNA ITS copies appear to be virtually identical and can sometimes be treated as a single gene. In this paper we examined ITS polymorphism within and among 13 diploid (A and C genomes), seven tetraploid (AB, AC and CC genomes) and four hexaploid (ACD genome) to infer the extent and direction of concerted evolution, and to reveal the phylogenetic and genome relationship among species of Avena. A total of 170 clones of the ITS1-5.8S-ITS2 fragment were sequenced to carry out haplotype and phylogenetic analysis. In addition, 111 Avena ITS sequences retrieved from GenBank were combined with 170 clones to construct a phylogeny and a network. We demonstrate the major divergence between the A and C genomes whereas the distinction among the A and B/D genomes was generally not possible. High affinity among the A(d) genome species A. damascena and the ACD genome species A. fatua was found, whereas the rest of the ACD genome hexaploids and the AACC tetraploids were highly affiliated with the A(l) genome diploid A. longiglumis. One of the AACC species A. murphyi showed the closest relationship with most of the hexaploid species. Both C(v) and C(p) genome species have been proposed as paternal donors of the C-genome carrying polyploids. Incomplete concerted evolution is responsible for the observed differences among different clones of a single Avena individual. The elimination of C-genome rRNA sequences and the resulting evolutionary inference of hexaploid species are discussed.

Characterization and expression analysis of waxy alleles in barley accessions

Granule Bound Starch Synthase I (GBSS I) encoded by the waxy gene plays an important role in accumulating amylose during the development of starch granules in barley. In this study, we isolated and characterized waxy alleles of three waxy (GSHO 908, GSHO 1828 and NA 40) and two non-waxy barley accessions (PI 483237 and CIho 15773), estimated the expression patterns of waxy genes via Real-time quantitative PCR (RT-qPCR), investigated promoter activity by analyzing promoter-GUS expression, and examined possible effects of waxy alleles on starch granule morphology in barley accessions by scanning electron microscopy (SEM). A 193-bp insertion in intron 1, a 15-bp insertion in the coding region, and some single nucleotide polymorphic sites were detected in the waxy barley accessions. In addition, a 397-bp deletion containing the TATA box, transcription starting point, exon 1 and partial intron 1 were also identified in the waxy barley accessions. RT-qPCR analysis showed that waxy accessions had lower waxy expression levels than those of non-waxy accessions. Transient expression assays showed that GUS activity driven by the 1,029-bp promoter of the non-waxy accessions was stronger than that driven by the 822-bp promoter of the waxy accessions. SEM revealed no apparent differences of starch granule morphology between waxy and non-waxy accessions. Our results showed that the 397-bp deletion identified in the waxy barley accessions is likely responsible for the reduction of waxy transcript, leading to lower concentrations of GBSS I protein thus lower amylose content.

Novel variants of HMW glutenin subunits from Aegilops section Sitopsis species in relation to evolution and wheat breeding

BackgroundHigh molecular weight glutenin subunits (HMW-GSs), encoded by the genes at Glu-1 loci in wheat and its related species, are significant in the determination of grain processing quality. However, the diversity and variations of HMW-GSs are relatively low in bread wheat. More interests are now focused on wheat wild relatives in Triticeae. The genus Aegilops represents an important germplasm for novel HWM-GSs and other useful genes for wheat genetic improvement.ResultsSix novel Glu-1 alleles and HMW-GSs were identified and characterized from three species of Aegilops section Sitopsis (S genome). Both open reading frames (ORFs) and promoter regions of these Glu-1 alleles were sequenced and characterized. The ORFs of Sitopsis Glu-1 genes are approximately 2.9 kb and 2.3 kb for x-type and y-type subunits, respectively. Although the primary structures of Sitopsis HMW-GSs are similar to those of previously reported ones, all six x-type or y-type subunits have the large fragment insertions. Our comparative analyses of the deduced amino acid sequences verified that Aegilops section Sitopsis species encode novel HMW-GSs with their molecular weights larger than almost all other known HMW-GSs. The Glu-1 promoter sequences share the high homology among S genome. Our phylogenetic analyses by both network and NJ tree indicated that there is a close phylogenetic evolutionary relationship of x-type and y-type subunit between S and D genome.ConclusionsThe large molecular weight of HMW-GSs from S genome is a unique feature identified in this study. Such large subunits are resulted from the duplications of repetitive domains in Sitopsis HMW-GSs. The unequal crossover events are the most likely mechanism of variations in glutenin subunits. The S genome-encoded subunits, 1Dx2.2 and 1Dx2.2* have independent origins, although they share similar evolutionary mechanism. As HMW-GSs play a key role in wheat baking quality, these large Sitopsis glutenin subunits can be used as special genetic resources for wheat quality improvement.

Structure and expression of barley starch phosphorylase genes

The function of starch phosphorylase has long been debated on the regulation of starch metabolism during the growth and development of plants. In this study, we isolated starch phosphorylase genes (Pho1 and Pho2) from barley, characterized their gene and protein structures, predicated their promoter’s cis-elements and analyzed expression patterns. Multiple alignments of these genes showed that (1) both Pho1 and Pho2 genes possess 15 exons and 14 introns in all but three of the species analyzed, Aegilops tauschii (for Pho1 which contains 16 exons and 15 introns), potato (for Pho1b which contains 14 exons and 13 introns), and Triticum uraru (for Pho2 which contains 15 exons and 14 introns); (2) the exon–intron junctions of Pho1 and Pho2 flanking the ligand-binding sites are more conservative than the other regions. Analysis of protein sequences revealed that Pho1 and Pho2 were highly homologous except for two regions, the N terminal domain and the L78 insertion region. The results of real-time quantitative PCR (RT-qPCR) indicated that Pho2 is mainly expressed in germinating seeds, and the expression of Pho1 is similar to that of starch synthesis genes during seed development in barley. Microarray-based analysis indicated that the accumulation of Pho1 or Pho2 transcripts exhibited uniform pattern both in various tissues and various stages of seed development among species of barley, rice, and Arabidopsis. Pho1 of barley was significantly down-regulated under cold and drought treatments, and up-regulated under stem rust infection. Pho2 exhibited similar expression to Pho1 in barley. However, significant difference in expression was not detected for either Pho1 or Pho2 under any of the investigated abiotic stresses. In Arabidopsis, significant down-regulation was detected for Pho1 (PHS1) under abscisic acid (ABA) and for Pho2 (PHS2) under cold, salt, and ABA. Our results provide valuable information to genetically manipulate phosphorylase genes and to further elucidate their regulatory mechanism in the starch biosynthetic pathway.

Characterization of barley Prp1 gene and its expression during seed development and under abiotic stress.

The pre-mRNA processing (Prp1) gene encodes a spliceosomal protein. It was firstly identified in fission yeast and plays a regular role during spliceosome activation and cell cycle. Plant Prp1 genes have only been identified from rice, Sorghum and Arabidopsisthaliana. In this study, we reported the identification and isolation of a novel Prp1 gene from barley, and further explored its expressional pattern by using real-time quantitative RT-PCR, promoter prediction and analysis of microarray data. The putative barley Prp1 protein has a similar primary structure features to those of other known Prp1 protein in this family. The results of amino acid comparison indicated that Prp1 protein of barley and other plant species has a highly conserved 3′ termnal region while their 5′ sequences greatly varied. The results of expressional analysis revealed that the expression level of barley Prp1 gene is always stable in different vegetative tissues, except it is up-regulated at the mid- and late stages of seed development or under the condition of cold stress. This kind of expressional pattern for barley Prp1 is also supported by our results of comparison of microarray data from barley, rice and Arabidopsis. For the molecular mechanism of its expressional pattern, we conclude that the expression of Prp1 gene may be up-regulated by the increase of pre-mRNAs and not be constitutive or ubiquitous.

Characterization of ω-secalin genes from rye, triticale, and a wheat 1BL/1RS translocation line.

Sixty-two DNA sequences for the coding regions of omega-secalin (ω-secalin) genes have been characterized from rye (Secale cereale L.), hexaploid and octoploid triticale (×Triticosecale Wittmack), and wheat (Triticum aestivum L.) 1BL/1RS translocation line. Only 19 out of the 62 ω-secalin gene sequences were full-length open reading frames (ORFs), which can be expressed into functional proteins. The other 43 DNA sequences were pseudogenes, as their ORFs were interrupted by one or a few stop codons or frameshift mutations. The 19 ω-secalin genes have a typical primary structure, which is different from wheat gliadins. There was no cysteine residue in ω-secalin proteins, and the potential celiac disease (CD) toxic epitope (PQQP) was identified to appear frequently in the repetitive domains. The ω-secalin genes from various cereal species shared high homology in their gene sequences. The ω-secalin gene family has involved fewer variations after the integration of the rye R chromosome or whole genome into the wheat or triticale genome. The higher Ka/Ks ratio (i.e. non-synonymous to synonymous substitutions per site) in ω-secalin pseudogenes than in ω-secalin ORFs indicate that the pseudogenes may be subject to a reduced selection pressure. Based on the conserved sequences of ω-secalin genes, it will be possible to manipulate the expression of this gene family in rye, triticale, or wheat 1BL/1RS translocation lines, to reduce its negative effects on grain quality.

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.

Structural variation and evolutionary relationship of novel HMW glutenin subunits from Elymus glaucus.

High molecular weight (HMW) glutenin subunits (GS) are important seed storage proteins relevant to the end-use quality of wheat and other cereal crops. Here we report the isolation and characterization of two novel HMW-GS alleles (1St 1.4 and 1St1.1) from the perennial Triticeae species Elymus glaucus. The amino acid (aa) sequences of E. glaucus 1St1.4 and 1St1.1 were predicted as 434 aa and 358 aa, respectively. Both subunits comprise a signal peptide with a conserved N-terminal domain, a central repetitive domain and a C-terminal domain. Elymus glaucus 1St 1.4 and 1St1.1 exhibit several distinct characteristics different from other known HMW-GSs. The lengths of repetitive domains in E. glaucus 1St 1.4 and 1St1.1 are substantially smaller than those of other known HMW-GSs, in which 1St1.1 (only 358 aa) is the smallest subunit identified so far. The N-terminal domains of E. glaucus 1St 1.4 and 1St1.1 are homologous to y-type subunits, whereas their C-terminal domains are similar to x-type subunits. Our results indicate that E. glaucus 1St 1.4 and 1St1.1 are novel HMW-GS variants or isoforms, and the characterization of both subunits can enhance our understanding on the structural differentiation and evolutionary relationship of HMW-GSs in Triticeae species.