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

Isocitrate lyase promotes Puccinia striiformis f. sp. tritici susceptibility in wheat (Triticum aestivum) by suppressing accumulation of glyoxylate cycle intermediates

Plant fungal parasites manipulate host metabolism to support their own survival. Among the many central metabolic pathways altered during infection, the glyoxylate cycle is frequently upregulated in both fungi and their host plants. Here, we examined the response of the glyoxylate cycle in bread wheat (Triticum aestivum) to infection by the obligate biotrophic fungal pathogen Puccinia striiformis f. sp. tritici (Pst). Gene expression analysis revealed that wheat genes encoding the two unique enzymes of the glyoxylate cycle, isocitrate lyase (TaICL) and malate synthase, diverged in their expression between susceptible and resistant Pst interactions. Focusing on TaICL, we determined that the TaICL B homoeolog is specifically upregulated during early stages of a successful Pst infection. Furthermore, disruption of the B homoeolog alone was sufficient to significantly perturb Pst disease progression. Indeed, Pst infection of the TaICL-B disruption mutant (TaICL-BY400*) was inhibited early during initial penetration, with the TaICL-BY400* line also accumulating high levels of malic acid, citric acid, and aconitic acid. Exogenous application of malic acid or aconitic acid also suppressed Pst infection, with trans-aconitic acid treatment having the most pronounced effect by decreasing fungal biomass 15-fold. Thus, enhanced TaICL-B expression during Pst infection may lower accumulation of malic acid and aconitic acid to promote Pst proliferation. As exogenous application of aconitic acid and malic acid has previously been shown to inhibit other critical pests and pathogens, we propose TaICL as a potential target for disruption in resistance breeding that could have wide-reaching protective benefits for wheat and beyond.

The translational landscape of bread wheat during grain development

The dynamics of gene expression in crop grains has typically been investigated at the transcriptional level. However, this approach neglects translational regulation, a widespread mechanism that rapidly modulates gene expression to increase the plasticity of organisms. Here, we performed ribosome profiling and polysome profiling to obtain a comprehensive translatome data set of developing bread wheat (Triticum aestivum) grains. We further investigated the genome-wide translational dynamics during grain development, revealing that the translation of many functional genes is modulated in a stage-specific manner. The unbalanced translation between subgenomes is pervasive, which increases the expression flexibility of allohexaploid wheat. In addition, we uncovered widespread previously unannotated translation events, including upstream open reading frames (uORFs), downstream open reading frames (dORFs), and open reading frames (ORFs) in long noncoding RNAs, and characterized the temporal expression dynamics of small ORFs. We demonstrated that uORFs act as cis-regulatory elements that can repress or even enhance the translation of mRNAs. Gene translation may be combinatorially modulated by uORFs, dORFs, and microRNAs. In summary, our study presents a translatomic resource that provides a comprehensive and detailed overview of the translational regulation in developing bread wheat grains. This resource will facilitate future crop improvements for optimal yield and quality.

Transcriptome meta-analysis of abiotic stresses-responsive genes and identification of candidate transcription factors for broad stress tolerance in wheat

Under field conditions, wheat is subjected to single or multiple stress conditions. The elucidation of the molecular mechanism of stress response is a key step to identify candidate genes for stress resistance in plants. In this study, RNA-seq data analysis identified 17.324, 10.562, 5.510, and 8.653 differentially expressed genes (DEGs) under salt, drought, heat, and cold stress, respectively. Moreover, the comparison of DEGs from each stress revealed 2374 shared genes from which 40% showed highly conserved expression patterns. Moreover, co-expression network analysis and GO enrichment revealed co-expression modules enriched with genes involved in transcription regulation, protein kinase pathway, and genes responding to phytohormones or modulating hormone levels. The expression of 15 selected transcription factor encoding genes was analyzed under abiotic stresses and ABA treatment in durum wheat. The identified transcription factor genes are excellent candidates for genetic engineering of stress tolerance in wheat.

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.

Genomic and Meiotic Changes Accompanying Polyploidization

