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

Upregulation of TaHSP90A transcripts enhances heat tolerance and increases grain yield in wheat under changing climate conditions

Plants have certain adaptation mechanisms to combat temperature extremes and fluctuations. The heat shock protein (HSP90A) plays a crucial role in plant defence mechanisms under heat stress. In silico analysis of the eight TaHSP90A transcripts showed diverse structural patterns in terms of intron/exons, domains, motifs and cis elements in the promoter region in wheat. These regions contained cis elements related to hormones, biotic and abiotic stress and development. To validate these findings, two contrasting wheat genotypes E-01 (thermo-tolerant) and SHP-52 (thermo-sensitive) were used to evaluate the expression pattern of three transcripts TraesCS2A02G033700.1, TraesCS5B02G258900.3 and TraesCS5D02G268000.2 in five different tissues at five different temperature regimes. Expression of TraesCS2A02G033700.1 was upregulated (2-fold) in flag leaf tissue after 1 and 4h of heat treatment in E-01. In contrast, SHP-52 showed downregulated expression after 1h of heat treatment. Additionally, it was shown that under heat stress, the increased expression of TaHSP90A led to an increase in grain production. As the molecular mechanism of genes involved in heat tolerance at the reproductive stage is mostly unknown, these results provide new insights into the role of TaHSP90A transcripts in developing phenotypic plasticity in wheat to develop heat-tolerant cultivars under the current changing climate scenario.

Transgenerational Paternal Inheritance of TaCKX GFMs Expression Patterns Indicate a Way to Select Wheat Lines with Better Parameters for Yield-Related Traits

Members of the TaCKX gene family (GFMs) encode the cytokinin oxygenase/dehydrogenase enzyme (CKX), which irreversibly degrades cytokinins in the organs of wheat plants; therefore, these genes perform a key role in the regulation of yield-related traits. The purpose of the investigation was to determine how expression patterns of these genes, together with the transcription factor-encoding gene TaNAC2-5A, and yield-related traits are inherited to apply this knowledge to speed up breeding processes. The traits were tested in 7 days after pollination (DAP) spikes and seedling roots of maternal and paternal parents and their F₂ progeny. The expression levels of most of them and the yield were inherited in F₂ from the paternal parent. Some pairs or groups of genes cooperated, and some showed opposite functions. Models of up- or down-regulation of TaCKX GFMs and TaNAC2-5A in low-yielding maternal plants crossed with higher-yielding paternal plants and their high-yielding F₂ progeny reproduced gene expression and yield of the paternal parent. The correlation coefficients between TaCKX GFMs, TaNAC2-5A, and yield-related traits in high-yielding F₂ progeny indicated which of these genes were specifically correlated with individual yield-related traits. The most common was expressed in 7 DAP spikes TaCKX2.1, which positively correlated with grain number, grain yield, spike number, and spike length, and seedling root mass. The expression levels of TaCKX1 or TaNAC2-5A in the seedling roots were negatively correlated with these traits. In contrast, the thousand grain weight (TGW) was negatively regulated by TaCKX2.2.2, TaCKX2.1, and TaCKX10 in 7 DAP spikes but positively correlated with TaCKX10 and TaNAC2-5A in seedling roots. Transmission of TaCKX GFMs and TaNAC2-5A expression patterns and yield-related traits from parents to the F₂ generation indicate their paternal imprinting. These newly shown data of nonmendelian epigenetic inheritance shed new light on crossing strategies to obtain a high-yielding F₂ generation.

Identification and Characterization of Phytocyanin Family Genes in Cotton Genomes

Phytocyanins (PCs) are a class of plant-specific blue copper proteins that have been demonstrated to play a role in electron transport and plant development. Through analysis of the copper ligand residues, spectroscopic properties, and domain architecture of the protein, PCs have been grouped into four subfamilies: uclacyanins (UCs), stellacyanins (SCs), plantacyanins (PLCs), and early nodulin-like proteins (ENODLs). The present study aimed to identify and characterise the PCs present in three distinct cotton species (Gossypium hirsutum, Gossyium arboreum, and Gossypium raimondii) through the identification of 98, 63, and 69 genes respectively. We grouped PCs into four clades by using bioinformatics analysis and sequence alignment, which exhibit variations in gene structure and motif distribution. PCs are distributed across all chromosomes in each of the three species, with varying numbers of exons per gene and multiple conserved motifs, and with a minimum of 1 and maximum of 11 exons found on one gene. Transcriptomic data and qRT-PCR analysis revealed that two highly differentiated PC genes were expressed at the fibre initiation stage, while three highly differentiated PCs were expressed at the fibre elongation stage. These findings serve as a foundation for further investigations aimed at understanding the contribution of this gene family in cotton fibre production.

Transcriptional gene silencing in bread wheat (Triticum aestivum L.) and its application to regulate male fertility for hybrid seed production

Transcriptional gene silencing (TGS) can offer a straightforward tool for functional analysis of plant genes, particularly in polyploid species such as wheat, where genetic redundancy poses a challenge in applying mutagenesis approaches, including CRISPR gene editing. In this study, we demonstrate efficient TGS in wheat, mediated by constitutive RNA expression of a promoter inverted repeat (pIR). pIR-mediated TGS of two anther-specific genes, TaMs45 and TaMs1, abolished their function resulting in male sterility. Whilst TGS of TaMs45 required transcriptional silencing of all three homoeologs, a B-genome-specific pIR for TaMs1 was sufficient to confer male sterility. We further show that the pIRs effect TGS of TaMs45 gene through DNA methylation of homologous promoter sequence, successfully suppressing transcription of all three homoeologs. Applying pIR-mediated TGS in wheat, we have generated a dominant male fertility system for production of hybrid seed and demonstrated the efficacy of this system under greenhouse and field conditions. This report describes the first successful TGS in wheat, whilst providing a dominant negative approach as alternative to gene knockout strategies for hybrid wheat breeding and seed production.

Exploiting the miniature inverted-repeat transposable elements insertion polymorphisms as an efficient DNA marker system for genome analysis and evolutionary studies in wheat and related species

Transposable elements (TEs) constitute ~80% of the complex bread wheat genome and contribute significantly to wheat evolution and environmental adaptation. We studied 52 TE insertion polymorphism markers to ascertain their efficiency as a robust DNA marker system for genetic studies in wheat and related species. Significant variation was found in miniature inverted-repeat transposable element (MITE) insertions in relation to ploidy with the highest number of "full site" insertions occurring in the hexaploids (32.6 ± 3.8), while the tetraploid and diploid progenitors had 22.3 ± 0.6 and 15.0 ± 3.5 "full sites," respectively, which suggested a recent rapid activation of these transposons after the formation of wheat. Constructed phylogenetic trees were consistent with the evolutionary history of these species which clustered mainly according to ploidy and genome types (SS, AA, DD, AABB, and AABBDD). The synthetic hexaploids sub-clustered near the tetraploid species from which they were re-synthesized. Preliminary genotyping in 104 recombinant inbred lines (RILs) showed predominantly 1:1 segregation for simplex markers, with four of these markers already integrated into our current DArT-and SNP-based linkage map. The MITE insertions also showed stability with no single excision observed. The MITE insertion site polymorphisms uncovered in this study are very promising as high-potential evolutionary markers for genomic studies in wheat.

N6-Methyladenosine dynamic changes and differential methylation in wheat grain development

More methylation changes occur in late interval than in early interval of wheat seed development with protein and the starch synthesis-related pathway enriched in the later stages. Wheat seed development is a critical process to determining wheat yield and quality, which is controlled by genetics, epigenetics and environments. The N6-methyladenosine (m⁶A) modification is a reversible and dynamic process and plays regulatory role in plant development and stress responses. To better understand the role of m⁶A in wheat grain development, we characterized the m⁶A modification at 10 day post-anthesis (DPA), 20 DPA and 30 DPA in wheat grain development. m⁶A-seq identified 30,615, 30,326, 27,676 high confidence m⁶A peaks from the 10DPA, 20DPA, and 30DPA, respectively, and enriched at 3'UTR. There were 29,964, 29,542 and 26,834 unique peaks identified in AN0942_10d, AN0942_20d and AN0942_30d. One hundred and forty-two genes were methylated by m⁶A throughout seed development, 940 genes methylated in early grain development (AN0942_20d vs AN0942_10d), 1542 genes in late grain development (AN0942_30d vs AN0942_20d), and 1190 genes between early and late development stage (AN0942_30d vs AN0942_10d). KEGG enrichment analysis found that protein-related pathways and the starch synthesis-related pathway were significantly enriched in the later stages of seed development. Our results provide novel knowledge on m⁶A dynamic changes and its roles in wheat grain development.

