Patterns of homoeologous gene expression shown by RNA sequencing in hexaploid bread wheat

The host-pathogen interaction between wheat and yellow rust induces temporally coordinated waves of gene expression

Background

Understanding how plants and pathogens modulate gene expression during the host-pathogen interaction is key to uncovering the molecular mechanisms that regulate disease progression. Recent advances in sequencing technologies have provided new opportunities to decode the complexity of such interactions. In this study, we used an RNA-based sequencing approach (RNA-seq) to assess the global expression profiles of the wheat yellow rust pathogen Puccinia striiformis f. sp. tritici (PST) and its host during infection.

Results

We performed a detailed RNA-seq time-course for a susceptible and a resistant wheat host infected with PST. This study (i) defined the global gene expression profiles for PST and its wheat host, (ii) substantially improved the gene models for PST, (iii) evaluated the utility of several programmes for quantification of global gene expression for PST and wheat, and (iv) identified clusters of differentially expressed genes in the host and pathogen. By focusing on components of the defence response in susceptible and resistant hosts, we were able to visualise the effect of PST infection on the expression of various defence components and host immune receptors.

Conclusions

Our data showed sequential, temporally coordinated activation and suppression of expression of a suite of immune-response regulators that varied between compatible and incompatible interactions. These findings provide the framework for a better understanding of how PST causes disease and support the idea that PST can suppress the expression of defence components in wheat to successfully colonize a susceptible host.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2684-4) contains supplementary material, which is available to authorized users.

Expression partitioning of homeologs and tandem duplications contribute to salt tolerance in wheat (Triticum aestivum L.)

Salt stress dramatically reduces crop yield and quality, but the molecular mechanisms underlying salt tolerance remain largely unknown. To explore the wheat transcriptional response to salt stress, we performed high-throughput transcriptome sequencing of 10-day old wheat roots under normal condition and 6, 12, 24 and 48 h after salt stress (HASS) in both a salt-tolerant cultivar and salt-sensitive cultivar. The results demonstrated global gene expression reprogramming with 36,804 genes that were up- or down-regulated in wheat roots under at least one stress condition compared with the controls and revealed the specificity and complexity of the functional pathways between the two cultivars. Further analysis showed that substantial expression partitioning of homeologous wheat genes occurs when the plants are subjected to salt stress, accounting for approximately 63.9% (2,537) and 66.1% (2,624) of the homeologous genes in ‘Chinese Spring’ (CS) and ‘Qing Mai 6’ (QM). Interestingly, 143 salt-responsive genes have been duplicated and tandemly arrayed on chromosomes during wheat evolution and polyploidization events, and the expression patterns of 122 (122/143, 85.3%) tandem duplications diverged dynamically over the time-course of salinity exposure. In addition, constitutive expression or silencing of target genes in Arabidopsis and wheat further confirmed our high-confidence salt stress-responsive candidates.

Transcriptome asymmetry in synthetic and natural allotetraploid wheats, revealed by RNA-sequencing

Allopolyploidization has occurred frequently within the Triticum-Aegilops complex which provides a suitable system to investigate how allopolyploidization shapes the expression patterns of duplicated homeologs. We have conducted transcriptome-profiling of leaves and young inflorescences in wild and domesticated tetraploid wheats, Triticum turgidum ssp. dicoccoides (BBAA) and ssp. durum (BBAA), an extracted tetraploid (BBAA), and a synthetic tetraploid (S(l) S(l) AA) wheat together with its diploid parents, Aegilops longissima (S(l) S(l) ) and Triticum urartu (AA). The two diploid species showed tissue-specific differences in genome-wide ortholog expression, which plays an important role in transcriptome shock-mediated homeolog expression rewiring and hence transcriptome asymmetry in the synthetic tetraploid. Further changes of homeolog expression apparently occurred in natural tetraploid wheats, which led to novel transcriptome asymmetry between the two subgenomes. In particular, our results showed that extremely biased homeolog expression can occur rapidly upon the allotetraploidzation and this trend is further enhanced in the course of domestication and evolution of polyploid wheats. Our results suggest that allopolyploidization is accompanied by distinct phases of homeolog expression changes, with parental legacy playing major roles in the immediate rewiring of homeolog expression upon allopolyploidization, while evolution and domestication under allotetraploidy drive further homeolog-expression changes toward re-established subgenome expression asymmetry.

