Genome architecture of the allotetraploid wild grass Aegilops ventricosa reveals its evolutionary history and contributions to wheat improvement.

The allotetraploid wild grass Aegilops ventricosa (2n=4X=28, genome DvDvNvNv) has been recognized as an important germplasm resource for wheat improvement due to its ability to tolerate biotic stresses. Especially 2NvS segment from Aegilops ventricosa, as a stable and effective resistance source, has greatly contributed to wheat improvement. The 2NvS/2AS translocation is a prevalent chromosomal translocation between common wheat and wild relatives, ranking just behind the 1B/1R translocation in importance for modern wheat breeding. Here, we assembled a high-quality chromosome-level reference genome of Ae. ventricosa RM271 with a total length of 8.67 Gb. Phylogenomic analyses revealed that the progenitor of the Dv subgenome of Ae. ventricosa was Ae. tauschii ssp. tauschii (genome DD); in contrast, the progenitor of the D subgenome of bread wheat (Triticum aestivum L.) was Ae. tauschii ssp. strangulata (genome DD). The oldest polyploidization time of Ae. ventricosa occurred ∼0.7 million years ago. The Dv subgenome of Ae. ventricosa was less conserved than the D subgenome of bread wheat. Construction of a graph-based pangenome of 2AS/6NvL (originally known as 2NvS) segments from Ae. ventricosa and other genomes in the Triticeae enables us identifying candidate resistance genes sourced from Ae. ventricosa. We identified 12 nonredundant introgressed segments from the Dv and Nv subgenomes using a large winter wheat collection representing the full diversity of the wheat European genetic pool, and 29.40% of European wheat varieties inherited at least one of these segments. The high-quality RM271 reference genome will provide a basis for cloning key genes, including the Yr17-Lr37-Sr38-Cre5 resistance gene cluster in Ae. ventricosa, and facilitate the full use of elite wild genetic resources to accelerate wheat improvement.

Genome resources for the elite bread wheat cultivar Aikang 58 and mining of elite homeologous haplotypes for accelerating wheat improvement.

Despite recent progress in crop genomics studies, the genomic changes brought about by modern breeding selection are still poorly understood, thus hampering genomics-assisted breeding, especially in polyploid crops with compound genomes such as common wheat (Triticum aestivum). In this work, we constructed genome resources for the modern elite common wheat variety Aikang 58 (AK58). Comparative genomics between AK58 and the landrace cultivar Chinese Spring (CS) shed light on genomic changes that occurred through recent varietal improvement. We also explored subgenome diploidization and divergence in common wheat and developed a homoeologous locus-based genome-wide association study (HGWAS) approach, which was more effective than single homoeolog-based GWAS in unraveling agronomic trait-associated loci. A total of 123 major HGWAS loci were detected using a genetic population derived from AK58 and CS. Elite homoeologous haplotypes (HHs), formed by combinations of subgenomic homoeologs of the associated loci, were found in both parents and progeny, and many could substantially improve wheat yield and related traits. We built a website where users can download genome assembly sequence and annotation data for AK58, perform blast analysis, and run JBrowse. Our work enriches genome resources for wheat, provides new insights into genomic changes during modern wheat improvement, and suggests that efficient mining of elite HHs can make a substantial contribution to genomics-assisted breeding in common wheat and other polyploid crops.

Intra-Varietal Diversity and Its Contribution to Wheat Evolution, Domestication, and Improvement in Wheat.