Wheat-Thinopyrum Substitution Lines Imprint Compensation Both From Recipients and Donors

Even frequently used in wheat breeding, we still have an insufficient understanding of the biology of the products via distant hybridization. In this study, a transcriptomic analysis was performed for six Triticum aestivum-Thinopyrum elongatum substitution lines in comparison with the host plants. All the six disomic substitution lines showed much stronger "transcriptomic-shock" occurred on alien genomes with 57.43-69.22% genes changed expression level but less on the recipient genome (2.19-8.97%). Genome-wide suppression of alien genes along chromosomes was observed with a high proportion of downregulated genes (39.69-48.21%). Oppositely, the wheat recipient showed genome-wide compensation with more upregulated genes, occurring on all chromosomes but not limited to the homeologous groups. Moreover, strong co-upregulation of the orthologs between wheat and Thinopyrum sub-genomes was enriched in photosynthesis with predicted chloroplastic localization, which indicates that the compensation happened not only on wheat host genomes but also on alien genomes.

Meta-analysis reveals consensus genomic regions associated with multiple disease resistance in wheat (Triticum aestivum L.)

In wheat, meta-QTLs (MQTLs) and candidate genes (CGs) were identified for multiple disease resistance (MDR). For this purpose, information was collected from 58 studies for mapping QTLs for resistance to one or more of the five diseases. As many as 493 QTLs were available from these studies, which were distributed in five diseases as follows: septoria tritici blotch (STB) 126 QTLs; septoria nodorum blotch (SNB), 103 QTLs; fusarium head blight (FHB), 184 QTLs; karnal bunt (KB), 66 QTLs; and loose smut (LS), 14 QTLs. Of these 493 QTLs, only 291 QTLs could be projected onto a consensus genetic map, giving 63 MQTLs. The CI of the MQTLs ranged from 0.04 to 15.31 cM with an average of 3.09 cM per MQTL. This is a ~ 4.39 fold reduction from the CI of QTLs, which ranged from 0 to 197.6 cM, with a mean of 13.57 cM. Of 63 MQTLs, 60 were anchored to the reference physical map of wheat (the physical interval of these MQTLs ranged from 0.30 to 726.01 Mb with an average of 74.09 Mb). Thirty-eight (38) of these MQTLs were verified using marker-trait associations (MTAs) derived from genome-wide association studies. As many as 874 CGs were also identified which were further investigated for differential expression using data from five transcriptome studies, resulting in 194 differentially expressed candidate genes (DECGs). Among the DECGs, 85 genes had functions previously reported to be associated with disease resistance. These results should prove useful for fine mapping and cloning of MDR genes and marker-assisted breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01282-z.

Transcriptome Analysis of Developing Wheat Grains at Rapid Expanding Phase Reveals Dynamic Gene Expression Patterns

Simple Summary

Understanding the regulatory mechanism underlying grain development is essential for wheat improvement. The early grain expanding phase boasts critical biological events like embryogenesis and initiation of grain filling. RNA sequencing analysis of this developmental stage revealed dynamic expressions of genes related to cell division, starch biosynthesis, and hormone biosynthesis. An unbalanced expression among triads may play critical roles as shown by multiple enriched metabolic pathways. Our work demonstrated complex regulation mechanisms in early grain development and provided useful information for future wheat improvement.

Abstract

Grain development, as a vital process in the crop’s life cycle, is crucial for determining crop quality and yield. The wheat grain expanding phase is the early process involving the rapid morphological changes and initiation of grain filling. However, little is known about the molecular basis of grain development at this stage. Here, we provide a time-series transcriptome profile of developing wheat grain at 0, 2, 4, 6, 8, and 10 days after pollination of the wheat landrace Chinese Spring. A total of 26,892 differentially expressed genes, including 1468 transcription factors, were found between adjacent time points. Co-expression cluster analysis and Gene Ontology enrichment revealed dynamic expressions of cell division and starch biosynthesis related structural genes and transcription factors. Moreover, diverse, differential and drastically varied expression trends of the key genes related to hormone metabolism were identified. Furthermore, ~30% of triads showed unbalanced expression patterns enriching for genes in multiple pivotal metabolic pathways. Hormone metabolism related genes, such as YUC10 (YUCCA flavin-containing monooxygenase 10), AOS2 (allene oxide synthase 2), CYP90D2 (cytochrome P450 90D2), and CKX1 (cytokinin dehydrogenase 1), were dominantly contributed by A or D homoeologs of the triads. Our study provided a systematic picture of transcriptional regulation of wheat grains at the early grain expanding phase which should deepen our understanding of wheat grain development and help in wheat yield improvement.

Gene expression profiling reveals transcription factor networks and subgenome bias during Brassica napus seed development

We profiled the global gene expression landscape across the reproductive lifecycle of Brassica napus. Comparative analysis of this nascent amphidiploid revealed the contribution of each subgenome to plant reproduction. Whole-genome transcription factor networks identified BZIP11 as a transcriptional regulator of early B. napus seed development. Knockdown of BZIP11 using RNA interference resulted in a similar reduction in gene activity of predicted gene targets, and a reproductive-lethal phenotype. Global mRNA profiling revealed lower accumulation of Cⁿ subgenome transcripts relative to the Aⁿ subgenome. Subgenome-specific transcription factor networks identified distinct transcription factor families enriched in each of the Aⁿ and Cⁿ subgenomes early in seed development. Analysis of laser-microdissected seed subregions further reveal subgenome expression dynamics in the embryo, endosperm and seed coat of early stage seeds. Transcription factors predicted to be regulators encoded by the Aⁿ subgenome are expressed primarily in the seed coat, whereas regulators encoded by the Cⁿ subgenome were expressed primarily in the embryo. Data suggest subgenome bias are characteristic features of the B. napus seed throughout development, and that such bias might not be universal across the embryo, endosperm and seed coat of the developing seed. Transcriptional networks spanning both the Aⁿ and Cⁿ genomes of the whole B. napus seed can identify valuable targets for seed development research and that -omics level approaches to studying gene regulation in B. napus can benefit from both broad and high-resolution analyses.

Genome-wide identification of seed storage protein gene regulators in wheat through coexpression analysis

Only a few transcriptional regulators of seed storage protein (SSP) genes have been identified in common wheat (Triticum aestivum L.). Coexpression analysis could be an efficient approach to characterize novel transcriptional regulators at the genome-scale considering the correlated expression between transcriptional regulators and target genes. As the A genome donor of common wheat, Triticum urartu is more suitable for coexpression analysis than common wheat considering the diploid genome and single gene copy. In this work, the transcriptome dynamics in endosperm of T. urartu throughout grain filling were revealed by RNA-Seq analysis. In the coexpression analysis, a total of 71 transcription factors (TFs) from 23 families were found to be coexpressed with SSP genes. Among these TFs, TuNAC77 enhanced the transcription of SSP genes by binding to cis-elements distributed in promoters. The homolog of TuNAC77 in common wheat, TaNAC77, shared an identical function, and the total SSPs were reduced by about 24% in common wheat when TaNAC77 was knocked down. This is the first genome-wide identification of transcriptional regulators of SSP genes in wheat, and the newly characterized transcriptional regulators will undoubtedly expand our knowledge of the transcriptional regulation of SSP synthesis.

Single-cell-type transcriptomic analysis reveals distinct gene expression profiles in wheat guard cells in response to abscisic acid

Stomatal closure, driven by shrinking guard cells in response to the accumulation of abscisic acid (ABA) under drought stress, has a great impact on plant growth and environmental acclimation. However, the molecular regulatory mechanism underlying the turgor alteration of guard cells remains elusive, especially in cereal grasses. Here, we develop a modified enzyme digestion-based approach for the isolation of wheat (Triticum aestivum L.) guard cells. With this approach, we can remove mesophyll, pavement cells and subsidiary cells successively from the epidermis of the trichomeless coleoptile in wheat and preserve guard cells on the cuticle layers in an intact and physiologically active conditions. Using a robust single-cell-type RNA sequencing analysis, we discovered 9829 differentially expressed genes (DEGs) as significantly up- or down-regulated in guard cells in response to ABA treatment. Transcriptome analysis revealed a large percent of DEGs encoding multiple phytohormone signalling pathways, transporters, calcium signalling components, protein kinases and other ABA signalling-related proteins, which are primarily involved in key signalling pathways in ABA-regulated stomatal control and stress response. Our findings provide valuable resource for investigating the transcriptional regulatory mechanism underlying wheat guard cells in response to ABA.

