Comprehensive functional analyses of expressed sequence tags in common wheat (Triticum aestivum)

Genome-wide association mapping of genomic regions associated with drought stress tolerance at seedling and reproductive stages in bread wheat

Understanding the genetic architecture of drought stress tolerance in bread wheat at seedling and reproductive stages is crucial for developing drought-tolerant varieties. In the present study, 192 diverse wheat genotypes, a subset from the Wheat Associated Mapping Initiative (WAMI) panel, were evaluated at the seedling stage in a hydroponics system for chlorophyll content (CL), shoot length (SLT), shoot weight (SWT), root length (RLT), and root weight (RWT) under both drought and optimum conditions. Following that, a genome-wide association study (GWAS) was carried out using the phenotypic data recorded during the hydroponics experiment as well as data available from previously conducted multi-location field trials under optimal and drought stress conditions. The panel had previously been genotyped using the Infinium iSelect 90K SNP array with 26,814 polymorphic markers. Using single as well as multi-locus models, GWAS identified 94 significant marker-trait associations (MTAs) or SNPs associated with traits recorded at the seedling stage and 451 for traits recorded at the reproductive stage. The significant SNPs included several novel, significant, and promising MTAs for different traits. The average LD decay distance for the whole genome was approximately 0.48 Mbp, ranging from 0.07 Mbp (chromosome 6D) to 4.14 Mbp (chromosome 2A). Furthermore, several promising SNPs revealed significant differences among haplotypes for traits such as RLT, RWT, SLT, SWT, and GY under drought stress. Functional annotation and in silico expression analysis revealed important putative candidate genes underlying the identified stable genomic regions such as protein kinases, O-methyltransferases, GroES-like superfamily proteins, NAD-dependent dehydratases, etc. The findings of the present study may be useful for improving yield potential, and stability under drought stress conditions.

Morpho-Physiological and Hormonal Response of Winter Wheat Varieties to Drought Stress at Stem Elongation and Anthesis Stages

Drought stress can significantly reduce wheat growth and development as well as grain yield. This study investigated morpho-physiological and hormonal (abscisic (ABA) and salicylic (SA) acids) responses of six winter wheat varieties during stem elongation and anthesis stage as well grain yield-related traits were measured after harvest. To examine drought response, plants were exposed to moderate non-lethal drought stress by withholding watering for 45 and 65% of the volumetric soil moisture content (VSMC) for 14 days at separate experiments for each of those two growth stages. During the stem elongation phase, ABA was increased, confirming the stress status of plants, and SA showed a tendency to increase, suggesting their role as stress hormones in the regulation of stress response, such as the increase in the number of leaves and tillers in drought stress conditions, and further keeping turgor pressure and osmotic adjustment in leaves. At the anthesis stage, heavier drought stress resulted in ABA accumulation in flag leaves that generated an integrated response of maturation, where ABA was not positively correlated with any of investigated traits. After harvest, the variety Bubnjar, followed by Pepeljuga and Anđelka, did not significantly decrease the number of grains per ear and 1000 kernel weight (except Anđelka) in drought treatments, thus, declaring them more tolerant to drought. On the other hand, Rujana, Fifi, and particularly Silvija experienced the highest reduction in grain yield-related traits, considering them drought-sensitive varieties.

Characteristics and Regulating Roles of Wheat TaHsfA2-13 in Abiotic Stresses

