Ancient hybridizations among the ancestral genomes of bread wheat

Heteroalleles in Common Wheat: Multiple Differences between Allelic Variants of the Gli-B1 Locus

The Gli-B1-encoded γ-gliadins and non-coding γ-gliadin DNA sequences for 15 different alleles of common wheat have been compared using seven tests: electrophoretic mobility (EM) and molecular weight (MW) of the encoded major γ-gliadin, restriction fragment length polymorphism patterns (RFLPs) (three different markers), Gli-B1-γ-gliadin-pseudogene known SNP markers (Single nucleotide polymorphisms) and sequencing the pseudogene GAG56B. It was discovered that encoded γ-gliadins, with contrasting EM, had similar MWs. However, seven allelic variants (designated from I to VII) differed among them in the other six tests: I (alleles Gli-B1i, k, m, o), II (Gli-B1n, q, s), III (Gli-B1b), IV (Gli-B1e, f, g), V (Gli-B1h), VI (Gli-B1d) and VII (Gli-B1a). Allele Gli-B1c (variant VIII) was identical to the alleles from group IV in four of the tests. Some tests might show a fine difference between alleles belonging to the same variant. Our results attest in favor of the independent origin of at least seven variants at the Gli-B1 locus that might originate from deeply diverged genotypes of the donor(s) of the B genome in hexaploid wheat and therefore might be called “heteroallelic”. The donor’s particularities at the Gli-B1 locus might be conserved since that time and decisively contribute to the current high genetic diversity of common wheat.

The unique disarticulation layer formed in the rachis of Aegilops longissima probably results from the spatial co-expression of Btr1 and Btr2

BACKGROUND AND AIMS

The brittle rachis trait is a feature of many wild grasses, particularly within the tribe Triticeae. Wild Hordeum and Triticum species form a disarticulation layer above the rachis node, resulting in the production of wedge-type dispersal units. In Aegilops longissima, only one or two of the nodes in the central portion of its rachis are brittle. In Triticeae species, the formation of a disarticulation layer above the rachis node requires the co-transcription of the two dominant and complementary genes Btr1 and Btr2. This study aims to establish whether homologues of Btr1 and/or Btr2 underlie the unusual brittle rachis phenotype observed in Ae. longissima.

METHODS

Scanning electron microscopy was used to examine the disarticulation surfaces. Quantitative RT-PCR and RNA in situ hybridization experiments were used to identify gene expression in the immature inflorescence.

KEY RESULTS

Analysis based on scanning electron microscopy was able to demonstrate that the disarticulation surfaces formed in the Ae. longissima rachis are morphologically indistinguishable from those formed in the rachises of wild Hordeum and Triticum species. RNA in situ hybridization showed that in the immature Ae. longissima inflorescence, the intensity of Btr1 transcription varied from high at the rachis base to low at its apex, while that of Btr2 was limited to the nodes in the central to distal portion of the rachis.

CONCLUSIONS

The disarticulation pattern shown by Ae. longissima results from the limitation of Btr1 and Btr2 co-expression to nodes lying in the centre of the rachis.

Molecular characterization revealed the role of catalases under abiotic and arsenic stress in bread wheat (Triticum aestivum L.)

Catalases are crucial antioxidant enzymes that reduce the excessive level of H₂O₂ caused by various environmental stresses and metal toxicity and hence protect the plant cells. In this study, a total of ten TaCAT genes, forming three homeologous groups, were identified in the genome of bread wheat (Triticum aestivum L.) and named as per the wheat gene symbolization guidelines. The identified catalases were characterized for various structural and physicochemical features. The proximal active-site (F(D/A)RERIPERVVHAKGASA) and heme-ligand (R(I/V)F(S/A)Y(A/S)DTQ) signature motifs, catalytic residues and peroxisomal targeting peptides were found conserved. Phylogenetic analysis clustered TaCATs into three classes, which showed conserved functional specialization based on their tissue specific expression. Modulated spatio-temporal expression of various TaCAT genes and alteration in total catalase enzyme activity during heat, drought, salt and arsenic (As^III and As^V) treatment suggested their roles in abiotic stress response and arsenic tolerance. Molecular cloning and overexpression of TaCAT3-B gene in Escherichia coli showed tolerance against heat, drought, salt and varied concentrations of As^III and As^V treatments. The results further confirmed their role in stress tolerance and recommended that these genes can be used in future stress management strategies for the development of abiotic and arsenic stress resistant transgenic crops.

Generation of Doubled Haploid Wheat-Triticum urartu Introgression Lines and Their Characterisation Using Chromosome-Specific KASP Markers

Wheat is one of the most important food and protein sources in the world and although, in recent years wheat breeders have achieved yield gains, they are not sufficient to meet the demands of an ever-growing population. Development of high yielding wheat varieties, resilient to abiotic and biotic stress resulting from climate change, has been limited by wheat's narrow genetic base. In contrast to wheat, the wild relatives of wheat provide a vast reservoir of genetic variation for most, if not all, agronomic traits. Previous studies by the authors have shown the transfer of genetic variation from T. urartu into bread wheat. However, before the introgression lines can be exploited for trait analysis, they are required to have stable transmission of the introgressions to the next generation. In this work, we describe the generation of 86 doubled haploid (DH) wheat-T. urartu introgression lines that carry homozygous introgressions which are stably inherited. The DH lines were characterised using the Axiom^® Wheat Relative Genotyping Array and 151 KASP markers to identify 65 unique T. urartu introgressions in a bread wheat background. DH production has helped accelerate the breeding process and facilitated the early release of homozygous wheat-T. urartu introgression lines. Together with the KASP markers, this valuable resource could greatly advance identification of beneficial alleles that can be used in wheat improvement.

Variation in Plastome Sizes Accompanied by Evolutionary History in Monogenomic Triticeae (Poaceae: Triticeae)

To investigate the pattern of chloroplast genome variation in Triticeae, we comprehensively analyzed the indels in protein-coding genes and intergenic sequence, gene loss/pseudonization, intron variation, expansion/contraction in inverted repeat regions, and the relationship between sequence characteristics and chloroplast genome size in 34 monogenomic Triticeae plants. Ancestral genome reconstruction suggests that major length variations occurred in four-stem branches of monogenomic Triticeae followed by independent changes in each genus. It was shown that the chloroplast genome sizes of monogenomic Triticeae were highly variable. The chloroplast genome of Pseudoroegneria, Dasypyrum, Lophopyrum, Thinopyrum, Eremopyrum, Agropyron, Australopyrum, and Henradia in Triticeae had evolved toward size reduction largely because of pseudogenes elimination events and length deletion fragments in intergenic. The Aegilops/Triticum complex, Taeniatherum, Secale, Crithopsis, Herteranthelium, and Hordeum in Triticeae had a larger chloroplast genome size. The large size variation in major lineages and their subclades are most likely consequences of adaptive processes since these variations were significantly correlated with divergence time and historical climatic changes. We also found that several intergenic regions, such as petN-trnC and psbE-petL containing unique genetic information, which can be used as important tools to identify the maternal relationship among Triticeae species. Our results contribute to the novel knowledge of plastid genome evolution in Triticeae.

A high-quality genome assembly highlights rye genomic characteristics and agronomically important genes

Rye is a valuable food and forage crop, an important genetic resource for wheat and triticale improvement and an indispensable material for efficient comparative genomic studies in grasses. Here, we sequenced the genome of Weining rye, an elite Chinese rye variety. The assembled contigs (7.74 Gb) accounted for 98.47% of the estimated genome size (7.86 Gb), with 93.67% of the contigs (7.25 Gb) assigned to seven chromosomes. Repetitive elements constituted 90.31% of the assembled genome. Compared to previously sequenced Triticeae genomes, Daniela, Sumaya and Sumana retrotransposons showed strong expansion in rye. Further analyses of the Weining assembly shed new light on genome-wide gene duplications and their impact on starch biosynthesis genes, physical organization of complex prolamin loci, gene expression features underlying early heading trait and putative domestication-associated chromosomal regions and loci in rye. This genome sequence promises to accelerate genomic and breeding studies in rye and related cereal crops.

Genome-Wide Identification, Evolution, and Expression Analysis of LBD Transcription Factor Family in Bread Wheat (Triticum aestivum L.)

