Triticum population sequencing provides insights into wheat adaptation

Chromosome painting reveals inter-chromosomal rearrangements and evolution of subgenome D of wheat

Aegilops species represent the most important gene pool for breeding bread wheat (Triticum aestivum). Thus, understanding the genome evolution, including chromosomal structural rearrangements and syntenic relationships among Aegilops species or between Aegilops and wheat, is important for both basic genome research and practical breeding applications. In the present study, we attempted to develop subgenome D-specific fluorescence in situ hybridization (FISH) probes by selecting D-specific oligonucleotides based on the reference genome of Chinese Spring. The oligo-based chromosome painting probes consisted of approximately 26 000 oligos per chromosome and their specificity was confirmed in both diploid and polyploid species containing the D subgenome. Two previously reported translocations involving two D chromosomes have been confirmed in wheat varieties and their derived lines. We demonstrate that the oligo painting probes can be used not only to identify the translocations involving D subgenome chromosomes, but also to determine the precise positions of chromosomal breakpoints. Chromosome painting of 56 accessions of Ae. tauschii from different origins led us to identify two novel translocations: a reciprocal 3D-7D translocation in two accessions and a complex 4D-5D-7D translocation in one accession. Painting probes were also used to analyze chromosomes from more diverse Aegilops species. These probes produced FISH signals in four different genomes. Chromosome rearrangements were identified in Aegilops umbellulata, Aegilops markgrafii, and Aegilops uniaristata, thus providing syntenic information that will be valuable for the application of these wild species in wheat breeding.

Population transcriptomic analysis identifies the comprehensive lncRNAs landscape of spike in wheat (Triticum aestivum L.)

BACKGROUND

Long noncoding RNAs (lncRNAs) are emerging as the important regulators involving in growth and development as well as stress response in plants. However, current lncRNA studies were mainly performed at the individual level and the significance of it is not well understood in wheat.

RESULTS

In this study, the lncRNA landscape of wheat spike was characterized through analysing a total of 186 spike RNA-seq datasets from 93 wheat genotypes. A total of 35,913 lncRNAs as well as 1,619 lncRNA-mRNA pairs comprised of 443 lncRNAs and 464 mRNAs were obtained. Compared to coding genes, these lncRNAs displayed rather low conservation among wheat and other gramineous species. Based on re-sequencing data, the genetic variations of these lncRNA were investigated and obvious genetic bottleneck were found on them during wheat domestication process. Furthermore, 122 lncRNAs were found to act as ceRNA to regulate endogenous competition. Finally, association and co-localization analysis of the candidate lncRNA-mRNA pairs identified 170 lncRNAs and 167 target mRNAs significantly associated with spike-related traits, including lncRNA.127690.1/TraesCS2A02G518500.1 (PMEI) and lncRNA.104854.1/TraesCS6A02G050300.1 (ATG5) associated with heading date and spike length, respectively.

CONCLUSIONS

This study reported the lncRNA landscape of wheat spike through the population transcriptome analysis, which not only contribute to better understand the wheat evolution from the perspective of lncRNA, but also lay the foundation for revealing roles of lncRNA playing in spike development.

Wheat genomic study for genetic improvement of traits in China

Bread wheat (Triticum aestivum L.) is a major crop that feeds 40% of the world's population. Over the past several decades, advances in genomics have led to tremendous achievements in understanding the origin and domestication of wheat, and the genetic basis of agronomically important traits, which promote the breeding of elite varieties. In this review, we focus on progress that has been made in genomic research and genetic improvement of traits such as grain yield, end-use traits, flowering regulation, nutrient use efficiency, and biotic and abiotic stress responses, and various breeding strategies that contributed mainly by Chinese scientists. Functional genomic research in wheat is entering a new era with the availability of multiple reference wheat genome assemblies and the development of cutting-edge technologies such as precise genome editing tools, high-throughput phenotyping platforms, sequencing-based cloning strategies, high-efficiency genetic transformation systems, and speed-breeding facilities. These insights will further extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture in China and throughout the world.

Evolution of the DEHYDRATION-RESPONSIVE ELEMENT-BINDING PROTEIN subfamily in green plants

Adapting to unfavorable environments is a necessary step in plant terrestrialization and radiation. The dehydration-responsive element-binding (DREB) protein subfamily plays a pivotal role in plant abiotic stress regulation. However, relationships between the origin and expansion of the DREB subfamily and adaptive evolution of land plants are still being elucidated. Here, we constructed the evolutionary history of the DREB subfamily by compiling APETALA2/ethylene-responsive element-binding protein superfamily genes from 169 representative species of green plants. Through extensive phylogenetic analyses and comparative genomic analysis, our results revealed that the DREB subfamily diverged from the ethylene-responsive factor (ERF) subfamily in the common ancestor of Zygnemophyceae and Embryophyta during the colonization of land by plants, followed by expansions to form three different ancient archetypal genes in Zygnemophyceae species, designated as groups archetype-I, archetype-II/III, and archetype-IV. Four large-scale expansions paralleling the evolution of land plants led to the nine-subgroup divergence of group archetype-II/III in angiosperms, and five whole-genome duplications during Brassicaceae and Poaceae radiation shaped the diversity of subgroup IIb-1. We identified a Poaceae-specific gene in subgroup IIb-1, ERF014, remaining in a Poaceae-specific microsynteny block and co-evolving with a small heat shock protein cluster. Expression analyses demonstrated that heat acclimation may have driven the neofunctionalization of ERF014s in Pooideae by engaging in the conserved heat-responsive module in Poaceae. This study provides insights into lineage-specific expansion and neofunctionalization in the DREB subfamily, together with evolutionary information valuable for future functional studies of plant stress biology.

Improving bread wheat yield through modulating an unselected AP2/ERF gene

Crop breeding heavily relies on natural genetic variation. However, additional new variations are desired to meet the increasing human demand. Inflorescence architecture determines grain number per spike, a major determinant of bread wheat (Triticum aestivum L.) yield. Here, using Brachypodium distachyon as a wheat proxy, we identified DUO-B1, encoding an APETALA2/ethylene response factor (AP2/ERF) transcription factor, regulating spike inflorescence architecture in bread wheat. Mutations of DUO-B1 lead to mild supernumerary spikelets, increased grain number per spike and, importantly, increased yield under field conditions without affecting other major agronomic traits. DUO-B1 suppresses cell division and promotes the expression of BH^t/WFZP, whose mutations could lead to branched 'miracle-wheat'. Pan-genome analysis indicated that DUO-B1 has not been utilized in breeding, and holds promise to increase wheat yield further.

Genome-Wide Association Study of Leaf Rust and Stem Rust Seedling and Adult Resistances in Tetraploid Wheat Accessions Harvested in Kazakhstan

Leaf rust (LR) and stem rust (SR) are diseases increasingly impacting wheat production worldwide. Fungal pathogens producing rust diseases in wheat may cause yield losses of up to 50−60%. One of the most effective methods for preventing such losses is the development of resistant cultivars with high yield potential. This goal can be achieved through complex breeding studies, including the identification of key genetic factors controlling rust disease resistance. The objective of this study was to identify sources of tetraploid wheat resistance to LR and SR races, both at the seedling growth stage in the greenhouse and at the adult plant stage in field experiments, under the conditions of the North Kazakhstan region. A panel consisting of 193 tetraploid wheat accessions was used in a genome-wide association study (GWAS) for the identification of quantitative trait loci (QTLs) associated with LR and SR resistance, using 16,425 polymorphic single-nucleotide polymorphism (SNP) markers in the seedling and adult stages of plant development. The investigated panel consisted of seven tetraploid subspecies (Triticum turgidum ssp. durum, ssp. turanicum, ssp. turgidum, ssp. polonicum, ssp. carthlicum, ssp. dicoccum, and ssp. dicoccoides). The GWAS, based on the phenotypic evaluation of the tetraploid collection’s reaction to the two rust species at the seedling (in the greenhouse) and adult (in the field) stages, revealed 38 QTLs (p < 0.001), comprising 17 for LR resistance and 21 for SR resistance. Ten QTLs were associated with the reaction to LR at the seedling stage, while six QTLs were at the adult plant stage and one QTL was at both the seedling and adult stages. Eleven QTLs were associated with SR response at the seedling stage, while nine QTLs were at the adult plant stage and one QTL was at both the seedling and adult stages. A comparison of these results with previous LR and SR studies indicated that 11 of the 38 QTLs are presumably novel loci. The QTLs identified in this work can potentially be used for marker-assisted selection of tetraploid and hexaploid wheat for the breeding of new LR- and SR-resistant cultivars.

