Genome evolution due to allopolyploidization in wheat

Evaluation and Identification of Powdery Mildew Resistance Genes in Aegilops tauschii and Emmer Wheat Accessions

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is a serious threat to wheat (Triticum aestivum L.) production. Narrow genetic basis of common wheat boosted the demand for diversified donors against powdery mildew. Aegilops tauschii Coss (2n = 2x = DD) and emmer wheat (2n = 4x = AABB), as the ancestor species of common wheat, are important gene donors for genetic improvement of common wheat. In this study, a total of 71 Ae. tauschii and 161 emmer wheat accessions were first evaluated for their powdery mildew resistance using the Bgt isolate E09. Thirty-three Ae. tauschii (46.5%) and 108 emmer wheat accessions (67.1%) were resistant. Then, all these accessions were tested by the diagnostic markers for 21 known Pm genes. The results showed that Pm2 alleles were detected in all the 71 Ae. tauschii and only Pm4 alleles were detected in 20 of 161 emmer wheat accessions. After haplotype analysis, we identified four Pm4 alleles (Pm4a, Pm4b, Pm4d, and Pm4f) in the emmer wheat accessions and three Pm2 alleles (Pm2d, Pm2e, and Pm2g) in the Ae. tauschii. Further resistance spectrum analysis indicated that these resistance accessions displayed different resistance reactions to different Bgt isolates, implying they may have other Pm genes apart from Pm2 and/or Pm4 alleles. Notably, a new Pm2 allele, Pm2S, was identified in Ae. tauschii, which contained a 64-bp deletion in the first exon and formed a new termination site at the 513th triplet of the shifted reading frame compared with reported Pm2 alleles. The phylogenetic tree of Pm2S showed that the kinship of Pm2S was close to Pm2h. To efficiently and accurately detect Pm2S and distinguish with other Pm2 alleles in Ae. tauschii background, a diagnostic marker, YTU-QS-3, was developed, and its effectiveness was verified. This study provided valuable Pm alleles and enriched the genetic diversity of the powdery mildew resistance in wheat improvement.

Major chromosome rearrangements in intergeneric wheat × rye hybrids in compatible and incompatible crosses detected by GBS read coverage analysis

The presence of incompatibility alleles in primary amphidiploids constitutes a reproductive barrier in newly synthesized wheat-rye hybrids. To overcome this barrier, the genome stabilization process includes large-scale chromosome rearrangements. In incompatible crosses resulting in fertile amphidiploids, the elimination of one of the incompatible alleles Eml-A1 or Eml-R1b can occur already in the somatic tissue of the wheat × rye hybrid embryo. We observed that the interaction of incompatible loci Eml-A1 of wheat and Eml-R1b of rye after overcoming embryo lethality leads to hybrid sterility in primary triticale. During subsequent seed reproductions (R1, R2 or R3) most of the chromosomes of A, B, D and R subgenomes undergo rearrangement or eliminations to increase the fertility of the amphidiploid by natural selection. Genotyping-by-sequencing (GBS) coverage analysis showed that improved fertility is associated with the elimination of entire and partial chromosomes carrying factors that either cause the disruption of plant development in hybrid plants or lead to the restoration of the euploid number of chromosomes (2n = 56) in the absence of one of the incompatible alleles. Highly fertile offspring obtained in compatible and incompatible crosses can be successfully adapted for the production of triticale pre-breeding stocks.

Dosage effect genes modulate grain development in synthesized Triticum durum-Haynaldia villosa allohexaploid

Polyploidization in plants often leads to increased cell size and grain size, which may be affected by the increased genome dosage and transcription abundance. The synthesized Triticum durum (AABB)-Haynaldia villosa (VV) amphiploid (AABBVV) has significantly increased grain size, especially grain length, than the tetraploid and diploid parents. To investigate how polyploidization affects grain development at the transcriptional level, we perform transcriptome analysis using the immature seeds of T. durum, H. villosa, and the amphiploid. The dosage effect genes are contributed more by differentially expressed genes from genome V of H. villosa. The dosage effect genes overrepresent grain development-related genes. Interestingly, the vernalization gene TaVRN1 is among the positive dosage effect genes in the T. durum‒H. villosa and T. turgidum‒Ae. tauschii amphiploids. The expression levels of TaVRN1 homologs are positively correlated with the grain size and weight. The TaVRN1-B1 or TaVRN1-D1 mutation shows delayed florescence, decreased cell size, grain size, and grain yield. These data indicate that dosage effect genes could be one of the important explanations for increased grain size by regulating grain development. The identification and functional validation of dosage effect genes may facilitate the finding of valuable genes for improving wheat yield.

The genetic consequences of range expansion and its influence on diploidization in polyploids

Despite newly formed polyploids being subjected to myriad fitness consequences, the relative prevalence of polyploidy both contemporarily and in ancestral branches of the tree of life suggests alternative advantages that outweigh these consequences. One proposed advantage is that polyploids may more easily colonize novel habitats such as deglaciated areas. However, previous research conducted in diploids suggests that range expansion comes with a fitness cost as deleterious mutations may fix rapidly on the expansion front. Here, we interrogate the potential consequences of expansion in polyploids by conducting spatially explicit forward-in-time simulations to investigate how ploidy and inheritance patterns impact the relative ability of polyploids to expand their range. We show that under realistic dominance models, autopolyploids suffer greater fitness reductions than diploids as a result of range expansion due to the fixation of increased mutational load that is masked in the range core. Alternatively, the disomic inheritance of allopolyploids provides a shield to this fixation resulting in minimal fitness consequences. In light of this advantage provided by disomy, we investigate how range expansion may influence cytogenetic diploidization through the reversion to disomy in autotetraploids. We show that under a wide range of parameters investigated for two models of diploidization, disomy frequently evolves more rapidly on the expansion front than in the range core, and that this dynamic inheritance model has additional effects on fitness. Together our results point to a complex interaction between dominance, ploidy, inheritance, and recombination on fitness as a population spreads across a geographic range.

Genomic and phenotypic characterization of finger millet indicates a complex diversification history

Advances in sequencing technologies mean that insights into crop diversification can now be explored in crops beyond major staples. We use a genome assembly of finger millet, an allotetraploid orphan crop, to analyze DArTseq single nucleotide polymorphisms (SNPs) at the whole and sub-genome level. A set of 8778 SNPs and 13 agronomic traits was used to characterize a diverse panel of 423 landraces from Africa and Asia. Through principal component analysis (PCA) and discriminant analysis of principal components, four distinct groups of accessions were identified that coincided with the primary geographic regions of finger millet cultivation. Notably, East Africa, presumed to be the crop's origin, exhibited the lowest genetic diversity. The PCA of phenotypic data also revealed geographic differentiation, albeit with differing relationships among geographic areas than indicated with genomic data. Further exploration of the sub-genomes A and B using neighbor-joining trees revealed distinct features that provide supporting evidence for the complex evolutionary history of finger millet. Although genome-wide association study found only a limited number of significant marker-trait associations, a clustering approach based on the distribution of marker effects obtained from a ridge regression genomic model was employed to investigate trait complexity. This analysis uncovered two distinct clusters. Overall, the findings suggest that finger millet has undergone complex and context-specific diversification, indicative of a lengthy domestication history. These analyses provide insights for the future development of finger millet.

Polyploid speciation in Zea (Poaceae): cytogenetic insights

MAIN CONCLUSION

The analysis of meiotic pairing affinities and genomic formulae in species and hybrids of Zea allowed us to speculate an evolutionary model to recreate the ancient polyploidization of maize and allied species. The meiotic pairing affinities and the genomic formulae analysis in Zea species and hybrids obtained in new and previous crosses, together with the molecular data known in the genus, allowed us to speculate an evolutionary model to attempt to recreate the ancient polyploidization process of Zea species. We propose that x = 5 semispecies are the ancestors of all modern species of the genus. The complex evolutionary process that originated the different taxa could be included hybridization between sympatric diploid ancestral semispecies (2n = 10) and recurrent duplication of the hybrid chromosome number, resulting in distinct auto- and allopolyploids. After the merger and doubling of independent genomes would have undergone cytological and genetical diploidization, implying revolutionary changes in genome organization and genic balance processes. Based on the meiotic behaviour of the 2n = 30 hybrids, that showed homoeology between the A subgenomes of all parental species, we propose that this subgenome A would be pivotal in all the species and would have conserved the rDNA sequences and the pairing regulator locus (PrZ). In the hypothetical model postulated here, the ancestral semispecies with the pivotal subgenome A would have had a wide geographic distribution, co-occurring and hybridizing with the semispecies harbouring B subgenomes, thus enabling sympatric speciation.