The lateral organ boundaries domain (LBD) genes, as the plant-specific transcription factor family, play a crucial role in controlling plant architecture and stress tolerance. Although it has been thoroughly characterized in many species, the LBD family was not well studied in wheat. Here, the wheat LBD family was systematically investigated through an in silico genome-wide search method. A total of 90 wheat LBD genes (TaLBDs) were identified, which were classified into class I containing seven subfamilies, and class II containing two subfamilies. Exon-intron structure, conserved protein motif, and cis-regulatory elements analysis showed that the members in the same subfamily shared similar gene structure organizations, supporting the classification. Furthermore, the expression patterns of these TaLBDs in different types of tissues and under diverse stresses were identified through public RNA-seq data analysis, and the regulation networks of TaLBDs involved were predicted. Finally, the expression levels of 12 TaLBDs were validated by quantitative PCR (qPCR) analysis and the homoeologous genes showed differential expression. Additionally, the genetic diversity of TaLBDs in the landrace population showed slightly higher than that of the genetically improved germplasm population while obvious asymmetry at the subgenome level. This study not only provided the potential targets for further functional analysis but also contributed to better understand the roles of LBD genes in regulating development and stress tolerance in wheat and beyond.

Molecular Characterization Revealed the Role of Thaumatin-Like Proteins of Bread Wheat in Stress Response

Thaumatin-like proteins (TLPs) are related to pathogenesis-related-5 (PR-5) family and involved in stress response. Herein, a total of 93 TLP genes were identified in the genome of Triticum aestivum. Further, we identified 26, 27, 39, and 37 TLP genes in the Brachypodium distachyon, Oryza sativa, Sorghum bicolor, and Zea mays genomes for comparative characterization, respectively. They could be grouped into small and long TLPs with conserved thaumatin signature motif. Tightly clustered genes exhibited conserved gene and protein structure. The physicochemical analyses suggested significant differences between small and long TLPs. Evolutionary analyses suggested the role of duplication events and purifying selection in the expansion of the TLP gene family. Expression analyses revealed the possible roles of TLPs in plant development and abiotic and fungal stress response. Recombinant expression of TaTLP2-B in Saccharomyces cerevisiae provided significant tolerance against cold, heat, osmotic, and salt stresses. The results depicted the importance of TLPs in cereal crops that would be highly useful in future crop improvement programs.

Genome-Wide Identification and Evolution of Receptor-Like Kinases (RLKs) and Receptor like Proteins (RLPs) in Brassica juncea

Brassica juncea, an allotetraploid species, is an important germplasm resource for canola improvement, due to its many beneficial agronomic traits, such as heat and drought tolerance and blackleg resistance. Receptor-like kinase (RLK) and receptor-like protein (RLP) genes are two types of resistance gene analogues (RGA) that play important roles in plant innate immunity, stress response and various development processes. In this study, genome wide analysis of RLKs and RLPs is performed in B. juncea. In total, 493 RLKs (LysM-RLKs and LRR-RLKs) and 228 RLPs (LysM-RLPs and LRR-RLPs) are identified in the genome of B. juncea, using RGAugury. Only 13.54% RLKs and 11.79% RLPs are observed to be grouped within gene clusters. The majority of RLKs (90.17%) and RLPs (52.83%) are identified as duplicates, indicating that gene duplications significantly contribute to the expansion of RLK and RLP families. Comparative analysis between B. juncea and its progenitor species, B. rapa and B. nigra, indicate that 83.62% RLKs and 41.98% RLPs are conserved in B. juncea, and RLPs are likely to have a faster evolution than RLKs. This study provides a valuable resource for the identification and characterisation of candidate RLK and RLP genes.

Coevolution in Hybrid Genomes: Nuclear-Encoded Rubisco Small Subunits and Their Plastid-Targeting Translocons Accompanying Sequential Allopolyploidy Events in Triticum

The Triticum/Aegilops complex includes hybrid species resulting from homoploid hybrid speciation and allopolyploid speciation. Sequential allotetra- and allohexaploidy events presumably result in two challenges for the hybrids, which involve 1) cytonuclear stoichiometric disruptions caused by combining two diverged nuclear genomes with the maternal inheritance of the cytoplasmic organellar donor; and 2) incompatibility of chimeric protein complexes with diverged subunits from nuclear and cytoplasmic genomes. Here, we describe coevolution of nuclear rbcS genes encoding the small subunits of Rubisco (ribulose 1,5-bisphosphate carboxylase/oxygenase) and nuclear genes encoding plastid translocons, which mediate recognition and translocation of nuclear-encoded proteins into plastids, in allopolyploid wheat species. We demonstrate that intergenomic paternal-to-maternal gene conversion specifically occurred in the genic region of the homoeologous rbcS3 gene from the D-genome progenitor of wheat (abbreviated as rbcS3D) such that it encodes a maternal-like or B-subgenome-like SSU3D transit peptide in allohexaploid wheat but not in allotetraploid wheat. Divergent and limited interaction between SSU3D and the D-subgenomic TOC90D translocon subunit is implicated to underpin SSU3D targeting into the chloroplast of hexaploid wheat. This implicates early selection favoring individuals harboring optimal maternal-like organellar SSU3D targeting in hexaploid wheat. These data represent a novel dimension of cytonuclear evolution mediated by organellar targeting and transportation of nuclear proteins.

Maximum Likelihood Estimation of Species Trees from Gene Trees in the Presence of Ancestral Population Structure

Though large multilocus genomic data sets have led to overall improvements in phylogenetic inference, they have posed the new challenge of addressing conflicting signals across the genome. In particular, ancestral population structure, which has been uncovered in a number of diverse species, can skew gene tree frequencies, thereby hindering the performance of species tree estimators. Here we develop a novel maximum likelihood method, termed TASTI (Taxa with Ancestral structure Species Tree Inference), that can infer phylogenies under such scenarios, and find that it has increasing accuracy with increasing numbers of input gene trees, contrasting with the relatively poor performances of methods not tailored for ancestral structure. Moreover, we propose a supertree approach that allows TASTI to scale computationally with increasing numbers of input taxa. We use genetic simulations to assess TASTI's performance in the three- and four-taxon settings and demonstrate the application of TASTI on a six-species Afrotropical mosquito data set. Finally, we have implemented TASTI in an open-source software package for ease of use by the scientific community.

Rhizospheric bacterial community structure of Triticum and Aegilops revealed by pyrosequencing analysis of the 16S rRNA gene: dominance of the A genome over the B and D genomes

This study examined the relationship between host plant and rhizospheric bacterial community structure, including composition and diversity, in Triticum and Aegilops species (12 and two accessions, respectively) as well as three closely related species, barley, rye and oat (four accessions), to explore the possibility that wheat root and rhizosphere interaction can be utilized for wheat breeding and biotechnology in the future. For this purpose, DNA was isolated from rhizospheric soil samples and one control non-rhizospheric soil sample, and the 16S rRNA gene region was amplified and subjected to DNA pyrosequencing. A total of 132,888 amplicons were analyzed. Bacterial composition at the phylum level was similar among the 18 rhizospheric samples; however, the proportion of Acidobacteria was much lower in these samples than in the control non-rhizospheric soil sample, indicating that rhizospheres influenced the bacterial composition even at the higher taxonomic level. Across host plant genome types (three levels of ploidy and three major genomes, A, B and D), there was no detectable difference in phylum composition or species diversity. Estimated bacterial species diversity was higher in the control soil sample than in plant rhizospheric soils, implying that bacterial species diversity was reduced in rhizospheres. A PCoA plot and UPGMA dendrogram based on the bacterial species composition showed that control soil was distantly located from the plant rhizospheric samples and that Triticum, Aegilops and related species were well separated. PERMANOVA analysis detected statistically significant differentiation among these four groups. Clustering of Triticum species suggested that the A genome was dominant over the B and D genomes, with respect to the influence on rhizospheric bacterial species composition. Although the cause was not investigated in this study, these results clearly indicated that the genetic constitution of the plant host exerted a strong influence on rhizospheric bacterial community structure.

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

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

Identification and Characterization of PLATZ Transcription Factors in Wheat