Identification of ankyrin-transmembrane-type subfamily genes in Triticeae species reveals TaANKTM2A-5 regulates powdery mildew resistance in wheat

The ankyrin-transmembrane (ANKTM) subfamily is the most abundant subgroup of the ANK superfamily, with critical roles in pathogen defense. However, the function of ANKTM proteins in wheat immunity remains largely unexplored. Here, a total of 381 ANKTMs were identified from five Triticeae species and Arabidopsis, constituting five classes. Among them, class a only contains proteins from Triticeae species and the number of ANKTM in class a of wheat is significantly larger than expected, even after consideration of the ploidy level. Tandem duplication analysis of ANKTM indicates that Triticum urartu, Triticum dicoccoides and wheat all had experienced tandem duplication events which in wheat-produced ANKTM genes all clustered in class a. The above suggests that not only did the genome polyploidization result in the increase of ANKTM gene number, but that tandem duplication is also a mechanism for the expansion of this subfamily. Micro-collinearity analysis of Triticeae ANKTMs indicates that some ANKTM type genes evolved into other types of ANKs in the evolution process. Public RNA-seq data showed that most of the genes in class d and class e are expressed, and some of them show differential responses to biotic stresses. Furthermore, qRT-PCR results showed that some ANKTMs in class d and class e responded to powdery mildew. Silencing of TaANKTM2A-5 by barley stripe mosaic virus-induced gene silencing compromised powdery mildew resistance in common wheat Bainongaikang58. Findings in this study not only help to understand the evolutionary process of ANKTM genes, but also form the basis for exploring disease resistance genes in the ANKTM gene family.

Alternative modes of introgression-mediated selection shaped crop adaptation to novel climates

Recent plant genomic studies provide fine-grained details on the evolutionary consequences of adaptive introgression during crop domestication. Modern genomic approaches and analytical methods now make it possible to better separate the introgression signal from the demographic signal thus providing a more comprehensive and complex picture of the role of introgression in local adaptation. Adaptive introgression has been fundamental for crop expansion and has involved complex patterns of gene flow. In addition to providing new and more favorable alleles of large effect, introgression during the early stages of domestication also increased allelic diversity at adaptive loci. Previous studies have largely underestimated the effect of such increased diversity following introgression. Recent genomic studies in wheat, potato, maize, grapevine, and ryegrass show that introgression of multiple genes, of as yet unknown effect, increased the effectiveness of purifying selection, and promoted disruptive or fluctuating selection in early cultivars and landraces. Historical selection processes associated with introgression from crop wild relatives provide an instructive analog for adaptation to current climate change and offer new avenues for crop breeding research that are expected to be instrumental for strengthening food security in the coming years.

Genome-Wide Identification and Characterization of RNA/DNA Differences Associated with Fusarium graminearum Infection in Wheat

RNA/DNA difference (RDD) is a post-transcriptional modification playing a crucial role in regulating diverse biological processes in eukaryotes. Although it has been extensively studied in plant chloroplast and mitochondria genomes, RDDs in plant nuclear genomes are not well studied at present. Here, we investigated the RDDs associated with fusarium head blight (FHB) through a novel method by comparing the RNA-seq data between Fusarium-infected and control samples of four wheat genotypes. A total of 187 high-confidence unique RDDs in 36 genes were identified, representing the first landscape of the FHB-responsive RDD in wheat. The majority (26) of these 36 RDD genes were correlated either positively or negatively with FHB levels. Effects of these RDDs on RNA and protein sequences have been identified, their editing frequency and the expression level of the corresponding genes provided, and the prediction of the effect on the minimum folding free energy of mRNA, miRNA binding, and colocation of RDDs with conserved domains presented. RDDs were predicted to induce modifications in the mRNA and protein structures of the corresponding genes. In two genes, TraesCS1B02G294300 and TraesCS3A02G263900, editing was predicted to enhance their affinity with tae-miR9661-5p and tae-miR9664-3p, respectively. To our knowledge, this study is the first report of the association between RDD and FHB in wheat; this will contribute to a better understanding of the molecular basis underlying FHB resistance, and potentially lead to novel strategies to improve wheat FHB resistance through epigenetic methods.

Population Genomics and Haplotype Analysis in Bread Wheat Identify a Gene Regulating Glume Pubescence

Glume hairiness or pubescence is an important morphological trait with high heritability to distinguish/characterize wheat and is related to the resistance to biotic and abiotic stresses. Hg1 (formerly named Hg) on chromosome arm 1AS controlled glume hairiness in wheat. Its genetic analysis and mapping have been widely studied, yet more useful and accurate information for fine mapping of Hg1 and identification of its candidate gene is lacking. The cloning of this gene has not yet been reported for the large complex wheat genome. Here, we performed a GWAS between SNP markers and glume pubescence (Gp) in a wheat population with 352 lines and further demonstrated the gene expression and haplotype analysis approach for isolating the Hg1 gene. One gene, TraesCSU02G143200 (TaELD1-1A), encoding glycosyltransferase-like ELD1/KOBITO 1, was identified as the most promising candidate gene of Hg1. The gene annotation, expression pattern, function SNP variation, haplotype analysis, and co-expression analysis in floral organ (spike) development indicated that it is likely to be involved in the regulation of glume pubescence. Our study demonstrates the importance of high-quality reference genomes and annotation information, as well as bioinformatics analysis, for gene cloning in wheat.

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

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

Genetic Diversity of Transcription Factor Genes in Triticum and Mining for Promising Haplotypes for Beneficial Agronomic Traits

Transcription factor (TF) is a class of the sequence-specific DNA-binding proteins that modulate the transcription of target genes, and thus regulate their expressions. Variations in TF are the crucial determinants for phenotypic traits. Although much progress has been made in the functions of TF genes in wheat, one of the most important staple crops globally, the diversity of TF genes in wheat and its progenitors are not well understood, especially the agronomically promising haplotypes have not yet been characterized. Here, we identified a total of 6,023 TF genes from hexaploid wheat through a genome-search method and classified them into 59 gene families based on the conserved domain. The characteristics and dN/dS values of these genes showed evidently selective effects. Based on re-sequencing data, we found a strong genetic bottleneck among these TF genes on A and D subgenomes while no found in B subgenome during wheat domestication. Combined with selective signals and known QTLs on the whole genome, 21 TF genes were preliminarily found to be associated with yield-related traits. The haplotype frequency of these TF genes was further investigated in bread wheat and its progenitors and 13 major haplotypes were the casual loci related to key traits. Finally, the tissue-specific TF genes were also identified using RNA-seq analysis. This study provided insights into the diversity and evolution of TF genes and the identified TF genes and excellent haplotypes associating with traits will contribute to wheat genetic improvement.

Dispersed emergence and protracted domestication of polyploid wheat uncovered by mosaic ancestral haploblock inference

Major crops are all survivors of domestication bottlenecks. Studies have focused on the genetic loci related to the domestication syndrome, while the contribution of ancient haplotypes remains largely unknown. Here, an ancestral genomic haploblock dissection method is developed and applied to a resequencing dataset of 386 tetraploid/hexaploid wheat accessions, generating a pan-ancestry haploblock map. Together with cytoplastic evidences, we reveal that domesticated polyploid wheat emerged from the admixture of six founder wild emmer lineages, which contributed the foundation of ancestral mosaics. The key domestication-related loci, originated over a wide geographical range, were gradually pyramided through a protracted process. Diverse stable-inheritance ancestral haplotype groups of the chromosome central zone are identified, revealing the expanding routes of wheat and the trends of modern wheat breeding. Finally, an evolution model of polyploid wheat is proposed, highlighting the key role of wild-to-crop and interploidy introgression, that increased genomic diversity following bottlenecks introduced by domestication and polyploidization.