Heat shock transcription factor (Hsf) exists widely in eukaryotes and responds to various abiotic stresses by regulating the expression of downstream transcription factors, functional enzymes, and molecular chaperones. In this study, TaHsfA2-13, a heat shock transcription factor belonging to A2 subclass, was cloned from wheat (Triticum aestivum) and its function was analyzed. TaHsfA2-13 encodes a protein containing 368 amino acids and has the basic characteristics of Hsfs. Multiple sequence alignment analysis showed that TaHsfA2-13 protein had the highest similarity with TdHsfA2c-like protein from Triticum dicoccoides, which reached 100%. The analysis of tissue expression characteristics revealed that TaHsfA2-13 was highly expressed in root, shoot, and leaf during the seedling stage of wheat. The expression of TaHsfA2-13 could be upregulated by heat stress, low temperature, H₂O₂, mannitol, salinity and multiple phytohormones. The TaHsfA2-13 protein was located in the nucleus under the normal growth conditions and showed a transcriptional activation activity in yeast. Further studies found that overexpression of TaHsfA2-13 in Arabidopsis thaliana Col-0 or athsfa2 mutant results in improved tolerance to heat stress, H₂O₂, SA and mannitol by regulating the expression of multiple heat shock protein (Hsp) genes. In summary, our study identified TaHsfA2-13 from wheat, revealed its regulatory function in varieties of abiotic stresses, and will provide a new target gene to improve stress tolerance for wheat breeding.

Transcriptome unveiled the gene expression patterns of root architecture in drought-tolerant and sensitive wheat genotypes

Drought is a big challenge for agricultural production. Root attributes are the important target traits for breeding high-yielding sustainable wheat varieties against ever changing climatic conditions. However, the transcriptomic of wheat concerning root architecture remained obscure. Here, we explored RNA-Seq based transcriptome to dissect putative genes involved in root system variations in naturally occurring six genotypes (drought-tolerant and sensitive) of wheat. Global RNA-Seq based root transcriptome analysis revealed single nucleotide polymorphisms (SNPs) variations and differentially expressed genes. Putative 56 SNPs were identified related to 15 genes involved in root architecture. Enrichment of these genes using GO terms demonstrated that differentially expressed genes (DEGs) are divided into sub-categories implicated in molecular functions, cellular components and biological processes. The KEGG analysis of DEGs in each comparison of genotype include metabolic, biosynthesis of secondary metabolites, microbial metabolism in diverse environments and biosynthesis of antibiotics. A deeper insight into DEGs unveiled various pathways involved in drought response and positive gravitropism. These genes belong to various transcription factor families such as DOF, C3H, MYB, and NAC involved in root developmental and stress-related pathways. Local White and UZ-11-CWA-8, which are drought-tolerant genotypes, harbor over-representation of most of DEGs or transcription factors. Notably, a microtubule-associated protein MAPRE1 belonging to RP/EB family recruited in positive gravitropism was enriched. Real-time PCR analysis revealed expression of MAPRE1 and PAL genes is consistent with RNA-seq data. The presented data and genetic resources seem valuable for providing genes involved in the root system architecture of drought-tolerant and susceptible genotypes.

RNA-Seq-based DNA marker analysis of the genetics and molecular evolution of Triticeae species

The release of high-quality chromosome-level genome sequences of members of the Triticeae tribe has greatly facilitated genetic and genomic analyses of important crops such as wheat (Triticum aestivum) and barley (Hordeum vulgare). Due to the large diploid genome size of Triticeae plants (ca. 5 Gbp), transcript analysis is an important method for identifying genetic and genomic differences among Triticeae species. In this review, we summarize our results of RNA-Seq analyses of diploid wheat accessions belonging to the genera Aegilops and Triticum. We also describe studies of the molecular relationships among these accessions and provide insight into the evolution of common hexaploid wheat. DNA markers based on polymorphisms within species can be used to map loci of interest. Even though the genome sequence of diploid Aegilops tauschii, the D-genome donor of common wheat, has been released, the diploid barley genome continues to provide key information about the physical structures of diploid wheat genomes. We describe how a series of RNA-Seq analyses of wheat relatives has helped uncover the structural and evolutionary features of genomic and genetic systems in wild and cultivated Triticeae species.

Development and Characterization of Near-Isogenic Lines Revealing Candidate Genes for a Major 7AL QTL Responsible for Heat Tolerance in Wheat