Crop genetic diversity is essential for adaptation and productivity in agriculture. A previous study revealed that poor allele diversity in wheat commercial cultivars is a major barrier to its further improvement. Homologs within a variety, including paralogs and orthologs in polyploid, account for a large part of the total genes of a species. Homolog diversity, intra-varietal diversity (IVD), and their functions have not been elucidated. Common wheat, an important food crop, is a hexaploid species with three subgenomes. This study analyzed the sequence, expression, and functional diversity of homologous genes in common wheat based on high-quality reference genomes of two representative varieties, a modern commercial variety Aikang 58 (AK58) and a landrace Chinese Spring (CS). A total of 85,908 homologous genes, accounting for 71.9% of all wheat genes, including inparalogs (IPs), outparalogs (OPs), and single-copy orthologs (SORs), were identified, suggesting that homologs are an important part of the wheat genome. The levels of sequence, expression, and functional variation in OPs and SORs were higher than that of IPs, which indicates that polyploids have more homologous diversity than diploids. Expansion genes, a specific type of OPs, made a great contribution to crop evolution and adaptation and endowed crop with special characteristics. Almost all agronomically important genes were from OPs and SORs, demonstrating their essential functions for polyploid evolution, domestication, and improvement. Our results suggest that IVD analysis is a novel approach for evaluating intra-genomic variations, and exploitation of IVD might be a new road for plant breeding, especially for polyploid crops, such as 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.

Low-affinity SPL binding sites contribute to subgenome expression divergence in allohexaploid wheat.

Expression divergence caused by genetic variation and crosstalks among subgenomes of the allohexaploid bread wheat (Triticum aestivum. L., BBAADD) is hypothesized to increase its adaptability and/or plasticity. However, the molecular basis of expression divergence remains unclear. Squamosa promoter-binding protein-like (SPL) transcription factors are critical for a wide array of biological processes. In this study, we constructed expression regulatory networks by combining DAP-seq for 40 SPLs, ATAC-seq, and RNA-seq. Our findings indicate that a group of low-affinity SPL binding regions (SBRs) were targeted by diverse SPLs and caused different sequence preferences around the core GTAC motif. The SBRs including the low-affinity ones are evolutionarily conserved, enriched GWAS signals related to important agricultural traits. However, those SBRs are highly diversified among the cis-regulatory regions (CREs) of syntenic genes, with less than 8% SBRs coexisting in triad genes, suggesting that CRE variations are critical for subgenome differentiations. Knocking out of TaSPL7A/B/D and TaSPL15A/B/D subfamily further proved that both high- and low-affinity SBRs played critical roles in the differential expression of genes regulating tiller number and spike sizes. Our results have provided baseline data for downstream networks of SPLs and wheat improvements and revealed that CRE variations are critical sources for subgenome divergence in the allohexaploid wheat.

Reference genome assemblies reveal the origin and evolution of allohexaploid oat.

Common oat (Avena sativa) is an important cereal crop serving as a valuable source of forage and human food. Although reference genomes of many important crops have been generated, such work in oat has lagged behind, primarily owing to its large, repeat-rich polyploid genome. Here, using Oxford Nanopore ultralong sequencing and Hi-C technologies, we have generated a reference-quality genome assembly of hulless common oat, comprising 21 pseudomolecules with a total length of 10.76 Gb and contig N50 of 75.27 Mb. We also produced genome assemblies for diploid and tetraploid Avena ancestors, which enabled the identification of oat subgenomes and provided insights into oat chromosomal evolution. The origin of hexaploid oat is inferred from whole-genome sequencing, chloroplast genomes and transcriptome assemblies of different Avena species. These findings and the high-quality reference genomes presented here will facilitate the full use of crop genetic resources to accelerate oat improvement.

A recent burst of gene duplications in Triticeae.

Gene duplication provides raw genetic materials for evolution and potentially novel genes for crop improvement. The two seminal genomic studies of Aegilops tauschii both mentioned the large number of genes independently duplicated in recent years, but the duplication mechanism and the evolutionary significance of these gene duplicates have not yet been investigated. Here, we found that a recent burst of gene duplications (hereafter abbreviated as the RBGD) has probably occurred in all sequenced Triticeae species. Further investigations of the characteristics of the gene duplicates and their flanking sequences suggested that transposable element (TE) activity may have been involved in generating the RBGD. We also characterized the duplication timing, retention pattern, diversification, and expression of the duplicates following the evolution of Triticeae. Multiple subgenome-specific comparisons of the duplicated gene pairs clearly supported extensive differential regulation and related functional diversity among such pairs in the three subgenomes of bread wheat. Moreover, several duplicated genes from the RBGD have evolved into key factors that influence important agronomic traits of wheat. Our results provide insights into a unique source of gene duplicates in Triticeae species, which has increased the gene dosage together with the two polyploidization events in the evolutionary history of wheat.