Evolution and origin of bread wheat

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

HAPPE: A Tool for Population Haplotype Analysis and Visualization in Editable Excel Tables

Haplotype identification, characterization and visualization are important for large-scale analysis and use in population genomics. Many tools have been developed to visualize haplotypes, but it is challenging to display both the pattern of haplotypes and the genotypes for each single SNP in the context of a large amount of genomic data. Here, we describe the tool HAPPE, which uses the agglomerative hierarchical clustering algorithm to characterize and visualize the genotypes and haplotypes in a phylogenetic context. The tool displays the plots by coloring the cells and/or their borders in Excel tables for any given gene and genomic region of interest. HAPPE facilitates informative displays wherein data in plots are easy to read and access. It allows parallel display of several lines of values, such as phylogenetic trees, P values of GWAS, the entry of genes or SNPs, and the sequencing depth at each position. These features are informative for the detection of insertion/deletions or copy number variations. Overall, HAPPE provides editable plots consisting of cells in Excel tables, which are user-friendly to non-programmers. This pipeline is coded in Python and is available at https://github.com/fengcong3/HAPPE.

The Triticum ispahanicum elongated glume locus P2 maps to chromosome 6A and is associated with the ectopic expression of SVP-A1

We propose the MADS-box transcription factor SVP-A1 as a promising candidate gene for the elongated glume locus P2, which maps to chromosome 6A instead of the previously proposed chromosome 7B. In rice and wheat, glume and floral organ length are positively correlated with grain size, making them an important target to increase grain size and potentially yield. The wheat subspecies Triticum ispahanicum is known to develop elongated glumes and floral organs as well as long grains. These multiple phenotypic effects are controlled by the P2 locus, which was previously mapped to wheat chromosome 7B. Using three mapping populations, we show that the long glume locus P2 does not map to chromosome 7B, but instead maps to a 1.68 Mbp interval on chromosome 6A. Within this interval, we identified SVP-A1, a MADS box transcription factor which is the direct ortholog of the maize gene underlying the 'pod corn' Tunicate locus and is a paralog to the T. polonicum elongated glume P1 gene. In T. ispahanicum, we identified a unique allele which has a 482-bp deletion in the SVP-A1 promoter and is associated with ectopic and higher expression of SVP-A1 in the elongated glumes and floral organs. We used near-isogenic lines (NILs) to show that P2 has a consistent positive effect on the length of glume, lemma, palea, spike and grain. Based on the mapping data, natural variation, biological function of SVP genes in cereals and expression analyses, we propose the MADS-box transcription factor SVP-A1 as a promising candidate for P2.

Pentaploidization Enriches the Genetic Diversity of Wheat by Enhancing the Recombination of AB Genomes

Allohexaploidization and continuous introgression play a key role in the origin and evolution of bread wheat. The genetic bottleneck of bread wheat resulting from limited germplasms involved in the origin and modern breeding may be compensated by gene flow from tetraploid wheat through introgressive hybridization. The inter-ploidy hybridization between hexaploid and tetraploid wheat generates pentaploid hybrids first, which absorbed genetic variations both from hexaploid and tetraploid wheat and have great potential for re-evolution and improvement in bread wheat. Therefore, understanding the effects of the pentaploid hybrid is of apparent significance in our understanding of the historic introgression and in informing breeding. In the current study, two sets of F2 populations of synthetic pentaploid wheat (SPW1 and SPW2) and synthetic hexaploid wheat (SHW1 and SHW2) were created to analyze differences in recombination frequency (RF) of AB genomes and distorted segregation of polymorphic SNP markers through SNP genotyping. Results suggested that (1) the recombination of AB genomes in the SPW populations was about 3- to 4-fold higher than that in the SHW populations, resulting from the significantly (P < 0.01) increased RF between adjacent and linked SNP loci, especially the variations that occurred in a pericentromeric region which would further enrich genetic diversity; (2) the crosses of hexaploid × tetraploid wheat could be an efficient way to produce pentaploid derivatives than the crosses of tetraploid × hexaploid wheat according to the higher germination rate found in the former crosses; (3) the high proportion of distorted segregation loci that skewed in favor of the female parent genotype/allele in the SPW populations might associate with the fitness and survival of the offspring. Based on the presented data, we propose that pentaploid hybrids should increasingly be used in wheat breeding. In addition, the contribution of gene flow from tetraploid wheat to bread wheat mediated by pentaploid introgressive hybridization also was discussed in the re-evolution of bread wheat.

Three-nucleotide periodicity of nucleotide diversity in a population enables the identification of open reading frames

Accurate prediction of open reading frames (ORFs) is important for studying and using genome sequences. Ribosomes move along mRNA strands with a step of three nucleotides and datasets carrying this information can be used to predict ORFs. The ribosome-protected footprints (RPFs) feature a significant 3-nt periodicity on mRNAs and are powerful in predicting translating ORFs, including small ORFs (sORFs), but the application of RPFs is limited because they are too short to be accurately mapped in complex genomes. In this study, we found a significant 3-nt periodicity in the datasets of populational genomic variants in coding sequences, in which the nucleotide diversity increases every three nucleotides. We suggest that this feature can be used to predict ORFs and develop the Python package ‘OrfPP’, which recovers ~83% of the annotated ORFs in the tested genomes on average, independent of the population sizes and the complexity of the genomes. The novel ORFs, including sORFs, identified from single-nucleotide polymorphisms are supported by protein mass spectrometry evidence comparable to that of the annotated ORFs. The application of OrfPP to tetraploid cotton and hexaploid wheat genomes successfully identified 76.17% and 87.43% of the annotated ORFs in the genomes, respectively, as well as 4704 sORFs, including 1182 upstream and 2110 downstream ORFs in cotton and 5025 sORFs, including 232 upstream and 234 downstream ORFs in wheat. Overall, we propose an alternative and supplementary approach for ORF prediction that can extend the studies of sORFs to more complex genomes.

Whole-genome resequencing of the wheat A subgenome progenitor Triticum urartu provides insights into its demographic history and geographic adaptation

Triticum urartu is the progenitor of the A subgenome in tetraploid and hexaploid wheat. Uncovering the landscape of genetic variations in T. urartu will help us understand the evolutionary and polyploid characteristics of wheat. Here, we investigated the population genomics of T. urartu by genome-wide sequencing of 59 representative accessions collected around the world. A total of 42.2 million high-quality single-nucleotide polymorphisms and 3 million insertions and deletions were obtained by mapping reads to the reference genome. The ancient T. urartu population experienced a significant reduction in effective population size (Ne) from ∼3 000 000 to ∼140 000 and subsequently split into eastern Mediterranean coastal and Mesopotamian-Transcaucasian populations during the Younger Dryas period. A map of allelic drift paths displayed splits and mixtures between different geographic groups, and a strong genetic drift towards hexaploid wheat was also observed, indicating that the direct donor of the A subgenome originated from northwestern Syria. Genetic changes were revealed between the eastern Mediterranean coastal and Mesopotamian-Transcaucasian populations in genes orthologous to those regulating plant development and stress responses. A genome-wide association study identified two single-nucleotide polymorphisms in the exonic regions of the SEMI-DWARF 37 ortholog that corresponded to the different T. urartu ecotype groups. Our study provides novel insights into the origin and genetic legacy of the A subgenome in polyploid wheat and contributes a gene repertoire for genomics-enabled improvements in wheat breeding.

Polyploidization affects the allelic variation of jasmonate-regulated protein Ta-JA1 belonging to the monocot chimeric jacalin (MCJ) family in wild emmer wheat