Wheat is one of the most important food crops in the world, but as a cool-season crop, it is more prone to heat stress, which severely affects crop production and grain quality. Heat tolerance in wheat is a quantitative trait, and the genes underlying reported quantitative trait loci (QTL) have rarely been identified. Near-isogenic lines (NILs) with a common genetic background but differing at a particular locus could turn quantitative traits into a Mendelian factor; therefore, they are suitable material for identifying candidate genes for targeted locus/loci. In this study, we developed and characterized NILs from two populations Cascades × Tevere and Cascades × W156 targeting a major 7AL QTL responsible for heat tolerance. Molecular marker screening and phenotyping for SPAD chlorophyll content and grain-yield-related traits confirmed four pairs of wheat NILs that contrasted for heat-stress responses. Genotyping the NILs using a 90K Infinium iSelect SNP array revealed five single nucleotide polymorphism (SNP) markers within the QTL interval that were distinguishable between the isolines. Seven candidate genes linked to the SNPs were identified as related to heat tolerance, and involved in important processes and pathways in response to heat stress. The confirmed multiple pairs of NILs and identified candidate genes in this study are valuable resources and information for further fine-mapping to clone major genes for heat tolerance.

Development of Drought-Tolerant Transgenic Wheat: Achievements and Limitations

Crop yield improvement is necessary to keep pace with increasing demand for food. Due to climatic variability, the incidence of drought stress at crop growth stages is becoming a major hindering factor to yield improvement. New techniques are required to increase drought tolerance along with improved yield. Genetic modification for increasing drought tolerance is highly desirable, and genetic engineering for drought tolerance requires the expression of certain stress-related genes. Genes have been identified which confer drought tolerance and improve plant growth and survival in transgenic wheat. However, less research has been conducted for the development of transgenic wheat as compared to rice, maize, and other staple food. Furthermore, enhanced tolerance to drought without any yield penalty is a major task of genetic engineering. In this review, we have focused on the progress in the development of transgenic wheat cultivars for improving drought tolerance and discussed the physiological mechanisms and testing of their tolerance in response to inserted genes under control or field conditions.

Comparative analysis of root transcriptome profiles between drought-tolerant and susceptible wheat genotypes in response to water stress

Water deficit is one of the major factors limiting crop productivity worldwide. Plant roots play a key role in uptaking water, perceiving and transducing of water deficit signals to shoot. Although the mechanisms of drought-tolerance have been reported recently, the transcriptional regulatory network of wheat root response to water stress has not been fully understood. In this study, drought-tolerant cultivar JM-262 and susceptible cultivar LM-2 are planted to characterize the root transcriptional changes and physiological responses to water deficit. A total of 8197 drought tolerance-associated differentially expressed genes (DEGs) are identified, these genes are mainly mapped to carbon metabolism, flavonoid biosynthesis, and phytohormone signal transduction. The number and expression level of DEGs involved in antioxidative and antiosmotic stresses are more enhanced in JM-262 under water stress. Furthermore, we find the DEGs related to root development are much more induced in JM-262 in phytohormone signal transduction and carbon metabolism pathway. In conclusion, JM-262 may alleviate the damage of drought by producing more osmoprotectants, ROS scavengers, biomass and energy. Interestingly, hormone signaling and cross-talk probably play an important role in promoting JM-262 greater root systems to take up more water, higher capabilities to induce more drought-related DEGs and higher resisitance to oxidative stresse.

Identification of Wheat Inflorescence Development-Related Genes Using a Comparative Transcriptomics Approach