Comparative transcriptome profiling of a resistant vs susceptible bread wheat (Triticum aestivum L.) cultivar in response to water deficit and cold stress

Bread wheat (Triticum aestivum L.) is one of the most important agricultural plants wearing abiotic stresses, such as water deficit and cold, that cause its productivity reduction. Since resistance to abiotic factors is a multigenic trait, therefore modern genome-wide approaches can help to involve various genetic material in breeding. One technique is full transcriptome analysis that reveals groups of stress response genes serving marker-assisted selection markers. Comparing transcriptome profiles of the same genetic material under several stresses is essential and makes the whole picture. Here, we addressed this by studying the transcriptomic response to water deficit and cold stress for two evolutionarily distant bread wheat varieties: stress-resistant cv. Saratovskaya 29 (S29) and stress-sensitive cv. Yanetzkis Probat (YP). For the first time, transcriptomes for these cultivars grown under abiotic stress conditions were obtained using Illumina based MACE technology. We identified groups of genes involved in response to cold and water deficiency stresses, including responses to each stress factor and both factors simultaneously that may be candidates for resistance genes. We discovered a core group of genes that have a similar pattern of stress-induced expression changes. The particular expression pattern was revealed not only for the studied varieties but also for the published transcriptomic data on cv. Jing 411 and cv. Fielder. Comparative transcriptome profiling of cv. S29 and cv. YP in response to water deficit and cold stress confirmed the hypothesis that stress-induced expression change is unequal within a homeologous gene group. As a rule, at least one changed significantly while the others had a relatively lower expression. Also, we found several SNPs distributed throughout the genomes of cv. S29 and cv. YP and distinguished the studied varieties from each other and the reference cv. Chinese Spring. Our results provide new data for genomics-assisted breeding of stress-tolerant wheat cultivars.

Transcriptomic profiling of wheat near-isogenic lines reveals candidate genes on chromosome 3A for pre-harvest sprouting resistance

BACKGROUND

Pre-harvest sprouting (PHS) in wheat can cause severe damage to both grain yield and quality. Resistance to PHS is a quantitative trait controlled by many genes located across all 21 wheat chromosomes. The study targeted a large-effect quantitative trait locus (QTL) QPhs.ccsu-3A.1 for PHS resistance using several sets previously developed near-isogenic lines (NILs). Two pairs of NILs with highly significant phenotypic differences between the isolines were examined by RNA sequencing for their transcriptomic profiles on developing seeds at 15, 25 and 35 days after pollination (DAP) to identify candidate genes underlying the QTL and elucidate gene effects on PHS resistance. At each DAP, differentially expressed genes (DEGs) between the isolines were investigated.

RESULTS

Gene ontology and KEGG pathway enrichment analyses of key DEGs suggested that six candidate genes underlie QPhs.ccsu-3A.1 responsible for PHS resistance in wheat. Candidate gene expression was further validated by quantitative RT-PCR. Within the targeted QTL interval, 16 genetic variants including five single nucleotide polymorphisms (SNPs) and 11 indels showed consistent polymorphism between resistant and susceptible isolines.

CONCLUSIONS

The targeted QTL is confirmed to harbor core genes related to hormone signaling pathways that can be exploited as a key genomic region for marker-assisted selection. The candidate genes and SNP/indel markers detected in this study are valuable resources for understanding the mechanism of PHS resistance and for marker-assisted breeding of the trait in wheat.

Multi-responses of O-methyltransferase genes to salt stress and fiber development of Gossypium species

BACKGROUND

O-methyltransferases (OMTs) are an important group of enzymes that catalyze the transfer of a methyl group from S-adenosyl-L-methionine to their acceptor substrates. OMTs are divided into several groups according to their structural features. In Gossypium species, they are involved in phenolics and flavonoid pathways. Phenolics defend the cellulose fiber from dreadful external conditions of biotic and abiotic stresses, promoting strength and growth of plant cell wall.

RESULTS

An OMT gene family, containing a total of 192 members, has been identified and characterized in three main Gossypium species, G. hirsutum, G. arboreum and G. raimondii. Cis-regulatory elements analysis suggested important roles of OMT genes in growth, development, and defense against stresses. Transcriptome data of different fiber developmental stages in Chromosome Substitution Segment Lines (CSSLs), Recombination Inbred Lines (RILs) with excellent fiber quality, and standard genetic cotton cultivar TM-1 demonstrate that up-regulation of OMT genes at different fiber developmental stages, and abiotic stress treatments have some significant correlations with fiber quality formation, and with salt stress response. Quantitative RT-PCR results revealed that GhOMT10_Dt and GhOMT70_At genes had a specific expression in response to salt stress while GhOMT49_At, GhOMT49_Dt, and GhOMT48_At in fiber elongation and secondary cell wall stages.

CONCLUSIONS

Our results indicate that O-methyltransferase genes have multi-responses to salt stress and fiber development in Gossypium species and that they may contribute to salt tolerance or fiber quality formation in Gossypium.

Alternative Splicing Diversified the Heat Response and Evolutionary Strategy of Conserved Heat Shock Protein 90s in Hexaploid Wheat (Triticum aestivum L.)

Crops are challenged by the increasing high temperature. Heat shock protein 90 (HSP90), a molecular chaperone, plays a critical role in the heat response in plants. However, the evolutionary conservation and divergence of HSP90s homeologs in polyploidy crops are largely unknown. Using the newly released hexaploid wheat reference sequence, we identified 18 TaHSP90s that are evenly distributed as homeologous genes among three wheat subgenomes, and were highly conserved in terms of sequence identity and gene structure among homeologs. Intensive time-course transcriptomes showed uniform expression and transcriptional response profiles among the three TaHSP90 homeologs. Based on the comprehensive isoforms generated by combining full-length single-molecule sequencing and Illumina short read sequencing, 126 isoforms, including 90 newly identified isoforms of TaHSP90s, were identified, and each TaHSP90 generated one to three major isoforms. Intriguingly, the numbers and the splicing modes of the major isoforms generated by three TaHSP90 homeologs were obviously different. Furthermore, the quantified expression profiles of the major isoforms generated by three TaHSP90 homeologs are also distinctly varied, exhibiting differential alternative splicing (AS) responses of homeologs. Our results showed that the AS diversified the heat response of the conserved TaHSP90s and provided a new perspective for understanding about functional conservation and divergence of homologous genes.

Transcriptome analysis of diploid and triploid Populus tomentosa

Triploid Chinese white poplar (Populus tomentosa Carr., Salicaceae) has stronger advantages in growth and better stress resistance and wood quality than diploid P. tomentosa. Using transcriptome sequencing technology to identify candidate transcriptome-based markers for growth vigor in young tree tissue is of great significance for the breeding of P. tomentosa varieties in the future. In this study, the cuttings of diploid and triploid P. tomentosa were used as plant materials, transcriptome sequencing was carried out, and their tissue culture materials were used for RT-qPCR verification of the expression of genes. The results showed that 12,240 differentially expressed genes in diploid and triploid P. tomentosa transcripts were annotated and enriched into 135 metabolic pathways. The top six pathways that enriched the most significantly different genes were plant-pathogen interaction, phenylpropanoid biosynthesis, MAPK signalling pathway-plant, ascorbate and aldarate metabolism, diterpenoid biosynthesis, and the betalain biosynthesis pathway. Ten growth-related genes were selected from pathways of plant hormone signal transduction and carbon fixation in photosynthetic organisms for RT-qPCR verification. The expression levels of MDH and CYCD3 in tissue-cultured and greenhouse planted triploid P. tomentosa were higher than those in tissue-cultured diploid P. tomentosa, which was consist ent with the TMM values calculated by transcriptome.

Development of RNA-seq-based molecular markers for characterizing Thinopyrum bessarabicum and Secale introgressions in wheat

Sequence-based markers have added a new dimension in the efficiency of identifying alien introgressions in wheat. Expressed sequence tag-sequence tagged sites (EST-STS) markers have proved useful in tracing alien chromatin. In this study, we report the development of Thinopyrum bessarabicum- and Secale anatolicum-specific EST-STS markers and their application in tracing respective alien chromatin introgressions in wheat. The parental lines, Chinese Spring (CS), ISR991.1 (CS/Th. bessarabicum amphidiploid), and ISR1049.2 (CS/Secale anatolicum amphidiploid), were used as core experimental materials. Using comparative analysis of RNA-Seq data, 10 903 and 10 660 candidate sequences specific to Th. bessarabicum and S. anatolicum, respectively, were assembled and identified. To validate the genome specificity of these candidate sequences, 68 and 64 EST-STS markers were developed from randomly selected candidate sequences of Th. bessarabicum and S. anatolicum, respectively, and tested on sets of alien addition lines. Fifty-five and 53 markers for Th. bessarabicum and S. anatolicum chromatin, respectively, were assigned to chromosomal location(s), covering all seven chromosomes. Approximately 83% of S. anatolicum-specific markers were transferable to S. cereale. The genome-specific candidate sequences identified and the EST-STS markers developed will be valuable resources for exploitation of Th. bessarabicum and Secale species diversity in wheat and triticale breeding.