Genome dosage alteration caused by chromosome pyramiding and shuffling effects on karyotypic heterogeneity, reproductive diversity, and phenotypic variation in Zea-Tripsacum allopolyploids

KEY MESSAGE

We developed an array of Zea-Tripsacum tri-hybrid allopolyploids with multiple ploidies. We unveiled that changes in genome dosage due to the chromosomes pyramiding and shuffling of three species effects karyotypic heterogeneity, reproductive diversity, and phenotypic variation in Zea-Tripsacum allopolyploids. Polyploidy, or whole genome duplication, has played a major role in evolution and speciation. The genomic consequences of polyploidy have been extensively studied in many plants; however, the extent of chromosomal variation, genome dosage, phenotypic diversity, and heterosis in allopolyploids derived from multiple species remains largely unknown. To address this question, we synthesized an allohexaploid involving Zea mays, Tripsacum dactyloides, and Z. perennis by chromosomal pyramiding. Subsequently, an allooctoploid and an allopentaploid were obtained by hybridization of the allohexaploid with Z. perennis. Moreover, we constructed three populations with different ploidy by chromosomal shuffling (allopentaploid × Z. perennis, allohexaploid × Z. perennis, and allooctoploid × Z. perennis). We have observed 3 types of sexual reproductive modes and 2 types of asexual reproduction modes in the tri-species hybrids, including 2n gamete fusion (2n + n), haploid gamete fusion (n + n), polyspermy fertilization (n + n + n) or 2n gamete fusion (n + 2n), haploid gametophyte apomixis, and asexual reproduction. The tri-hybrids library presents extremely rich karyotype heterogeneity. Chromosomal compensation appears to exist between maize and Z. perennis. A rise in the ploidy of the trihybrids was linked to a higher frequency of chromosomal translocation. Variation in the degree of phenotypic diversity observed in different segregating populations suggested that genome dosage effects phenotypic manifestation. These findings not only broaden our understanding of the mechanisms of polyploid formation and reproductive diversity but also provide a novel insight into genome pyramiding and shuffling driven genome dosage effects and phenotypic diversity.

The Calcium-Dependent Protein Kinase TaCDPK7 Positively Regulates Wheat Resistance to Puccinia striiformis f. sp. tritici

Ca2+ plays a crucial role as a secondary messenger in plant development and response to abiotic/biotic stressors. Calcium-dependent protein kinases (CDPKs/CPKs) are essential Ca2+ sensors that can convert Ca2+ signals into downstream phosphorylation signals. However, there is limited research on the function of CDPKs in the context of wheat-Puccinia striiformis f. sp. tritici (Pst) interaction. In this study, we aimed to address this gap by identifying putative CDPK genes from the wheat reference genome and organizing them into four phylogenetic clusters (I-IV). To investigate the expression patterns of the TaCDPK family during the wheat-Pst interaction, we analyzed time series RNA-seq data and further validated the results through qRT-PCR assays. Among the TaCDPK genes, TaCDPK7 exhibited a significant induction during the wheat-Pst interaction, suggesting that it has a potential role in wheat resistance to Pst. To gain further insights into the function of TaCDPK7, we employed virus-induced gene silencing (VIGS) to knock down its expression which resulted in impaired wheat resistance to Pst, accompanied by decreased accumulation of hydrogen peroxide (H2O2), increased fungal biomass ratio, reduced expression of defense-related genes, and enhanced pathogen hyphal growth. These findings collectively suggest that TaCDPK7 plays an important role in wheat resistance to Pst. In summary, this study expands our understanding of wheat CDPKs and provides novel insights into their involvement in the wheat-Pst interaction.

Hybridization in agricultural weeds: A review from ecological, evolutionary, and management perspectives

Agricultural weeds frequently hybridize with each other or with related crop species. Some hybrid weeds exhibit heterosis (hybrid vigor), which may be stabilized through mechanisms like genome duplication or vegetative reproduction. Even when heterosis is not stabilized, hybridization events diversify weed gene pools and often enable adaptive introgression. Consequently, hybridization may promote weed evolution and exacerbate weed-crop competition. However, hybridization does not always increase weediness. Even when viable and fertile, hybrid weeds sometimes prove unsuccessful in crop fields. This review provides an overview of weed hybridization and its management implications. We describe intrinsic and extrinsic factors that influence hybrid fitness in agroecosystems. We also survey the rapidly growing literature on crop-weed hybridization and the link between hybridization and invasiveness. These topics are increasingly relevant in this era of genetic tools for crop improvement, intensive and simplified cropping systems, and globalized trade. The review concludes with suggested research priorities, including hybridization in the context of climate change, plant-insect interactions, and redesigned weed management programs. From a weed management perspective, hybridization is one of many reasons that researchers and land managers must diversify their weed control toolkits.

Molecular cytological analysis of alien introgressions in common wheat lines created by crossing of Triticum aestivum with T. dicoccoides and T. dicoccum

Wild and domesticated emmer (ВВАА, 2n = 28) are of significant interest for expanding the genetic diversity of common wheat as sources of a high protein and microelement grain content, resistance to many biotic and abiotic factors. Particular interest in these species is also determined by their close relationship with Triticum aestivum L., which facilitates interspecific hybridization. The objective of this work was to analyze the nature of alien introgressions in hybrid lines from crossing common wheat varieties with T. dicoccoides and T. dicoccum, and to assess the effect of their genome fragments on the cytological stability of introgression lines. A C-banding technique and genotyping with SNP and SSR markers were used to determine localization and length of introgression fragments. Assessment of cytological stability was carried out on the basis of chromosome behavior in microsporogenesis. A molecular cytogenetic analysis of introgression wheat lines indicated that the inclusion of the genetic material of wild and domesticated emmer was carried out mainly in the form of whole arms or large fragments in the chromosomes of the B genome and less extended inserts in the A genome. At the same time, the highest frequency of introgressions of the emmer genome was observed in chromosomes 1A, 1B, 2B, and 3B. The analysis of the final stage of meiosis showed a high level of cytological stability in the vast majority of introgression wheat lines (meiotic index was 83.0-99.0 %), which ensures the formation of functional gametes in an amount sufficient for successful reproduction. These lines are of interest for the selection of promising material with agronomically valuable traits and their subsequent inclusion in the breeding process.

Einkorn genomics reveals ancient roots of domesticated wheat

Wheat is an important staple food crop that underwent complex genome duplications. During domestication, genetic changes occurred, improving modern wheat, but understanding its phylogenetic history has been lacking. Mahmood and Mansoor discuss a recent publication demonstrating the evolutionary history of domesticated wheat (Triticum monococcum), providing opportunities for advancements in cereal improvement.

Comparative Analysis Reveals Different Evolutionary Fates and Biological Functions in Wheat Duplicated Genes (Triticum aestivum L.)

Wheat (Triticum aestivum L.) is a staple food crop that provides 20% of total human calorie consumption. Gene duplication has been considered to play an important role in evolution by providing new genetic resources. However, the evolutionary fates and biological functions of the duplicated genes in wheat remain to be elucidated. In this study, the resulting data showed that the duplicated genes evolved faster with shorter gene lengths, higher codon usage bias, lower expression levels, and higher tissue specificity when compared to non-duplicated genes. Our analysis further revealed functions of duplicated genes in various biological processes with significant enrichment to environmental stresses. In addition, duplicated genes derived from dispersed, proximal, tandem, transposed, and whole-genome duplication differed in abundance, evolutionary rate, gene compactness, expression pattern, and genetic diversity. Tandem and proximal duplicates experienced stronger selective pressure and showed a more compact gene structure with diverse expression profiles than other duplication modes. Moreover, genes derived from different duplication modes showed an asymmetrical evolutionary pattern for wheat A, B, and D subgenomes. Several candidate duplication hotspots associated with wheat domestication or polyploidization were characterized as potential targets for wheat molecular breeding. Our comprehensive analysis revealed the evolutionary trajectory of duplicated genes and laid the foundation for future functional studies on wheat.