Hybridization and polyploidy have been considered as significant evolutionary forces in adaptation and speciation, especially among plants. Interspecific gene flow generates novel genetic variants adaptable to different environments, but it is also a gene introgression mechanism in crops to increase their agronomical yield. An estimate of 9% of interspecific hybridization has been reported although the frequency varies among taxa. Homoploid hybrid speciation is rare compared to allopolyploidy. Chromosome doubling after hybridization is the result of cellular defects produced mainly during meiosis. Unreduced gametes, which are formed at an average frequency of 2.52% across species, are the result of altered spindle organization or orientation, disturbed kinetochore functioning, abnormal cytokinesis, or loss of any meiotic division. Meiotic changes and their genetic basis, leading to the cytological diploidization of allopolyploids, are just beginning to be understood especially in wheat. However, the nature and mode of action of homoeologous recombination suppressor genes are poorly understood in other allopolyploids. The merger of two independent genomes causes a deep modification of their architecture, gene expression, and molecular interactions leading to the phenotype. We provide an overview of genomic changes and transcriptomic modifications that particularly occur at the early stages of allopolyploid formation.

Spatial and Temporal Transcriptomic Heredity and Asymmetry in an Artificially Constructed Allotetraploid Wheat (AADD)

Polyploidy, or whole-genome duplication (WGD), often induces dramatic changes in gene expression due to "transcriptome shock. " However, questions remain about how allopolyploidy (the merging of multiple nuclear genomes in the same nucleus) affects gene expression within and across multiple tissues and developmental stages during the initial foundation of allopolyploid plants. Here, we systematically investigated the immediate effect of allopolyploidy on gene expression variation in an artificial allopolyploidy system consisting of a constructed allotetraploid wheat (AADD genome, accession AT2) and its diploid progenitors Triticum urartu and Aegilops tauschii. We performed comprehensive RNA sequencing of 81 samples from different genotypes, tissues, and developmental stages. First, we found that intrinsic interspecific differences between the diploid parents played a major role in establishing the expression architecture of the allopolyploid. Nonetheless, allopolyploidy per se also induced dramatic and asymmetric patterns of differential gene expression between the subgenomes, and genes from the D subgenome exhibited a more drastic response. Second, analysis of homoeolog expression bias (HEB) revealed that the D subgenome exhibited significant expression bias and that de novo-generated HEB was attributed mainly to asymmetrical differential gene expression. Homoeolog-specific expression (HSE) analyses showed that the cis-only regulatory pattern was predominant in AT2, reflecting significant divergence between the parents. Co-expression network analysis revealed that homoeolog expression connectivity (HEC) was significantly correlated with sequence divergence in cis elements between subgenomes. Interestingly, allopolyploidy-induced reconstruction of network modules was also associated with different HSE patterns. Finally, a transcriptome atlas of spike development demonstrated that the phenotypic similarity of AT2 to T. urartu may be attributed to the combination of relatively stable expression of A-subgenome genes and drastic downregulation of their D-subgenome homoeologs. These findings provide a broad, multidimensional characterization of allopolyploidy-induced transcriptomic responses and suggest that allopolyploidy can have immediate and complex regulatory effects on the expression of nuclear genes.

Multi-tissue transcriptome analysis of two Begonia species reveals dynamic patterns of evolution in the chalcone synthase gene family