The jasmonate-regulated protein Ta-JA1 belongs to the monocot chimeric jacalin (MCJ) family and plays a vital role in stress resistance in wheat. However, the impact of wheat polyploidization on Ta-JA1 remains unclear. In this study, 149 members of the MCJ family were identified among members of Triticeae using a genome-wide approach. The genes were resolved into three clades; MCJ genes in each clade were derived from different donor genes during evolution. Segmental duplication may have been the primary driver, compared with tandem duplication, of expansion in the MCJ family of wheat. Gene loss and acquisition occurred during tetraploidization, and the core expansion of the family occurred after tetraploidization. Sequencing data for 2104 accessions of T. aestivum and 99 accessions of T. dicoccoides showed that Ta-JA1-2A and Ta-JA1 were highly conserved in common wheat, and four alleles (TdJA1-Ax2, TdJA1-Ay2, TdJA1-Ax3, and TdJA1-Ay3) were detected in T. dicoccoides. Using gene-specific markers, one AsJA1-B allele was detected in 11 Ae. speltoides accessions and one TuJA1-Ax1 allele was detected in 70 T. urartu accessions. Six alleles were detected on chromosome 2A: TdJA1-Ax1 (13 accessions), TdJA1-Ay1 (57 accessions), TdJA1-Ax2 (23 accessions), TdJA1-Ay2 (42 accessions), TdJA1-Ax3 (29 accessions), and TdJA1-Ay3 (251 accessions). Only one allele (TdJA1-B) on chromosome 2B was detected in 415 T. dicoccoides accessions. A geographical distribution analysis revealed that Israel hosted higher allelic variation than other regions. Quantitative reverse transcription PCR analysis indicated that divergence in expression has occurred among Ta-JA1 alleles and, notably, TdJA1-Ax1 and TdJA1-Ay1 showed significantly higher expression levels than the other four allelic types in T. dicoccoides. The present results contribute to an improved understanding of the effects of polyploidization on the MCJ gene family and the functions of Ta-JA1, and may be useful to enrich common wheat germplasm resources.

The integrated genomics of crop domestication and breeding

As a major event in human civilization, wild plants were successfully domesticated to be crops, largely owing to continuing artificial selection. Here, we summarize new discoveries made during the past decade in crop domestication and breeding. The construction of crop genome maps and the functional characterization of numerous trait genes provide foundational information. Approaches to read, interpret, and write complex genetic information are being leveraged in many plants for highly efficient de novo or re-domestication. Understanding the underlying mechanisms of crop microevolution and applying the knowledge to agricultural productions will give possible solutions for future challenges in food security.

Enhancer transcription detected in the nascent transcriptomic landscape of bread wheat

The precise spatiotemporal gene expression is orchestrated by enhancers that lack general sequence features and thus are difficult to be computationally identified. By nascent RNA sequencing combined with epigenome profiling, we detect active transcription of enhancers from the complex bread wheat genome. We find that genes associated with transcriptional enhancers are expressed at significantly higher levels, and enhancer RNA is more precise and robust in predicting enhancer activity compared to chromatin features. We demonstrate that sub-genome-biased enhancer transcription could drive sub-genome-biased gene expression. This study highlights enhancer transcription as a hallmark in regulating gene expression in wheat.

New insights into the evolution of CAF1 family and utilization of TaCAF1Ia1 specificity to reveal the origin of the maternal progenitor for common wheat

INTRODUCTION

Until now, the most likely direct maternal progenitor (AABB) for common wheat (AABBDD) has yet to be identified. Here, we try to solve this particular problem with the specificity of a novel gene family in wheat and by using large population of rare germplasm resources.

OBJECTIVES

Dissect the novelty of TaCAF1Ia subfamily in wheat. Exploit the conservative and specific characteristics of TaCAF1Ia1 to reveal the origin of the maternal progenitor for common wheat.

METHODS

Phylogenetic and collinear analysis of TaCAF1 genes were performed to identify the evolutionary specificity of TaCAF1Ia subfamily. The large-scale expression patterns and interaction patterns analysis of CCR4-NOT complex were used to clarify the expressed and structural specificity of TaCAF1Ia subfamily in wheat. The population resequencing and phylogeny analysis of the TaCAF1Ia1 were utilized for the traceability analysis to understand gene-pool exchanges during the transferring and subsequent development from tetraploid to hexaploidy wheat.

RESULTS

TaCAF1Ia is a novel non-typical CAF1 subfamily without DEDD (Asp-Glu-Asp-Asp) domain, whose members were extensively duplicated in wheat genome. The replication events had started and constantly evolved from ancestor species. Specifically, it was found that a key member CAF1Ia1 was highly specialized and only existed in the subB genome and S genome. Unlike CAF1s reported in other plants, TaCAF1Ia genes may be new factors for anther development. These atypical TaCAF1s could also form CCR4-NOT complex in wheat but with new interaction sites. Utilizing the particular but conserved characteristics of the TaCAF1Ia1 gene, the comparative analysis of haplotypes composition for TaCAF1Ia1 were identified among wheat populations with different ploidy levels. Based on this, the dual-lineages origin model of maternal progenitor for common wheat and potential three-lineages domestication model for cultivated tetraploid wheat were proposed.

CONCLUSION

This study brings fresh insights for revealing the origin of wheat and the function of CAF1 in plants.

Genome-wide identification, expression analysis and evolutionary relationships of the IQ67-domain gene family in common wheat (Triticum aestivum L.) and its progenitors

BACKGROUND

The plant-specific IQ67-domain (IQD) gene family plays an important role in plant development and stress responses. However, little is known about the IQD family in common wheat (Triticum aestivum L), an agriculturally important crop that provides more than 20% of the calories and protein consumed in the modern human diet.

RESULTS

We identified 125 IQDs in the wheat genome and divided them into four subgroups by phylogenetic analysis. The IQDs belonging to the same subgroup had similar exon-intron structure and conserved motif composition. Polyploidization contributed significantly to the expansion of IQD genes in wheat. Characterization of the expression profile of these genes revealed that a few T. aestivum (Ta)IQDs showed high tissue-specificity. The stress-induced expression pattern also revealed a potential role of TaIQDs in environmental adaptation, as TaIQD-2A-2, TaIQD-3A-9 and TaIQD-1A-7 were significantly induced by cold, drought and heat stresses, and could be candidates for future functional characterization. In addition, IQD genes in the A, B and D subgenomes displayed an asymmetric evolutionary pattern, as evidenced by their different gain or loss of member genes, expression levels and nucleotide diversity.

CONCLUSIONS

This study elucidated the potential biological functions and evolutionary relationships of the IQD gene family in wheat and revealed the divergent fates of IQD genes during polyploidization.

Fine mapping and distribution analysis of hybrid necrosis genes Ne1 and Ne2 in wheat in China

Flanking markers useful for identifying hybrid necrosis alleles were identified by fine mapping Ne1 and Ne2 and the distribution of the two necrosis genes was investigated in Chinese elite wheat varieties. Hybrid necrosis of wheat is caused by the interaction of two dominant complementary genes Ne1 and Ne2 present separately in normal parents and is regarded as a barrier to gene transfer in wheat breeding. However, the necrosis alleles still occur at a high frequency in modern wheat varieties. In this study, we constructed two high-density genetic maps of Ne1 and Ne2 in winter wheat. In these cultivars, Ne1 was found to be located in a span interval of 0.50 centimorgan (cM) on chromosome 5BL delimited by markers Nwu_5B_4137 and Nwu_5B_5114, while Ne2 co-segregated with markers Lseq102 and TC67744 on 2BS. Statistical analysis confirmed that the dosage effect of Ne1 and Ne2 also existed in moderate and severe hybrid necrosis systems, and the symptoms of necrosis can also be affected by the genetic background. Furthermore, we clarified the discrete distribution and proportion of the Ne1 and Ne2 in the 10 China's agro-ecological production zones. We concluded that 26.2% and 33.2% of the 1364 cultivars (lines) were genotyped with Ne1Ne1ne2ne2 and ne1ne1Ne2Ne2, respectively and introduced modern cultivars should directly affect the frequencies of necrosis genes in modern Chinese cultivars (lines), especially that of Ne2. Taking investigations in spring wheat together, we proposed that hybrid necrosis alleles could positively affect breeding owing to their linked excellent genes such as Lr13. Additionally, based on the pedigrees and hybridization tests, we speculated that the Ne1 and Ne2 in winter wheat may directly originate from wild emmer and introduced cultivars or hexaploid triticale, respectively.