Contribution of Duplicated Nucleotide-Binding Leucine-Rich Repeat (NLR) Genes to Wheat Disease Resistance

Wheat has a large and diverse repertoire of NLRs involved in disease resistance, with over 1500 NLRs detected in some studies. These NLR genes occur as singletons or clusters containing copies of NLRs from different phylogenetic clades. The number of NLRs and cluster size can differ drastically among ecotypes and cultivars. Primarily, duplication has led to the evolution and diversification of NLR genes. Among the various mechanisms, whole genome duplication (WGD) is the most intense and leading cause, contributing to the complex evolutionary history and abundant gene set of hexaploid wheat. Tandem duplication or recombination is another major mechanism of NLR gene expansion in wheat. The diversity and divergence of duplicate NLR genes are responsible for the broad-spectrum resistance of most plant species with limited R genes. Understanding the mechanisms underlying the rapid evolution and diversification of wheat NLR genes will help improve disease resistance in crops. The present review focuses on the diversity and divergence of duplicate NLR genes and their contribution to wheat disease resistance. Moreover, we provide an overview of disease resistance-associated gene duplication and the underlying strategies in wheat.

Genome-wide analysis of heavy metal ATPases (HMAs) in Poaceae species and their potential role against copper stress in Triticum aestivum

Plants require copper for normal growth and development and have evolved an efficient system for copper management based on transport proteins such as P_1B-ATPases, also known as heavy metal ATPases (HMAs). Here, we report HMAs in eleven different Poaceae species, including wheat. Furthermore, the possible role of wheat HMAs in copper stress was investigated. BlastP searches identified 27 HMAs in wheat, and phylogenetic analysis based on the Maximum Likelihood method demonstrated a separation into four distinct clades. Conserved motif analysis, domain identification, gene structure, and transmembrane helices number were also identified for wheat HMAs using computational tools. Wheat seedlings grown hydroponically were subjected to elevated copper and demonstrated toxicity symptoms with effects on fresh weight and changes in expression of selected HMAs TaHMA7, TaHMA8, and TaHMA9 were upregulated in response to elevated copper, suggesting a role in wheat copper homeostasis. Further investigations on these heavy metal pumps can provide insight into strategies for enhancing crop heavy metal tolerance in the face of heavy metal pollution.

Differential alternative polyadenylation of homoeologous genes of allohexaploid wheat ABD subgenomes during drought stress response

Because allohexaploid wheat genome contains ABD subgenomes, how the expression of homoeologous genes is coordinated remains largely unknown, particularly at the co-transcriptional level. Alternative polyadenylation (APA) is an important part of co-transcriptional regulation, which is crucial in developmental processes and stress responses. Drought stress is a major threat to the stable yield of wheat. Focusing on APA, we used poly(A) tag sequencing to track poly(A) site dynamics in wheat under drought stress. The results showed that drought stress led to extensive APA involving 37-47% of differentially expressed genes in wheat. Significant poly(A) site switching was found in stress-responsive genes. Interestingly, homoeologous genes exhibit unequal numbers of poly(A) sites, divergent APA patterns with tissue specificity and time-course dynamics, and distinct 3'-UTR length changes. Moreover, differentially expressed transcripts in leaves and roots used different poly(A) signals, the up- and downregulated isoforms had distinct preferences for non-canonical poly(A) sites. Genes that encode key polyadenylation factors showed differential expression patterns under drought stress. In summary, poly(A) signals and the changes in core poly(A) factors may widely affect the selection of poly(A) sites and gene expression levels during the response to drought stress, and divergent APA patterns among homoeologous genes add extensive plasticity to this responsive network. These results not only reveal the significant role of APA in drought stress response, but also provide a fresh perspective on how homoeologous genes contribute to adaptability through transcriptome diversity. In addition, this work provides information about the ends of transcripts for a better annotation of the wheat genome.

Global transcriptome analysis of allopolyploidization reveals large-scale repression of the D-subgenome in synthetic hexaploid wheat

Synthetic hexaploid wheat (SHW) lines are created as pre-breeding germplasm to diversify the D subgenome of hexaploid wheat and capitalize upon the untapped genetic diversity of the Aegilops tauschii gene pool. However, the phenotypes observed in the Ae. tauschii parents are not always recovered in the SHW lines, possibly due to inter-subgenome interactions. To elucidate this post-polyploidization genome reprogramming phenomenon, we performed RNA-seq of four SHW lines and their corresponding tetraploid and diploid parents, across ten tissues and three biological replicates. Homoeologue expression bias (HEB) analysis using more than 18,000 triads suggests massive suppression of homoeoalleles of the D subgenome in SHWs. Comparative transcriptome analysis of the whole-genome gene set further corroborated this finding. Alternative splicing analysis of the high-confidence genes indicates an additional layer of complexity where all five splice events are identified, and retained intron is predominant. Homoeologue expression upon resynthesis of hexaploid wheat has implications to the usage and handling of this germplasm in breeding as it relates to capturing the effects of epistatic interaction across subgenomes upon polyploidization. Special considerations must be given to this germplasm in pre-breeding activities to consider the extent of the inter-subgenome interactions on gene expression and their impact on traits for crop improvement.

Asymmetric gene expression and cell-type-specific regulatory networks in the root of bread wheat revealed by single-cell multiomics analysis

BACKGROUND

Homoeologs are defined as homologous genes resulting from allopolyploidy. Bread wheat, Triticum aestivum, is an allohexaploid species with many homoeologs. Homoeolog expression bias, referring to the relative contribution of homoeologs to the transcriptome, is critical for determining the traits that influence wheat growth and development. Asymmetric transcription of homoeologs has been so far investigated in a tissue or organ-specific manner, which could be misleading due to a mixture of cell types.

RESULTS

Here, we perform single nuclei RNA sequencing and ATAC sequencing of wheat root to study the asymmetric gene transcription, reconstruct cell differentiation trajectories and cell-type-specific gene regulatory networks. We identify 22 cell types. We then reconstruct cell differentiation trajectories that suggest different origins between epidermis/cortex and endodermis, distinguishing bread wheat from Arabidopsis. We show that the ratio of asymmetrically transcribed triads varies greatly when analyzing at the single-cell level. Hub transcription factors determining cell type identity are also identified. In particular, we demonstrate that TaSPL14 participates in vasculature development by regulating the expression of BAM1. Combining single-cell transcription and chromatin accessibility data, we construct the pseudo-time regulatory network driving root hair differentiation. We find MYB3R4, REF6, HDG1, and GATAs as key regulators in this process.

CONCLUSIONS

Our findings reveal the transcriptional landscape of root organization and asymmetric gene transcription at single-cell resolution in polyploid wheat.

WheatCENet: A Database for Comparative Co-expression Networks Analysis of Allohexaploid Wheat and Its Progenitors

Genetic and epigenetic changes after polyploidization events could result in variable gene expression and modified regulatory networks. Here, using large-scale transcriptome data, we constructed co-expression networks for diploid, tetraploid, and hexaploid wheat species, and built a platform for comparing co-expression networks of allohexaploid wheat and its progenitors, named WheatCENet. WheatCENet is a platform for searching and comparing specific functional co-expression networks, as well as identifying the related functions of the genes clustered therein. Functional annotations like pathways, gene families, protein-protein interactions, microRNAs (miRNAs), and several lines of epigenome data are integrated into this platform, and Gene Ontology (GO) annotation, gene set enrichment analysis (GSEA), motif identification, and other useful tools are also included. Using WheatCENet, we found that the network of WHEAT ABERRANT PANICLE ORGANIZATION 1 (WAPO1) has more co-expressed genes related to spike development in hexaploid wheat than its progenitors. We also found a novel motif of CCWWWWWWGG (CArG) specifically in the promoter region of WAPO-A1, suggesting that neofunctionalization of the WAPO-A1 gene affects spikelet development in hexaploid wheat. WheatCENet is useful for investigating co-expression networks and conducting other analyses, and thus facilitates comparative and functional genomic studies in wheat. WheatCENet is freely available at http://bioinformatics.cpolar.cn/WheatCENet and http://bioinformatics.cau.edu.cn/WheatCENet.