The PLATZ (plant AT-rich protein and zinc-binding protein) transcription factor family is a class of plant-specific zinc-dependent DNA-binding proteins. PLATZ has essential roles in seed endosperm development, as well as promoting cell proliferation duration in the earlier stages of the crops. In the present study, 62 TaPLATZ genes were identified from the wheat genome, and they were unequally distributed on 15 chromosomes. According to the phylogenetic analysis, 62 TaPLATZ genes were classified into six groups, including two groups that were unique in wheat. Members in the same groups shared similar exon-intron structures. The polyploidization, together with genome duplication of wheat, plays a crucial role in the expansion of the TaPLATZs family. Transcriptome data indicated a distinct divergence expression pattern of TaPLATZ genes that could be clustered into four modules. The TaPLATZs in Module b possessed a seed-specific expression pattern and displayed obvious high expression in the earlier development stage of seeds. Subcellular localization data of TaPLATZs suggesting that they likely perform a function as a conventional transcription factor. This study provides insight into understanding the structure divergence, evolutionary features, expression profiles, and potential function of PLATZ in wheat.

Evolution of Homeologous Gene Expression in Polyploid Wheat

Polyploidization has played a prominent role in the evolutionary history of plants. Two recent and sequential allopolyploidization events have resulted in the formation of wheat species with different ploidies, and which provide a model to study the effects of polyploidization on the evolution of gene expression. In this study, we identified differentially expressed genes (DEGs) between four BBAA tetraploid wheats of three different ploidy backgrounds. DEGs were found to be unevenly distributed among functional categories and duplication modes. We observed more DEGs in the extracted tetraploid wheat (ETW) than in natural tetraploid wheats (TD and TTR13) as compared to a synthetic tetraploid (AT2). Furthermore, DEGs showed higher Ka/Ks ratios than those that did not show expression changes (non-DEGs) between genotypes, indicating DEGs and non-DEGs experienced different selection pressures. For A-B homeolog pairs with DEGs, most of them had only one differentially expressed copy, however, when both copies of a homeolog pair were DEGs, the A and B copies were more likely to be regulated to the same direction. Our results suggest that both cis- and inter-subgenome trans-regulatory changes are important drivers in the evolution of homeologous gene expression in polyploid wheat, with ploidy playing a significant role in the process.

Multiple wheat genomes reveal global variation in modern breeding

Advances in genomics have expedited the improvement of several agriculturally important crops but similar efforts in wheat (Triticum spp.) have been more challenging. This is largely owing to the size and complexity of the wheat genome¹, and the lack of genome-assembly data for multiple wheat lines^2,3. Here we generated ten chromosome pseudomolecule and five scaffold assemblies of hexaploid wheat to explore the genomic diversity among wheat lines from global breeding programs. Comparative analysis revealed extensive structural rearrangements, introgressions from wild relatives and differences in gene content resulting from complex breeding histories aimed at improving adaptation to diverse environments, grain yield and quality, and resistance to stresses^4,5. We provide examples outlining the utility of these genomes, including a detailed multi-genome-derived nucleotide-binding leucine-rich repeat protein repertoire involved in disease resistance and the characterization of Sm1⁶, a gene associated with insect resistance. These genome assemblies will provide a basis for functional gene discovery and breeding to deliver the next generation of modern wheat cultivars.

Comparison of multiple genome assemblies from wheat reveals extensive diversity that results from the complex breeding history of wheat and provides a basis for further potential improvements to this important food crop.

Seed-Derived Microbial Colonization of Wild Emmer and Domesticated Bread Wheat (Triticum dicoccoides and T. aestivum) Seedlings Shows Pronounced Differences in Overall Diversity and Composition

The composition of the plant microbiota may be altered by ecological and evolutionary changes in the host population. Seed-associated microbiota, expected to be largely vertically transferred, have the potential to coadapt with their host over generations. Strong directional selection and changes in the genetic composition of plants during domestication and cultivation may have impacted the assembly and transmission of seed-associated microbiota. Nonetheless, the effect of plant speciation and domestication on the composition of these microbes is poorly understood. Here, we have investigated the composition of bacteria and fungi associated with the wild emmer wheat (Triticum dicoccoides) and domesticated bread wheat (Triticum aestivum). We show that vertically transmitted bacteria, but not fungi, of domesticated bread wheat species T. aestivum are less diverse and more inconsistent among individual plants compared to those of the wild emmer wheat species T. dicoccoides. We propagated wheat seeds under sterile conditions to characterize the colonization of seedlings by seed-associated microbes. Hereby, we show markedly different community compositions and diversities of leaf and root colonizers of the domesticated bread wheat compared to the wild emmer wheat. By propagating the wild emmer wheat and domesticated bread wheat in two different soils, we furthermore reveal a small effect of plant genotype on microbiota assembly. Our results suggest that domestication and prolonged breeding have impacted the vertically transferred bacteria, but only to a lesser extent have affected the soil-derived microbiota of bread wheat.IMPORTANCE Genetic and physiological changes associated with plant domestication have been studied for many crop species. Still little is known about the impact of domestication on the plant-associated microbiota. In this study, we analyze the seed-associated and soil-derived bacterial and fungal microbiota of domesticated bread wheat and wild emmer wheat. We show a significant difference in the seed-associated, but not soil-derived, bacterial communities of the wheat species. Interestingly, we find less pronounced effects on the fungal communities. Overall, this study provides novel insight into the diversity of vertically transmitted microbiota of wheat and thereby contributes to our understanding of wheat as a "metaorganism." Insight into the wheat microbiota is of fundamental importance for the development of improved crops.

Production of synthetic wheat lines to exploit the genetic diversity of emmer wheat and D genome containing Aegilops species in wheat breeding

Due to the accumulation of various useful traits over evolutionary time, emmer wheat (Triticum turgidum subsp. dicoccum and dicoccoides, 2n = 4x = 28; AABB), durum wheat (T. turgidum subsp. durum, 2n = 4x = 28; AABB), T. timopheevii (2n = 4x = 28; AAGG) and D genome containing Aegilops species offer excellent sources of novel variation for the improvement of bread wheat (T. aestivum L., AABBDD). Here, we made 192 different cross combinations between diverse genotypes of wheat and Aegilops species including emmer wheat × Ae. tauschii (2n = DD or DDDD), durum wheat × Ae. tauschii, T. timopheevii × Ae. tauschii, Ae. crassa × durum wheat, Ae. cylindrica × durum wheat and Ae. ventricosa × durum wheat in the field over three successive years. We successfully recovered 56 different synthetic hexaploid and octaploid F₂ lines with AABBDD, AABBDDDD, AAGGDD, D¹D¹X^crX^crAABB, D^cD^cC^cC^cAABB and D^vD^vN^vN^vAABB genomes via in vitro rescue of F₁ embryos and spontaneous production of F₂ seeds on the F_l plants. Cytogenetic analysis of F₂ lines showed that the produced synthetic wheat lines were generally promising stable amphiploids. Contribution of D genome bearing Aegilops and the less-investigated emmer wheat genotypes as parents in the crosses resulted in synthetic amphiploids which are a valuable resource for bread wheat breeding.

Extensive Genetic Variation at the Sr22 Wheat Stem Rust Resistance Gene Locus in the Grasses Revealed Through Evolutionary Genomics and Functional Analyses

In the last 20 years, severe wheat stem rust outbreaks have been recorded in Africa, Europe, and Central Asia. This previously well controlled disease, caused by the fungus Puccinia graminis f. sp. tritici, has reemerged as a major threat to wheat cultivation. The stem rust (Sr) resistance gene Sr22 encodes a nucleotide-binding and leucine-rich repeat receptor which confers resistance to the highly virulent African stem rust isolate Ug99. Here, we show that the Sr22 gene is conserved among grasses in the Triticeae and Poeae lineages. Triticeae species contain syntenic loci with single-copy orthologs of Sr22 on chromosome 7, except Hordeum vulgare, which has experienced major expansions and rearrangements at the locus. We also describe 14 Sr22 sequence variants obtained from both Triticum boeoticum and the domesticated form of this species, T. monococcum, which have been postulated to encode both functional and nonfunctional Sr22 alleles. The nucleotide sequence analysis of these alleles identified historical sequence exchange resulting from recombination or gene conversion, including breakpoints within codons, which expanded the coding potential at these positions by introduction of nonsynonymous substitutions. Three Sr22 alleles were transformed into wheat cultivar Fielder and two postulated resistant alleles from Schomburgk (hexaploid wheat introgressed with T. boeoticum segment carrying Sr22) and T. monococcum accession PI190945, respectively, conferred resistance to P. graminis f. sp. tritici race TTKSK, thereby unequivocally confirming Sr22 effectiveness against Ug99. The third allele from accession PI573523, previously believed to confer susceptibility, was confirmed as nonfunctional against Australian P. graminis f. sp. tritici race 98-1,2,3,5,6.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Triticum population sequencing provides insights into wheat adaptation