Homology-mediated inter-chromosomal interactions in hexaploid wheat lead to specific subgenome territories following polyploidization and introgression.

BACKGROUND: Polyploidization and introgression are major events driving plant genome evolution and influencing crop breeding. However, the mechanisms underlying the higher-order chromatin organization of subgenomes and alien chromosomes are largely unknown. RESULTS: We probe the three-dimensional chromatin architecture of Aikang 58 (AK58), a widely cultivated allohexaploid wheat variety in China carrying the 1RS/1BL translocation chromosome. The regions involved in inter-chromosomal interactions, both within and between subgenomes, have highly similar sequences. Subgenome-specific territories tend to be connected by subgenome-dominant homologous transposable elements (TEs). The alien 1RS chromosomal arm, which was introgressed from rye and differs from its wheat counterpart, has relatively few inter-chromosome interactions with wheat chromosomes. An analysis of local chromatin structures reveals topologically associating domain (TAD)-like regions covering 52% of the AK58 genome, the boundaries of which are enriched with active genes, zinc-finger factor-binding motifs, CHH methylation, and 24-nt small RNAs. The chromatin loops are mostly localized around TAD boundaries, and the number of gene loops is positively associated with gene activity. CONCLUSIONS: The present study reveals the impact of the genetic sequence context on the higher-order chromatin structure and subgenome stability in hexaploid wheat. Specifically, we characterized the sequence homology-mediated inter-chromosome interactions and the non-canonical role of subgenome-biased TEs. Our findings may have profound implications for future investigations of the interplay between genetic sequences and higher-order structures and their consequences on polyploid genome evolution and introgression-based breeding of crop plants.

A Novel Wheat C-bZIP Gene, TabZIP14-B, Participates in Salt and Freezing Tolerance in Transgenic Plants.

The group C-bZIP transcription factors (TFs) are involved in diverse biological processes, such as the regulation of seed storage protein (SSP) production and the responses to pathogen challenge and abiotic stress. However, our knowledge of the abiotic functions of group C-bZIP genes in wheat remains limited. Here, we present the function of a novel TabZIP14-B gene in wheat. This gene belongs to the group C-bZIP TFs and contains six exons and five introns; three haplotypes were identified among accessions of tetraploid and hexaploid wheat. A subcellular localization analysis indicated that TabZIP14-B was targeted to the nucleus of tobacco epidermal cells. A transactivation assay demonstrated that TabZIP14-B showed transcriptional activation ability and was capable of binding the abscisic acid (ABA) responsive element (ABRE) in yeast. RT-qPCR revealed that TabZIP14-B was expressed in the roots, stems, leaves, and young spikes and was up-regulated by exogenous ABA, salt, low-temperature, and polyethylene glycol (PEG) stress treatments. Furthermore, Arabidopsis plants overexpressing TabZIP14-B exhibited enhanced tolerance to salt, freezing stresses and ABA sensitivity. Overexpression of TabZIP14-B resulted in increased expression of the AtRD29A, AtCOR47, AtRD20, AtGSTF6, and AtRAB18 genes and changes in several physiological characteristics. These results suggest that TabZIP14-B could function as a positive regulator in mediating the abiotic stress response.

Candidate loci involved in domestication and improvement detected by a published 90K wheat SNP array.