Genome–Transcriptome Transition Approaches to Characterize Anthocyanin Biosynthesis Pathway Genes in Blue, Black and Purple Wheat

Colored wheat has gained enormous attention from the scientific community, but the information available on the anthocyanin biosynthetic genes is very minimal. The study involved their genome-wide identification, in silico characterization and differential expression analysis among purple, blue, black and white wheat lines. The recently released wheat genome mining putatively identified eight structural genes in the anthocyanin biosynthesis pathway with a total of 1194 isoforms. Genes showed distinct exon architecture, domain profile, regulatory elements, chromosome emplacement, tissue localization, phylogeny and synteny, indicative of their unique function. RNA sequencing of developing seeds from colored (black, blue and purple) and white wheats identified differential expressions in 97 isoforms. The F3H on group two chromosomes and F3′5′H on 1D chromosomes could be significant influencers in purple and blue color development, respectively. Apart from a role in anthocyanin biosynthesis, these putative structural genes also played an important role in light, drought, low temperature and other defense responses. The information can assist in targeted anthocyanin production in the wheat seed endosperm.

The TaTCP4/10-B1 cascade regulates awn elongation in wheat (Triticum aestivum L.)

Awns are important morphological markers for wheat and exert a strong physiological effect on wheat yield. The awn elongation suppressor B1 has recently been cloned through association and linkage analysis in wheat. However, the mechanism of awn inhibition centered around B1 remains to be clarified. Here, we identified an allelic variant in the coding region of B1 through analysis of re-sequencing data; this variant causes an amino acid substitution and premature termination, resulting in a long-awn phenotype. Transcriptome analysis indicated that B1 inhibited awn elongation by impeding cytokinin- and auxin-promoted cell division. Moreover, B1 directly repressed the expression of TaRAE2 and TaLks2, whose orthologs have been reported to promote awn development in rice or barley. More importantly, we found that TaTCP4 and TaTCP10 synergistically inhibited the expression of B1, and a G-to-A mutation in the B1 promoter attenuated its inhibition by TaTCP4/10. Taken together, our results reveal novel mechanisms of awn development and provide genetic resources for trait improvement in wheat.

Genome sequencing of Sitopsis species provides insights into their contribution to the B subgenome of bread wheat

Wheat (Triticum aestivum, BBAADD) is an allohexaploid species that originated from two polyploidization events. The progenitors of the A and D subgenomes have been identified as Triticum urartu and Aegilops tauschii, respectively. Current research suggests that Aegilops speltoides is the closest but not the direct ancestor of the B subgenome. However, whether Ae. speltoides has contributed genomically to the wheat B subgenome and which chromosome regions are conserved between Ae. speltoides and the B subgenome remain unclear. Here, we assembled a high-quality reference genome for Ae. speltoides, resequenced 53 accessions from seven species (Aegilops bicornis, Aegilops longissima, Aegilops searsii, Aegilops sharonensis, Ae. speltoides, Aegilops mutica [syn. Amblyopyrum muticum], and Triticum dicoccoides) and revealed their genomic contributions to the wheat B subgenome. Our results showed that centromeric regions were particularly conserved between Aegilops and Triticum and revealed 0.17 Gb of conserved blocks between Ae. speltoides and the B subgenome. We classified five groups of conserved and non-conserved genes between Aegilops and Triticum, revealing their biological characteristics, differentiation in gene expression patterns, and collinear relationships between Ae. speltoides and the wheat B subgenome. We also identified gene families that expanded in Ae. speltoides during its evolution and 789 genes specific to Ae. speltoides. These genes can serve as genetic resources for improvement of adaptability to biotic and abiotic stress. The newly constructed reference genome and large-scale resequencing data for Sitopsis species will provide a valuable genomic resource for wheat genetic improvement and genomic studies.

Single-gene FISH maps and major chromosomal rearrangements in Elymus sibiricus and E. nutans

BACKGROUND

Chromosomal variations have been revealed in both E. sibiricus and E. nutans, but chromosomal structural variations, such as intra-genome translocations and inversions, are still not recognized due to the cytological limitations of previous studies. Furthermore, the syntenic relationship between both species and wheat chromosomes remains unknown.

RESULTS

Fifty-nine single-gene fluorescence in situ hybridization (FISH) probes, including 22 single-gene probes previously mapped on wheat chromosomes and other newly developed probes from the cDNA of Elymus species, were used to characterize the chromosome homoeologous relationship and collinearity of both E. sibiricus and E. nutans with those of wheat. Eight species-specific chromosomal rearrangements (CRs) were exclusively identified in E. sibiricus, including five pericentric inversions in 1H, 2H, 3H, 6H and 2St; one possible pericentric inversion in 5St; one paracentric inversion in 4St; and one reciprocal 4H/6H translocation. Five species-specific CRs were identified in E. nutans, including one possible pericentric inversion in 2Y, three possible pericentric multiple-inversions in 1H, 2H and 4Y, and one reciprocal 4Y/5Y translocation. Polymorphic CRs were detected in three of the six materials in E. sibiricus, which were mainly represented by inter-genomic translocations. More polymorphic CRs were identified in E. nutans, including duplication and insertion, deletion, pericentric inversion, paracentric inversion, and intra- or inter-genomic translocation in different chromosomes.

CONCLUSIONS

The study first identified the cross-species homoeology and the syntenic relationship between E. sibiricus, E. nutans and wheat chromosomes. There are distinct different species-specific CRs between E. sibiricus and E. nutans, which may be due to their different polyploidy processes. The frequencies of intra-species polymorphic CRs in E. nutans were higher than that in E. sibiricus. To conclude, the results provide new insights into genome structure and evolution and will facilitate the utilization of germplasm diversity in both E. sibiricus and E. nutans.

Does one subgenome become dominant in the formation and evolution of a polyploid?

BACKGROUND

Polyploids are common in flowering plants and they tend to have more expanded ranges of distributions than their diploid progenitors. Possible mechanisms underlying polyploid success have been intensively investigated. Previous studies showed that polyploidy generates novel changes and that subgenomes in allopolyploid species often differ in gene number, gene expression levels and levels of epigenetic alteration. It is widely believed that such differences are the results of conflicts among the subgenomes. These differences have been treated by some as subgenome dominance, and it is claimed that the magnitude of subgenome dominance increases in polyploid evolution.

SCOPE

In addition to changes which occurred during evolution, differences between subgenomes of a polyploid species may also be affected by differences between the diploid donors and changes which occurred during polyploidization. The variable genome components in many plant species are extensive, which would result in exaggerated differences between a subgenome and its progenitor when a single genotype or a small number of genotypes are used to represent a polyploid or its donors. When artificially resynthesized polyploids are used as surrogates for newly formed genotypes which have not been exposed to evolutionary selection, differences between diploid genotypes available today and those involved in the formation of the natural polyploid genotypes must also be considered.

CONCLUSIONS

Contrary to the now widely held views that subgenome biases in polyploids are the results of conflicts among the subgenomes and that one of the parental subgenomes generally retains more genes which are more highly expressed, available results show that subgenome biases mainly reflect legacy from the progenitors and that they can be detected before the completion of polyploidization events. Further, there is no convincing evidence that the magnitudes of subgenome biases have significantly changed during evolution for any of the allopolyploid species assessed.

Novel Tetraploid Triticale (Einkorn Wheat × Rye)-A Source of Stem Rust Resistance