Genome sequences of three Aegilops species of the section Sitopsis reveal phylogenetic relationships and provide resources for wheat improvement

Aegilops is a close relative of wheat (Triticum spp.), and Aegilops species in the section Sitopsis represent a rich reservoir of genetic diversity for the improvement of wheat. To understand their diversity and advance their utilization, we produced whole-genome assemblies of Aegilops longissima and Aegilops speltoides. Whole-genome comparative analysis, along with the recently sequenced Aegilops sharonensis genome, showed that the Ae. longissima and Ae. sharonensis genomes are highly similar and are most closely related to the wheat D subgenome. By contrast, the Ae. speltoides genome is more closely related to the B subgenome. Haplotype block analysis supported the idea that Ae. speltoides genome is closest to the wheat B subgenome, and highlighted variable and similar genomic regions between the three Aegilops species and wheat. Genome-wide analysis of nucleotide-binding leucine-rich repeat (NLR) genes revealed species-specific and lineage-specific NLR genes and variants, demonstrating the potential of Aegilops genomes for wheat improvement.

Genomic sequencing of Thinopyrum elongatum chromosome arm 7EL, carrying fusarium head blight resistance, and characterization of its impact on the transcriptome of the introgressed line CS-7EL

Background

The tall wheatgrass species Thinopyrum elongatum carries a strong fusarium head blight (FHB) resistance locus located on the long arm of chromosome 7 (7EL) as well as resistance to leaf and stem rusts, all diseases with a significant impact on wheat production. Towards understanding the contribution of Th. elongatum 7EL to improvement of disease resistance in wheat, the genomic sequence of the 7EL fragment present in the wheat Chinese Spring (CS) telosomic addition line CS-7EL was determined and the contribution and impact of 7EL on the rachis transcriptome during FHB infection was compared between CS and CS-7EL.

Results

We assembled the Th. elongatum 7EL chromosome arm using a reference-guided approach. Combining this assembly with the available reference sequence for CS hexaploid wheat provided a reliable reference for interrogating the transcriptomic differences in response to infection conferred by the 7EL fragment. Comparison of the transcriptomes of rachis tissues from CS and CS-7EL showed expression of Th. elongatum transcripts as well as modulation of wheat transcript expression profiles in the CS-7EL line. Expression profiles at 4 days after infection with Fusarium graminearum, the causal agent of FHB, showed an increased in expression of genes associated with an effective defense response, in particular glucan endo-1,3-beta-glucosidases and chitinases, in the FHB-resistant line CS-7EL while there was a larger increase in differential expression for genes associated with the level of fungal infection in the FHB-susceptible line CS. One hundred and seven 7EL transcripts were expressed in the smallest 7EL region defined to carry FHB resistance.

Conclusion

7EL contributed to CS-7EL transcriptome by direct expression and through alteration of wheat transcript profiles. FHB resistance in CS-7EL was associated with transcriptome changes suggesting a more effective defense response. A list of candidate genes for the FHB resistance locus on 7EL has been established.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08433-8.

Population genomic analysis reveals domestication of cultivated rye from weedy rye

Rye (Secale cereale) is an important crop with multiple uses and a valuable genetic resource for wheat breeding. However, due to its complex genome and outcrossing nature, the origin of cultivated rye remains elusive. The geneticist N.I. Vavilov proposed that cultivated rye had been domesticated from weedy rye, rather than directly from wild species like other crops. Unraveling the domestication history of rye will extend our understanding of crop evolution and upend our inherent understanding of agricultural weeds. To this end, in this study we generated the 8.5 Tb of whole-genome resequencing data from 116 worldwide accessions of wild, weedy, and cultivated rye, and demonstrated that cultivated rye was domesticated directly from weedy relatives with a similar but enhanced genomic selection by humans. We found that a repertoire of genes that experienced artificial selection is associated with important agronomic traits, including shattering, grain yield, and disease resistance. Furthermore, we identified a composite introgression in cultivated rye from the wild perennial Secale strictum and detected a 2-Mb introgressed fragment containing a candidate ammonium transporter gene with potential effect on the grain yield and plant growth of rye. Taken together, our findings unravel the domestication history of cultivated rye, suggest that interspecific introgression serves as one of the likely causes of obscure species taxonomy of the genus Secale, and provide an important resource for future rye and wheat breeding.

Genomic variants affecting homoeologous gene expression dosage contribute to agronomic trait variation in allopolyploid wheat

Allopolyploidy greatly expands the range of possible regulatory interactions among functionally redundant homoeologous genes. However, connection between the emerging regulatory complexity and expression and phenotypic diversity in polyploid crops remains elusive. Here, we use diverse wheat accessions to map expression quantitative trait loci (eQTL) and evaluate their effects on the population-scale variation in homoeolog expression dosage. The relative contribution of cis- and trans-eQTL to homoeolog expression variation is strongly affected by both selection and demographic events. Though trans-acting effects play major role in expression regulation, the expression dosage of homoeologs is largely influenced by cis-acting variants, which appear to be subjected to selection. The frequency and expression of homoeologous gene alleles showing strong expression dosage bias are predictive of variation in yield-related traits, and have likely been impacted by breeding for increased productivity. Our study highlights the importance of genomic variants affecting homoeolog expression dosage in shaping agronomic phenotypes and points at their potential utility for improving yield in polyploid crops.

Genome-wide identification of the NLR gene family in Haynaldia villosa by SMRT-RenSeq

Background

Nucleotide-binding and leucine-rich repeat (NLR) genes have attracted wide attention due to their crucial role in protecting plants from pathogens. SMRT-RenSeq, combining PacBio sequencing after resistance gene enrichment sequencing (RenSeq), is a powerful method for selectively capturing and sequencing full-length NLRs. Haynaldia villosa, a wild grass species with a proven potential for wheat improvement, confers resistance to multiple diseases. So, genome-wide identification of the NLR gene family in Haynaldia villosa by SMRT-RenSeq can facilitate disease resistance genes exploration.

Results

In this study, SMRT-RenSeq was performed to identify the genome-wide NLR complement of H. villosa. In total, 1320 NLRs were annotated in 1169 contigs, including 772 complete NLRs. All the complete NLRs were phylogenetically analyzed and 11 main clades with special characteristics were derived. NLRs could be captured with high efficiency when aligned with cloned R genes, and cluster expansion in some specific gene loci was observed. The physical location of NLRs to individual chromosomes in H. villosa showed a perfect homoeologous relationship with group 1, 2, 3, 5 and 6 of other Triticeae species, however, NLRs physically located on 4VL were largely in silico predicted to be located on the homoeologous group 7. Fifteen types of integrated domains (IDs) were integrated in 52 NLRs, and Kelch and B3 NLR-IDs were found to have expanded in H. villosa, while DUF948, NAM-associated and PRT_C were detected as unique integrated domains implying the new emergence of NLR-IDs after H. villosa diverged from other species.

Conclusion

SMRT-RenSeq is a powerful tool to identify NLR genes from wild species using the baits of the evolutionary related species with reference sequences. The availability of the NLRs from H. villosa provide a valuable library for R gene mining and transfer of disease resistance into wheat.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08334-w.

Genome sequencing-based coverage analyses facilitate high-resolution detection of deletions linked to phenotypes of gamma-irradiated wheat mutants

BACKGROUND