Bread wheat expanded its habitat from a core area of the Fertile Crescent to global environments within ~10,000 years. The genetic mechanisms of this remarkable evolutionary success are not well understood. By whole-genome sequencing of populations from 25 subspecies within the genera Triticum and Aegilops, we identified composite introgression from wild populations contributing to a substantial portion (4-32%) of the bread wheat genome, which increased the genetic diversity of bread wheat and allowed its divergent adaptation. Meanwhile, convergent adaptation to human selection showed 2- to 16-fold enrichment relative to random expectation-a certain set of genes were repeatedly selected in Triticum species despite their drastic differences in ploidy levels and growing zones, indicating the important role of evolutionary constraints in shaping the adaptive landscape of bread wheat. These results showed the genetic necessities of wheat as a global crop and provided new perspectives on transferring adaptive success across species for crop improvement.

Origin and adaptation to high altitude of Tibetan semi-wild wheat

Tibetan wheat is grown under environmental constraints at high-altitude conditions, but its underlying adaptation mechanism remains unknown. Here, we present a draft genome sequence of a Tibetan semi-wild wheat (Triticum aestivum ssp. tibetanum Shao) accession Zang1817 and re-sequence 245 wheat accessions, including world-wide wheat landraces, cultivars as well as Tibetan landraces. We demonstrate that high-altitude environments can trigger extensive reshaping of wheat genomes, and also uncover that Tibetan wheat accessions accumulate high-altitude adapted haplotypes of related genes in response to harsh environmental constraints. Moreover, we find that Tibetan semi-wild wheat is a feral form of Tibetan landrace, and identify two associated loci, including a 0.8-Mb deletion region containing Brt1/2 homologs and a genomic region with TaQ-5A gene, responsible for rachis brittleness during the de-domestication episode. Our study provides confident evidence to support the hypothesis that Tibetan semi-wild wheat is de-domesticated from local landraces, in response to high-altitude extremes.

Mechanism of high altitude adaptation of wheat remains unknown. Here, the authors assemble the draft genome of a Tibetan semi-wild wheat accession and resequence 245 wheat accessions to reveal that Tibetan semi-wild wheat has been de-domesticated from local landraces to adapt to high altitude.

Genome-wide characterization and expression analysis suggested diverse functions of the mechanosensitive channel of small conductance-like (MSL) genes in cereal crops

Mechanosensitive ion channels are pore-forming transmembrane proteins that allow ions to move down their electrochemical gradient in response to mechanical stimuli. They participate in many plant developmental processes including the maintenance of plastid shape, pollen tube growth, etc. Herein, a total of 11, 10, 6, 30, 9, and 8 MSL genes were identified in Aegilops tauschii, Hordeum vulgare, Sorghum bicolor, Triticum aestivum, Triticum urartu, and Zea mays, respectively. These genes were located on various chromosomes of their respective cereal, while MSLs of T. urartu were found on scaffolds. The phylogenetic analysis, subcellular localization, and sequence homology suggested clustering of MSLs into two classes. These genes consisted of cis-regulatory elements related to growth and development, responsive to light, hormone, and stress. Differential expression of various MSL genes in tissue developmental stages and stress conditions revealed their precise role in development and stress responses. Altered expression during CaCl2 stress suggested their role in Ca2+ homeostasis and signaling. The co-expression analysis suggested their interactions with other genes involved in growth, defense responses etc. A comparative expression profiling of paralogous genes revealed either retention of function or pseudo-functionalization. The present study unfolded various characteristics of MSLs in cereals, which will facilitate their in-depth functional characterization in future studies.

Genetic Diversity and Population Structure of Asian and European Common Wheat Accessions Based on Genotyping-By-Sequencing

Obtaining information on the genetic diversity and population structure of germplasm facilitates its use in wheat breeding programs. Recently, with the development of next-generation sequencing technology, genotyping-by-sequencing (GBS) has been used as a high-throughput and cost-effective molecular tool for examination of the genetic diversity of wheat breeding lines. In this study, GBS was used to characterize a population of 180 accessions of common wheat originating from Asia and Europe between the latitudes 30° and 45°N. In total, 24,767 high-quality single-nucleotide polymorphism (SNP) markers were used for analysis of genetic diversity and population structure. The B genome contained the highest number of SNPs, followed by the A and D genomes. The polymorphism information content was in the range of 0.1 to 0.4, with a mean of 0.26. The distribution of SNPs markers on the 21 chromosomes ranged from 243 on chromosome 4D to 2,337 on chromosome 3B. Structure and cluster analyses divided the panel of accessions into two subgroups (G1 and G2). G1 principally consisted of European and partial Asian accessions, and G2 comprised mainly accessions from the Middle East and partial Asia. Molecular analysis of variance showed that the genetic variation was greater within groups (99%) than between groups (1%). Comparison of the two subgroups indicated that G1 and G2 contained a high level of genetic diversity. The genetic diversity of G2 was slightly higher as indicated by the observed heterozygosity (H o) = 0.23, and unbiased diversity index (uh) = 0.34. The present results will not only help breeders to understand the genetic diversity of wheat germplasm on the Eurasian continent between the latitudes of 30° and 45°N, but also provide valuable information for wheat genetic improvement through introgression of novel genetic variation in this region.

Application of Genomics Tools in Wheat Breeding to Attain Durable Rust Resistance

Wheat is an important source of dietary protein and calories for the majority of the world's population. It is one of the largest grown cereal in the world occupying over 215 M ha. Wheat production globally is challenged by biotic stresses such as pests and diseases. Of the 50 diseases of wheat that are of economic importance, the three rust diseases are the most ubiquitous causing significant yield losses in the majority of wheat production environments. Under severe epidemics they can lead to food insecurity threats amid the continuous evolution of new races of the pathogens, shifts in population dynamics and their virulence patterns, thereby rendering several effective resistance genes deployed in wheat breeding programs vulnerable. This emphasizes the need to identify, characterize, and deploy effective rust-resistant genes from diverse sources into pre-breeding lines and future wheat varieties. The use of genetic resistance has been marked as eco-friendly and to curb the further evolution of rust pathogens. Deployment of multiple rust resistance genes including major and minor genes in wheat lines could enhance the durability of resistance thereby reducing pathogen evolution. Advances in next-generation sequencing (NGS) platforms and associated bioinformatics tools have revolutionized wheat genomics. The sequence alignment of the wheat genome is the most important landmark which will enable genomics to identify marker-trait associations, candidate genes and enhanced breeding values in genomic selection (GS) studies. High throughput genotyping platforms have demonstrated their role in the estimation of genetic diversity, construction of the high-density genetic maps, dissecting polygenic traits, and better understanding their interactions through GWAS (genome-wide association studies) and QTL mapping, and isolation of R genes. Application of breeder's friendly KASP assays in the wheat breeding program has expedited the identification and pyramiding of rust resistance alleles/genes in elite lines. The present review covers the evolutionary trends of the rust pathogen and contemporary wheat varieties, and how these research strategies galvanized to control the wheat killer genus Puccinia. It will also highlight the outcome and research impact of cost-effective NGS technologies and cloning of rust resistance genes amid the public availability of common and tetraploid wheat reference genomes.

Genome-wide identification and expression profiling of cytosine-5 DNA methyltransferases during drought and heat stress in wheat (Triticum aestivum)