Selection is one of the most important forces in crop evolution. Common wheat is a major world food crop and a typical allopolyploid with a huge and complex genome. We applied four approaches to detect loci selected in wheat during domestication and improvement. A total of 7,984 candidate loci were detected, accounting for 23.3% of all 34,317 SNPs analysed, a much higher proportion than estimated in previous reports. We constructed a first generation wheat selection map which revealed the following new insights on genome-wide selection: (1) diversifying selection acted by increasing, decreasing or not affecting gene frequencies; (2) the number of loci under selection during domestication was much higher than that during improvement; (3) the contribution to wheat improvement by the D sub-genome was relatively small due to the bottleneck of hexaploidisation and diversity can be expanded by using synthetic wheat and introgression lines; and (4) clustered selection regions occur throughout the wheat genome, including the centromere regions. This study will not only help future wheat breeding and evolutionary studies, but will also accelerate study of other crops, especially polyploids.

An assessment of heavy ion irradiation mutagenesis for reverse genetics in wheat (Triticum aestivum L.).

Reverse genetic techniques harnessing mutational approaches are powerful tools that can provide substantial insight into gene function in plants. However, as compared to diploid species, reverse genetic analyses in polyploid plants such as bread wheat can present substantial challenges associated with high levels of sequence and functional similarity amongst homoeologous loci. We previously developed a high-throughput method to identify deletions of genes within a physically mutagenized wheat population. Here we describe our efforts to combine multiple homoeologous deletions of three candidate disease susceptibility genes (TaWRKY11, TaPFT1 and TaPLDß1). We were able to produce lines featuring homozygous deletions at two of the three homoeoloci for all genes, but this was dependent on the individual mutants used in crossing. Intriguingly, despite extensive efforts, viable lines possessing homozygous deletions at all three homoeoloci could not be produced for any of the candidate genes. To investigate deletion size as a possible reason for this phenomenon, we developed an amplicon sequencing approach based on synteny to Brachypodium distachyon to assess the size of the deletions removing one candidate gene (TaPFT1) in our mutants. These analyses revealed that genomic deletions removing the locus are relatively large, resulting in the loss of multiple additional genes. The implications of this work for the use of heavy ion mutagenesis for reverse genetic analyses in wheat are discussed.

Dynamic evolution of rht-1 homologous regions in grass genomes.

Hexaploid bread wheat contains A, B, and D three subgenomes with its well-characterized ancestral genomes existed at diploid and tetraploid levels, making the wheat act as a good model species for studying evolutionary genomic dynamics. Here, we performed intra- and inter-species comparative analyses of wheat and related grass genomes to examine the dynamics of homologous regions surrounding Rht-1, a well-known "green revolution" gene. Our results showed that the divergence of the two A genomes in the Rht-1 region from the diploid and tetraploid species is greater than that from the tetraploid and hexaploid wheat. The divergence of D genome between diploid and hexaploid is lower than those of A genome, suggesting that D genome diverged latter than others. The divergence among the A, B and D subgenomes was larger than that among different ploidy levels for each subgenome which mainly resulted from genomic structural variation of insertions and, perhaps deletions, of the repetitive sequences. Meanwhile, the repetitive sequences caused genome expansion further after the divergence of the three subgenomes. However, several conserved non-coding sequences were identified to be shared among the three subgenomes of wheat, suggesting that they may have played an important role to maintain the homolog of three subgenomes. This is a pilot study on evolutionary dynamics across the wheat ploids, subgenomes and differently related grasses. Our results gained new insights into evolutionary dynamics of Rht-1 region at sequence level as well as the evolution of wheat during the plolyploidization process.

Aegilops tauschii draft genome sequence reveals a gene repertoire for wheat adaptation.