Begonia is an important horticultural plant group, as well as one of the most speciose Angiosperm genera, with over 2000 described species. Genus wide studies of genome size have shown that Begonia has a highly variable genome size, and analysis of paralog pairs has previously suggested that Begonia underwent a whole genome duplication. We address the contribution of gene duplication to the generation of diversity in Begonia using a multi-tissue RNA-seq approach. We chose to focus on chalcone synthase (CHS), a gene family having been shown to be involved in biotic and abiotic stress responses in other plant species, in particular its importance in maximising the use of variable light levels in tropical plants. We used RNA-seq to sample six tissues across two closely related but ecologically and morphologically divergent species, Begonia conchifolia and B. plebeja, yielding 17,012 and 19,969 annotated unigenes respectively. We identified the chalcone synthase gene family members in our Begonia study species, as well as in Hillebrandia sandwicensis, the monotypic sister genus to Begonia, Cucumis sativus, Arabidopsis thaliana, and Zea mays. Phylogenetic analysis suggested the CHS gene family has high duplicate turnover, all members of CHS identified in Begonia arising recently, after the divergence of Begonia and Cucumis. Expression profiles were similar within orthologous pairs, but we saw high inter-ortholog expression variation. Sequence analysis showed relaxed selective constraints on some ortholog pairs, with substitutions at conserved sites. Evidence of pseudogenisation and species specific duplication indicate that lineage specific differences are already beginning to accumulate since the divergence of our study species. We conclude that there is evidence for a role of gene duplication in generating diversity through sequence and expression divergence in Begonia.

Homoeolog Inference Methods Requiring Bidirectional Best Hits or Synteny Miss Many Pairs

Homoeologs are pairs of genes or chromosomes in the same species that originated by speciation and were brought back together in the same genome by allopolyploidization. Bioinformatic methods for accurate homoeology inference are crucial for studying the evolutionary consequences of polyploidization, and homoeology is typically inferred on the basis of bidirectional best hit (BBH) and/or positional conservation (synteny). However, these methods neglect the fact that genes can duplicate and move, both prior to and after the allopolyploidization event. These duplications and movements can result in many-to-many and/or nonsyntenic homoeologs-which thus remain undetected and unstudied. Here, using the allotetraploid upland cotton (Gossypium hirsutum) as a case study, we show that conventional approaches indeed miss a substantial proportion of homoeologs. Additionally, we found that many of the missed pairs of homoeologs are broadly and highly expressed. A gene ontology analysis revealed a high proportion of the nonsyntenic and non-BBH homoeologs to be involved in protein translation and are likely to contribute to the functional repertoire of cotton. Thus, from an evolutionary and functional genomics standpoint, choosing a homoeolog inference method which does not solely rely on 1:1 relationship cardinality or synteny is crucial for not missing these potentially important homoeolog pairs.

Genomics of Evolutionary Novelty in Hybrids and Polyploids

It has long been recognized that hybridization and polyploidy are prominent processes in plant evolution. Although classically recognized as significant in speciation and adaptation, recognition of the importance of interspecific gene flow has dramatically increased during the genomics era, concomitant with an unending flood of empirical examples, with or without genome doubling. Interspecific gene flow is thus increasingly thought to lead to evolutionary innovation and diversification, via adaptive introgression, homoploid hybrid speciation and allopolyploid speciation. Less well understood, however, are the suite of genetic and genomic mechanisms set in motion by the merger of differentiated genomes, and the temporal scale over which recombinational complexity mediated by gene flow might be expressed and exposed to natural selection. We focus on these issues here, considering the types of molecular genetic and genomic processes that might be set in motion by the saltational event of genome merger between two diverged species, either with or without genome doubling, and how these various processes can contribute to novel phenotypes. Genetic mechanisms include the infusion of new alleles and the genesis of novel structural variation including translocations and inversions, homoeologous exchanges, transposable element mobilization and novel insertional effects, presence-absence variation and copy number variation. Polyploidy generates massive transcriptomic and regulatory alteration, presumably set in motion by disrupted stoichiometries of regulatory factors, small RNAs and other genome interactions that cascade from single-gene expression change up through entire networks of transformed regulatory modules. We highlight both these novel combinatorial possibilities and the range of temporal scales over which such complexity might be generated, and thus exposed to natural selection and drift.

Homoeologous gene expression and co-expression network analyses and evolutionary inference in allopolyploids