DNA methylation is a potential epigenetic mechanism that regulates genome stability, development, and stress mitigation in plants. It is mediated by cytosine-5 DNA methyltransferases (C5-MTases). We identified 52 wheat C5-MTases; and based on domain structure and phylogenetics, these 52 C5-MTases were classified into four sub-families including MET, CMT, DRM and DNMT2; and were distributed on 18 chromosomes. Cis-acting regulatory elements analysis identified abiotic stress-responsive, phytohormone-responsive, development-related and light-related elements in the promoters of TaC5-MTases. We also examined the transcript abundance of TaC5-MTases in different tissues, developmental stages and under abiotic stresses. Notably, most of the TaC5-MTases (TaCMT2, TaCMT3b, TaCMT3c, TaMET1, TaDRM10, TaDNMT2) showed differential regulation of their transcript abundance during drought and heat stress. Overall, the above results provide significant insights into the expression and the probable functions of TaC5-MTases and will also expedite future research programs to explore the mechanisms of epigenetic regulation in wheat.

Bread wheat TaSPO11-1 exhibits evolutionarily conserved function in meiotic recombination across distant plant species

The manipulation of meiotic recombination in crops is essential to develop new plant varieties rapidly, helping to produce more cultivars in a sustainable manner. One option is to control the formation and repair of the meiosis-specific DNA double-strand breaks (DSBs) that initiate recombination between the homologous chromosomes and ultimately lead to crossovers. These DSBs are introduced by the evolutionarily conserved topoisomerase-like protein SPO11 and associated proteins. Here, we characterized the homoeologous copies of the SPO11-1 protein in hexaploid bread wheat (Triticum aestivum). The genome contains three SPO11-1 gene copies that exhibit 93-95% identity at the nucleotide level, and clearly the A and D copies originated from the diploid ancestors Triticum urartu and Aegilops tauschii, respectively. Furthermore, phylogenetic analysis of 105 plant genomes revealed a clear partitioning between monocots and dicots, with the seven main motifs being almost fully conserved, even between clades. The functional similarity of the proteins among monocots was confirmed through complementation analysis of the Oryza sativa (rice) spo11-1 mutant by the wheat TaSPO11-1-5D coding sequence. Also, remarkably, although the wheat and Arabidopsis SPO11-1 proteins share only 55% identity and the partner proteins also differ, the TaSPO11-1-5D cDNA significantly restored the fertility of the Arabidopsis spo11-1 mutant, indicating a robust functional conservation of the SPO11-1 protein activity across distant plants. These successful heterologous complementation assays, using both Arabidopsis and rice hosts, are good surrogates to validate the functionality of candidate genes and cDNA, as well as variant constructs, when the transformation and mutant production in wheat is much longer and more tedious.

SPO11.2 is essential for programmed double-strand break formation during meiosis in bread wheat (Triticum aestivum L.)

Meiotic recombination is initiated by formation of DNA double-strand breaks (DSBs). This involves a protein complex that includes in plants the two similar proteins, SPO11-1 and SPO11-2. We analysed the sequences of SPO11-2 in hexaploid bread wheat (Triticum aestivum), as well as in its diploid and tetraploid progenitors. We investigated its role during meiosis using single, double and triple mutants. The three homoeologous SPO11-2 copies of hexaploid wheat exhibit high nucleotide and amino acid similarities with those of the diploids, tetraploids and Arabidopsis. Interestingly, however, two nucleotides deleted in exon-2 of the A copy lead to a premature stop codon and suggest that it encodes a non-functional protein. Remarkably, the mutation was absent from the diploid A-relative Triticum urartu, but present in the tetraploid Triticum dicoccoides and in different wheat cultivars indicating that the mutation occurred after the first polyploidy event and has since been conserved. We further show that triple mutants with all three copies (A, B, D) inactivated are sterile. Cytological analyses of these mutants show synapsis defects, accompanied by severe reductions in bivalent formation and numbers of DMC1 foci, thus confirming the essential role of TaSPO11-2 in meiotic recombination in wheat. In accordance with its 2-nucleotide deletion in exon-2, double mutants for which only the A copy remained are also sterile. Notwithstanding, some DMC1 foci remain visible in this mutant, suggesting a residual activity of the A copy, albeit not sufficient to restore fertility.

Structural Variations Affecting Genes and Transposable Elements of Chromosome 3B in Wheats

Structural variations (SVs) such as copy number and presence–absence variations are polymorphisms that are known to impact genome composition at the species level and are associated with phenotypic variations. In the absence of a reference genome sequence, their study has long been hampered in wheat. The recent production of new wheat genomic resources has led to a paradigm shift, making possible to investigate the extent of SVs among cultivated and wild accessions. We assessed SVs affecting genes and transposable elements (TEs) in a Triticeae diversity panel of 45 accessions from seven tetraploid and hexaploid species using high-coverage shotgun sequencing of sorted chromosome 3B DNA and dedicated bioinformatics approaches. We showed that 23% of the genes are variable within this panel, and we also identified 330 genes absent from the reference accession Chinese Spring. In addition, 60% of the TE-derived reference markers were absent in at least one accession, revealing a high level of intraspecific and interspecific variability affecting the TE space. Chromosome extremities are the regions where we observed most of the variability, confirming previous hypotheses made when comparing wheat with the other grasses. This study provides deeper insights into the genomic variability affecting the complex Triticeae genomes at the intraspecific and interspecific levels and suggests a phylogeny with independent hybridization events leading to different hexaploid species.

Wheat Disease Resistance Genes and Their Diversification Through Integrated Domain Fusions

Plants are in a constant evolutionary arms race with their pathogens. At the molecular level, the plant nucleotide-binding leucine-rich repeat receptors (NLRs) family has coevolved with rapidly evolving pathogen effectors. While many NLRs utilize variable leucine-rich repeats (LRRs) to detect effectors, some have gained integrated domains (IDs) that may be involved in receptor activation or downstream signaling. The major objectives of this project were to identify NLR genes in wheat (Triticum aestivum L.) and assess IDs associated with immune signaling (e.g., kinase and transcription factor domains). We identified 2,151 NLR-like genes in wheat, of which 1,298 formed 547 gene clusters. Among the non-toll/interleukin-1 receptor NLR (non-TNL)-like genes, 1,552 encode LRRs, 802 are coiled-coil (CC) domain-encoding (CC-NBS-LRR or CNL) genes, and three encode resistance to powdery mildew 8 (RPW8) domains (RPW8-NBS-LRR or RNL). The expansion of the NLR gene family in wheat is attributable to its origin by recent polyploidy events. Gene clusters were likely formed by tandem duplications, and wheat NLR phylogenetic relationships were similar to those in barley and Aegilops. We also identified wheat NLR-ID fusion proteins as candidates for NLR functional diversification, often as kinase and transcription factor domains. Comparative analyses of the IDs revealed evolutionary conservation of more than 80% amino acid sequence similarity. Homology assessment indicates that these domains originated as functional non-NLR-encoding genes that were incorporated into NLR-encoding genes through duplication events. We also found that many of the NLR-ID genes encode alternative transcripts that include or exclude IDs, a phenomenon that seems to be conserved among species. To verify this, we have analyzed the alternative transcripts that include or exclude an ID of an NLR-ID from another monocotyledon species, rice (Oryza sativa). This indicates that plants employ alternative splicing to regulate IDs, possibly using them as baits, decoys, and functional signaling components. Genomic and expression data support the hypothesis that wheat uses alternative splicing to include and exclude IDs from NLR proteins.

Recent Progress in Germplasm Evaluation and Gene Mapping to Enable Breeding of Drought-Tolerant Wheat

There is a need to increase wheat productivity to meet the food demands of the ever-growing human population. However, accelerated development of high yielding varieties is hindered by drought, which is worsening due to climate change. In this context, germplasm diversity is central to the development of drought-tolerant wheat. Extensive collections of these genetic resources are conserved in national and international genebanks. In addition to phenotypic assessments, the use of advanced molecular techniques (e.g., genotype by sequencing) to identify quantitative trait loci (QTLs) for drought tolerance related traits is useful for genome- and marker-assisted selection based approaches. Therefore, to assist wheat breeders at a critical time, we searched the recent peer-reviewed literature (2011-current), first, to identify wheat germplasm observed to be useful genetic sources for drought tolerance, and second, to report QTLs associated with drought tolerance. Though many breeders limit the parents used in breeding programs to a familiar core collection, the results of this review show that larger germplasm collections have been sources of useful genes for drought tolerance in wheat. The review also demonstrates that QTLs for drought tolerance in wheat are associated with diverse physio-morphological traits, at different growth stages. Here, we also briefly discuss the potential of genome engineering/editing to improve drought tolerance in wheat. The use of CRISPR-Cas9 and other gene-editing technologies can be used to fine-tune the expression of genes controlling drought adaptive traits, while high throughput phenotyping (HTP) techniques can potentially accelerate the selection process. These efforts are empowered by wheat researcher consortia.