Inflorescence represents the highly specialized plant tissue producing the grains. Although key genes regulating flower initiation and development are conserved, the mechanism regulating fertility is still not well explained. To identify genes and gene network underlying inflorescence morphology and fertility of bread wheat, expressed sequence tags (ESTs) from different tissues were analyzed using a comparative transcriptomics approach. Based on statistical comparison of EST frequencies of individual genes in EST pools representing different tissues and verification with RT-PCR and RNA-seq data, 170 genes of 59 gene sets predominantly expressed in the inflorescence were obtained. Nearly one-third of the gene sets displayed differentiated expression profiles in terms of their subgenome orthologs. The identified genes, most of which were predominantly expressed in anthers, encode proteins involved in wheat floral identity determination, anther and pollen development, pollen-pistil interaction, and others. Particularly, 25 annotated gene sets are associated with pollen wall formation, of which 18 encode enzymes or proteins participating in lipid metabolic pathway, including fatty acid ω-hydroxylation, alkane and fatty alcohol biosynthesis, and glycerophospholipid metabolism. We showed that the comparative transcriptomics approach was effective in identifying genes for reproductive development and found that lipid metabolism was particularly active in wheat anthers.

Drought Response in Wheat: Key Genes and Regulatory Mechanisms Controlling Root System Architecture and Transpiration Efficiency

Abiotic stresses such as, drought, heat, salinity, and flooding threaten global food security. Crop genetic improvement with increased resilience to abiotic stresses is a critical component of crop breeding strategies. Wheat is an important cereal crop and a staple food source globally. Enhanced drought tolerance in wheat is critical for sustainable food production and global food security. Recent advances in drought tolerance research have uncovered many key genes and transcription regulators governing morpho-physiological traits. Genes controlling root architecture and stomatal development play an important role in soil moisture extraction and its retention, and therefore have been targets of molecular breeding strategies for improving drought tolerance. In this systematic review, we have summarized evidence of beneficial contributions of root and stomatal traits to plant adaptation to drought stress. Specifically, we discuss a few key genes such as, DRO1 in rice and ERECTA in Arabidopsis and rice that were identified to be the enhancers of drought tolerance via regulation of root traits and transpiration efficiency. Additionally, we highlight several transcription factor families, such as, ERF (ethylene response factors), DREB (dehydration responsive element binding), ZFP (zinc finger proteins), WRKY, and MYB that were identified to be both positive and negative regulators of drought responses in wheat, rice, maize, and/or Arabidopsis. The overall aim of this review is to provide an overview of candidate genes that have been identified as regulators of drought response in plants. The lack of a reference genome sequence for wheat and non-transgenic approaches for manipulation of gene functions in wheat in the past had impeded high-resolution interrogation of functional elements, including genes and QTLs, and their application in cultivar improvement. The recent developments in wheat genomics and reverse genetics, including the availability of a gold-standard reference genome sequence and advent of genome editing technologies, are expected to aid in deciphering of the functional roles of genes and regulatory networks underlying adaptive phenological traits, and utilizing the outcomes of such studies in developing drought tolerant cultivars.

De Novo Assembled Wheat Transcriptomes Delineate Differentially Expressed Host Genes in Response to Leaf Rust Infection

Pathogens like Puccinia triticina, the causal organism for leaf rust, extensively damages wheat production. The interaction at molecular level between wheat and the pathogen is complex and less explored. The pathogen induced response was characterized using mock- or pathogen inoculated near-isogenic wheat lines (with or without seedling leaf rust resistance gene Lr28). Four Serial Analysis of Gene Expression libraries were prepared from mock- and pathogen inoculated plants and were subjected to Sequencing by Oligonucleotide Ligation and Detection, which generated a total of 165,767,777 reads, each 35 bases long. The reads were processed and multiple k-mers were attempted for de novo transcript assembly; 22 k-mers showed the best results. Altogether 21,345 contigs were generated and functionally characterized by gene ontology annotation, mining for transcription factors and resistance genes. Expression analysis among the four libraries showed extensive alterations in the transcriptome in response to pathogen infection, reflecting reorganizations in major biological processes and metabolic pathways. Role of auxin in determining pathogenesis in susceptible and resistant lines were imperative. The qPCR expression study of four LRR-RLK (Leucine-rich repeat receptor-like protein kinases) genes showed higher expression at 24 hrs after inoculation with pathogen. In summary, the conceptual model of induced resistance in wheat contributes insights on defense responses and imparts knowledge of Puccinia triticina-induced defense transcripts in wheat plants.