Polyploidy is a widespread phenomenon throughout eukaryotes. Due to the coexistence of duplicated genomes, polyploids offer unique challenges for estimating gene expression levels, which is essential for understanding the massive and various forms of transcriptomic responses accompanying polyploidy. Although previous studies have explored the bioinformatics of polyploid transcriptomic profiling, the causes and consequences of inaccurate quantification of transcripts from duplicated gene copies have not been addressed. Using transcriptomic data from the cotton genus (Gossypium) as an example, we present an analytical workflow to evaluate a variety of bioinformatic method choices at different stages of RNA-seq analysis, from homoeolog expression quantification to downstream analysis used to infer key phenomena of polyploid expression evolution. In general, EAGLE-RC and GSNAP-PolyCat outperform other quantification pipelines tested, and their derived expression dataset best represents the expected homoeolog expression and co-expression divergence. The performance of co-expression network analysis was less affected by homoeolog quantification than by network construction methods, where weighted networks outperformed binary networks. By examining the extent and consequences of homoeolog read ambiguity, we illuminate the potential artifacts that may affect our understanding of duplicate gene expression, including an overestimation of homoeolog co-regulation and the incorrect inference of subgenome asymmetry in network topology. Taken together, our work points to a set of reasonable practices that we hope are broadly applicable to the evolutionary exploration of polyploids.

Cis-trans controls and regulatory novelty accompanying allopolyploidization

Allopolyploidy is a prevalent process in plants, having important physiological, ecological and evolutionary consequences. Transcriptomic responses to genomic merger and doubling have been demonstrated in many allopolyploid systems, encompassing a diversity of phenomena including homoeolog expression bias, genome dominance, expression-level dominance and revamping of co-expression networks. Notwithstanding the foregoing, there remains a need to develop a conceptual framework that will stimulate a deeper understanding of these diverse phenomena and their mechanistic interrelationships. Here we introduce considerations relevant to this framework with a focus on cis-trans interactions among duplicated genes and alleles in hybrids and allopolyploids. By extending classic allele-specific expression analysis to the allopolyploid level, we distinguish the distinct effects of progenitor regulatory interactions from the novel intergenomic interactions that arise from genome merger and allopolyploidization. This perspective informs experiments designed to reveal the molecular genetic basis of gene regulatory control, and will facilitate the disentangling of genetic from epigenetic and higher-order effects that impact gene expression. Finally, we suggest that the extended cis-trans model may help conceptually unify several presently disparate hallmarks of allopolyploid evolution, including genome-wide expression dominance and biased fractionation, and lead to a new level of understanding of phenotypic novelty accompanying polyploidy.

Boosting Genetic Gain in Allogamous Crops via Speed Breeding and Genomic Selection

Breeding schemes that utilize modern breeding methods like genomic selection (GS) and speed breeding (SB) have the potential to accelerate genetic gain for different crops. We investigated through stochastic computer simulation the advantages and disadvantages of adopting both GS and SB (SpeedGS) into commercial breeding programs for allogamous crops. In addition, we studied the effect of omitting one or two selection stages from the conventional phenotypic scheme on GS accuracy, genetic gain, and inbreeding. As an example, we simulated GS and SB for five traits (heading date, forage yield, seed yield, persistency, and quality) with different genetic architectures and heritabilities (0.7, 0.3, 0.4, 0.1, and 0.3; respectively) for a tall fescue breeding program. We developed a new method to simulate correlated traits with complex architectures of which effects can be sampled from multiple distributions, e.g. to simulate the presence of both minor and major genes. The phenotypic selection scheme required 11 years, while the proposed SpeedGS schemes required four to nine years per cycle. Generally, SpeedGS schemes resulted in higher genetic gain per year for all traits especially for traits with low heritability such as persistency. Our results showed that running more SB rounds resulted in higher genetic gain per cycle when compared to phenotypic or GS only schemes and this increase was more pronounced per year when cycle time was shortened by omitting cycle stages. While GS accuracy declined with additional SB rounds, the decline was less in round three than in round two, and it stabilized after the fourth SB round. However, more SB rounds resulted in higher inbreeding rate, which could limit long-term genetic gain. The inbreeding rate was reduced by approximately 30% when generating the initial population for each cycle through random crosses instead of generating half-sib families. Our study demonstrated a large potential for additional genetic gain from combining GS and SB. Nevertheless, methods to mitigate inbreeding should be considered for optimal utilization of these highly accelerated breeding programs.