Among cereals, triticale (×Trititcoseale Wittmack ex A. Camus) represents a number of advantages such as high grain yield even in marginal environments, tolerance to drought, cold and acid soils, as well as lower production costs. Together with high biomass of grain and straw, triticale is also considered as an industrial energy crop. As an artificial hybrid, it has not evolved naturally, which is reflected in narrow genetic diversity causing a resistance collapse in recent years. Here, we describe a novel, synthetic tetraploid triticale, which was developed by the crossing of rye (Secale cereale L.) with einkorn wheat (Triticum monococcum spp. monococcum), which possess Sr35 stem rust resistance gene. Three subsequent generations of alloploids were obtained by chromosome doubling followed by self-pollination. The cytogenetic analyses revealed that the amphiploids possess a set of 28 chromosomes (14 of Am-genome and 14 of R-genome). The values of the most important yield-shaping traits for these tetraploid triticale form, including thousand-grain weight, plant height and stem length were higher compared to parental genotypes, as well as standard hexaploid triticale cultivars. This study shows that this tetraploid triticale genetic stock can be an interesting pre-breeding germplasm for triticale improvement or can be developed as a new alternative crop.

Evolutionary reconstruction, nomenclature and functional meta-analysis of the Kiwellin protein family

Crop diseases caused by pathogens critically affect global food security and plant ecology. Pathogens are well adapted to their host plants and have developed sophisticated mechanisms allowing successful colonization. Plants in turn have taken measures to counteract pathogen attacks resulting in an evolutionary arms race. Recent studies provided mechanistic insights into how two plant Kiwellin proteins from Zea mays mitigate the activity of the chorismate mutase Cmu1, a virulence factor secreted by the fungal pathogen Ustilago maydis during maize infection. Formerly identified as human allergens in kiwifruit, the biological function of Kiwellins is apparently linked to plant defense. We combined the analysis of proteome data with structural predictions to obtain a holistic overview of the Kiwellin protein family, that is subdivided into proteins with and without a N-terminal kissper domain. We found that Kiwellins are evolutionarily conserved in various plant species. At median five Kiwellin paralogs are encoded in each plant genome. Structural predictions revealed that Barwin-like proteins and Kiwellins cannot be discriminated purely at the sequence level. Our data shows that Kiwellins emerged in land plants (embryophyta) and are not present in fungi as suggested earlier. They evolved via three major duplication events that lead to clearly distinguishable subfamilies. We introduce a systematic Kiwellin nomenclature based on a detailed evolutionary reconstruction of this protein family. A meta-analysis of publicly available transcriptome data demonstrated that Kiwellins can be differentially regulated upon the interaction of plants with pathogens but also with symbionts. Furthermore, significant differences in Kiwellin expression levels dependent on tissues and cultivars were observed. In summary, our study sheds light on the evolution and regulation of a large protein family and provides a framework for a more detailed understanding of the molecular functions of Kiwellins.

Genome-wide association study reveals structural chromosome variations with phenotypic effects in wheat (Triticum aestivum L.)

Structural chromosome variations (SCVs) are large-scale genomic variations that can be detected by fluorescence in situ hybridization (FISH). SCVs have played important roles in the genome evolution of wheat (Triticum aestivum L.), but little is known about their genetic effects. In this study, a total of 543 wheat accessions from the Chinese wheat mini-core collection and the Shanxi Province wheat collection were used for chromosome analysis using oligonucleotide probe multiplex FISH. A total of 139 SCVs including translocations, pericentric inversions, presence/absence variations (PAVs), and copy number variations (CNVs) in heterochromatin were identified at 230 loci. The distribution frequency of SCVs varied between ecological regions and between landraces and modern cultivars. Structural analysis using SCVs as markers clearly divided the landraces and modern cultivars into different groups. There are very clear instances illustrating alien introgression and wide application of foreign germplasms improved the chromosome diversity of Chinese modern wheat cultivars. A genome-wide association study (GWAS) identified 29 SCVs associated with 12 phenotypic traits, and five (RT4AS•4AL-1DS/1DL•1DS-4AL, Mg2A-3, Mr3B-10, Mr7B-13, and Mr4A-7) of them were further validated using a doubled haploid population and advanced sib-lines, implying the potential value of these SCVs. Importantly, the number of favored SCVs that were associated with agronomic trait improvement was significantly higher in modern cultivars compared to landraces, indicating positive selection in wheat breeding. This study demonstrates the significant effects of SCVs during wheat breeding and provides an efficient method of mining favored SCVs in wheat and other crops.

Alteration of synonymous codon usage bias accompanies polyploidization in wheat

The diploidization of polyploid genomes is accompanied by genomic variation, including synonymous nucleotide substitutions that may lead to synonymous codon usage bias (SCUB). SCUB can mirror the evolutionary specialization of plants, but its effect on the formation of polyploidies is not well documented. We explored this issue here with hexaploid wheat and its progenitors. Synonymous codons (SCs) ending in either cytosine (NNC) or guanidine (NNG) were more frequent than those ending in either adenosine (NNA) or thymine (NNT), and the preference for NNC/G codons followed the increase in genome ploidy. The ratios between NNC/G and NNA/T codons gradually decreased in genes with more introns, and the difference in these ratios between wheat and its progenitors diminished with increasing ploidy. SCUB frequencies were heterogeneous among exons, and the bias preferred to NNA/T in more internal exons, especially for genes with more exons; while the preference did not appear to associate with ploidy. The SCUB alteration of the progenitors was different during the formation of hexaploid wheat, so that SCUB was the homogeneous among A, B and D subgenomes. DNA methylation-mediated conversion from cytosine to thymine weakened following the increase of genome ploidy, coinciding with the stronger bias for NNC/G SCs in the genome as a function of ploidy, suggesting that SCUB contribute to the epigenetic variation in hexaploid wheat. The patterns in SCUB mirrored the formation of hexaploid wheat, which provides new insight into genome shock-induced genetic variation during polyploidization. SCs representing non-neutral synonymous mutations can be used for genetic dissection and improvement of agricultural traits of wheat and other polyploidies.

Isolation and Sequencing of Chromosome Arm 7RS of Rye, Secale cereale

Rye (Secale cereale) is a climate-resilient cereal grown extensively as grain or forage crop in Northern and Eastern Europe. In addition to being an important crop, it has been used to improve wheat through introgression of genomic regions for improved yield and disease resistance. Understanding the genomic diversity of rye will assist both the improvement of this crop and facilitate the introgression of more valuable traits into wheat. Here, we isolated and sequenced the short arm of rye chromosome 7 (7RS) from Triticale 380SD using flow cytometry and compared it to the public Lo7 rye whole genome reference assembly. We identify 2747 Lo7 genes present on the isolated chromosome arm and two clusters containing seven and sixty-five genes that are present on Triticale 380SD 7RS, but absent from Lo7 7RS. We identified 29 genes that are not assigned to chromosomal locations in the Lo7 assembly but are present on Triticale 380SD 7RS, suggesting a chromosome arm location for these genes. Our study supports the Lo7 reference assembly and provides a repertoire of genes on Triticale 7RS.

Genome-Wide Association Study in Bread Wheat Identifies Genomic Regions Associated with Grain Yield and Quality under Contrasting Water Availability

The objectives of this study were to identify genetic loci in the bread wheat genome that would influence yield stability and quality under water stress, and to identify accessions that can be recommended for cultivation in dry and hot regions. We performed a genome-wide association study (GWAS) using a panel of 232 wheat accessions spanning diverse ecogeographic regions. Plants were evaluated in the Israeli Northern Negev, under two environments: water-limited (D; 250 mm) and well-watered (W; 450 mm) conditions; they were genotyped with ~71,500 SNPs derived from exome capture sequencing. Of the 14 phenotypic traits evaluated, 12 had significantly lower values under D compared to W conditions, while the values for two traits were higher under D. High heritability (H² = 0.5–0.9) was observed for grain yield, spike weight, number of grains per spike, peduncle length, and plant height. Days to heading and grain yield could be partitioned based on accession origins. GWAS identified 154 marker-trait associations (MTAs) for yield and quality-related traits, 82 under D and 72 under W, and identified potential candidate genes. We identified 24 accessions showing high and/or stable yields under D conditions that can be recommended for cultivation in regions under the threat of global climate change.

Retrieving a disrupted gene encoding phospholipase A for fibre enhancement in allotetraploid cultivated cotton