Transcriptome analysis of the Chinese bread wheat cultivar Yunong 201 and its ethyl methanesulfonate mutant line

Roche 454 next-generation sequencing was applied to obtain extensive information about the transcriptomes of the bread wheat cultivar Yunong 201 and its EMS mutant line Yunong 3114. Totals of 1.43 million and 1.44 million raw reads were generated, 14,432, 17,845 and 27,867 isotigs were constructed using the reads in Yunong 201, Yunong 3114 and their combination, respectively. Moreover, 29,042, 34,722, and 48,486 unigenes were generated in Yunong 201, Yunong 3114, and combined cultivars, respectively. A total of 50,382 and 59,891 unigenes from the Yunong 201 and Yunong 3114 were mapped on different chromosomes. Of all unigenes, 1363 DEGs were identified in Yunong 201 and Yunong 3114. qRT-PCR analysis confirmed the expression profiles of 40 candidate unigenes possibly related to abiotic stresses. The expression patterns of four annotated DEGs were also verified in the two wheat cultivars under abiotic stresses. This study provided useful information for further analysis of wheat functional genomics.

Genome-wide fungal stress responsive miRNA expression in wheat

MicroRNAs (miRNAs) are small non-coding class of RNAs. They were identified in many plants with their diverse regulatory roles in several cellular and metabolic processes. A number of miRNAs were involved in biotic and abiotic stress responses. Here, fungal stress responsive wheat miRNAs were analyzed by using miRNA-microarray strategy. Two different fungi (Fusarium culmorum and Bipolaris sorokiniana) were inoculated on resistant and sensitive wheat cultivars. A total of 87 differentially regulated miRNAs were detected in the 8 × 15 K array including all of the available plant miRNAs. Using bioinformatics tools, the target transcripts of responsive miRNAs were predicted, and related biological processes and mechanisms were assessed. A number of the miRNAs such as miR2592s, miR869.1, miR169b were highly differentially regulated showing more than 200-fold change upon fungal-inoculation. Some of the miRNAs were identified as fungal-inoculation responsive for the first time. The analyses showed that some of the differentially regulated miRNAs targeted resistance-related genes such as LRR, glucuronosyl transferase, peroxidase and Pto kinase. The comparison of the two miRNA-microarray analyses indicated that fungal-responsive wheat miRNAs were differentially regulated in pathogen- and cultivar-specific manners.

Transcriptome-wide analysis of WRKY transcription factors in wheat and their leaf rust responsive expression profiling

WRKY, a plant-specific transcription factor family, has important roles in pathogen defense, abiotic cues and phytohormone signaling, yet little is known about their roles and molecular mechanism of function in response to rust diseases in wheat. We identified 100 TaWRKY sequences using wheat Expressed Sequence Tag database of which 22 WRKY sequences were novel. Identified proteins were characterized based on their zinc finger motifs and phylogenetic analysis clustered them into six clades consisting of class IIc and class III WRKY proteins. Functional annotation revealed major functions in metabolic and cellular processes in control plants; whereas response to stimuli, signaling and defense in pathogen inoculated plants, their major molecular function being binding to DNA. Tag-based expression analysis of the identified genes revealed differential expression between mock and Puccinia triticina inoculated wheat near isogenic lines. Gene expression was also performed with six rust-related microarray experiments at Gene Expression Omnibus database. TaWRKY10, 15, 17 and 56 were common in both tag-based and microarray-based differential expression analysis and could be representing rust specific WRKY genes. The obtained results will bestow insight into the functional characterization of WRKY transcription factors responsive to leaf rust pathogenesis that can be used as candidate genes in molecular breeding programs to improve biotic stress tolerance in wheat.