Genome-wide analysis of the HSP101/CLPB gene family for heat tolerance in hexaploid wheat

Heat Shock Protein 101 (HSP101), the homolog of Caseinolytic Protease B (CLPB) proteins, has functional conservation across species to play roles in heat acclimation and plant development. In wheat, several TaHSP101/CLPB genes were identified, but have not been comprehensively characterized. Given the complexity of a polyploid genome with its phenomena of homoeologous expression bias, detailed analysis on the whole TaCLPB family members is important to understand the genetic basis of heat tolerance in hexaploid wheat. In this study, a genome-wide analysis revealed thirteen members of TaCLPB gene family and their expression patterns in various tissues, developmental stages, and stress conditions. Detailed characterization of TaCLPB gene and protein structures suggested potential variations of the sub-cellular localization and their functional regulations. We revealed homoeologous specific variations among TaCLPB gene copies that have not been reported earlier. A study of the Chromosome 1 TaCLPB in four wheat genotypes demonstrated unique patterns of the homoeologous gene expression under moderate and extreme heat treatments. The results give insight into the strategies to improve heat tolerance by targeting one or some of the TaCLPB genes in wheat.

Integrative analysis of hexaploid wheat roots identifies signature components during iron starvation

Iron (Fe) is an essential micronutrient for all organisms. In crop plants, Fe deficiency can decrease crop yield significantly; however, our current understanding of how major crops respond to Fe deficiency remains limited. Herein, the effect of Fe deprivation at both the transcriptomic and metabolic level in hexaploid wheat was investigated. Genome-wide gene expression reprogramming was observed in wheat roots subjected to Fe starvation, with a total of 5854 genes differentially expressed. Homoeologue and subgenome-specific analysis unveiled the induction-biased contribution from the A and B genomes. In general, the predominance of genes coding for nicotianamine synthase, yellow stripe-like transporters, metal transporters, ABC transporters, and zinc-induced facilitator-like protein was noted. Expression of genes related to the Strategy II mode of Fe uptake was also predominant. Our transcriptomic data were in agreement with the GC-MS analysis that showed the enhanced accumulation of various metabolites such as fumarate, malonate, succinate, and xylofuranose, which could be contributing to Fe mobilization. Interestingly, Fe starvation leads to a significant temporal increase of glutathione S-transferase at both the transcriptional level and enzymatic activity level, which indicates the involvement of glutathione in response to Fe stress in wheat roots. Taken together, our result provides new insight into the wheat response to Fe starvation at the molecular level and lays the foundation to design new strategies for the improvement of Fe nutrition in crops.

Breaking Free: The Genomics of Allopolyploidy-Facilitated Niche Expansion in White Clover

The merging of distinct genomes, allopolyploidization, is a widespread phenomenon in plants. It generates adaptive potential through increased genetic diversity, but examples demonstrating its exploitation remain scarce. White clover (Trifolium repens) is a ubiquitous temperate allotetraploid forage crop derived from two European diploid progenitors confined to extreme coastal or alpine habitats. We sequenced and assembled the genomes and transcriptomes of this species complex to gain insight into the genesis of white clover and the consequences of allopolyploidization. Based on these data, we estimate that white clover originated ∼15,000 to 28,000 years ago during the last glaciation when alpine and coastal progenitors were likely colocated in glacial refugia. We found evidence of progenitor diversity carryover through multiple hybridization events and show that the progenitor subgenomes have retained integrity and gene expression activity as they traveled within white clover from their original confined habitats to a global presence. At the transcriptional level, we observed remarkably stable subgenome expression ratios across tissues. Among the few genes that show tissue-specific switching between homeologous gene copies, we found flavonoid biosynthesis genes strongly overrepresented, suggesting an adaptive role of some allopolyploidy-associated transcriptional changes. Our results highlight white clover as an example of allopolyploidy-facilitated niche expansion, where two progenitor genomes, adapted and confined to disparate and highly specialized habitats, expanded to a ubiquitous global presence after glaciation-associated allopolyploidization.

Genome-Wide Identification and Homoeologous Expression Analysis of PP2C Genes in Wheat (Triticum aestivum L.)

Plant protein phosphatase 2Cs (PP2Cs) play crucial roles in phytohormone signaling, developmental processes, and both biotic and abiotic stress responses. However, little research has been conducted on the PP2C gene family in hexaploid wheat (Triticum aestivum L.), which is an important cereal crop. In this study, a genome-wide investigation of TaPP2C gene family was performed. A total of 257 homoeologs of 95 TaPP2C genes were identified, of which 80% of genes had all the three homoeologs across A, B, and D subgenomes. Domain analysis indicated that all the TaPP2C homoeologs harbored the type 2C phosphatase domains. Based on the phylogenetic analysis, TaPP2Cs were divided into 13 groups (A-M) and 4 single branches, which corresponded to the results of gene structure and protein motif analyses. Results of chromosomal location and synteny relationship analysis of TaPP2C homoeologs revealed that known chromosome translocation events and pericentromeric inversions were responsible for the formation of TaPP2C gene family. Expression patterns of TaPP2C homoeologs in various tissues and under diverse stress conditions were analyzed using publicly available RNA-seq data. The results suggested that TaPP2C genes regulate wheat developmental processes and stress responses. Homoeologous expression patterns of TaPP2C triad homoeologs from A, B, and D subgenomes, revealed expression bias within triads under the normal condition, and variability in expression under different stress treatments. Quantitative real-time PCR (qRT-PCR) analysis of eight TaPP2C genes in group A revealed that they were all up-regulated after abscisic acid treatment. Some genes in group A also responded to other phytohormones such as methyl jasmonate and gibberellin. Yeast two-hybrid assays showed that group A TaPP2Cs also interacted with TaSnRK2.1 and TaSnRK2.2 from subclass II, besides with subclass III TaSnRK2s. TaPP2C135 in group A was transformed into Arabidopsis and germination assay revealed that ectopic expression of TaPP2C135 in Arabidopsis enhanced its tolerance to ABA. Overall, these results enhance our understanding of the function of TaPP2Cs in wheat, and provide novel insights into the roles of group A TaPP2Cs. This information will be useful for in-depth functional analysis of TaPP2Cs in future studies and for wheat breeding.

Insights into transcriptional characteristics and homoeolog expression bias of embryo and de-embryonated kernels in developing grain through RNA-Seq and Iso-Seq

Bread wheat (Triticum aestivum L.) is an allohexaploid, and the transcriptional characteristics of the wheat embryo and endosperm during grain development remain unclear. To analyze the transcriptome, we performed isoform sequencing (Iso-Seq) for wheat grain and RNA sequencing (RNA-Seq) for the embryo and de-embryonated kernels. The differential regulation between the embryo and de-embryonated kernels was found to be greater than the difference between the two time points for each tissue. Exactly 2264 and 4790 tissue-specific genes were found at 14 days post-anthesis (DPA), while 5166 and 3784 genes were found at 25 DPA in the embryo and de-embryonated kernels, respectively. Genes expressed in the embryo were more likely to be related to nucleic acid and enzyme regulation. In de-embryonated kernels, genes were rich in substance metabolism and enzyme activity functions. Moreover, 4351, 4641, 4516, and 4453 genes with the A, B, and D homoeoloci were detected for each of the four tissues. Expression characteristics suggested that the D genome may be the largest contributor to the transcriptome in developing grain. Among these, 48, 66, and 38 silenced genes emerged in the A, B, and D genomes, respectively. Gene ontology analysis showed that silenced genes could be inclined to different functions in different genomes. Our study provided specific gene pools of the embryo and de-embryonated kernels and a homoeolog expression bias model on a large scale. This is helpful for providing new insights into the molecular physiology of wheat.

An extracted tetraploid wheat harbouring the BBAA component of common wheat shows anomalous shikimate and sucrose metabolism

BACKGROUND

The BBAA subgenomes of hexaploid common wheat are structurally intact, which makes it possible to extract the BBAA subgenomes to constitute a novel plant type, namely, extracted tetraploid wheat (ETW). ETW displays multiple abnormal phenotypes such as massively reduced biomass and abnormal spike development, compared to extant tetraploid wheat with a BBAA genome. The genetic, biochemical and physiological basis underlying the phenotypic abnormality of ETW remains unknown.