After polyploidization originated from one interspecific hybridization event in Gossypium, Gossypium barbadense evolved to produce extra-long staple fibres than Gossypium hirsutum (Upland cotton), which produces a higher fibre yield. The genomic diversity between G. barbadense and G. hirsutum thus provides a genetic basis for fibre trait variation. Recently, rapid accumulation of gene disruption or deleterious mutation was reported in allotetraploid cotton genomes, with unknown impacts on fibre traits. Here, we identified gene disruptions in allotetraploid G. hirsutum (18.14%) and G. barbadense (17.38%) through comparison with their presumed diploid progenitors. Relative to conserved genes, these disrupted genes exhibited faster evolution rate, lower expression level and altered gene co-expression networks. Within a module regulating fibre elongation, a hub gene experienced gene disruption in G. hirsutum after polyploidization, with a 2-bp deletion in the coding region of GhNPLA1D introducing early termination of translation. This deletion was observed in all of the 34 G. hirsutum landraces and 36 G. hirsutum cultivars, but not in 96% of 57 G. barbadense accessions. Retrieving the disrupted gene GhNPLA1D using its homoeolog GhNPLA1A achieved longer fibre length in G. hirsutum. Further enzyme activity and lipids analysis confirmed that GhNPLA1A encodes a typical phospholipase A and promotes cotton fibre elongation via elevating intracellular levels of linolenic acid and 34:3 phosphatidylinositol. Our work opens a strategy for identifying disrupted genes and retrieving their functions in ways that can provide valuable resources for accelerating fibre trait enhancement in cotton breeding.

Effects of homoeologous exchange on gene expression and alternative splicing in a newly formed allotetraploid wheat

Homoeologous exchange (HE) is a major mechanism generating post-polyploidization genetic variation with important evolutionary consequences. However, the direct impacts of HE on gene expression and transcript diversity in allopolyploids without the intertwined evolutionary processes remain to be fully understood. Here, we analyzed high-throughput RNA-seq data of young leaves from plant groups of a synthetic allotetraploid wheat (AADD), which contained variable numbers of HEs. We aimed to investigate if and to which extent HE directly impacts gene expression and alternative splicing (AS). We found that HE impacts expression of genes located within HE regions primarily via a cis-acting dosage effect, which led to significant changes in the total expression level of homoeologous gene pairs, especially for homoeologs whose original expression was biased. In parallel, HE also influences expression of a large number of genes residing in non-HE regions by trans-regulation leading to convergent expression of homoeologs. Intriguingly, when taking the original relative homoeolog expression states into account, homoeolog pairs under trans-effect are more prone to manifesting a convergent response to the HEs whereas those under cis-regulation tended to show further exacerbated subgenome-biased expression. Moreover, HE-induced quantitative, largely individual-specific, changes of AS events were detected. Similar to homoeologous expression, homoeo-AS events under trans-effect were more responsive to HE. HE therefore exerts multifaceted immediate effects on gene expression and, to a less extent, on individualized transcript diversity in nascent allopolyploidy.

Challenges and opportunities for proteomics and the improvement of bread wheat quality

Wheat remains a critical global food source, pressured by climate change and the need to maximize yield, improve processing and nutritional quality and ensure safety. An enormous amount of research has been conducted to understand gluten protein composition and structure in relation to end-use quality, yet progress has become stagnant. This is mainly due to the need and inability to biochemically characterize the intact functional glutenin polymer in order to correlate to quality, necessitating reduction to monomeric subunits and a loss of contextual information. While some individual gluten proteins might have a positive or negative influence on gluten quality, it is the sum total of these proteins, their relative and absolute expression, their sub-cellular trafficking, the amount and size of glutenin polymers, and ratios between gluten protein classes that define viscoelasticity of gluten. The sub-cellular trafficking of gluten proteins during seed maturation is still not completely clear and there is evidence of dual pathways and therefore different destinations for proteins, either constitutively or temporally. The trafficking of proteins is also unclear in endosperm cells as they undergo programmed cell death; Golgi disappear around 12 DPA but protein filling continues at least to 25 DPA. Modulation of the timing of cellular events will invariably affect protein deposition and therefore gluten strength and function. Existing and emerging proteomics technologies such as proteoform profiling and top-down proteomics offer new tools to study gluten protein composition as a whole system and identify compositional patterns that can modify gluten structure with improved functionality.

Transcription factor retention through multiple polyploidization steps in wheat

Whole-genome duplication is widespread in plant evolutionary history and is followed by nonrandom gene loss to return to a diploid state. Across multiple angiosperm species, the retained genes tend to be dosage-sensitive regulatory genes such as transcription factors, yet data for younger polyploid species is sparse. Here, we analyzed the retention, expression, and genetic variation in transcription factors in the recent allohexaploid bread wheat (Triticum aestivum L.). By comparing diploid, tetraploid, and hexaploid wheat, we found that, following each of two hybridization and whole-genome duplication events, the proportion of transcription factors in the genome increased. Transcription factors were preferentially retained over other genes as homoeologous groups in tetraploid and hexaploid wheat. Across cultivars, transcription factor homoeologs contained fewer deleterious missense mutations than nontranscription factors, suggesting that transcription factors are maintained as three functional homoeologs in hexaploid wheat populations. Transcription factor homoeologs were more strongly coexpressed than nontranscription factors, indicating conservation of function between homoeologs. We found that the B3, MADS-M-type, and NAC transcription factor families were less likely to have three homoeologs present than other families, which was associated with low expression levels and high levels of tandem duplication. Together, our results show that transcription factors are preferentially retained in polyploid wheat genomes although there is variation between families. Knocking out one transcription factor homoeolog to alter gene dosage, using TILLING or CRISPR, could generate new phenotypes for wheat breeding.

Utility of Triti-Map for bulk-segregated mapping of causal genes and regulatory elements in Triticeae

Triticeae species, including wheat, barley, and rye, are critical for global food security. Mapping agronomically important genes is crucial for elucidating molecular mechanisms and improving crops. However, Triticeae includes many wild relatives with desirable agronomic traits, and frequent introgressions occurred during Triticeae evolution and domestication. Thus, Triticeae genomes are generally large and complex, making the localization of genes or functional elements that control agronomic traits challenging. Here, we developed Triti-Map, which contains a suite of user-friendly computational packages specifically designed and optimized to overcome the obstacles of gene mapping in Triticeae, as well as a web interface integrating multi-omics data from Triticeae for the efficient mining of genes or functional elements that control particular traits. The Triti-Map pipeline accepts both DNA and RNA bulk-segregated sequencing data as well as traditional QTL data as inputs for locating genes and elucidating their functions. We illustrate the usage of Triti-Map with a combination of bulk-segregated ChIP-seq data to detect a wheat disease-resistance gene with its promoter sequence that is absent from the reference genome and clarify its evolutionary process. We hope that Triti-Map will facilitate gene isolation and accelerate Triticeae breeding.

Evolution and origin of bread wheat

Bread wheat (Triticum aestivum, genome BBAADD) is a young hexaploid species formed only 8,500-9,000 years ago through hybridization between a domesticated free-threshing tetraploid progenitor, genome BBAA, and Aegilops tauschii, the diploid donor of the D subgenome. Very soon after its formation, it spread globally from its cradle in the fertile crescent into new habitats and climates, to become a staple food of humanity. This extraordinary global expansion was probably enabled by allopolyploidy that accelerated genetic novelty through the acquisition of new traits, new intergenomic interactions, and buffering of mutations, and by the attractiveness of bread wheat's large, tasty, and nutritious grain with high baking quality. New genome sequences suggest that the elusive donor of the B subgenome is a distinct (unknown or extinct) species rather than a mosaic genome. We discuss the origin of the diploid and tetraploid progenitors of bread wheat and the conflicting genetic and archaeological evidence on where it was formed and which species was its free-threshing tetraploid progenitor. Wheat experienced many environmental changes throughout its evolution, therefore, while it might adapt to current climatic changes, efforts are needed to better use and conserve the vast gene pool of wheat biodiversity on which our food security depends.

Transfer of the ph1b Deletion Chromosome 5B From Chinese Spring Wheat Into a Winter Wheat Line and Induction of Chromosome Rearrangements in Wheat-Aegilops biuncialis Hybrids