About 8,000 years ago in the Fertile Crescent, a spontaneous hybridization of the wild diploid grass Aegilops tauschii (2n = 14; DD) with the cultivated tetraploid wheat Triticum turgidum (2n = 4x = 28; AABB) resulted in hexaploid wheat (T. aestivum; 2n = 6x = 42; AABBDD). Wheat has since become a primary staple crop worldwide as a result of its enhanced adaptability to a wide range of climates and improved grain quality for the production of baker's flour. Here we describe sequencing the Ae. tauschii genome and obtaining a roughly 90-fold depth of short reads from libraries with various insert sizes, to gain a better understanding of this genetically complex plant. The assembled scaffolds represented 83.4% of the genome, of which 65.9% comprised transposable elements. We generated comprehensive RNA-Seq data and used it to identify 43,150 protein-coding genes, of which 30,697 (71.1%) were uniquely anchored to chromosomes with an integrated high-density genetic map. Whole-genome analysis revealed gene family expansion in Ae. tauschii of agronomically relevant gene families that were associated with disease resistance, abiotic stress tolerance and grain quality. This draft genome sequence provides insight into the environmental adaptation of bread wheat and can aid in defining the large and complicated genomes of wheat species.

TaCKX6-D1, the ortholog of rice OsCKX2, is associated with grain weight in hexaploid wheat.

The cytokinin oxidase/dehydrogenase (CKX) gene plays a principal role in controlling cytokinin levels and has been shown to be a major quantitative trait locus (QTL) affecting grain number in rice. However, the function and evaluation of the haplotypes of the wheat CKX gene have yet to be illustrated. In this study, TaCKX6-D1, a wheat ortholog of rice OsCKX2, was cloned and its haplotype variants were determined to be significantly associated with the 1000-grain weight on the basis of linkage mapping, association analysis and gene expression analysis. Five TaCKX6-D1 haplotypes, designated a-e, were identified. An indel marker was developed to identify haplotype a, which was associated with higher grain weight. Haplotype a showed decreased expression relative to haplotype b in seeds at 8 d after pollination. Sequence variations among modern cultivars, landraces and wild species suggest a significant domestication signature at the TaCKX6-D1 locus in Chinese wheat germplasm. TaCKX6-D1 may serve as a useful gene for the breeding of high-yielding wheat. A strategy for allele mining and utilization of TaCKX6-D1 was proposed. Our study also sheds light on the mechanisms of grain development and domestication of wheat, as well as the functional divergence of orthologs in comparative genomics.

Molecular characterization of 60 isolated wheat MYB genes and analysis of their expression during abiotic stress.

The proteins of the MYB superfamily play central roles in developmental processes and defence responses in plants. Sixty unique wheat MYB genes that contain full-length cDNA sequences were isolated. These 60 genes were grouped into three categories, namely one R1R2R3-MYB, 22 R2R3-MYBs, and 37 MYB-related members. The sequence composition of the R2 and R3 repeats was conserved among the 22 wheat R2R3-MYB proteins. Phylogenetic comparison of the members of this superfamily among wheat, rice, and Arabidopsis revealed that the putative functions of some wheat MYB proteins were clustered into the Arabidopsis functional clades. Tissue-specific expression profiles showed that most of the wheat MYB genes were expressed in all of the tissues examined, suggesting that wheat MYB genes take part in multiple cellular processes. The expression analysis during abiotic stress identified a group of MYB genes that respond to one or more stress treatments. The overexpression of a salt-inducible gene, TaMYB32, enhanced the tolerance to salt stress in transgenic Arabidopsis. This study is the first comprehensive study of the MYB gene family in Triticeae.

Discovery of quantitative trait loci for crossability from a synthetic wheat genotype.