The Brittle Rachis Trait in Species Belonging to the Triticeae and Its Controlling Genes Btr1 and Btr2

In many non-cultivated angiosperm species, seed dispersal is facilitated by the shattering of the seed head at maturity; in the Triticeae tribe, to which several of the world's most important cereals belong, shattering takes the form of a disarticulation of the rachis. The products of the genes Btr1 and Btr2 are both required for disarticulation to occur above the rachis nodes within the genera Hordeum (barley) and Triticum/Aegilops (wheat). Here, it has been shown that both Btr1 and Btr2 are specific to the Triticeae tribe, although likely paralogs (Btr1-like and Btr2-like) are carried by the family Poaceae including Triticeae. Aegilops tauschii (the donor of the bread wheat D genome) lacks a copy of Btr1 and disarticulation in this species occurs below, rather than above the rachis node; thus, the product of Btr1 appears to be required for disarticulation to occur above the rachis node.

Illegitimate Recombination Between Homeologous Genes in Wheat Genome

Polyploidies produce a large number of duplicated regions and genes in genomes, which have a long-term impact and stimulate genetic innovation. The high similarity between homeologous chromosomes, forming different subgenomes, or homologous regions after genome repatterning, may permit illegitimate DNA recombination. Here, based on gene colinearity, we aligned the (sub)genomes of common wheat (Triticum aestivum, AABBDD genotype) and its relatives, including Triticum urartu (AA), Aegilops tauschii (DD), and T. turgidum ssp. dicoccoides (AABB) to detect the homeologous (paralogous or orthologous) colinear genes within and between (sub)genomes. Besides, we inferred more ancient paralogous regions produced by a much ancient grass-common tetraploidization. By comparing the sequence similarity between paralogous and orthologous genes, we assumed abnormality in the topology of constructed gene trees, which could be explained by gene conversion as a result of illegitimate recombination. We found large numbers of inferred converted genes (>2,000 gene pairs) suggested long-lasting genome instability of the hexaploid plant, and preferential donor roles by DD genes. Though illegitimate recombination was much restricted, duplicated genes produced by an ancient whole-genome duplication, which occurred millions of years ago, also showed evidence of likely gene conversion. As to biological function, we found that ~40% catalytic genes in colinearity, including those involved in starch biosynthesis, were likely affected by gene conversion. The present study will contribute to understanding the functional and structural innovation of the common wheat genome.

Differential chromatin accessibility landscape reveals structural and functional features of the allopolyploid wheat chromosomes

Background

Our understanding of how the complexity of the wheat genome influences the distribution of chromatin states along the homoeologous chromosomes is limited. Using a differential nuclease sensitivity assay, we investigate the chromatin states of the coding and repetitive regions of the allopolyploid wheat genome.

Results

Although open chromatin is found to be significantly enriched around genes, the majority of MNase-sensitive regions are located within transposable elements (TEs). Chromatin of the smaller D genome is more accessible than that of the larger A and B genomes. Chromatin states of different TEs vary among families and are influenced by the TEs’ chromosomal position and proximity to genes. While the chromatin accessibility of genes is influenced by proximity to TEs, and not by their position on the chromosomes, we observe a negative chromatin accessibility gradient along the telomere-centromere axis in the intergenic regions, positively correlated with the distance between genes. Both gene expression levels and homoeologous gene expression bias are correlated with chromatin accessibility in promoter regions. The differential nuclease sensitivity assay accurately predicts previously detected centromere locations. SNPs located within more accessible chromatin explain a higher proportion of genetic variance for a number of agronomic traits than SNPs located within more closed chromatin.

Conclusions

Chromatin states in the wheat genome are shaped by the interplay of repetitive and gene-encoding regions that are predictive of the functional and structural organization of chromosomes, providing a powerful framework for detecting genomic features involved in gene regulation and prioritizing genomic variation to explain phenotypes.

Molecular characterization of ascorbate peroxidase (APX) and APX-related (APX-R) genes in Triticum aestivum L

Ascorbate peroxidases (APXs) are heme-dependent H₂O₂ scavenging enzymes involved in myriad biological processes. Herein, a total of 21 TaAPX and six TaAPX-R genes were identified from the A, B and D sub-genomes of Triticum aestivum L. The occurrence of three paralogous gene pairs with unequal evolutionary rate suggested functional divergence. The phylogenetic analysis formed four distinct clades having conserved gene and protein architecture, and sub-cellular localization. The tertiary structure analysis revealed the presence of helices and coils and residues involved in ligand binding. Transcriptional profiling of each TaAPX and TaAPX-R gene suggested their specific role during development and stress response. Modulated transcript expression and APX enzyme activity during various stress conditions indicated their role in stress response. Interaction analyses suggested their association with other genes, miRNAs and various legends. The present study reported numerous features of these genes, and may provide a platform for their detailed functional characterization in future studies.

Evolutionary Insights into the Nature of Plant Domestication

Domestication is a co-evolutionary process that occurs when wild plants are brought into cultivation by humans, leading to origin of new species and/or differentiated populations that are critical for human survival. Darwin used domesticated species as early models for evolution, highlighting their variation and the key role of selection in species differentiation. Over the last two decades, a growing synthesis of plant genetics, genomics, and archaeobotany has led to challenges to old orthodoxies and the advent of fresh perspectives on how crop domestication and diversification proceed. I discuss four new insights into plant domestication - that in general domestication is a protracted process, that unconscious (natural) selection plays a prominent role, that interspecific hybridization may be an important mechanism for crop species diversification and range expansion, and that similar genes across multiple species underlies parallel/convergent phenotypic evolution between domesticated taxa. Insights into the evolutionary origin and diversification of crop species can help us in developing new varieties (and possibly even new species) to deal with current and future environmental challenges in a sustainable manner.

Chromosome Pairing in Polyploid Grasses

Polyploids are species in which three or more sets of chromosomes coexist. Polyploidy frequently occurs in plants and plays a major role in their evolution. Based on their origin, polyploid species can be divided into two groups: autopolyploids and allopolyploids. The autopolyploids arise by multiplication of the chromosome sets from a single species, whereas allopolyploids emerge from the hybridization between distinct species followed or preceded by whole genome duplication, leading to the combination of divergent genomes. Having a polyploid constitution offers some fitness advantages, which could become evolutionarily successful. Nevertheless, polyploid species must develop mechanism(s) that control proper segregation of genetic material during meiosis, and hence, genome stability. Otherwise, the coexistence of more than two copies of the same or similar chromosome sets may lead to multivalent formation during the first meiotic division and subsequent production of aneuploid gametes. In this review, we aim to discuss the pathways leading to the formation of polyploids, the occurrence of polyploidy in the grass family (Poaceae), and mechanisms controlling chromosome associations during meiosis, with special emphasis on wheat.

Reduced chromatin accessibility underlies gene expression differences in homologous chromosome arms of diploid Aegilops tauschii and hexaploid wheat

Background

Polyploidy is centrally important in the evolution and domestication of plants because it leads to major genomic changes, such as altered patterns of gene expression, which are thought to underlie the emergence of new traits. Despite the common occurrence of these globally altered patterns of gene expression in polyploids, the mechanisms involved are not well understood.

Results

Using a precisely defined framework of highly conserved syntenic genes on hexaploid wheat chromosome 3DL and its progenitor 3 L chromosome arm of diploid Aegilops tauschii, we show that 70% of these gene pairs exhibited proportionately reduced gene expression, in which expression in the hexaploid context of the 3DL genes was ∼40% of the levels observed in diploid Ae tauschii. Several genes showed elevated expression during the later stages of grain development in wheat compared with Ae tauschii. Gene sequence and methylation differences probably accounted for only a few cases of differences in gene expression. In contrast, chromosome-wide patterns of reduced chromatin accessibility of genes in the hexaploid chromosome arm compared with its diploid progenitor were correlated with both reduced gene expression and the imposition of new patterns of gene expression.

Conclusions