Effective utilization of genetic diversity in wild relatives to improve wheat requires recombination between wheat and alien chromosomes. However, this is suppressed by the Pairing homoeologous gene, Ph1, on the long arm of wheat chromosome 5B. A deletion mutant of the Ph1 locus (ph1b) has been used widely to induce homoeologous recombination in wheat × alien hybrids. However, the original ph1b mutation, developed in Chinese Spring (CS) background has poor agronomic performance. Hence, alien introgression lines are first backcrossed with adapted wheat genotypes and after this step, alien chromosome segments are introduced into breeding lines. In this work, the ph1b mutation was transferred from two CSph1b mutants into winter wheat line Mv9kr1. Homozygous genotypes Mv9kr1 ph1b/ph1b exhibited improved plant and spike morphology compared to Chinese Spring. Flow cytometric chromosome analysis confirmed reduced DNA content of the mutant 5B chromosome in both wheat genotype relative to the wild type chromosome. The ph1b mutation in the Mv9kr1 genotype allowed wheat-alien chromosome pairing in meiosis of Mv9kr1ph1b_K × Aegilops biuncialis F₁ hybrids, predominantly with the M^b-genome chromosomes of Aegilops relative to those of the U^b genome. High frequency of wheat-Aegilops chromosome interactions resulted in rearranged chromosomes identified in the new Mv9kr1ph1b × Ae. Biuncialis amphiploids, making these lines valuable sources for alien introgressions. The new Mv9kr1ph1b mutant genotype is a unique resource to support alien introgression breeding of hexaploid wheat.

Comprehensive molecular evaluation of the histone methyltransferase gene family and their important roles in two-line hybrid wheat

BACKGROUND

Histone methylation usually plays important roles in plant development through post-translational regulation and may provide a new visual field for heterosis. The histone methyltransferase gene family has been identified in various plants, but its members and functions in hybrid wheat related in heterosis is poorly studied.

RESULTS

In this study, 175 histone methyltransferase (HMT) genes were identified in wheat, including 152 histone lysine methyltransferase (HKMT) genes and 23 protein arginine methyltransferase (PRMT) genes. Gene structure analysis, physicochemical properties and subcellular localization predictions of the proteins, exhibited the adequate complexity of this gene family. As an allohexaploid species, the number of the genes (seven HKMTs orthologous groups and four PRMTs orthologous groups) in wheat were about three times than those in diploids and showed certain degrees of conservation, while only a small number of subfamilies such as ASH-like and Su-(var) subfamilies have expanded their members. Transcriptome analysis showed that HMT genes were mainly expressed in the reproductive organs. Expression analysis showed that some TaHMT genes with different trends in various hybrid combinations may be regulated by lncRNAs with similar expression trends. Pearson correlation analysis of the expression of TaHMT genes and two yield traits indicated that four DEGs may participate in the yield heterosis of two-line hybrid wheat. ChIP-qPCR results showed that the histone modifications (H3K4me3, H3K36me3 and H3K9ac) enriched in promoter regions of three TaCCA1 genes which are homologous to Arabidopsis heterosis-related CCA1/LHY genes. The higher expression levels of TaCCA1 in F₁ than its parents are positive with these histone modifications. These results showed that histone modifications may play important roles in wheat heterosis.

CONCLUSIONS

Our study identified characteristics of the histone methyltransferase gene family and enhances the understanding of the evolution and function of these members in allohexaploid wheat. The causes of heterosis of two-line hybrid wheat were partially explained from the perspective of histone modifications.

Genetic Diversity of Rhanterium eppaposum Oliv. Populations in Kuwait as Revealed by GBS

Natural populations of Rhanterium eppaposum Oliv. (Arfaj), a perennial forage shrub, have depleted due to unethical human interventions and climate change in Kuwait. Therefore, there is an urgent need to conserve this native plant through the assessment of its genetic diversity and population structure. Genotyping by sequencing (GBS) has recently emerged as a powerful tool for the molecular diversity analysis of higher plants without prior knowledge of their genome. This study represents the first effort in using GBS to discover genome-wide single nucleotide polymorphisms (SNPs) of local Rhanterium plants to assess the genetic diversity present in landraces collected from six different locations in Kuwait. The study generated a novel set of 11,231 single nucleotide polymorphisms (SNPs) and indels (insertions and deletions) in 98 genotypes of Rhanterium. The analysis of molecular variance (AMOVA) revealed ~1.5% variation residing among the six populations, ~5% among the individuals within the population and 93% variation present within the populations (F_ST = 0.029; p = 0.0). Bayesian and UPGMA analyses identified two admixed clusters of the tested samples; however, the principal coordinates analysis returned the complete population as a single group. Mantel's test returned a very weak correlation coefficient of r² = 0.101 (p = 0.00) between the geographic and genetic distance. These findings are useful for the native species to formulate conservation strategies for its sustainable management and desert rehabilitation.

Spatial Features and Functional Implications of Plant 3D Genome Organization

The advent of high-throughput sequencing-based methods for chromatin conformation, accessibility, and immunoprecipitation assays has been a turning point in 3D genomics. Altogether, these new tools have been pushing upward the interpretation of pioneer cytogenetic evidence for a higher order in chromatin packing. Here, we review the latest development in our understanding of plant spatial genome structures and different levels of organization and discuss their functional implications. Then, we spotlight the complexity of organellar (i.e., mitochondria and plastids) genomes and discuss their 3D packing into nucleoids. Finally, we propose unaddressed research axes to investigate functional links between chromatin-like dynamics and transcriptional regulation within organellar nucleoids.

Dissection of Structural Reorganization of Wheat 5B Chromosome Associated With Interspecies Recombination Suppression

Chromosomal rearrangements that lead to recombination suppression can have a significant impact on speciation, and they are also important for breeding. The regions of recombination suppression in wheat chromosome 5B were identified based on comparisons of the 5B map of a cross between the Chinese Spring (CS) variety of hexaploid wheat and CS-5Bdic (genotype CS with 5B substituted with its homologue from tetraploid Triticum dicoccoides) with several 5B maps of tetraploid and hexaploid wheat. In total, two regions were selected in which recombination suppression occurred in cross CS × CS-5Bdic when compared with other maps: one on the short arm, 5BS_RS, limited by markers BS00009810/BS00022336, and the second on the long arm, 5BL_RS, between markers Ra_c10633_2155 and BS00087043. The regions marked as 5BS_RS and 5BL_RS, with lengths of 5 Mb and 3.6 Mb, respectively, were mined from the 5B pseudomolecule of CS and compared to the homoeologous regions (7.6 and 3.8 Mb, respectively) of the 5B pseudomolecule of Zavitan (T. dicoccoides). It was shown that, in the case of 5BS_RS, the local heterochromatin islands determined by the satellite DNA (119.2) and transposable element arrays, as well as the dissimilarity caused by large insertions/deletions (chromosome rearrangements) between 5BSs aestivum/dicoccoides, are likely the key determinants of recombination suppression in the region. Two major and two minor segments with significant loss of similarity were recognized within the 5BL_RS region. It was shown that the loss of similarity, which can lead to suppression of recombination in the 5BL_RS region, is caused by chromosomal rearrangements, driven by the activity of mobile genetic elements (both DNA transposons and long terminal repeat retrotransposons) and their divergence during evolution. It was noted that the regions marked as 5BS_RS and 5BL_RS are associated with chromosomal rearrangements identified earlier by С-banding analysis of intraspecific polymorphism of tetraploid emmer wheat. The revealed divergence in 5BS_RS and 5BL_RS may be a consequence of interspecific hybridization, plant genetic adaptation, or both.

Progenitor species hold untapped diversity for potential climate-responsive traits for use in wheat breeding and crop improvement

Climate change will have numerous impacts on crop production worldwide necessitating a broadening of the germplasm base required to source and incorporate novel traits. Major variation exists in crop progenitor species for seasonal adaptation, photosynthetic characteristics, and root system architecture. Wheat is crucial for securing future food and nutrition security and its evolutionary history and progenitor diversity offer opportunities to mine favourable functional variation in the primary gene pool. Here we provide a review of the status of characterisation of wheat progenitor variation and the potential to use this knowledge to inform the use of variation in other cereal crops. Although significant knowledge of progenitor variation has been generated, we make recommendations for further work required to systematically characterise underlying genetics and physiological mechanisms and propose steps for effective use in breeding. This will enable targeted exploitation of useful variation, supported by the growing portfolio of genomics and accelerated breeding approaches. The knowledge and approaches generated are also likely to be useful across wider crop improvement.