dbWFA: a web-based database for functional annotation of Triticum aestivum transcripts

The functional annotation of genes based on sequence homology with genes from model species genomes is time-consuming because it is necessary to mine several unrelated databases. The aim of the present work was to develop a functional annotation database for common wheat Triticum aestivum (L.). The database, named dbWFA, is based on the reference NCBI UniGene set, an expressed gene catalogue built by expressed sequence tag clustering, and on full-length coding sequences retrieved from the TriFLDB database. Information from good-quality heterogeneous sources, including annotations for model plant species Arabidopsis thaliana (L.) Heynh. and Oryza sativa L., was gathered and linked to T. aestivum sequences through BLAST-based homology searches. Even though the complexity of the transcriptome cannot yet be fully appreciated, we developed a tool to easily and promptly obtain information from multiple functional annotation systems (Gene Ontology, MapMan bin codes, MIPS Functional Categories, PlantCyc pathway reactions and TAIR gene families). The use of dbWFA is illustrated here with several query examples. We were able to assign a putative function to 45% of the UniGenes and 81% of the full-length coding sequences from TriFLDB. Moreover, comparison of the annotation of the whole T. aestivum UniGene set along with curated annotations of the two model species assessed the accuracy of the annotation provided by dbWFA. To further illustrate the use of dbWFA, genes specifically expressed during the early cell division or late storage polymer accumulation phases of T. aestivum grain development were identified using a clustering analysis and then annotated using dbWFA. The annotation of these two sets of genes was consistent with previous analyses of T. aestivum grain transcriptomes and proteomes.

Database URL:urgi.versailles.inra.fr/dbWFA/

A hybrid distance measure for clustering expressed sequence tags originating from the same gene family

BACKGROUND

Clustering is a key step in the processing of Expressed Sequence Tags (ESTs). The primary goal of clustering is to put ESTs from the same transcript of a single gene into a unique cluster. Recent EST clustering algorithms mostly adopt the alignment-free distance measures, where they tend to yield acceptable clustering accuracies with reasonable computational time. Despite the fact that these clustering methods work satisfactorily on a majority of the EST datasets, they have a common weakness. They are prone to deliver unsatisfactory clustering results when dealing with ESTs from the genes derived from the same family. The root cause is the distance measures applied on them are not sensitive enough to separate these closely related genes.

METHODOLOGY/PRINCIPAL FINDINGS

We propose a hybrid distance measure that combines the global and local features extracted from ESTs, with the aim to address the clustering problem faced by ESTs derived from the same gene family. The clustering process is implemented using the DBSCAN algorithm. We test the hybrid distance measure on the ten EST datasets, and the clustering results are compared with the two alignment-free EST clustering tools, i.e. wcd and PEACE. The clustering results indicate that the proposed hybrid distance measure performs relatively better (in terms of clustering accuracy) than both EST clustering tools.

CONCLUSIONS/SIGNIFICANCE

The clustering results provide support for the effectiveness of the proposed hybrid distance measure in solving the clustering problem for ESTs that originate from the same gene family. The improvement of clustering accuracies on the experimental datasets has supported the claim that the sensitivity of the hybrid distance measure is sufficient to solve the clustering problem.

Co-expression and co-responses: within and beyond transcription

Whole genome sequencing, the relative ease of transcript profiling by the use of microarrays and latterly RNA sequencing approaches have facilitated the capture of vast amounts of transcript data. However, despite the enormous progress made in gene annotation a substantial proportion of genes remain to be annotated at the functional level. Considerable progress has, however, been made by searching for transcriptional coordination between genes of known function and non-annotated genes on the premise that such co-expressed genes tend to be functionally related. Here we review progress made following this approach as well as its expansion to include phenotypic information from other levels of cellular organization such as proteomic and metabolomic data as well as physiological and developmental phenotypes.