Altered expression of TaRSL4 gene by genome interplay shapes root hair length in allopolyploid wheat

Polyploidy is a major driving force in plant evolution and speciation. Phenotypic changes often arise with the formation, natural selection and domestication of polyploid plants. However, little is known about the consequence of hybridization and polyploidization on root hair development. Here, we report that root hair length of synthetic and natural allopolyploid wheats is significantly longer than those of their diploid progenitors, whereas no difference is observed between allohexaploid and allotetraploid wheats. The expression of wheat gene TaRSL4, an orthologue of AtRSL4 controlling the root hair development in Arabidopsis, was positively correlated with the root hair length in diploid and allotetraploid wheats. Moreover, transcript abundance of TaRSL4 homoeologue from A genome (TaRSL4-A) was much higher than those of other genomes in natural allopolyploid wheat. Notably, increased root hair length by overexpression of the TaRSL4-A in wheat led to enhanced shoot fresh biomass under nutrient-poor conditions. Our observations indicate that increased root hair length in allohexaploid wheat originated in the allotetraploid progenitors and altered expression of TaRSL4 gene by genome interplay shapes root hair length in allopolyploid wheat.

A genome-wide survey of DNA methylation in hexaploid wheat

BACKGROUND

DNA methylation is an important mechanism of epigenetic gene expression control that can be passed between generations. Here, we use sodium bisulfite treatment and targeted gene enrichment to study genome-wide methylation across the three sub-genomes of allohexaploid wheat.

RESULTS

While the majority of methylation is conserved across all three genomes we demonstrate that differential methylation exists between the sub-genomes in approximately equal proportions. We correlate sub-genome-specific promoter methylation with decreased expression levels and show that altered growing temperature has a small effect on methylation state, identifying a small but functionally relevant set of methylated genes. Finally, we demonstrate long-term methylation maintenance using a comparison between the D sub-genome of hexaploid wheat and its progenitor Aegilops tauschii.

CONCLUSIONS

We show that tri-genome methylation is highly conserved with the diploid wheat progenitor while sub-genome-specific methylation shows more variation.

Homoeologous copy-specific expression patterns of MADS-box genes for floral formation in allopolyploid wheat

The consensus model for floral organ formation in higher plants, the so-called ABCDE model, proposes that floral whorl-specific combinations of class A, B, C, D, and E genes specify floral organ identity. Class A, B, C, D and E genes encode MADS-box transcription factors; the single exception being the class A gene APETALA2. Bread wheat (Triticum aestivum) is a hexaploid species with a genome constitution AABBDD; the hexaploid originated from a cross between tetraploid T. turgidum (AABB) and diploid Aegilops tauschii (DD). Tetraploid wheat is thought to have originated from a cross between the diploid species T. urartu (AA) and Ae. speltoides (BB). Consequently, the hexaploid wheat genome contains triplicated homoeologous copies (homoeologs) of each gene derived from the different ancestral diploid species. In this study, we examined the expression patterns of homoeologs of class B, C and D MADS-box genes during floral development. For the class B gene wheat PISTILLATA2 (WPI2), the homoeologs from the A and D genomes were expressed, while expression of the B genome homoeolog was suppressed. For the class C gene wheat AGAMOUS1 (WAG1), the homoeologs on the A and B genomes were expressed, while expression of the D genome homoeolog was suppressed. For the class D gene wheat SEEDSTICK (WSTK), the B genome homoeolog was preferentially expressed. These differential patterns of homoeolog expression were consistently observed among different hexaploid wheat varieties and synthetic hexaploid wheat lines developed by artificial crosses between tetraploid wheat and Ae. tauschii. These results suggest that homoeolog-specific regulation of the floral MADS-box genes occurs in allopolyploid wheat.

Molecular and phylogenetic characterization of the homoeologous EPSP Synthase genes of allohexaploid wheat, Triticum aestivum (L.)

Background