Rubiscosome gene expression is balanced across the hexaploid wheat genome

Functional and active Rubisco is essential for CO₂ fixation and is a primary target for engineering approaches to increasing crop yields. However, the assembly and maintenance of active Rubisco are dependent on the coordinated biosynthesis of at least 11 nuclear-encoded proteins, termed the 'Rubiscosome'. Using publicly available gene expression data for wheat (Triticum aestivum L.), we show that the expression of Rubiscosome genes is balanced across the three closely related subgenomes that form the allohexaploid genome. Each subgenome contains a near complete set of homoeologous genes and contributes equally to overall expression, both under optimal and under heat stress conditions. The expression of the wheat thermo-tolerant Rubisco activase isoform 1β increases under heat stress and remains balanced across the subgenomes, albeit with a slight shift towards greater contribution from the D subgenome. The findings show that the gene copies in all three subgenomes need to be accounted for when designing strategies for crop improvement.

Genetic characterization and curation of diploid A-genome wheat species

A-genome diploid wheats represent the earliest domesticated and cultivated wheat species in the Fertile Crescent and include the donor of the wheat A sub-genome. The A-genome species encompass the cultivated einkorn (Triticum monococcum L. subsp. monococcum), wild einkorn (T. monococcum L. subsp. aegilopoides (Link) Thell.), and Triticum urartu. We evaluated the collection of 930 accessions in the Wheat Genetics Resource Center (WGRC) using genotyping by sequencing and identified 13,860 curated single-nucleotide polymorphisms. Genomic analysis detected misclassified and genetically identical (>99%) accessions, with most of the identical accessions originating from the same or nearby locations. About 56% (n = 520) of the WGRC A-genome species collections were genetically identical, supporting the need for genomic characterization for effective curation and maintenance of these collections. Population structure analysis confirmed the morphology-based classifications of the accessions and reflected the species geographic distributions. We also showed that T. urartu is the closest A-genome diploid to the A-subgenome in common wheat (Triticum aestivum L.) through phylogenetic analysis. Population analysis within the wild einkorn group showed three genetically distinct clusters, which corresponded with wild einkorn races α, β, and γ described previously. The T. monococcum genome-wide FST scan identified candidate genomic regions harboring a domestication selection signature at the Non-brittle rachis 1 (Btr1) locus on the short arm of chromosome 3Am at ∼70 Mb. We established an A-genome core set (79 accessions) based on allelic diversity, geographical distribution, and available phenotypic data. The individual species core set maintained at least 79% of allelic variants in the A-genome collection and constituted a valuable genetic resource to improve wheat and domesticated einkorn in breeding programs.

Tracing the Evolution of the Angiosperm Genome from the Cytogenetic Point of View

Cytogenetics constitutes a branch of genetics that is focused on the cellular components, especially chromosomes, in relation to heredity and genome structure, function and evolution. The use of modern cytogenetic approaches and the latest microscopes with image acquisition and processing systems enables the simultaneous two- or three-dimensional, multicolour visualisation of both single-copy and highly-repetitive sequences in the plant genome. The data that is gathered using the cytogenetic methods in the phylogenetic background enable tracing the evolution of the plant genome that involve changes in: (i) genome sizes; (ii) chromosome numbers and morphology; (iii) the content of repetitive sequences and (iv) ploidy level. Modern cytogenetic approaches such as FISH using chromosome- and genome-specific probes have been widely used in studies of the evolution of diploids and the consequences of polyploidy. Nowadays, modern cytogenetics complements analyses in other fields of cell biology and constitutes the linkage between genetics, molecular biology and genomics.

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.

Transposable Element Populations Shed Light on the Evolutionary History of Wheat and the Complex Co-Evolution of Autonomous and Non-Autonomous Retrotransposons

Wheat has one of the largest and most repetitive genomes among major crop plants, containing over 85% transposable elements (TEs). TEs populate genomes much in the way that individuals populate ecosystems, diversifying into different lineages, sub-families and sub-populations. The recent availability of high-quality, chromosome-scale genome sequences from ten wheat lines enables a detailed analysis how TEs evolved in allohexaploid wheat, its diploids progenitors, and in various chromosomal haplotype segments. LTR retrotransposon families evolved into distinct sub-populations and sub-families that were active in waves lasting several hundred thousand years. Furthermore, It is shown that different retrotransposon sub-families were active in the three wheat sub-genomes, making them useful markers to study and date polyploidization events and chromosomal rearrangements. Additionally, haplotype-specific TE sub-families are used to characterize chromosomal introgressions in different wheat lines. Additionally, populations of non-autonomous TEs co-evolved over millions of years with their autonomous partners, leading to complex systems with multiple types of autonomous, semi-autonomous and non-autonomous elements. Phylogenetic and TE population analyses revealed the relationships between non-autonomous elements and their mobilizing autonomous partners. TE population analysis provided insights into genome evolution of allohexaploid wheat and genetic diversity of species, and may have implication for future crop breeding.

Wheat in vivo RNA structure landscape reveals a prevalent role of RNA structure in modulating translational subgenome expression asymmetry

BACKGROUND

Polyploidy, especially allopolyploidy, which entails merging divergent genomes via hybridization and whole-genome duplication (WGD), is a major route to speciation in plants. The duplication among the parental genomes (subgenomes) often leads to one subgenome becoming dominant over the other(s), resulting in subgenome asymmetry in gene content and expression. Polyploid wheats are allopolyploids with most genes present in two (tetraploid) or three (hexaploid) functional copies, which commonly show subgenome expression asymmetry. It is unknown whether a similar subgenome asymmetry exists during translation. We aim to address this key biological question and explore the major contributing factors to subgenome translation asymmetry.

RESULTS

Here, we obtain the first tetraploid wheat translatome and reveal that subgenome expression asymmetry exists at the translational level. We further perform in vivo RNA structure profiling to obtain the wheat RNA structure landscape and find that mRNA structure has a strong impact on translation, independent of GC content. We discover a previously uncharacterized contribution of RNA structure in subgenome translation asymmetry. We identify 3564 single-nucleotide variations (SNVs) across the transcriptomes between the two tetraploid wheat subgenomes, which induce large RNA structure disparities. These SNVs are highly conserved within durum wheat cultivars but are divergent in both domesticated and wild emmer wheat.

CONCLUSIONS

We successfully determine both the translatome and in vivo RNA structurome in tetraploid wheat. We reveal that RNA structure serves as an important modulator of translational subgenome expression asymmetry in polyploids. Our work provides a new perspective for molecular breeding of major polyploid crops.

Karyotype Reorganization in Wheat-Rye Hybrids Obtained via Unreduced Gametes: Is There a Limit to the Chromosome Number in Triticale?

To date, few data have been accumulated on the contribution of meiotic restitution to the formation of Triticum aestivum hybrid karyotypes. In this study, based on FISH and C-banding, karyotype reorganization was observed in three groups of F₅ wheat–rye hybrids 1R(1A) × R. Aberrations, including aneuploidy, telocentrics, and Robertsonian translocations, were detected in all groups. Some of the Group 1 plants and all of the Group 2 plants only had a 4R4R pair (in addition to 1R1R), which was either added or substituted for its homeolog in ABD subgenomes. In about 82% of meiocytes, 4R4R formed bivalents, which indicates its competitiveness. The rest of the Group 1 plants had 2R and 7R chromosomes in addition to 1R1R. Group 3 retained all their rye chromosomes, with a small aneuploidy on the wheat chromosomes. A feature of the meiosis in the Group 3 plants was asynchronous cell division and omission of the second division. Diploid gametes did not form because of the significant disturbances during gametogenesis. As a result, the frequency of occurrence of the formed dyads was negatively correlated (r = −0.73) with the seed sets. Thus, meiotic restitution in the 8n triticale does not contribute to fertility or increased ploidy in subsequent generations.