Gamma-irradiated mutants of Triticum aestivum L., hexaploid wheat, provide novel and agriculturally important traits and are used as breeding materials. However, the identification of causative genomic regions of mutant phenotypes is challenging because of the large and complicated genome of hexaploid wheat. Recently, the combined use of high-quality reference genome sequences of common wheat and cost-effective resequencing technologies has made it possible to evaluate genome-wide polymorphisms, even in complex genomes.

RESULTS

To investigate whether the genome sequencing approach can effectively detect structural variations, such as deletions, frequently caused by gamma irradiation, we selected a grain-hardness mutant from the gamma-irradiated population of Japanese elite wheat cultivar "Kitahonami." The Hardness (Ha) locus, including the puroindoline protein-encoding genes Pina-D1 and Pinb-D1 on the short arm of chromosome 5D, primarily regulates the grain hardness variation in common wheat. We performed short-read genome sequencing of wild-type and grain-hardness mutant plants, and subsequently aligned their short reads to the reference genome of the wheat cultivar "Chinese Spring." Genome-wide comparisons of depth-of-coverage between wild-type and mutant strains detected ~ 130 Mbp deletion on the short arm of chromosome 5D in the mutant genome. Molecular markers for this deletion were applied to the progeny populations generated by a cross between the wild-type and the mutant. A large deletion in the region including the Ha locus was associated with the mutant phenotype, indicating that the genome sequencing is a powerful and efficient approach for detecting a deletion marker of a gamma-irradiated mutant phenotype. In addition, we investigated a pre-harvest sprouting tolerance mutant and identified a 67.8 Mbp deletion on chromosome 3B where Viviparous-B1 and GRAS family transcription factors are located. Co-dominant markers designed to detect the deletion-polymorphism confirmed the association with low germination rate, leading to pre-harvest sprouting tolerance.

CONCLUSIONS

Short read-based genome sequencing of gamma-irradiated mutants facilitates the identification of large deletions linked to mutant phenotypes when combined with segregation analyses in progeny populations. This method allows effective application of mutants with agriculturally important traits in breeding using marker-assisted selection.

Genomic insights into the evolution of Echinochloa species as weed and orphan crop

As one of the great survivors of the plant kingdom, barnyard grasses (Echinochloa spp.) are the most noxious and common weeds in paddy ecosystems. Meanwhile, at least two Echinochloa species have been domesticated and cultivated as millets. In order to better understand the genomic forces driving the evolution of Echinochloa species toward weed and crop characteristics, we assemble genomes of three Echinochloa species (allohexaploid E. crus-galli and E. colona, and allotetraploid E. oryzicola) and re-sequence 737 accessions of barnyard grasses and millets from 16 rice-producing countries. Phylogenomic and comparative genomic analyses reveal the complex and reticulate evolution in the speciation of Echinochloa polyploids and provide evidence of constrained disease-related gene copy numbers in Echinochloa. A population-level investigation uncovers deep population differentiation for local adaptation, multiple target-site herbicide resistance mutations of barnyard grasses, and limited domestication of barnyard millets. Our results provide genomic insights into the dual roles of Echinochloa species as weeds and crops as well as essential resources for studying plant polyploidization, adaptation, precision weed control and millet improvements.

Characterization of TaSPP-5A gene associated with sucrose content in wheat (Triticum aestivum L.)

BACKGROUND

Sucrose, the major product of photosynthesis and the primary sugar transported as a soluble carbohydrate via the phloem, is a critical determinant for harvest yield in wheat crops. Sucrose-phosphatase (SPP) catalyzes the final step in the sucrose biosynthesis pathway, implying its essential role in the plant.

RESULT

In this study, wheat SPP homologs genes were isolated from chromosomes 5A, 5B, and 5D, designated as TaSPP-5A, TaSPP-5B, and TaSPP-5D, respectively. Sequence alignment showed one 1-bp Insertion-deletion (InDel) and three single nucleotide polymorphisms (SNPs) at TaSPP-5A coding region, forming two haplotypes, TaSPP-5Aa and TaSPP-5Ab, respectively. A derived cleaved amplified polymorphism sequence (dCAPS) marker, TaSPP-5A-dCAPS, was developed to discriminate allelic variation based on the polymorphism at position 1242 (C-T). A total of 158 varieties were used to perform a TaSPP-5A marker-trait association analysis, where two haplotypes were significantly associated with sucrose content in two environments and with thousand-grain weight (TGW) and grain length (GL) in three environments. Quantitative real-time PCR further revealed that TaSPP-5Aa showed relatively higher expression than TaSPP-5Ab in wheat seedling leaves, generally associating with increased sucrose content and TGW. The expression of TaSPP-5A and sucrose content in TaSPP-5Aa haplotypes were also higher than those in TaSPP-5Ab haplotypes under both 20% PEG-6000 and 100 μM ABA treatment. Sequence alignment showed that the two TaSPP-5A haplotypes comprised 11 SNPs from -395 to -1962 bp at TaSPP-5A promoter locus, participating in the formation of several conserved sequences, may account for the high expression of TaSPP-5A in TaSPP-5Aa haplotypes. In addition, the distribution analysis of TaSPP-5A haplotypes revealed that TaSPP-5Aa was preferred in the natural wheat population, being strongly positively selected in breeding programs.

CONCLUSION

According to the SNPs detected in the TaSPP-5A sequence, two haplotypes, TaSPP-5Aa and TaSPP-5Ab, were identified among wheat accessions, which potential value for sucrose content selection was validated by association analysis. Our results indicate that the favorable allelic variation TaSPP-5Aa should be valuable in enhancing grain yield by improving the sucrose content. Furthermore, a functional marker, TaSPP-5A-dCAPS, can be used for marker-assisted selection to improve grain weight in wheat and provides insights into the biological function of TaSPP-5A gene.

The Landscape of Autophagy-Related (ATG) Genes and Functional Characterization of TaVAMP727 to Autophagy in Wheat

Autophagy is an indispensable biological process and plays crucial roles in plant growth and plant responses to both biotic and abiotic stresses. This study systematically identified autophagy-related proteins (ATGs) in wheat and its diploid and tetraploid progenitors and investigated their genomic organization, structure characteristics, expression patterns, genetic variation, and regulation network. We identified a total of 77, 51, 29, and 30 ATGs in wheat, wild emmer, T. urartu and A. tauschii, respectively, and grouped them into 19 subfamilies. We found that these autophagy-related genes (ATGs) suffered various degrees of selection during the wheat’s domestication and breeding processes. The genetic variations in the promoter region of Ta2A_ATG8a were associated with differences in seed size, which might be artificially selected for during the domestication process of tetraploid wheat. Overexpression of TaVAMP727 improved the cold, drought, and salt stresses resistance of the transgenic Arabidopsis and wheat. It also promoted wheat heading by regulating the expression of most ATGs. Our findings demonstrate how ATGs regulate wheat plant development and improve abiotic stress resistance. The results presented here provide the basis for wheat breeding programs for selecting varieties of higher yield which are capable of growing in colder, drier, and saltier areas.

Modified expression of TaCYP78A5 enhances grain weight with yield potential by accumulating auxin in wheat (Triticum aestivum L.)

Increasing grain yield has always been the primary goal of crop breeding. KLUH/CYP78A5 has been shown to affect seed size in several plant species, but the relevant molecular mechanism is still unclear and there are no reports of this gene contributing to yield. Here, we demonstrate that modified expression of TaCYP78A5 can enhance wheat grain weight and grain yield per plant by accumulating auxin. TaCYP78A5 is highly expressed in maternal tissues, including ovary and seed coat during wheat development. The constitutive overexpression of TaCYP78A5 leads to significantly increased seed size and weight but not grain yield per plant due to the strengthening of apical dominance. However, localized overexpression of TaCYP78A5 in maternal integument enhances grain weight and grain yield per plant by 4.3%-18.8% and 9.6%-14.7%, respectively, in field trials. Transcriptome and hormone metabolome analyses reveal that TaCYP78A5 participates in auxin synthesis pathway and promotes auxin accumulation and cell wall remodelling in ovary. Phenotype investigation and cytological observation show that localized overexpression of TaCYP78A5 in ovary results in delayed flowering and prolonged proliferation of maternal integument cells, which promote grain enlargement. Moreover, naturally occurring variations in the promoter of TaCYP78A5-2A contribute to thousand-grain weight (TGW) and grain yield per plant of wheat;TaCYP78A5-2A haplotype Ap-HapII with higher activity is favourable for improving grain weight and grain yield per plant and has been positively selected in wheat breeding. Then, a functional marker of TaCYP78A5 haplotype Ap-HapII is developed for marker-assisted selection in wheat grain and yield improvement.