RESULTS

To explore the biochemical basis of these phenotypic abnormalities, we analysed the metabolomic and proteomic profiles and quantified 46 physiological traits of ETW in comparison with its common wheat donor (genome BBAADD), and a durum tetraploid wheat cultivar (genome BBAA). Among these three types of wheat, ETW showed a saliently different pattern of nutrient accumulation and seed quality, markedly lower concentrations of many metabolites involved in carbohydrate metabolism, and higher concentrations of many metabolites related to amino acids. Among the metabolites, changes in shikimate and sucrose were the most conspicuous. Higher levels of shikimate and lower levels of sucrose influence many metabolic processes including carbohydrate and amino acid metabolism, which may contribute to the phenotypic abnormalities. Gene expression assay showed downregulation of a shikimate degradation enzyme (5-enolpyruvylshikimate-3-phosphate synthase) coding gene and upregulation of several genes coding for the sucrose hydrolysis enzyme, which could explain the higher levels of shikimate and lower levels of sucrose, respectively.

CONCLUSIONS

Our results suggest that significant and irreversible biochemical changes have occurred in the BBAA subgenomes of common wheat during the course of its co-evolution with the DD subgenome at the hexaploid level.

Broadening the bread wheat D genome

Although Ae. tauschii has been extensively utilised for wheat breeding, the D-genome-containing allopolyploids have largely remained unexploited. In this review, we discuss approaches that can be used to exploit the D genomes of the different Aegilops species for the improvement of bread wheat. The D genome of allohexaploid bread wheat (Triticum aestivum, 2n = AABBDD) is the least diverse of the three wheat genomes and is unarguably less diverse than that of diploid progenitor Aegilops tauschii (2n = DD). Useful genetic variation and phenotypic traits also exist within each of the wheat group species containing a copy of the D genome: allopolyploid Aegilops species Ae. cylindrica (2n = D^cD^cC^cC^c), Ae. crassa 4x (2n = D¹D¹X^crX^cr), Ae. crassa 6x (2n = D¹D¹X^crX^crD^crD^cr), Ae. ventricosa (2n = D^vD^vN^vN^v), Ae. vavilovii (2n = D¹D¹X^crX^crS^vS^v) and Ae. juvenalis (2n = D¹D¹X^crX^crU^jU^j). Although Ae. tauschii has been extensively utilised for wheat breeding, the D-genome-containing allopolyploids have largely remained unexploited. Some of these D genomes appear to be modified relative to the bread wheat and Ae. tauschii D genomes, and others present in the allopolyploids may also contain useful variation as a result of adaptation to an allopolyploid, multi-genome environment. We summarise the genetic relationships, karyotypic variation and phenotypic traits known to be present in each of the D genome species that could be of relevance for bread wheat improvement and discuss approaches that can be used to exploit the D genomes of the different Aegilops species for the improvement of bread wheat. Better understanding of factors controlling chromosome inheritance and recombination in wheat group interspecific hybrids, as well as effective utilisation of new and developing genetics and genomics technologies, have great potential to improve the agronomic potential of the bread wheat D genome.

Generation of herbicide tolerance traits and a new selectable marker in wheat using base editing

Developing herbicide-tolerant varieties by genome editing holds great promise for addressing the worsening weed problems in wheat cultivation¹. Here, we generated transgene-free wheat germplasms harbouring herbicide tolerance mutations that confer tolerance to sulfonylurea-, imidazolinone- and aryloxyphenoxy propionate-type herbicides by base editing the acetolactate synthase (ALS) and acetyl-coenzyme A carboxylase genes. These stackable herbicide tolerance traits provide a potentially powerful tool for weed management. In addition, we found that base editing at the wheat ALS Pro-174 codon (TaALS-P174) endowed wheat with sufficient resistance to nicosulfuron herbicide in MS growth medium to allow selection. When the TaALS-P174 editor was coupled with editors for other targets of interest, co-editing occurred in the nicosulfuron-resistant plants, and selection for resistance in growth medium enriched the frequency of coupled targets by several-fold. This selectable co-editing system has the potential to greatly bolster adoption of base editing for crop improvement applications.

Genome-wide identification and abiotic stress-responsive pattern of heat shock transcription factor family in Triticum aestivum L

Background

Enhancement of crop productivity under various abiotic stresses is a major objective of agronomic research. Wheat (Triticum aestivum L.) as one of the world’s staple crops is highly sensitive to heat stress, which can adversely affect both yield and quality. Plant heat shock factors (Hsfs) play a crucial role in abiotic and biotic stress response and conferring stress tolerance. Thus, multifunctional Hsfs may be potentially targets in generating novel strains that have the ability to survive environments that feature a combination of stresses.

Result

In this study, using the released genome sequence of wheat and the novel Hsf protein HMM (Hidden Markov Model) model constructed with the Hsf protein sequence of model monocot (Oryza sativa) and dicot (Arabidopsis thaliana) plants, genome-wide TaHsfs identification was performed. Eighty-two non-redundant and full-length TaHsfs were randomly located on 21 chromosomes. The structural characteristics and phylogenetic analysis with Arabidopsis thaliana, Oryza sativa and Zea mays were used to classify these genes into three major classes and further into 13 subclasses. A novel subclass, TaHsfC3 was found which had not been documented in wheat or other plants, and did not show any orthologous genes in A. thaliana, O. sativa, or Z. mays Hsf families. The observation of a high proportion of homeologous TaHsf gene groups suggests that the allopolyploid process, which occurred after the fusion of genomes, contributed to the expansion of the TaHsf family. Furthermore, TaHsfs expression profiling by RNA-seq revealed that the TaHsfs could be responsive not only to abiotic stresses but also to phytohormones. Additionally, the TaHsf family genes exhibited class-, subclass- and organ-specific expression patterns in response to various treatments.

Conclusions

A comprehensive analysis of Hsf genes was performed in wheat, which is useful for better understanding one of the most complex Hsf gene families. Variations in the expression patterns under different abiotic stress and phytohormone treatments provide clues for further analysis of the TaHsfs functions and corresponding signal transduction pathways in wheat.

Electronic supplementary material

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

Analysis of transcriptional responses in root tissue of bread wheat landrace (Triticum aestivum L.) reveals drought avoidance mechanisms under water scarcity

In this study, high-throughput sequencing (RNA-Seq) was utilized to evaluate differential expression of transcripts and their related genes involved in response to terminal drought in root tissues of bread wheat landrace (L-82) and drought-sensitive genotype (Marvdasht). Subsets of 460 differentially expressed genes (DEGs) in drought-tolerant genotype and 236 in drought-sensitive genotype were distinguished and functionally annotated with 105 gene ontology (GO) terms and 77 metabolic pathways. Transcriptome profiling of drought-resistant genotype “L-82” showed up-regulation of genes mostly involved in Oxidation-reduction process, secondary metabolite biosynthesis, abiotic stress response, transferase activity and heat shock proteins. On the other hand, down-regulated genes mostly involved in signaling, oxidation-reduction process, secondary metabolite biosynthesis, auxin-responsive protein and lipid metabolism. We hypothesized that the drought tolerance in “L-82” was a result of avoidance strategies. Up-regulation of genes related to the deeper root system and adequate hydraulic characteristics to allow water uptake under water scarcity confirms our hypothesis. The transcriptomic sequences generated in this study provide information about mechanisms of acclimation to drought in the selected bread wheat landrace, “L-82”, and will help us to unravel the mechanisms underlying the ability of crops to reproduce and keep its productivity even under drought stress.

De novo assembly of wheat root transcriptomes and transcriptional signature of longitudinal differentiation

Hidden underground, root systems constitute an important part of the plant for its development, nourishment and sensing the soil environment around it, but we know very little about its genetic regulation in crop plants like wheat. In the present study, we de novo assembled the root transcriptomes in reference cultivar Chinese Spring from RNA-seq reads generated by the 454-GS-FLX and HiSeq platforms. The FLX reads were assembled into 24,986 transcripts with completeness of 54.84%, and the HiSeq reads were assembled into 91,543 high-confidence protein-coding transcripts, 2,404 low-confidence protein-coding transcripts, and 13,181 non-coding transcripts with the completeness of >90%. Combining the FLX and HiSeq assemblies, we assembled a root transcriptome of 92,335 ORF-containing transcripts. Approximately 7% of the coding transcripts and ~2% non-coding transcripts are not present in the current wheat genome assembly. Functional annotation of both assemblies showed similar gene ontology patterns and that ~7% coding and >5% non-coding transcripts are root-specific. Transcription quantification identified 1,728 differentially expressed transcripts between root tips and maturation zone, and functional annotation of these transcripts captured a transcriptional signature of longitudinal development of wheat root. With the transcriptomic resources developed, this study provided the first view of wheat root transcriptome under different developmental zones and laid a foundation for molecular studies of wheat root development and growth using a reverse genetic approach.