Our pilot-scale analyses show that chromatin compaction may orchestrate reduced gene expression levels in the hexaploid chromosome arm of wheat compared to its diploid progenitor chromosome arm.

Homoeologous exchanges occur through intragenic recombination generating novel transcripts and proteins in wheat and other polyploids

Recombination between homeologous chromosomes, also known as homeologous exchange (HE), plays a significant role in shaping genome structure and gene expression in interspecific hybrids and allopolyploids of several plant species. However, the molecular mechanisms that govern HEs are not well understood. Here, we studied HE events in the progeny of a nascent allotetraploid (genome AADD) derived from two diploid progenitors of hexaploid bread wheat using cytological and whole-genome sequence analyses. In total, 37 HEs were identified and HE junctions were mapped precisely. HEs exhibit typical patterns of homologous recombination hotspots, being biased toward low-copy, subtelomeric regions of chromosome arms and showing association with known recombination hotspot motifs. But, strikingly, while homologous recombination preferentially takes place upstream and downstream of coding regions, HEs are highly enriched within gene bodies, giving rise to novel recombinant transcripts, which in turn are predicted to generate new protein fusion variants. To test whether this is a widespread phenomenon, a dataset of high-resolution HE junctions was analyzed for allopolyploid Brassica, rice, Arabidopsis suecica, banana, and peanut. Intragenic recombination and formation of chimeric genes was detected in HEs of all species and was prominent in most of them. HE thus provides a mechanism for evolutionary novelty in transcript and protein sequences in nascent allopolyploids.

Supernumerary B chromosomes of Aegilops speltoides undergo precise elimination in roots early in embryo development

Not necessarily all cells of an organism contain the same genome. Some eukaryotes exhibit dramatic differences between cells of different organs, resulting from programmed elimination of chromosomes or their fragments. Here, we present a detailed analysis of programmed B chromosome elimination in plants. Using goatgrass Aegilops speltoides as a model, we demonstrate that the elimination of B chromosomes is a strictly controlled and highly efficient root-specific process. At the onset of embryo differentiation B chromosomes undergo elimination in proto-root cells. Independent of centromere activity, B chromosomes demonstrate nondisjunction of chromatids and lagging in anaphase, leading to micronucleation. Chromatin structure and DNA replication differ between micronuclei and primary nuclei and degradation of micronucleated DNA is the final step of B chromosome elimination. This process might allow root tissues to survive the detrimental expression, or overexpression of B chromosome-located root-specific genes with paralogs located on standard chromosomes.

B chromosomes are supernumerary chromosomes exhibiting dramatic differences between different organs in same species. Here, the authors show programmed B chromosome elimination in goatgrass starts at the onset of embryo differentiation by nondisjunction of chromatids, anaphase lagging, and ends with the degradation of micronucleated DNA.

Preferential Subgenome Elimination and Chromosomal Structural Changes Occurring in Newly Formed Tetraploid Wheat-Aegilops ventricosa Amphiploid (AABBDvDvNvNv)

Artificial allopolyploids derived from the genera Triticum and Aegilops have been used as genetic resources for wheat improvement and are a classic example of evolution via allopolyploidization. In this study, we investigated chromosomes and subgenome transmission behavior in the newly formed allopolyploid of wheat group via multicolor Fluorescence in situ hybridization (mc-FISH), using pSc119.2, pTa535, and (GAA)₇ as probe combinations, to enabled us to precisely identify individual chromosomes in 381 S₃ and S₄ generations plants derived from reciprocal crosses between Ae. ventricosa (D^vD^vN^vN^v) and T. turgidum (AABB). A higher rate of aneuploidy, constituting 66.04–86.41% individuals, was observed in these two early generations. Of the four constituent subgenomes, D^v showed the highest frequency of elimination, followed by N^v and B, while A was the most stable. In addition, structural chromosomal changes occurred ubiquitously in the selfed progenies of allopolyploids. Among the constituent subgenomes, B showed the highest number of aberrations. In terms of chromosomal dynamics, there was no significant association between the chromosomal behavior model and the cytoplasm, with the exception of chromosomal loss in the D^v subgenome. The chromosome loss frequency in the D^v subgenome was significantly higher in the T. turgidum × Ae. ventricosa cross than in the Ae. ventricosa × T. turgidum cross. This result indicates that, although the D subgenome showed great instability, allopolyploids containing D subgenome could probably be maintained after a certain hybridization in which the D subgenome donor was used as the maternal parent at its onset stage. Our findings provide valuable information pertaining to the behavior patterns of subgenomes during allopolyploidization. Moreover, the allopolyploids developed here could be used as potential resources for the genetic improvement of wheat.

Chromosome restructuring among hybridizing wild wheats

The wheat group offers an outstanding system to address the interplay between hybridization, chromosomal evolution and biological diversification. Most diploid wild wheats originated following hybridization between the A-genome lineage and the B-genome lineage some 4 Myr ago, resulting in an admixed D-genome lineage that presented dramatic radiation accompanied by considerable changes in genome size and chromosomal rearrangements. Comparative profiling of low-copy genes, repeated sequences and transposable elements among those divergent species characterized by different karyotypes highlights high genome dynamics and sheds new light on the processes underlying chromosomal evolution in wild wheats. One of the hybrid clades presents upsizing of metacentric chromosomes going along with the proliferation of specific repeats (i.e. 'genomic obesity'), whereas other species show stable genome size associated with increasing chromosomal asymmetry. Genetic and ecological variation in those specialized species suggest that genome restructuring was coupled with adaptive processes to support the evolution of a majority of acrocentric chromosomes. This synthesis of current knowledge on genome restructuring across the diversity of wild wheats paves the way towards surveys based on latest sequencing technologies to characterize valuable resources and address the significance of chromosomal evolution in species with complex genomes.

Genome-Wide Discovery of G-Quadruplexes in Wheat: Distribution and Putative Functional Roles

G-quadruplexes are nucleic acid secondary structures formed by a stack of square planar G-quartets. G-quadruplexes were implicated in many biological functions including telomere maintenance, replication, transcription, and translation, in many species including humans and plants. For wheat, however, though it is one of the world's most important staple food, no G-quadruplex studies have been reported to date. Here, we computationally identify putative G4 structures (G4s) in wheat genome for the first time and compare its distribution across the genome against five other genomes (human, maize, Arabidopsis, rice, and sorghum). We identified close to 1 million G4 motifs with a density of 76 G4s/Mb across the whole genome and 93 G4s/Mb over genic regions. Remarkably, G4s were enriched around three regions, two located on the antisense and one on the sense strand at the following positions: 1) the transcription start site (TSS) (antisense), 2) the first coding domain sequence (CDS) (antisense), and 3) the start codon (sense). Functional enrichment analysis revealed that the gene models containing G4 motifs within these peaks were associated with specific gene ontology (GO) terms, such as developmental process, localization, and cellular component organization or biogenesis. We investigated genes encoding MADS-box transcription factors and showed examples of G4 motifs within critical regulatory regions in the VRN-1 genes in wheat. Furthermore, comparison with other plants showed that monocots share a similar distribution of G4s, but Arabidopsis shows a unique G4 distribution. Our study shows for the first time the prevalence and possible functional roles of G4s in wheat.

Genome-wide sequence information reveals recurrent hybridization among diploid wheat wild relatives

Many conflicting hypotheses regarding the relationships among crops and wild species closely related to wheat (the genera Aegilops, Amblyopyrum, and Triticum) have been postulated. The contribution of hybridization to the evolution of these taxa is intensely discussed. To determine possible causes for this, and provide a phylogeny of the diploid taxa based on genome-wide sequence information, independent data were obtained from genotyping-by-sequencing and a target-enrichment experiment that returned 244 low-copy nuclear loci. The data were analyzed using Bayesian, likelihood and coalescent-based methods. D statistics were used to test if incomplete lineage sorting alone or together with hybridization is the source for incongruent gene trees. Here we present the phylogeny of all diploid species of the wheat wild relatives. We hypothesize that most of the wheat-group species were shaped by a primordial homoploid hybrid speciation event involving the ancestral Triticum and Am. muticum lineages to form all other species except Ae. speltoides. This hybridization event was followed by multiple introgressions affecting all taxa except Triticum. Mostly progenitors of the extant species were involved in these processes, while recent interspecific gene flow seems insignificant. The composite nature of many genomes of wheat-group taxa results in complicated patterns of diploid contributions when these lineages are involved in polyploid formation, which is, for example, the case for tetraploid and hexaploid wheats. Our analysis provides phylogenetic relationships and a testable hypothesis for the genome compositions in the basic evolutionary units within the wheat group of Triticeae.