Crossability between wheat and rye is an important trait for wheat improvement. No quantitative trait loci (QTLs) were detected from wheat ancestors previously. The objectives of this study were to dissect the QTLs for crossability using 111 introgression lines (ILs) derived from synthetic hexaploid wheat. A total of 1275 SSR markers were screened for polymorphism between the two parents, and 552 markers of them displayed polymorphism, of which 64 were selected for genotyping the 111 BC(5)F(6) ILs. Field trials were performed in a Latinized α-lattice design in Luoyang and Jiaozuo of Henan Province of China in 2007-2008 and 2008-2009 cropping seasons. One-way ANOVA and interval mapping (IM) analysis were used to detect QTL for crossability between wheat and rye. A total of 13 putative QTLs were detected. Five of them, QCa.caas.1A, QCa.caas.2D, QCa.caas.4B, QCa.caas.5B and QCa.caas.6A, were detected in both trials and three of them, QCa.caas.2D, QCa.caas.4B and QCa.caas.6A, were novel. The positive effect allele of the four QTLs came from the donor parent Am3 except QCa.caas.6A that came from the recurrent parent Laizhou953. ILs with both higher positive effect alleles and favorable agronomic traits developed in present study are elite germplasm for wide crossing in wheat. Results from the current study suggest that wheat ancestors can be rich in new sources of crossability genes.

Rapid evolution and complex structural organization in genomic regions harboring multiple prolamin genes in the polyploid wheat genome.

Genes encoding wheat prolamins belong to complicated multi-gene families in the wheat genome. To understand the structural complexity of storage protein loci, we sequenced and analyzed orthologous regions containing both gliadin and LMW-glutenin genes from the A and B genomes of a tetraploid wheat species, Triticum turgidum ssp. durum. Despite their physical proximity to one another, the gliadin genes and LMW-glutenin genes are organized quite differently. The gliadin genes are found to be more clustered than the LMW-glutenin genes which are separated from each other by much larger distances. The separation of the LMW-glutenin genes is the result of both the insertion of large blocks of repetitive DNA owing to the rapid amplification of retrotransposons and the presence of genetic loci interspersed between them. Sequence comparisons of the orthologous regions reveal that gene movement could be one of the major factors contributing to the violation of microcolinearity between the homoeologous A and B genomes in wheat. The rapid sequence rearrangements and differential insertion of repetitive DNA has caused the gene islands to be not conserved in compared regions. In addition, we demonstrated that the i-type LMW-glutenin originated from a deletion of 33-bps in the 5' coding region of the m-type gene. Our results show that multiple rounds of segmental duplication of prolamin genes have driven the amplification of the omega-gliadin genes in the region; such segmental duplication could greatly increase the repetitive DNA content in the genome depending on the amount of repetitive DNA present in the original duplicate region.

Development, utilization of introgression lines using a synthetic wheat as donor.

A series of introgression lines (ILs) were generated from repeated backcrossing between the exotic hexaploid wheat genotype Am3 and the common wheat genotype Laizhou953. Am3 was synthesized by crossing Triticum carthlicum with Aegilops tauschii and was used as the donor parent in this study, and Laizhou953 was used as the recurrent parent. Two hundred and five SSR markers showing polymorphism between the two parents were used to identify the introgressed Am3 chromosome segments in 97 BC4F3 ILs. The introgressed segments in each line and the length of the introgressed segments were estimated according to the wheat SSR consensus map. The introgressed segments from Am3 in the 97 lines covers 37.7% of the donor genome. The introgressed segments were most found on 2D, 3B, 6B, and 1D with coverage of 59.8, 59.5, 59.1, and 59% of the chromosomes, respectively. None of the 97 lines tested contained chromosome 4D segments introgressed from Am3. Introgressed segments for each of the chromosomes were mapped using the consensus wheat linkage map. Nine agronomic traits from BC4F3 lines were evaluated and the phenotype showed most lines have the tendency to be more similar to the recurrent parent. There were lines showing better agronomic traits than the recurrent parent, which indicated the introgression of favorable alleles from the exotic hexaploid wheat into the elite cultivar Laizhou953. Marker and phenotype data were used to identify quantitative trait loci (QTLs) controlling these nine traits. In total, 38, 33, and 28 putative QTLs were detected for seven of the nine traits in 2003, 2004, and 2005, respectively. Some of these agronomic important QTLs were detected in more than one season.