Population genomic analysis of Aegilops tauschii identifies targets for bread wheat improvement

Aegilops tauschii, the diploid wild progenitor of the D subgenome of bread wheat, is a reservoir of genetic diversity for improving bread wheat performance and environmental resilience. Here we sequenced 242 Ae. tauschii accessions and compared them to the wheat D subgenome to characterize genomic diversity. We found that a rare lineage of Ae. tauschii geographically restricted to present-day Georgia contributed to the wheat D subgenome in the independent hybridizations that gave rise to modern bread wheat. Through k-mer-based association mapping, we identified discrete genomic regions with candidate genes for disease and pest resistance and demonstrated their functional transfer into wheat by transgenesis and wide crossing, including the generation of a library of hexaploids incorporating diverse Ae. tauschii genomes. Exploiting the genomic diversity of the Ae. tauschii ancestral diploid genome permits rapid trait discovery and functional genetic validation in a hexaploid background amenable to breeding.

Boosting the power of transcriptomics by developing an efficient gene expression profiling approach

Recent advances in plant genomics are scaling up gene expression profiling from the individual level to the population level, making transcriptomics a more powerful tool while deciphering the genome function. This study developed an efficient 3'RNA-seq method, Simplified Poly(A) Anchored Sequencing (SiPAS), to perform large-scale experiments of gene expression quantification. Aside from being cost-effective, by conducting a comprehensive performance assessment of SiPAS in hexaploid wheat, we demonstrated that SiPAS is highly sensitive, accurate, and reproducible while quantifying gene expression. Our method is anticipated to boost studies of population transcriptomics in plants and improve our understanding of genome biology.

Genome-wide association mapping of Fusarium crown rot resistance in Aegilops tauschii

Fusarium crown rot (FCR), caused by various Fusarium species, is a primary fungal disease in most wheat-growing regions worldwide. A. tauschii, the diploid wild progenitor of the D-genome of common wheat, is a reservoir of genetic diversity for improving bread wheat biotic and abiotic resistance/tolerance. A worldwide collection of 286 A. tauschii accessions was used to evaluate FCR resistance. Population structure analysis revealed that 115 belonged to the A. tauschii ssp. strangulata subspecies, and 171 belonged to the A. tauschii ssp. tauschii subspecies. Five accessions with disease index values lower than 20 showed moderate resistance to FCR. These five originated from Afghanistan, China, Iran, Uzbekistan, and Turkey, all belonging to the tauschii subspecies. Genome-wide association mapping using 6,739 single nucleotide polymorphisms (SNPs) revealed that two SNPs on chromosome 2D and four SNPs on chromosome 7D were significantly associated with FCR resistance. Almost all FCR resistance alleles were presented in accessions from the tauschii subspecies, and only 4, 11, and 19 resistance alleles were presented in accessions from the strangulata subspecies. Combining phenotypic correlation analysis and genome-wide association mapping confirmed that FCR resistance loci were independent of flowering time, heading date, and plant height in this association panel. Six genes encoding disease resistance-related proteins were selected as candidates for further validation. The identified resistant A. tauschii accessions will provide robust resistance gene sources for breeding FCR-resistant cultivars. The associated loci/genes will accelerate and improve FCR in breeding programs by deploying marker-assisted selection.

Whole-genome resequencing of worldwide wild and domestic sheep elucidates genetic diversity, introgression and agronomically important loci

Domestic sheep and their wild relatives harbor substantial genetic variants that can form the backbone of molecular breeding, but their genome landscapes remain understudied. Here, we present a comprehensive genome resource for wild ovine species, landraces and improved breeds of domestic sheep, comprising high-coverage (∼16.10×) whole genomes of 810 samples from 7 wild species and 158 diverse domestic populations. We detected, in total, ∼121.2 million single nucleotide polymorphisms, ∼61 million of which are novel. Some display significant (P < 0.001) differences in frequency between wild and domestic species, or are private to continent-wide or individual sheep populations. Retained or introgressed wild gene variants in domestic populations have contributed to local adaptation, such as the variation in the HBB associated with plateau adaptation. We identified novel and previously reported targets of selection on morphological and agronomic traits such as stature, horn, tail configuration, and wool fineness. We explored the genetic basis of wool fineness and unveiled a novel mutation (chr25: T7,068,586C) in the 3′-UTR of IRF2BP2 as plausible causal variant for fleece fiber diameter. We reconstructed prehistorical migrations from the Near Eastern domestication center to South-and-Southeast Asia and found two main waves of migrations across the Eurasian Steppe and the Iranian Plateau in the Early and Late Bronze Ages. Our findings refine our understanding of genome variation as shaped by continental migrations, introgression, adaptation, and selection of sheep.

Wheat speciation and adaptation: perspectives from reticulate evolution

Reticulate evolution through the interchanging of genetic components across organisms can impact significantly on the fitness and adaptation of species. Bread wheat (Triticum aestivum subsp. aestivum) is one of the most important crops in the world. Allopolyploid speciation, frequent hybridization, extensive introgression, and occasional horizontal gene transfer (HGT) have been shaping a typical paradigm of reticulate evolution in bread wheat and its wild relatives, which is likely to have a substantial influence on phenotypic traits and environmental adaptability of bread wheat. In this review, we outlined the evolutionary history of bread wheat and its wild relatives with a highlight on the interspecific hybridization events, demonstrating the reticulate relationship between species/subspecies in the genera Triticum and Aegilops. Furthermore, we discussed the genetic mechanisms and evolutionary significance underlying the introgression of bread wheat and its wild relatives. An in-depth understanding of the evolutionary process of Triticum species should be beneficial to future genetic study and breeding of bread wheat.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

New insights into the dispersion history and adaptive evolution of taxon Aegilops tauschii in China

Aegilops tauschii, the wild progenitor of wheat D-genome and a valuable germplasm for wheat improvement, has a wide natural distribution from eastern Turkey to China. However, the phylogenetic relationship and dispersion history of Ae. tauschii in China has not been scientifically clarified. In this study, we genotyped 208 accessions (with 104 in China) using ddRAD sequencing and 55K SNP array, and classified the population into six sublineages. Three possible spreading routes or events were identified, resulting in specific distribution patterns, with four sublineages found in Xinjiang, one in Qinghai, two in Shaanxi and one in Henan. We also established the correlation of SNP-based, karyotype-based and spike-morphology-based techniques to demonstrate the internal classification of Ae. tauschii, and developed consensus dataset with 1245 putative accessions by merging data previously published. Our analysis suggested that eight inter-lineage accessions could be assigned to the putative Lineage 3 and these accessions would help to conserve the genetic diversity of the species. By developing the consensus phylogenetic relationships of Ae. tauschii, our work validated the hypothesis on the dispersal history of Ae. tauschii in China, and contributed to the efficient and comprehensive germplasm-mining of the species.

High molecular weight glutenin gene diversity in Aegilops tauschii demonstrates unique origin of superior wheat quality