5-Enolpyruvylshikimate-3-phosphate synthase (EPSPS) is the sixth and penultimate enzyme in the shikimate biosynthesis pathway, and is the target of the herbicide glyphosate.  The EPSPS genes of allohexaploid wheat (Triticum aestivum, AABBDD) have not been well characterized. Herein, the three homoeologous copies of the allohexaploid wheat EPSPS gene were cloned and characterized.

Methods

Genomic and coding DNA sequences of EPSPS from the three related genomes of allohexaploid wheat were isolated using PCR and inverse PCR approaches from soft white spring “Louise’.  Development of genome-specific primers allowed the mapping and expression analysis of TaEPSPS-7A1, TaEPSPS-7D1, and TaEPSPS-4A1 on chromosomes 7A, 7D, and 4A, respectively. Sequence alignments of cDNA sequences from wheat and wheat relatives served as a basis for phylogenetic analysis.

Results

The three genomic copies of wheat EPSPS differed by insertion/deletion and single nucleotide polymorphisms (SNPs), largely in intron sequences. RT-PCR analysis and cDNA cloning revealed that EPSPS is expressed from all three genomic copies. However, TaEPSPS-4A1 is expressed at much lower levels than TaEPSPS-7A1 and TaEPSPS-7D1 in wheat seedlings. Phylogenetic analysis of 1190-bp cDNA clones from wheat and wheat relatives revealed that: 1) TaEPSPS-7A1 is most similar to EPSPS from the tetraploid AB genome donor, T. turgidum (99.7 % identity); 2) TaEPSPS-7D1 most resembles EPSPS from the diploid D genome donor, Aegilops tauschii (100 % identity); and 3) TaEPSPS-4A1 resembles EPSPS from the diploid B genome relative, Ae. speltoides (97.7 % identity). Thus, EPSPS sequences in allohexaploid wheat are preserved from the most two recent ancestors. The wheat EPSPS genes are more closely related to Lolium multiflorum and Brachypodium distachyon than to Oryza sativa (rice).

Conclusions

The three related EPSPS homoeologues of wheat exhibited conservation of the exon/intron structure and of coding region sequence, but contained significant sequence variation within intron regions. The genome-specific primers developed will enable future characterization of natural and induced variation in EPSPS sequence and expression.  This can be useful in investigating new causes of glyphosate herbicide resistance.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2084-1) contains supplementary material, which is available to authorized users.

Genome-Wide Comparative Analysis of Flowering-Related Genes in Arabidopsis, Wheat, and Barley

Early flowering is an important trait influencing grain yield and quality in wheat (Triticum aestivum L.) and barley (Hordeum vulgare L.) in short-season cropping regions. However, due to large and complex genomes of these species, direct identification of flowering genes and their molecular characterization remain challenging. Here, we used a bioinformatic approach to predict flowering-related genes in wheat and barley from 190 known Arabidopsis (Arabidopsis thaliana (L.) Heynh.) flowering genes. We identified 900 and 275 putative orthologs in wheat and barley, respectively. The annotated flowering-related genes were clustered into 144 orthologous groups with one-to-one, one-to-many, many-to-one, and many-to-many orthology relationships. Our approach was further validated by domain and phylogenetic analyses of flowering-related proteins and comparative analysis of publicly available microarray data sets for in silico expression profiling of flowering-related genes in 13 different developmental stages of wheat and barley. These further analyses showed that orthologous gene pairs in three critical flowering gene families (PEBP, MADS, and BBX) exhibited similar expression patterns among 13 developmental stages in wheat and barley, suggesting similar functions among the orthologous genes with sequence and expression similarities. The predicted candidate flowering genes can be confirmed and incorporated into molecular breeding for early flowering wheat and barley in short-season cropping regions.

Nonadditive gene expression in polyploids

Allopolyploidy involves hybridization and duplication of divergent parental genomes and provides new avenues for gene expression. The expression levels of duplicated genes in polyploids can show deviation from parental additivity (the arithmetic average of the parental expression levels). Nonadditive expression has been widely observed in diverse polyploids and comprises at least three possible scenarios: (a) The total gene expression level in a polyploid is similar to that of one of its parents (expression-level dominance); (b) total gene expression is lower or higher than in both parents (transgressive expression); and (c) the relative contribution of the parental copies (homeologs) to the total gene expression is unequal (homeolog expression bias). Several factors may result in expression nonadditivity in polyploids, including maternal-paternal influence, gene dosage balance, cis- and/or trans-regulatory networks, and epigenetic regulation. As our understanding of nonadditive gene expression in polyploids remains limited, a new generation of investigators should explore additional phenomena (i.e., alternative splicing) and use other high-throughput "omics" technologies to measure the impact of nonadditive expression on phenotype, proteome, and metabolome.