Cloning and characterization of the homoeologous genes for the Rec8-like meiotic cohesin in polyploid wheat

BACKGROUND

Meiosis is a specialized cell division critical for gamete production in the sexual reproduction of eukaryotes. It ensures genome integrity and generates genetic variability as well. The Rec8-like cohesin is a cohesion protein essential for orderly chromosome segregation in meiotic cell division. The Rec8-like genes and cohesins have been cloned and characterized in diploid models, but not in polyploids. The present study aimed to clone the homoeologous genes (homoeoalleles) for Rec8-like cohesin in polyploid wheat, an important food crop for humans, and to characterize their structure and function under a polyploid condition.

RESULTS

We cloned two Rec8-like homoeoalleles from tetraploid wheat (TtRec8-A1 and TtRec8-B1) and one from hexaploid wheat (TaRec8-D1), and performed expression and functional analyses of the homoeoalleles. Also, we identified other two Rec8 homoeoalleles in hexaploid wheat (TaRec8-A1 and TaRec8-B1) and the one in Aegilops tauschii (AetRec8-D1) by referencing the DNA sequences of the Rec8 homoeoalleles cloned in this study. The coding DNA sequences (CDS) of these six Rec8 homoeoalleles are all 1,827 bp in length, encoding 608 amino acids. They differed from each other primarily in introns although single nucleotide polymorphisms were detected in CDS. Substantial difference was observed between the homoeoalleles from the subgenome B (TtRec8-B1 and TaRec8-B1) and those from the subgenomes A and D (TtRec8-A1, TaRec8-A1, and TaRec8-D1). TtRec8-A1 expressed dominantly over TtRec8-B1, but comparably to TaRec8-D1, in polyploid wheat. In addition, we developed the antibody against wheat Rec8 and used the antibody to detect Rec8 cohesin in the Western blotting and subcellular localization analyses.

CONCLUSIONS

The Rec8 homoeoalleles from the subgenomes A and D are transcriptionally more active than the one from the subgenome B in polyploid wheat. The structural variation and differential expression of the Rec8 homoeoalleles indicate a unique cross-genome coordination of the homoeologous genes in polyploid wheat, and imply the distinction of the wheat subgenome B from the subgenomes A and D in the origin and evolution.

Three genomes differentially contribute to the seedling lateral root number in allohexaploid wheat: evidence from phenotype evolution and gene expression

Common wheat is an allohexaploid (BBAADD) that originated from the hybridization and polyploidization of the diploid Aegilops tauschii (DD) with the allotetraploid Triticum turgidum (BBAA). Phenotypic changes often arise with the formation and evolution of allopolyploid wheat, but little is known about the evolution of root traits in different wheat species with varying ploidy levels. Here, we reported that the lateral root number on the primary root (LRNPR) of synthetic and natural allohexaploid wheats (BBAADD) is significantly higher than that of their allotetraploid (BBAA) and diploid (AA and SS) progenitors, but is much lower than that of their diploid (DD) progenitors. The expression of the wheat gene TaLBD16, an ortholog of the Arabidopsis LATERAL ORGAN BOUNDARIES-DOMAIN16/ASYMMETRIC LEAVES2-LIKE18 (LBD16), which is involved in lateral root development in Arabidopsis, was positively correlated with the LRNPR in diploid and allopolyploid wheats. In natural and synthetic allohexaploid wheats, the transcript of the TaLBD16 from the D genome (TaLBD16-D) was relatively more abundant compared with TaLBD16-A and TaLBD16-B. Consistent with the observed variation in LRNPR, the divergence in the expression of TaLBD16 homoeologous genes occurred before the formation of polyploidy wheat. Collectively, our observations indicate that the D genome played a crucial role in the increased lateral root number of allohexaploid wheats compared with their allotetraploid progenitors, and that TaLBD16-D was one of the key genes involved in the formation of lateral root number during wheat evolution.

The transcriptional landscape of polyploid wheat

The coordinated expression of highly related homoeologous genes in polyploid species underlies the phenotypes of many of the world's major crops. Here we combine extensive gene expression datasets to produce a comprehensive, genome-wide analysis of homoeolog expression patterns in hexaploid bread wheat. Bias in homoeolog expression varies between tissues, with ~30% of wheat homoeologs showing nonbalanced expression. We found expression asymmetries along wheat chromosomes, with homoeologs showing the largest inter-tissue, inter-cultivar, and coding sequence variation, most often located in high-recombination distal ends of chromosomes. These transcriptionally dynamic genes potentially represent the first steps toward neo- or subfunctionalization of wheat homoeologs. Coexpression networks reveal extensive coordination of homoeologs throughout development and, alongside a detailed expression atlas, provide a framework to target candidate genes underpinning agronomic traits in wheat.

Phylogenetic analyses and in-seedling expression of ammonium and nitrate transporters in wheat

Plants deploy several ammonium transporter (AMT) and nitrate transporter (NRT) genes to acquire NH₄⁺ and NO₃⁻ from the soil into the roots and then transport them to other plant organs. Coding sequences of wheat genes obtained from ENSEMBL were aligned to known AMT and NRT sequences of Arabidopsis, barley, maize, rice, and wheat to retrieve homologous genes. Bayesian phylogenetic relationships among these genes showed distinct classification of sequences with significant homology to NRT1, NRT2, and NRT3 (NAR2). Inter-species gene duplication analysis showed that eight AMT and 77 NRT genes were orthologous to the AMT and NRT genes of aforementioned plant species. Expression patterns of these genes were studied via whole transcriptome sequencing of 21-day old seedlings of five spring wheat lines. Eight AMT and 52 NRT genes were differentially expressed between root and shoot; and 131 genes did not express neither in root nor in shoot of 21-day old seedlings. Homeologous genes in the A, B, and D genomes, characterized by high sequence homology, revealed that their counterparts exhibited different expression patterns. This complement and evolutionary relationship of wheat AMT and NRT genes is expected to help in development of wheat germplasm with increased efficiency in nitrogen uptake and usage.

Functional Divergence between Subgenomes and Gene Pairs after Whole Genome Duplications

Gene loss following whole genome duplication (WGD) is often biased, with one subgenome retaining more ancestral genes and the other sustaining more gene deletions. While bias toward the greater expression of gene copies on one subgenome can explain bias in gene loss, this raises the question to what drives differences in gene expression levels between subgenomes. Differences in chromatin modifications and epigenetic markers between subgenomes in several model species are now being identified, providing an explanation for bias in gene expression between subgenomes. WGDs can be classified into duplications with higher, biased gene loss and bias in gene expression between subgenomes versus those with lower, unbiased rates of gene loss and an absence of detectable bias between subgenomes; however, the originally proposed link between these two classes and whether WGD results from an allo- or autopolyploid event is inconsistent with recent data from the allopolyploid Capsella bursa-pastoris. The gene balance hypothesis can explain bias in the functional categories of genes retained following WGD, the difference in gene loss rates between unbiased and biased WGDs, and how plant genomes have avoided being overrun with genes encoding dose-sensitive subunits of multiprotein complexes. Comparisons of gene expression patterns between retained transcription factor pairs in maize suggest the high degree of retention for WGD-derived pairs of transcription factors may instead be explained by the older duplication-degeneration-complementation model.

Genomic Imprinting Was Evolutionarily Conserved during Wheat Polyploidization

Genomic imprinting is an epigenetic phenomenon that causes genes to be differentially expressed depending on their parent of origin. To evaluate the evolutionary conservation of genomic imprinting and the effects of ploidy on this process, we investigated parent-of-origin-specific gene expression patterns in the endosperm of diploid (Aegilops spp), tetraploid, and hexaploid wheat (Triticum spp) at various stages of development via high-throughput transcriptome sequencing. We identified 91, 135, and 146 maternally or paternally expressed genes (MEGs or PEGs, respectively) in diploid, tetraploid, and hexaploid wheat, respectively, 52.7% of which exhibited dynamic expression patterns at different developmental stages. Gene Ontology enrichment analysis suggested that MEGs and PEGs were involved in metabolic processes and DNA-dependent transcription, respectively. Nearly half of the imprinted genes exhibited conserved expression patterns during wheat hexaploidization. In addition, 40% of the homoeolog pairs originating from whole-genome duplication were consistently maternally or paternally biased in the different subgenomes of hexaploid wheat. Furthermore, imprinted expression was found for 41.2% and 50.0% of homolog pairs that evolved by tandem duplication after genome duplication in tetraploid and hexaploid wheat, respectively. These results suggest that genomic imprinting was evolutionarily conserved between closely related Triticum and Aegilops species and in the face of polyploid hybridization between species in these genera.