Central to the diversity of wheat products was the origin of hexaploid bread wheat, which added the D-genome of Aegilops tauschii to tetraploid wheat giving rise to superior dough properties in leavened breads. The polyploidization, however, imposed a genetic bottleneck, with only limited diversity introduced in the wheat D-subgenome. To understand genetic variants for quality, we sequenced 273 accessions spanning the known diversity of Ae. tauschii. We discovered 45 haplotypes in Glu-D1, a major determinant of quality, relative to the two predominant haplotypes in wheat. The wheat allele 2 + 12 was found in Ae. tauschii Lineage 2, the donor of the wheat D-subgenome. Conversely, the superior quality wheat allele 5 + 10 allele originated in Lineage 3, a recently characterized lineage of Ae. tauschii, showing a unique origin of this important allele. These two wheat alleles were also quite similar relative to the total observed molecular diversity in Ae. tauschii at Glu-D1. Ae. tauschii is thus a reservoir for unique Glu-D1 alleles and provides the genomic resource to begin utilizing new alleles for end-use quality improvement in wheat breeding programs.

Genetic identification and characterization of chromosomal regions for kernel length and width increase from tetraploid wheat

BACKGROUND

Improvement of wheat gercTriticum aestivum L.) yield could relieve global food shortages. Kernel size, as an important component of 1000-kernel weight (TKW), is always a significant consideration to improve yield for wheat breeders. Wheat related species possesses numerous elite genes that can be introduced into wheat breeding. It is thus vital to explore, identify, and introduce new genetic resources for kernel size from wheat wild relatives to increase wheat yield.

RESULTS

In the present study, quantitative trait loci (QTL) for kernel length (KL) and width (KW) were detected in a recombinant inbred line (RIL) population derived from a cross between a wild emmer accession 'LM001' and a Sichuan endemic tetraploid wheat 'Ailanmai' using the Wheat 55 K single nucleotide polymorphism (SNP) array-based constructed linkage map and phenotype from six different environments. We identified eleven QTL for KL and KW including two major ones QKL.sicau-AM-3B and QKW.sicau-AM-4B, the positive alleles of which were from LM001 and Ailanmai, respectively. They explained 17.57 to 44.28% and 13.91 to 39.01% of the phenotypic variance, respectively. For these two major QTL, Kompetitive allele-specific PCR (KASP) markers were developed and used to successfully validate their effects in three F3 populations and two natural populations containing a panel of 272 Chinese wheat landraces and that of 300 Chinese wheat cultivars, respectively. QKL.sicau-AM-3B was located at 675.6-695.4 Mb on chromosome arm 3BL. QKW.sicau-AM-4B was located at 444.2-474.0 Mb on chromosome arm 4BL. Comparison with previous studies suggested that these two major QTL were likely new loci. Further analysis indicated that the positive alleles of QKL.sicau-AM-3B and QKW.sicau-AM-4B had a great additive effect increasing TKW by 6.01%. Correlation analysis between KL and other agronomic traits showed that KL was significantly correlated to spike length, length of uppermost internode, TKW, and flag leaf length. KW was also significantly correlated with TKW. Four genes, TRIDC3BG062390, TRIDC3BG062400, TRIDC4BG037810, and TRIDC4BG037830, associated with kernel development were predicted in physical intervals harboring these two major QTL on wild emmer and Chinese Spring reference genomes.

CONCLUSIONS

Two stable and major QTL for KL and KW across six environments were detected and verified in three biparental populations and two natural populations. Significant relationships between kernel size and yield-related traits were identified. KASP markers tightly linked the two major QTL could contribute greatly to subsequent fine mapping. These results suggested the application potential of wheat related species in wheat genetic improvement.

Rewilding crops for climate resilience: economic analysis and de novo domestication strategies

To match predicted population growth, annual food production should be doubled by 2050. This is not achievable by current agronomical and breeding practices, due to the impact of climate changes and associated abiotic stresses on agricultural production systems. Here, we analyze the impact of global climate trends on crop productivity and show that the overall loss in crop production from climate-driven abiotic stresses may exceed US$170 billion year-1 and represents a major threat to global food security. We also show that abiotic stress tolerance had been present in wild progenitors of modern crops but was lost during their domestication. We argue for a major shift in our paradigm of crop breeding, focusing on climate resilience, and call for a broader use of wild relatives as a major tool in this process. We argue that, while molecular tools are currently in place to harness the potential of climate-resilient genes present in wild relatives, the complex polygenic nature of tolerance traits remains a major bottleneck in this process. Future research efforts should be focused not only on finding appropriate wild relatives but also on development of efficient cell-based high-throughput phenotyping platforms allowing assessment of the in planta operation of key genes.

TdPm60 identified in wild emmer wheat is an ortholog of Pm60 and constitutes a strong candidate for PmG16 powdery mildew resistance

We identified TdPm60 alleles from wild emmer wheat (WEW), an ortholog of Pm60 from T. urartu, which constitutes a strong candidate for PmG16 mildew resistance. Deployment of PmG16 in Israeli modern bread wheat cultivar Ruta improved the resistance to several local Bgt isolates. Wild emmer wheat (WEW), the tetraploid progenitor of durum and bread wheat, is a valuable genetic resource for resistance to powdery mildew fungal disease caused by Blumeria graminis f. sp. tritici (Bgt). PmG16 gene, derived from WEW, confers high resistance to most tested Bgt isolates. We mapped PmG16 to a 1.4-cM interval between the flanking markers uhw386 and uhw390 on Chromosome 7AL. Based on gene annotation of WEW reference genome Zavitan_V1, 34 predicted genes were identified within the ~ 3.48-Mb target region. Six genes were annotated as associated with disease resistance, of which TRIDC7AG077150.1 was found to be highly similar to Pm60, previously cloned from Triticum urartu, and resides in the same syntenic region. The functional molecular marker (FMM) for Pm60 (M-Pm60-S1) co-segregated with PmG16, suggesting the Pm60 ortholog from WEW (designated here as TdPm60) as a strong candidate for PmG16. Sequence alignment identified only eight SNPs that differentiate between TdPm60 and TuPm60. Furthermore, TdPm60 was found to be present also in the WEW donor lines of the powdery mildew resistance genes MlIW172 and MlIW72, mapped to the same region of Chromosome 7AL as PmG16, suggesting that TdPm60 constitutes a candidate also for these genes. Furthermore, screening of additional 230 WEW accessions with Pm60 specific markers revealed 58 resistant accessions from the Southern Levant that harbored TdPm60, while none of the susceptible accessions showed the presence of this gene. Deployment of PmG16 in Israeli modern bread wheat cultivar Ruta conferred resistance against several local Bgt isolates.

Differentially Amplified Repetitive Sequences Among Aegilops tauschii Subspecies and Genotypes

Genomic repetitive sequences commonly show species-specific sequence type, abundance, and distribution patterns, however, their intraspecific characteristics have been poorly described. We quantified the genomic repetitive sequences and performed single nucleotide polymorphism (SNP) analysis between 29 Ae. tauschii genotypes and subspecies using publicly available raw genomic Illumina sequence reads and used fluorescence in situ hybridization (FISH) to experimentally analyze some repeats. The majority of the identified repetitive sequences had similar contents and proportions between anathera, meyeri, and strangulata subspecies. However, two Ty3/gypsy retrotransposons (CL62 and CL87) showed significantly higher abundances, and CL1, CL119, CL213, CL217 tandem repeats, and CL142 retrotransposon (Ty1/copia type) showed significantly lower abundances in subspecies strangulata compared with the subspecies anathera and meyeri. One tandem repeat and 45S ribosomal DNA (45S rDNA) abundances showed a high variation between genotypes but their abundances were not subspecies specific. Phylogenetic analysis using the repeat abundances of the aforementioned clusters placed the strangulata subsp. in a distinct clade but could not discriminate anathera and meyeri. A near complete differentiation of anathera and strangulata subspecies was observed using SNP analysis; however, var. meyeri showed higher genetic diversity. FISH using major tandem repeats couldn't detect differences between subspecies, although (GAA)₁₀ signal patterns generated two different karyotype groups. Taken together, the different classes of repetitive DNA sequences have differentially accumulated between strangulata and the other two subspecies of Ae. tauschii that is generally in agreement with spike morphology, implying that factors affecting repeatome evolution are variable even among highly closely related lineages.