Genome-Scale Analysis of Homologous Genes among Subgenomes of Bread Wheat (Triticum aestivum L.)

Determining the distribution and correspondence of genome-scale homologous genes in wheat are effective ways to uncover chromosome rearrangement that has occurred during crop evolution and domestication, which can contribute to improvements in crop breeding. High-resolution and comprehensive analysis of the wheat genome by the International Wheat Genome Sequencing Consortium (IWGSC) revealed a total of 88,733 high-confidence homologous genes of four major types (1:1:1, 1:1:0, 0:1:1 and 1:0:1) among the A, B and D subgenomes of wheat. This data was used to compare homologous gene densities among chromosomes, clarify their distribution and correspondence relationship, and compare their functional enrichment. The average density of 1:1:1 homologous genes was about 10 times more than the density of the other three types of homologous genes, although the homologous gene densities of the various chromosomes were similar within each homologous type. Three regions of exceptional density were detected in 1:1:1 homologous genes, the isolate peak on the tail of chromosome 4A, and the desert regions at the start of chromosome 7A and 7D. The correspondence between homologous genes of the wheat subgenomes demonstrated translocation between the tail segments of chromosome 4A and 5A, and the inversion of the segment of original 5A and 7B into the tail of 4A. The homologous genes on the inserting segments of 5A and 7B to 4A were highly enriched in nitrogen, primary metabolite and small molecular metabolism processes, compared with genes on other regions of the original 4A chromosome. This study provides a refined genome-scale reference of homologous genes for wheat molecular research and breeding, which will help to broaden the application of the wheat genome and can be used as a template for research on other polyploid plants.

Partitioning and physical mapping of wheat chromosome 3B and its homoeologue 3E in Thinopyrum elongatum by inducing homoeologous recombination

We performed homoeologous recombination-based partitioning and physical mapping of wheat chromosome 3B and Th. elongatum chromosome 3E, providing a unique physical framework of this homoeologous pair for genome studies. The wheat (Triticum aestivum, 2n = 6x = 42, AABBDD) and Thinopyrum elongatum (2n = 2x = 14, EE) genomes can be differentiated from each other by fluorescent genomic in situ hybridization (FGISH) as well as molecular markers. This has facilitated homoeologous recombination-based partitioning and engineering of their genomes for physical mapping and alien introgression. Here, we constructed a special wheat genotype, which was double monosomic for wheat chromosome 3B and Th. elongatum chromosome 3E and homozygous for the ph1b mutant, to induce 3B-3E homoeologous recombination. Totally, 81 3B-3E recombinants were recovered and detected in the primary, secondary, and tertiary homoeologous recombination cycles by FGISH. Comparing to the primary recombination, the secondary and tertiary recombination shifted toward the proximal regions due to the increase in homology between the pairing partners. The 3B-3E recombinants were genotyped by high-throughput wheat 90-K single nucleotide polymorphism (SNP) arrays and their recombination breakpoints physically mapped based on the FGISH patterns and SNP results. The 3B-3E recombination physically partitioned chromosome 3B into 38 bins, and 429 SNPs were assigned to the distinct bins. Integrative analysis of FGISH and SNP results led to the construction of a composite bin map for chromosome 3B. Additionally, we developed 22 SNP-derived semi-thermal asymmetric reverse PCR markers specific for chromosome 3E and constructed a comparative map of homoeologous chromosomes 3E, 3B, 3A, and 3D. In summary, this work provides a unique physical framework for further studies of the 3B-3E homoeologous pair and diversifies the wheat genome for wheat improvement.

Pathway-specific enzymes from bamboo and crop leaves biosynthesize anti-nociceptive C-glycosylated flavones

C-glycosylated flavones (CGFs) are promising candidates as anti-nociceptive compounds. The leaves of bamboo and related crops in the grass family are a largely unexploited bioresource with a wide array of CGFs. We report here pathway-specific enzymes including C-glycosyltransferases (CGTs) and P450 hydroxylases from cereal crops and bamboo species accumulating abundant CGFs. Mining of CGTs and engineering of P450s that decorate the flavonoid skeleton allowed the production of desired CGFs (with yield of 20–40 mg/L) in an Escherichia coli cell factory. We further explored the antinociceptive activity of major CGFs in mice models and identified isoorientin as the most potent, with both neuroanalgesic and anti-inflammatory effects superior to clinical drugs such as rotundine and aspirin. Our discovery of the pain-alleviating flavonoids elicited from bamboo and crop leaves establishes this previously underutilized source, and sheds light on the pathway and pharmacological mechanisms of the compounds.

Yuwei Sun, Zhuo Chen, Jingya Yang et al. identify bamboo as a rich source of C-glycosylated flavonoids that reduces pain and inflammation. They identify isoorientin as the most potent C-glycosylated flavonoid, superior to aspirin, and report new enzymes that synthesize pain-alleviating C-glycosylated flavonoids.

The complete chloroplast genomes of seventeen Aegilops tauschii: genome comparative analysis and phylogenetic inference

The D genome progenitor of bread wheat, Aegilops tauschii Cosson (DD, 2n = 2x = 14), which is naturally distributed in Central Eurasia, ranging from northern Syria and Turkey to western China, is considered a potential genetic resource for improving bread wheat. In this study, the chloroplast (cp) genomes of 17 Ae. tauschii accessions were reconstructed. The cp genome sizes ranged from 135,551 bp to 136,009 bp and contained a typical quadripartite structure of angiosperms. Within these genomes, we identified a total of 124 functional genes, including 82 protein-coding genes, 34 transfer RNA genes and eight ribosomal RNA genes, with 17 duplicated genes in the IRs. Although the comparative analysis revealed that the genomic structure (gene order, gene number and IR/SC boundary regions) is conserved, a few variant loci were detected, predominantly in the non-coding regions (intergenic spacer regions). The phylogenetic relationships determined based on the complete genome sequences were consistent with the hypothesis that Ae. tauschii populations in the Yellow River region of China originated in South Asia not Xinjiang province or Iran, which could contribute to more effective utilization of wild germplasm resources. Furthermore, we confirmed that Ae. tauschii was derived from monophyletic speciation rather than hybrid speciation at the cp genome level. We also identified four variable genomic regions, rpl32-trnL-UAG, ccsA-ndhD, rbcL-psaI and rps18-rpl20, showing high levels of nucleotide polymorphisms, which may accordingly prove useful as cpDNA markers in studying the intraspecific genetic structure and diversity of Ae. tauschii.

Heteromorphic seeds of wheat wild relatives show germination niche differentiation

Crop wild relatives are fundamental genetic resources for crop improvement. Wheat wild relatives often produce heteromorphic seeds that differ in morphological and physiological traits. Several Aegilops and Triticum species possess, within the same spikelet, a dimorphic seed pair, with one seed being larger than the other. A comprehensive analysis is needed to understand which traits are involved in seed dimorphism and if these aspects of variation in dimorphic pairs are functionally related. To this end, dispersal units of Triticum urartu and five Aegilops species were X-rayed and the different seed morphs weighed. Germination tests were carried out on seeds, both dehulled and left in their dispersal units. Controlled ageing tests were performed to detect differences in seed longevity among seed morphs, and the antioxidant profile was assessed in terms of antioxidant compounds equipment and expression of selected antioxidant genes. We used PCA to group seed morphs sharing similar patterns of germination traits, longevity estimates and antioxidant profile. Different seed morphs differed significantly in terms of mass, final germination, germination timing, longevity estimates and antioxidant profile in most of the tested species. Small seeds germinated slower, had lower germination when left in their dispersal units, a higher antioxidant potential and were longer-lived than large seeds. The antioxidant gene expression varied between morphs, with different patterns across species but not clearly reflecting the phenotypic observations. The results highlight different trait trade-offs in dimorphic seeds of Aegilops and T. urartu, affecting their germination phenology and longevity, thereby resulting in recruitment niche differentiation.