Gene Co-expression Network Analysis Suggests the Existence of Transcriptional Modules Containing a High Proportion of Transcriptionally Differentiated Homoeologs in Hexaploid Wheat

Genome duplications aid in the formation of novel molecular networks through regulatory differentiation of the duplicated genes and facilitate adaptation to environmental change. Hexaploid wheat, Triticum aestivum, contains three homoeologous chromosome sets, the A-, B-, and D-subgenomes, which evolved through interspecific hybridization and subsequent whole-genome duplication. The divergent expression patterns of the homoeologs in hexaploid wheat suggest that they have undergone transcriptional and/or functional differentiation during wheat evolution. However, the distribution of transcriptionally differentiated homoeologs in gene regulatory networks and their related biological functions in hexaploid wheat are still largely unexplored. Therefore, we retrieved 727 publicly available wheat RNA-sequencing (RNA-seq) datasets from various tissues, developmental stages, and conditions, and identified 10,415 expressed homoeologous triplets. Examining the co-expression modules in the wheat transcriptome, we found that 66% of the expressed homoeologous triplets possess all three homoeologs grouped in the same co-expression modules. Among these, 15 triplets contain co-expressed homoeologs with differential expression levels between homoeoalleles across ≥ 95% of the 727 RNA-seq datasets, suggesting a consistent trend of homoeolog expression bias. In addition, we identified 2,831 differentiated homoeologs that showed gene expression patterns that deviated from those of the other two homoeologs. We found that seven co-expression modules contained a high proportion of such differentiated homoeologs, which accounted for ≥ 20% of the genes in each module. We also found that five of the co-expression modules are abundantly composed of genes involved in biological processes such as chloroplast biogenesis, RNA metabolism, putative defense response, putative posttranscriptional modification, and lipid metabolism, thereby suggesting that, the differentiated homoeologs might highly contribute to these biological functions in the gene network of hexaploid wheat.

POTAGE: A Visualisation Tool for Speeding up Gene Discovery in Wheat

POPSEQ Ordered Triticum aestivum Gene Expression (POTAGE) is a web application which accelerates the process of identifying candidate genes for quantitative trait loci (QTL) in hexaploid wheat. This is achieved by leveraging several of the most commonly used data sets in wheat research. These include the Chromosome Survey Sequences, their order along the chromosomes determined by the population sequencing (POPSEQ) approach, the gene predictions and RNA-Seq expression data. POTAGE aggregates those data sets and provides an intuitive interface for biologists to explore the expression of the predicted genes and their functional annotation in a chromosomal context. The interface accelerates some of the laborious and repetitive tasks commonly undertaken in the process of identifying and prioritising genes which may underlie QTL. We illustrate the utility of POTAGE by showing how a short-list of candidate genes can quickly be identified for a QTL linked to pre-harvest sprouting - a major cause of quality and yield loss in wheat production. The candidate genes identified using POTAGE included TaMKK3, which was recently reported as a causal gene for seed dormancy in wheat, and a mutation in its barley ortholog has been shown to reduce pre-harvest sprouting. POTAGE is available at http://crobiad.agwine.adelaide.edu.au/potage .

Conclusive evidence for hexasomic inheritance in chrysanthemum based on analysis of a 183 k SNP array

BACKGROUND

Cultivated chrysanthemum is an outcrossing hexaploid (2n = 6× = 54) with a disputed mode of inheritance. In this paper, we present a single nucleotide polymorphism (SNP) selection pipeline that was used to design an Affymetrix Axiom array with 183 k SNPs from RNA sequencing data (1). With this array, we genotyped four bi-parental populations (with sizes of 405, 53, 76 and 37 offspring plants respectively), and a cultivar panel of 63 genotypes. Further, we present a method for dosage scoring in hexaploids from signal intensities of the array based on mixture models (2) and validation of selection steps in the SNP selection pipeline (3). The resulting genotypic data is used to draw conclusions on the mode of inheritance in chrysanthemum (4), and to make an inference on allelic expression bias (5).

RESULTS

With use of the mixture model approach, we successfully called the dosage of 73,936 out of 183,130 SNPs (40.4%) that segregated in any of the bi-parental populations. To investigate the mode of inheritance, we analysed markers that segregated in the large bi-parental population (n = 405). Analysis of segregation of duplex x nulliplex SNPs resulted in evidence for genome-wide hexasomic inheritance. This evidence was substantiated by the absence of strong linkage between markers in repulsion, which indicated absence of full disomic inheritance. We present the success rate of SNP discovery out of RNA sequencing data as affected by different selection steps, among which SNP coverage over genotypes and use of different types of sequence read mapping software. Genomic dosage highly correlated with relative allele coverage from the RNA sequencing data, indicating that most alleles are expressed according to their genomic dosage.

CONCLUSIONS

The large population, genotyped with a very large number of markers, is a unique framework for extensive genetic analyses in hexaploid chrysanthemum. As starting point, we show conclusive evidence for genome-wide hexasomic inheritance.

Revisiting Pivotal-Differential Genome Evolution in Wheat

An interesting and possibly unique pattern of genome evolution following polyploidy can be observed among allopolyploids of the Triticum and Aegilops genera (wheat group). Most polyploids in this group are presumed to share a common unaltered (pivotal) subgenome (U, D, or A) together with one or two modified (differential) subgenomes, a status that has been referred to as 'pivotal-differential' genome evolution. In this review we discuss various mechanisms that could be responsible for this evolutionary pattern, as well as evidence for and against the putative evolutionary mechanisms involved. We suggest that, in light of recent advances in genome sequencing and related technologies in the wheat group, the time has come to reopen the investigation into pivotal-differential genome evolution.

Spatiotemporal expression patterns of wheat amino acid transporters reveal their putative roles in nitrogen transport and responses to abiotic stress

Amino acid transporters have roles in amino acid uptake from soil, long-distance transport, remobilization from vegetative tissues and accumulation in grain. Critically, the majority of wheat grain nitrogen is derived from amino acids remobilized from vegetative organs. However, no systematic analysis of wheat AAT genes has been reported to date. Here, 283 full length wheat AAT genes representing 100 distinct groups of homeologs were identified and curated by selectively consolidating IWGSC CSSv2 and TGACv1 Triticum aestivum genome assemblies and reassembling or mapping of IWGSC CSS chromosome sorted reads to fill any gaps. Gene expression profiling was performed using public RNA-seq data from root, leaf, stem, spike, grain and grain cells (transfer cell (TC), aleurone cell (AL), and starchy endosperm (SE)). AATs highly expressed in roots are good candidates for amino acid uptake from soil whilst AATs highly expressed in senescing leaves and stems may be involved in translocation to grain. AATs in TC (TaAAP2 and TaAAP19) and SE (TaAAP13) may play important roles in determining grain protein content and grain yield. The expression levels of AAT homeologs showed unequal contributions in response to abiotic stresses and development, which may aid wheat adaptation to a wide range of environments.

Spatiotemporal expression patterns of wheat amino acid transporters reveal their putative roles in nitrogen transport and responses to abiotic stress

Amino acid transporters have roles in amino acid uptake from soil, long-distance transport, remobilization from vegetative tissues and accumulation in grain. Critically, the majority of wheat grain nitrogen is derived from amino acids remobilized from vegetative organs. However, no systematic analysis of wheat AAT genes has been reported to date. Here, 283 full length wheat AAT genes representing 100 distinct groups of homeologs were identified and curated by selectively consolidating IWGSC CSSv2 and TGACv1 Triticum aestivum genome assemblies and reassembling or mapping of IWGSC CSS chromosome sorted reads to fill any gaps. Gene expression profiling was performed using public RNA-seq data from root, leaf, stem, spike, grain and grain cells (transfer cell (TC), aleurone cell (AL), and starchy endosperm (SE)). AATs highly expressed in roots are good candidates for amino acid uptake from soil whilst AATs highly expressed in senescing leaves and stems may be involved in translocation to grain. AATs in TC (TaAAP2 and TaAAP19) and SE (TaAAP13) may play important roles in determining grain protein content and grain yield. The expression levels of AAT homeologs showed unequal contributions in response to abiotic stresses and development, which may aid wheat adaptation to a wide range of environments.