Temporal transcriptome profiling reveals expression partitioning of homeologous genes contributing to heat and drought acclimation in wheat (Triticum aestivum L.)

BACKGROUND

Hexaploid wheat (Triticum aestivum) is a globally important crop. Heat, drought and their combination dramatically reduce wheat yield and quality, but the molecular mechanisms underlying wheat tolerance to extreme environments, especially stress combination, are largely unknown. As an allohexaploid, wheat consists of three closely related subgenomes (A, B, and D), and was reported to show improved tolerance to stress conditions compared to tetraploid. But so far very little is known about how wheat coordinates the expression of homeologous genes to cope with various environmental constraints on the whole-genome level.

RESULTS

To explore the transcriptional response of wheat to the individual and combined stress, we performed high-throughput transcriptome sequencing of seedlings under normal condition and subjected to drought stress (DS), heat stress (HS) and their combination (HD) for 1 h and 6 h, and presented global gene expression reprograms in response to these three stresses. Gene Ontology (GO) enrichment analysis of DS, HS and HD responsive genes revealed an overlap and complexity of functional pathways between each other. Moreover, 4,375 wheat transcription factors were identified on a whole-genome scale based on the released scaffold information by IWGSC, and 1,328 were responsive to stress treatments. Then, the regulatory network analysis of HSFs and DREBs implicated they were both involved in the regulation of DS, HS and HD response and indicated a cross-talk between heat and drought stress. Finally, approximately 68.4 % of homeologous genes were found to exhibit expression partitioning in response to DS, HS or HD, which was further confirmed by using quantitative RT-PCR and Nullisomic-Tetrasomic lines.

CONCLUSIONS

A large proportion of wheat homeologs exhibited expression partitioning under normal and abiotic stresses, which possibly contributes to the wide adaptability and distribution of hexaploid wheat in response to various environmental constraints.

Evolution of gene expression after gene amplification

We took a rather unique approach to investigate the conservation of gene expression of prolamin storage protein genes across two different subfamilies of the Poaceae. We took advantage of oat plants carrying single maize chromosomes in different cultivars, called oat–maize addition (OMA) lines, which permitted us to determine whether regulation of gene expression was conserved between the two species. We found that γ-zeins are expressed in OMA7.06, which carries maize chromosome 7 even in the absence of the trans-acting maize prolamin-box-binding factor (PBF), which regulates their expression. This is likely because oat PBF can substitute for the function of maize PBF as shown in our transient expression data, using a γ-zein promoter fused to green fluorescent protein (GFP). Despite this conservation, the younger, recently amplified prolamin genes in maize, absent in oat, are not expressed in the corresponding OMAs. However, maize can express the oldest prolamin gene, the wheat high-molecular weight glutenin Dx5 gene, even when maize Pbf is knocked down (through PbfRNAi), and/or another maize transcription factor, Opaque-2 (O2) is knocked out (in maize o2 mutant). Therefore, older genes are conserved in their regulation, whereas younger ones diverged during evolution and eventually acquired a new repertoire of suitable transcriptional activators.

Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement

Upland cotton is a model for polyploid crop domestication and transgenic improvement. Here we sequenced the allotetraploid Gossypium hirsutum L. acc. TM-1 genome by integrating whole-genome shotgun reads, bacterial artificial chromosome (BAC)-end sequences and genotype-by-sequencing genetic maps. We assembled and annotated 32,032 A-subgenome genes and 34,402 D-subgenome genes. Structural rearrangements, gene loss, disrupted genes and sequence divergence were more common in the A subgenome than in the D subgenome, suggesting asymmetric evolution. However, no genome-wide expression dominance was found between the subgenomes. Genomic signatures of selection and domestication are associated with positively selected genes (PSGs) for fiber improvement in the A subgenome and for stress tolerance in the D subgenome. This draft genome sequence provides a resource for engineering superior cotton lines.