Characterisation of the nicotianamine aminotransferase and deoxymugineic acid synthase genes essential to Strategy II iron uptake in bread wheat (Triticum aestivum L.)

Iron (Fe) uptake in graminaceous plant species occurs via the release and uptake of Fe-chelating compounds known as mugineic acid family phytosiderophores (MAs). In the MAs biosynthetic pathway, nicotianamine aminotransferase (NAAT) and deoxymugineic acid synthase (DMAS) enzymes catalyse the formation of 2'-deoxymugineic acid (DMA) from nicotianamine (NA). Here we describe the identification and characterisation of six TaNAAT and three TaDMAS1 genes in bread wheat (Triticum aestivum L.). The coding sequences of all six TaNAAT homeologs consist of seven exons with ≥88.0% nucleotide sequence identity and most sequence variation present in the first exon. The coding sequences of the three TaDMAS1 homeologs consist of three exons with ≥97.8% nucleotide sequence identity. Phylogenetic analysis revealed that the TaNAAT and TaDMAS1 proteins are most closely related to the HvNAAT and HvDMAS1 proteins of barley and that there are two distinct groups of TaNAAT proteins-TaNAAT1 and TaNAAT2 -that correspond to the HvNAATA and HvNAATB proteins, respectively. Quantitative reverse transcription-PCR analysis revealed that the TaNAAT2 genes are expressed at highest levels in anther tissues whilst the TaNAAT1 and TaDMAS1 genes are expressed at highest levels in root tissues of bread wheat. Furthermore, the TaNAAT1, TaNAAT2 and TaDMAS1 genes were differentially regulated by plant Fe status and their expression was significantly upregulated in root tissues from day five onwards during a seven-day Fe deficiency treatment. The identification and characterization of the TaNAAT1, TaNAAT2 and TaDMAS1 genes provides a valuable genetic resource for improving bread wheat growth on Fe deficient soils and enhancing grain Fe nutrition.

The defence-associated transcriptome of hexaploid wheat displays homoeolog expression and induction bias

Bread wheat (Triticum aestivum L.) is an allopolyploid species containing three ancestral genomes. Therefore, three homoeologous copies exist for the majority of genes in the wheat genome. Whether different homoeologs are differentially expressed (homoeolog expression bias) in response to biotic and abiotic stresses is poorly understood. In this study, we applied a RNA-seq approach to analyse homoeolog-specific global gene expression patterns in wheat during infection by the fungal pathogen Fusarium pseudograminearum, which causes crown rot disease in cereals. To ensure specific detection of homoeologs, we first optimized read alignment methods and validated the results experimentally on genes with known patterns of subgenome-specific expression. Our global analysis identified widespread patterns of differential expression among homoeologs, indicating homoeolog expression bias underpins a large proportion of the wheat transcriptome. In particular, genes differentially expressed in response to Fusarium infection were found to be disproportionately contributed from B and D subgenomes. In addition, we found differences in the degree of responsiveness to pathogen infection among homoeologous genes with B and D homoeologs exhibiting stronger responses to pathogen infection than A genome copies. We call this latter phenomenon as 'homoeolog induction bias'. Understanding how homoeolog expression and induction biases operate may assist the improvement of biotic stress tolerance in wheat and other polyploid crop species.

Proteome Profiling of Wheat Shoots from Different Cultivars

Wheat is a cereal grain and one of the world's major food crops. Recent advances in wheat genome sequencing are by now facilitating its genomic and proteomic analyses. However, little is known about possible differences in total protein levels of hexaploid versus tetraploid wheat cultivars, and also knowledge of phosphorylated wheat proteins is still limited. Here, we performed a detailed analysis of the proteome of seedling leaves from two hexaploid wheat cultivars (Triticum aestivum L. Pavon 76 and USU-Apogee) and one tetraploid wheat (T. turgidum ssp. durum cv. Senatore Cappelli). Our shotgun proteomics data revealed that, whereas we observed some significant differences, overall a high similarity between hexaploid and tetraploid varieties with respect to protein abundance was observed. In addition, already at the seedling stage, a small set of proteins was differential between the small (USU-Apogee) and larger hexaploid wheat cultivars (Pavon 76), which could potentially act as growth predictors. Finally, the phosphosites identified in this study can be retrieved from the in-house developed plant PTM-Viewer (bioinformatics.psb.ugent.be/webtools/ptm_viewer/), making this the first searchable repository for phosphorylated wheat proteins. This paves the way for further in depth, quantitative (phospho)proteome-wide differential analyses upon a specific trigger or environmental change.

Identification and molecular characterization of the nicotianamine synthase gene family in bread wheat

Nicotianamine (NA) is a non-protein amino acid involved in fundamental aspects of metal uptake, transport and homeostasis in all plants and constitutes the biosynthetic precursor of mugineic acid family phytosiderophores (MAs) in graminaceous plant species. Nicotianamine synthase (NAS) genes, which encode enzymes that synthesize NA from S-adenosyl-L-methionine (SAM), are differentially regulated by iron (Fe) status in most plant species and plant genomes have been found to contain anywhere from 1 to 9 NAS genes. This study describes the identification of 21 NAS genes in the hexaploid bread wheat (Triticum aestivum L.) genome and their phylogenetic classification into two distinct clades. The TaNAS genes are highly expressed during germination, seedling growth and reproductive development. Fourteen of the clade I NAS genes were up-regulated in root tissues under conditions of Fe deficiency. Protein sequence analyses revealed the presence of endocytosis motifs in all of the wheat NAS proteins as well as chloroplast, mitochondrial and secretory transit peptide signals in four proteins. These results greatly expand our knowledge of NAS gene families in graminaceous plant species as well as the genetics underlying Fe nutrition in bread wheat.

A large-scale chromosome-specific SNP discovery guideline

Single-nucleotide polymorphisms (SNPs) are the most prevalent type of variation in genomes that are increasingly being used as molecular markers in diversity analyses, mapping and cloning of genes, and germplasm characterization. However, only a few studies reported large-scale SNP discovery in Aegilops tauschii, restricting their potential use as markers for the low-polymorphic D genome. Here, we report 68,592 SNPs found on the gene-related sequences of the 5D chromosome of Ae. tauschii genotype MvGB589 using genomic and transcriptomic sequences from seven Ae. tauschii accessions, including AL8/78, the only genotype for which a draft genome sequence is available at present. We also suggest a workflow to compare SNP positions in homologous regions on the 5D chromosome of Triticum aestivum, bread wheat, to mark single nucleotide variations between these closely related species. Overall, the identified SNPs define a density of 4.49 SNPs per kilobyte, among the highest reported for the genic regions of Ae. tauschii so far. To our knowledge, this study also presents the first chromosome-specific SNP catalog in Ae. tauschii that should facilitate the association of these SNPs with morphological traits on chromosome 5D to be ultimately targeted for wheat improvement.

Gene expression analysis reveals important pathways for drought response in leaves and roots of a wheat cultivar adapted to rainfed cropping in the Cerrado biome

Drought limits wheat production in the Brazilian Cerrado biome. In order to search for candidate genes associated to the response to water deficit, we analyzed the gene expression profiles, under severe drought stress, in roots and leaves of the cultivar MGS1 Aliança, a well-adapted cultivar to the Cerrado. A set of 4,422 candidate genes was found in roots and leaves. The number of down-regulated transcripts in roots was higher than the up-regulated transcripts, while the opposite occurred in leaves. The number of common transcripts between the two tissues was 1,249, while 2,124 were specific to roots and 1,049 specific to leaves. Quantitative RT-PCR analysis revealed a 0.78 correlation with the expression data. The candidate genes were distributed across all chromosomes and component genomes, but a greater number was mapped on the B genome, particularly on chromosomes 3B, 5B and 2B. When considering both tissues, 116 different pathways were induced. One common pathway, among the top three activated pathways in both tissues, was starch and sucrose metabolism. These results pave the way for future marker development and selection of important genes and are useful for understanding the metabolic pathways involved in wheat drought response.

HANDS2: accurate assignment of homoeallelic base-identity in allopolyploids despite missing data

Characterization of homoeallelic base-identity in allopolyploids is difficult since homeologous subgenomes are closely related and becomes further challenging if diploid-progenitor data is missing. We present HANDS2, a next-generation sequencing-based tool that enables highly accurate (>90%) genome-wide discovery of homeolog-specific base-identity in allopolyploids even in the absence of a diploid-progenitor. We applied HANDS2 to the transcriptomes of various cruciferous plants belonging to genus Brassica. Our results suggest that the three C genomes in Brassica are more similar to each other than the three A genomes, and provide important insights into the relationships between various Brassica tetraploids and their diploid-progenitors at a single-base resolution.