The Aegilops tauschii genome reveals multiple impacts of transposons

Wheat genomics: genomes, pangenomes, and beyond

There is an urgent need to improve wheat for upcoming challenges, including biotic and abiotic stresses. Sustainable wheat improvement requires the introduction of new genes and alleles in high-yielding wheat cultivars. Using new approaches, tools, and technologies to identify and introduce new genes in wheat cultivars is critical. High-quality genomes, transcriptomes, and pangenomes provide essential resources and tools to examine wheat closely to identify and manipulate new and targeted genes and alleles. Wheat genomics has improved excellently in the past 5 years, generating multiple genomes, pangenomes, and transcriptomes. Leveraging these resources allows us to accelerate our crop improvement pipelines. This review summarizes the progress made in wheat genomics and trait discovery in the past 5 years.

Integrative gene duplication and genome-wide analysis as an approach to facilitate wheat reverse genetics: An example in the TaCIPK family

INTRODUCTION

Reverse genetic studies conducted in the plant with a complex or polyploidy genome enriched with large gene families (like wheat) often meet challenges in identifying the key candidate genes related to important traits and prioritizing the genes for functional experiments.

OBJECTIVE

To overcome the above-mentioned challenges of reverse genetics, this work aims to establish an efficient multi-species strategy for genome-wide gene identification and prioritization of the key candidate genes.

METHODS

We established the integrative gene duplication and genome-wide analysis (iGG analysis) as a strategy for pinpointing key candidate genes deserving functional research. The iGG captures the evolution, and the expansion/contraction of large gene families across phylogeny-related species and integrates spatial-temporal expression information for gene function inference. Transgenic approaches were also employed to functional validation.

RESULTS

As a proof-of-concept for the iGG analysis, we took the wheat calcineurin B-like protein-interacting protein kinases (CIPKs) family as an example. We identified CIPKs from seven monocot species, established the orthologous relationship of CIPKs between rice and wheat, and characterized Triticeae-specific CIPK duplicates (e.g., CIPK4 and CIPK17). Integrated with our analysis of CBLs and CBL-CIPK interaction, we revealed that divergent expressions of TaCBLs and TaCIPKs could play an important role in keeping the stoichiometric balance of CBL-CIPK. Furthermore, we validated the function of TaCIPK17-A2 in the regulation of drought tolerance by using transgenic approaches. Overexpression of TaCIPK17 enhanced antioxidant capacity and improved drought tolerance in wheat.

CONCLUSION

The iGG analysis leverages evolutionary and comparative genomics of crops with large genomes to rapidly highlight the duplicated genes potentially associated with speciation, domestication and/or particular traits that deserve reverse-genetic functional studies. Through the identification of Triticeae-specific TaCIPK17 duplicates and functional validation, we demonstrated the effectiveness of the iGG analysis and provided a new target gene for improving drought tolerance in wheat.

Weed biology and management in the multi-omics era: Progress and perspectives

Weeds pose a significant threat to crop production, resulting in substantial yield reduction. In addition, they possess robust weedy traits that enable them to survive in extreme environments and evade human control. In recent years, the application of multi-omics biotechnologies has helped to reveal the molecular mechanisms underlying these weedy traits. In this review, we systematically describe diverse applications of multi-omics platforms for characterizing key aspects of weed biology, including the origins of weed species, weed classification, and the underlying genetic and molecular bases of important weedy traits such as crop-weed interactions, adaptability to different environments, photoperiodic flowering responses, and herbicide resistance. In addition, we discuss limitations to the application of multi-omics techniques in weed science, particularly compared with their extensive use in model plants and crops. In this regard, we provide a forward-looking perspective on the future application of multi-omics technologies to weed science research. These powerful tools hold great promise for comprehensively and efficiently unraveling the intricate molecular genetic mechanisms that underlie weedy traits. The resulting advances will facilitate the development of sustainable and highly effective weed management strategies, promoting greener practices in agriculture.

Characterization of sucrose nonfermenting-1-related protein kinase 2 (SnRK2) gene family in Haynaldia villosa demonstrated SnRK2.9-V enhances drought and salt stress tolerance of common wheat

BACKGROUND

The sucrose nonfermenting-1-related protein kinase 2 (SnRK2) plays a crucial role in responses to diverse biotic/abiotic stresses. Currently, there are reports on these genes in Haynaldia villosa, a diploid wild relative of wheat.

RESULTS

To understand the evolution of SnRK2-V family genes and their roles in various stress conditions, we performed genome-wide identification of the SnRK2-V gene family in H. villosa. Ten SnRK2-V genes were identified and characterized for their structures, functions and spatial expressions. Analysis of gene exon/intron structure further revealed the presence of evolutionary paths and replication events of SnRK2-V gene family in the H. villosa. In addition, the features of gene structure, the chromosomal location, subcellular localization of the gene family were investigated and the phylogenetic relationship were determined using computational approaches. Analysis of cis-regulatory elements of SnRK2-V gene members revealed their close correlation with different phytohormone signals. The expression profiling revealed that ten SnRK2-V genes expressed at least one tissue (leave, stem, root, or grain), or in response to at least one of the biotic (stripe rust or powdery mildew) or abiotic (drought or salt) stresses. Moreover, SnRK2.9-V was up-regulated in H. villosa under the drought and salt stress and overexpressing of SnRK2.9-V in wheat enhanced drought and salt tolerances via enhancing the genes expression of antioxidant enzymes, revealing a potential value of SnRK2.9-V in wheat improvement for salt tolerance.

CONCLUSION

Our present study provides a basic genome-wide overview of SnRK2-V genes in H. villosa and demonstrates the potential use of SnRK2.9-V in enhancing the drought and salt tolerances in common wheat.

Long-read sequencing reveals the structural complexity of genomic integration of HPV DNA in cervical cancer cell lines

BACKGROUND

Cervical cancer (CC) causes more than 311,000 deaths annually worldwide. The integration of human papillomavirus (HPV) is a crucial genetic event that contributes to cervical carcinogenesis. Despite HPV DNA integration is known to disrupt the genomic architecture of both the host and viral genomes in CC, the complexity of this process remains largely unexplored.

RESULTS

In this study, we conducted whole-genome sequencing (WGS) at 55-65X coverage utilizing the PacBio long-read sequencing platform in SiHa and HeLa cells, followed by comprehensive analyses of the sequence data to elucidate the complexity of HPV integration. Firstly, our results demonstrated that PacBio long-read sequencing effectively identifies HPV integration breakpoints with comparable accuracy to targeted-capture Next-generation sequencing (NGS) methods. Secondly, we constructed detailed models of complex integrated genome structures that included both the HPV genome and nearby regions of the human genome by utilizing PacBio long-read WGS. Thirdly, our sequencing results revealed the occurrence of a wide variety of genome-wide structural variations (SVs) in SiHa and HeLa cells. Additionally, our analysis further revealed a potential correlation between changes in gene expression levels and SVs on chromosome 13 in the genome of SiHa cells.

CONCLUSIONS

Using PacBio long-read sequencing, we have successfully constructed complex models illustrating HPV integrated genome structures in SiHa and HeLa cells. This accomplishment serves as a compelling demonstration of the valuable capabilities of long-read sequencing in detecting and characterizing HPV genomic integration structures within human cells. Furthermore, these findings offer critical insights into the complex process of HPV16 and HPV18 integration and their potential contribution to the development of cervical cancer.

Transcriptome analysis in Aegilops tauschii unravels further insights into genetic control of stripe rust resistance

MAIN CONCLUSION

The Aegilops tauschii resistant accession prevented the pathogen colonization by controlling the sugar flow and triggering the hypersensitive reaction. This study suggested that NBS-LRRs probably induce resistance through bHLH by controlling JA- and SA-dependent pathways. Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst) is one of wheat's most destructive fungal diseases that causes a severe yield reduction worldwide. The most effective and economically-friendly strategy to manage this disease is genetic resistance which can be achieved through deploying new and effective resistance genes. Aegilops tauschii, due to its small genome and co-evolution with Pst, can provide detailed information about underlying resistance mechanisms. Hence, we used RNA-sequencing approach to identify the transcriptome variations of two contrasting resistant and susceptible Ae. tauschii accessions in interaction with Pst and differentially expressed genes (DEGs) for resistance to stripe rust. Gene ontology, pathway analysis, and search for functional domains, transcription regulators, resistance genes, and protein-protein interactions were used to interpret the results. The genes encoding NBS-LRR, CC-NBS-kinase, and phenylalanine ammonia-lyase, basic helix-loop-helix (bHLH)-, basic-leucine zipper (bZIP)-, APETALA2 (AP2)-, auxin response factor (ARF)-, GATA-, and LSD-like transcription factors were up-regulated exclusively in the resistant accession. The key genes involved in response to salicylic acid, amino sugar and nucleotide sugar metabolism, and hypersensitive response contributed to plant defense against stripe rust. The activation of jasmonic acid biosynthesis and starch and sucrose metabolism pathways under Pst infection in the susceptible accession explained the colonization of the host. Overall, this study can fill the gaps in the literature on host-pathogen interaction and enrich the Ae. tauschii transcriptome sequence information. It also suggests candidate genes that could guide future breeding programs attempting to develop rust-resistant cultivars.

A platform for whole-genome speed introgression from Aegilops tauschii to wheat for breeding future crops

Breeding new and sustainable crop cultivars of high yields and desirable traits has been a major challenge for ensuring food security for the growing global human population. For polyploid crops such as wheat, introducing genetic variation from wild relatives of its subgenomes is a key strategy to improve the quality of their breeding pools. Over the past decades, considerable progress has been made in speed breeding, genome sequencing, high-throughput phenotyping and genomics-assisted breeding, which now allows us to realize whole-genome introgression from wild relatives to modern crops. Here, we present a standardized protocol to rapidly introgress the entire genome of Aegilops tauschii, the progenitor of the D subgenome of bread wheat, into elite wheat backgrounds. This protocol integrates multiple modern high-throughput technologies and includes three major phases: development of synthetic octaploid wheat, generation of hexaploid A. tauschii-wheat introgression lines (A-WIs) and homozygosis of the generated A-WIs. Our approach readily generates stable introgression lines in 2 y, thus greatly accelerating the generation of A-WIs and the introduction of desirable genes from A. tauschii to wheat cultivars. These A-WIs are valuable for wheat-breeding programs and functional gene discovery. The current protocol can be easily modified and used for introgressing the genomes of wild relatives to other polyploid crops.

Cytological analysis of the diploid-like inheritance of newly synthesized allotetraploid wheat

Polyploidization is a process which is related to species hybridization and whole genome duplication. It is widespread among angiosperm evolution and is essential for speciation and diversification. Allopolyploidization is mainly derived from interspecific hybridization and is believed to pose chromosome imbalances and genome instability caused by meiotic irregularity. However, the self-compatible allopolyploid in wild nature is cytogenetically and genetically stable. Whether this stabilization form was achieved in initial generation or a consequence of long term of evolution was largely unknown. Here, we synthesized a series of nascent allotetraploid wheat derived from three diploid genomes of A, S*, and D. The chromosome numbers of the majority of the progeny derived from these newly formed allotetraploid wheat plants were found to be relatively consistent, with each genome containing 14 chromosomes. In meiosis, bivalent was the majority of the chromosome configuration in metaphase I which supports the stable chromosome number inheritance in the nascent allotetraploid. These findings suggest that diploidization occurred in the newly formed synthetic allotetraploid wheat. However, we still detected aneuploids in a proportion of newly formed allotetraploid wheat, and meiosis of these materials present more irregular chromosome behavior than the euploid. We found that centromere pairing and centromere clustering in meiosis was affected in the aneuploids, which suggest that aneuploidy may trigger the irregular interactions of centromere in early meiosis which may take participate in promoting meiosis stabilization in newly formed allotetraploid wheat.

Deciphering the evolution and complexity of wheat germplasm from a genomic perspective

Bread wheat provides an essential fraction of the daily calorific intake for humanity. Due to its huge and complex genome, progress in studying on the wheat genome is substantially trailed behind those of the other two major crops, rice and maize, for at least a decade. With rapid advances in genome assembling and reduced cost of high-throughput sequencing, emerging de novo genome assemblies of wheat and whole-genome sequencing data are leading to a paradigm shift in wheat research. Here, we review recent progress in dissecting the complex genome and germplasm evolution of wheat since the release of the first high-quality wheat genome. New insights have been gained in the evolution of wheat germplasm during domestication and modern breeding progress, genomic variations at multiple scales contributing to the diversity of wheat germplasm, and complex transcriptional and epigenetic regulations of functional genes in polyploid wheat. Genomics databases and bioinformatics tools meeting the urgent needs of wheat genomics research are also summarized. The ever-increasing omics data, along with advanced tools and well-structured databases, are expected to accelerate deciphering the germplasm and gene resources in wheat for future breeding advances.

Identification and introgression of a novel leaf rust resistance gene from Thinopyrum intermedium chromosome 7Js into wheat

KEY MESSAGE

A novel leaf rust resistance locus located on a terminal segment (0-69.29 Mb) of Thinopyrum intermedium chromosome arm 7JsS has been introduced into wheat genome for disease resistance breeding. Xiaoyan 78829, a wheat-Thinopyrum intermedium partial amphiploid, exhibits excellent resistance to fungal diseases in wheat. To transfer its disease resistance to common wheat (Triticum aestivum), we previously developed a translocation line WTT26 using chromosome engineering. Disease evaluation showed that WTT26 was nearly immune to 14 common races of leaf rust pathogen (Puccinia triticina) and highly resistant to Ug99 race PTKST of stem rust pathogen (P. graminis f. sp. tritici) at the seedling stage. It also displayed high adult plant resistance to powdery mildew (caused by Blumeria graminis f. sp. tritici). Cytogenetic and molecular marker analysis revealed that WTT26 carried a T4BS·7JsS chromosome translocation. Once transferred into the susceptible wheat genetic background, chromosome 7JsS exhibited its resistance to leaf rust, indicating that the resistance locus was located on this alien chromosome. To enhance the usefulness of this locus in wheat breeding, we further developed several new translocation lines with small Th. intermedium segments using irradiation and developed 124 specific markers using specific-locus amplified fragment sequencing, which increased the marker density of chromosome 7JsS. Furthermore, a refined physical map of chromosome 7JsS was constructed with 74 specific markers, and six bins were thus arranged according to the co-occurrence of markers and alien chromosome segments. Combining data from specific marker amplification and resistance evaluation, we mapped a new leaf rust resistance locus in the 0-69.29 Mb region on chromosome 7JsS. The translocation lines carrying the new leaf rust resistance locus and its linked markers will contribute to wheat disease-resistance breeding.

Unravelling the due importance of pseudogenes and their resurrection in plants

The complexities of a genome are underpinned to the vast expanses of the intergenic region, which constitutes ∼97-98% of the genome. This region is essentially composed of what is colloquially referred to as the "junk DNA" and is composed of various elements like transposons, repeats, pseudogenes, etc. The latter have long been considered as dead elements merely contributing to transcriptional noise in the genome. Many studies now describe the previously unknown regulatory functions of these genes. Recent advances in the Next-generation sequencing (NGS) technologies have allowed unprecedented access to these regions. With the availability of whole genome sequences of more than 788 different plant species in past 20 years, genome annotation has become feasible like never before. Different bioinformatic pipelines are available for the identification of pseudogenes. However, still little is known about their biological functions. The functional validation of these genes remains challenging and research in this area is still in infancy, particularly in plants. CRISPR/Cas-based genome editing could provide solutions to understand the biological roles of these genes by allowing creation of precise edits within these genes. The possibility of pseudogene reactivation or resurrection as has been demonstrated in a few studies might open new avenues of genetic manipulation to yield a desirable phenotype. This review aims at comprehensively summarizing the progress made with regards to the identification of pseudogenes and understanding their biological functions in plants.

Fine mapping of two recessive powdery mildew resistance genes from Aegilops tauschii accession CIae8

KEY MESSAGE

Two recessive powdery mildew resistance loci pmAeCIae8_2DS and pmAeCIae8_7DS from Aegilops tauschii were mapped and two synthesized hexaploid wheat lines were developed by distant hybridization. Wheat powdery mildew (Pm), one of the worldwide destructive fungal diseases, causes significant yield loss up to 30%. The identification of new Pm resistance genes will enrich the genetic diversity of wheat breeding for Pm resistance. Aegilops tauschii is the ancestor donor of sub-genome D of hexaploid wheat. It provides beneficial genes that can be easily transferred into wheat by producing synthetic hexaploid wheat followed by genetic recombination. We assessed the Pm resistance level of 35 Ae. tauschii accessions from different origins. Accession CIae8 exhibited high Pm resistance. Inheritance analysis and gene mapping were performed using F2 and F2:3 populations derived from the cross between CIae8 and a Pm susceptible accession PI574467. The Pm resistance of CIae8 was controlled by two independent recessive genes. Bulked segregate analysis using a 55 K SNP array revealed the SNPs were mainly enriched into genome regions, i.e. 2DS (13.5-20 Mb) and 7DS (4.0-15.5 Mb). The Pm resistance loci were named as pmAeCIae8_2DS and pmAeCIae8_7DS, respectively. By recombinant screening, we narrowed the pmAeCIae8_2DS into a 370-kb interval flanked by markers CINAU-AE7800 (14.89 Mb) and CINAU-AE20 (15.26 Mb), and narrowed the pmAeCIae8_7DS into a 260-kb interval flanked by markers CINAU-AE58 (4.72 Mb) and CINAU-AE25 (4.98 Mb). The molecular markers closely linked with the resistance loci were developed, and two synthesized hexaploid wheat (SHW) lines were produced. These laid the foundation for cloning of the two resistance loci and for transferring the resistance into common wheat.

A high-resolution genotype-phenotype map identifies the TaSPL17 controlling grain number and size in wheat

BACKGROUND

Large-scale genotype-phenotype association studies of crop germplasm are important for identifying alleles associated with favorable traits. The limited number of single-nucleotide polymorphisms (SNPs) in most wheat genome-wide association studies (GWASs) restricts their power to detect marker-trait associations. Additionally, only a few genes regulating grain number per spikelet have been reported due to sensitivity of this trait to variable environments.

RESULTS

We perform a large-scale GWAS using approximately 40 million filtered SNPs for 27 spike morphology traits. We detect 132,086 significant marker-trait associations and the associated SNP markers are located within 590 associated peaks. We detect additional and stronger peaks by dividing spike morphology into sub-traits relative to GWAS results of spike morphology traits. We propose that the genetic dissection of spike morphology is a powerful strategy to detect signals for grain yield traits in wheat. The GWAS results reveal that TaSPL17 positively controls grain size and number by regulating spikelet and floret meristem development, which in turn leads to enhanced grain yield per plant. The haplotypes at TaSPL17 indicate geographical differentiation, domestication effects, and breeding selection.

CONCLUSION

Our study provides valuable resources for genetic improvement of spike morphology and a fast-forward genetic solution for candidate gene detection and cloning in wheat.

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

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

Whole genome doubling-induced the enrichment of H3K27me3 in genes carrying specific TEs in Aegilops tauschii

Polyploidization plays important roles in the evolution and breeding of the common wheat. Aegilops tauschii, the D-genome progenitor of the common wheat, provides a valuable pool of resistance genes to multiple diseases. Extensive studies focus on the exploration of these genes for wheat improvement. However, few studies have unveiled alternations on genome-wide expression pattern and histone modifications induced by whole-genome doubling (WGD) process. In this study, we conducted transcriptome analysis for the diploid and tetraploid Ae. taushcii lines using the leaf and root tissues. Both lines tend to display similar tissue-specific pattern. Interestingly, we found that TEs located in genic regions were depleted of the repressive histone mark H3K27me3, whereas their adjacent chromatin was enriched with H3K27me3. The tetraploid line exhibited higher levels of H3K27me3 in those regions than the diploid line, particularly for genic regions associated with TEs of the long interspersed nuclear elements (LINEs), CACTA, PIF/Harbinger, Tc1/Mariner and unclassed DNA transposon. Surprisingly, the expression levels of these TEs cognate genes were negatively associated with the levels of H3K27me3 between the tetraploid and diploid lines, suggesting the five types of TEs located within genic regions might be involved in the regulation of the ploidy-related gene expression, possibly through differential enrichment of H3K27me3 in the genic regions. These findings will help to understand the potential role of specific types of TEs on transcription in response to WGD.

Triticeae crop genome biology: an endless frontier

Triticeae, the wheatgrass tribe, includes several major cereal crops and their wild relatives. Major crops within the Triticeae are wheat, barley, rye, and oat, which are important for human consumption, animal feed, and rangeland protection. Species within this tribe are known for their large genomes and complex genetic histories. Powered by recent advances in sequencing technology, researchers worldwide have made progress in elucidating the genomes of Triticeae crops. In addition to assemblies of high-quality reference genomes, pan-genome studies have just started to capture the genomic diversities of these species, shedding light on our understanding of the genetic basis of domestication and environmental adaptation of Triticeae crops. In this review, we focus on recent signs of progress in genome sequencing, pan-genome analyses, and resequencing analysis of Triticeae crops. We also propose future research avenues in Triticeae crop genomes, including identifying genome structure variations, the association of genomic regions with desired traits, mining functions of the non-coding area, introgression of high-quality genes from wild Triticeae resources, genome editing, and integration of genomic resources.

Duplicate Genes Contribute to Variability in Abiotic Stress Resistance in Allopolyploid Wheat

Gene duplication is a universal biological phenomenon that drives genomic variation and diversity, plays a crucial role in plant evolution, and contributes to innovations in genetic engineering and crop development. Duplicated genes participate in the emergence of novel functionality, such as adaptability to new or more severe abiotic stress resistance. Future crop research will benefit from advanced, mechanistic understanding of the effects of gene duplication, especially in the development and deployment of high-performance, stress-resistant, elite wheat lines. In this review, we summarize the current knowledge of gene duplication in wheat, including the principle of gene duplication and its effects on gene function, the diversity of duplicated genes, and how they have functionally diverged. Then, we discuss how duplicated genes contribute to abiotic stress response and the mechanisms of duplication. Finally, we have a future prospects section that discusses the direction of future efforts in the short term regarding the elucidation of replication and retention mechanisms of repetitive genes related to abiotic stress response in wheat, excellent gene function research, and practical applications.

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.

Epigenetic Landscape Is Largely Shaped by Diversiform Transposons in Aegilops tauschii

Transposons (TEs) account for more than 80% of the wheat genome, the highest among all known crop species. They play an important role in shaping the elaborate genomic landscape, which is the key to the speciation of wheat. In this study, we analyzed the association between TEs, chromatin states, and chromatin accessibility in Aegilops tauschii, the D genome donor of bread wheat. We found that TEs contributed to the complex but orderly epigenetic landscape as chromatin states showed diverse distributions on TEs of different orders or superfamilies. TEs also contributed to the chromatin state and openness of potential regulatory elements, affecting the expression of TE-related genes. Some TE superfamilies, such as hAT-Ac, carry active/open chromatin regions. In addition, the histone mark H3K9ac was found to be associated with the accessibility shaped by TEs. These results suggest the role of diversiform TEs in shaping the epigenetic landscape and in gene expression regulation in Aegilops tauschii. This has positive implications for understanding the transposon roles in Aegilops tauschii or the wheat D genome.

A multi-omic resource of wheat seed tissues for nutrient deposition and improvement for human health

As a globally important staple crop, wheat seeds provide us with nutrients and proteins. The trend of healthy dietary has become popular recently, emphasizing the consumption of whole-grain wheat products and the dietary benefits. However, the dynamic changes in nutritional profiles of different wheat seed regions (i.e., the embryo, endosperm and outer layers) during developmental stages and the molecular regulation have not been well studied. Here, we provide this multi-omic resource of wheat seeds and describe the generation, technical assessment and preliminary analyses. This resource includes a time-series RNA-seq dataset of the embryo, endosperm and outer layers of wheat seeds and their corresponding metabolomic dataset, covering the middle and late stages of seed development. Our RNA-seq experiments profile the expression of 63,708 genes, while the metabolomic data includes the abundance of 984 metabolites. We believe that this was the first reported transcriptome and metabolome dataset of wheat seeds that helps understand the molecular regulation of the deposition of beneficial nutrients and hence improvements for nutritional and processing quality traits.

Increasingly amplified stimulation mediated by TaNAC69-B is crucial for the leaf senescence in wheat

Leaf senescence involves massive multidimensional alterations, such as nutrient redistribution, and is closely related to crop yield and quality. No apical meristem, Arabidopsis transcription activation factor, and Cup-shaped cotyledon (NAC)-type transcription factors integrate various signals and modulate an enormous number of target genes to ensure the appropriate progression of leaf senescence. However, few leaf senescence-related NACs have been functionally characterized in wheat. Based on our previous RNA-sequencing (RNA-seq) data, we focused on a NAC family member, TaNAC69-B, which is increasingly expressed during leaf senescence in wheat. Overexpression of TaNAC69-B led to precocious leaf senescence in wheat and Arabidopsis, and affected several agricultural traits in transgenic wheat. Moreover, impaired expression of TaNAC69-B by virus-induced gene silencing retarded the leaf senescence in wheat. By RNA-seq and quantitative real-time polymerase chain reaction analysis, we confirmed that some abscisic acid (ABA) biosynthesis genes, including AAO3 and its ortholog in wheat, TraesCS2B02G270600 (TaAO3-B), were elevated by the overexpression of TaNAC69-B. Consistently, we observed more severe ABA-induced leaf senescence in TaNAC69-B-OE wheat and Arabidopsis plants. Furthermore, we determined that TaNAC69-B bound to the NAC binding site core (CGT) on the promoter regions of AAO3 and TaAO3-B. Moreover, we confirmed elevated ABA levels in TaNAC69-B-OE wheat lines. Although TaNAC69-B shares 39.83% identity (amino acid) with AtNAP, TaNAC69-B did not completely restore the delayed leaf senescence in the atnap mutant. Collectively, our results revealed a positive feedback loop, consisting of TaNAC69-B, ABA biosynthesis and leaf senescence, that is essential for the regulation of leaf senescence in wheat.

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

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

Novel function of a putative TaCOBL ortholog associated with cold response

The plant COBRA protein family plays an important role in secondary cell wall biosynthesis and the orientation of cell expansion. The COBRA gene family has been well studied in Arabidopsis thaliana, maize, rice, etc., but no systematic studies were conducted in wheat. In this study, the full-length sequence of TaCOBLs was obtained by homology cloning from wheat, and a conserved motif analysis confirmed that TaCOBLs belonged to the COBRA protein family. qRT-PCR results showed that the TaCOBL transcripts were induced by abiotic stresses, including cold, drought, salinity, and abscisic acid (ABA). Two haplotypes of TaCOBL-5B (Hap5B-a and Hap5B-b), harboring one indel (----/TATA) in the 5' flanking region (- 550 bp), were found on chromosome 5BS. A co-dominant marker, Ta5BF/Ta5BR, was developed based on the polymorphism of the two TaCOBL-5B haplotypes. Significant correlations between the two TaCOBL-5B haplotypes and cold resistance were observed under four environmental conditions. Hap5B-a, a favored haplotype acquired during wheat polyploidization, may positively contribute to enhanced cold resistance in wheat. Based on the promoter activity analysis, the Hap5B-a promoter containing a TATA-box was more active than that of Hap5B-b without the TATA-box under low temperature. Our study provides valuable information indicating that the TaCOBL genes are associated with cold response in wheat.

Orthologous genes Pm12 and Pm21 from two wild relatives of wheat show evolutionary conservation but divergent powdery mildew resistance

Wheat powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is a devastating disease that threatens wheat production worldwide. Pm12, which originated from Aegilops speltoides, a wild relative of wheat, confers strong resistance to powdery mildew and therefore has potential use in wheat breeding. Using susceptible mutants induced by gamma irradiation, we physically mapped and isolated Pm12 and showed it to be orthologous to Pm21 from Dasypyrum villosum, also a wild relative of wheat. The resistance function of Pm12 was validated via ethyl methanesulfonate mutagenesis, virus-induced gene silencing, and stable genetic transformation. Evolutionary analysis indicates that the Pm12/Pm21 loci in wheat species are relatively conserved but dynamic. Here, we demonstrated that the two orthologous genes, Pm12 and Pm21, possess differential resistance against the same set of Bgt isolates. Overexpression of the coiled-coil domains of both PM12 and PM21 induces cell death in Nicotiana benthamiana leaves. However, their full-length forms display different cell death-inducing activities caused by their distinct intramolecular interactions. Cloning of Pm12 will facilitate its application in wheat breeding programs. This study also gives new insight into two orthologous resistance genes, Pm12 and Pm21, which show different race specificities and intramolecular interaction patterns.

Weed genomics: yielding insights into the genetics of weedy traits for crop improvement

Weeds cause tremendous economic and ecological damage worldwide. The number of genomes established for weed species has sharply increased during the recent decade, with some 26 weed species having been sequenced and de novo genomes assembled. These genomes range from 270 Mb (Barbarea vulgaris) to almost 4.4 Gb (Aegilops tauschii). Importantly, chromosome-level assemblies are now available for 17 of these 26 species, and genomic investigations on weed populations have been conducted in at least 12 species. The resulting genomic data have greatly facilitated studies of weed management and biology, especially origin and evolution. Available weed genomes have indeed revealed valuable weed-derived genetic materials for crop improvement. In this review, we summarize the recent progress made in weed genomics and provide a perspective for further exploitation in this emerging field.

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.

Advances in the Mining of Disease Resistance Genes from Aegilops tauschii and the Utilization in Wheat

Aegilops tauschii is one of the malignant weeds that affect wheat production and is also the wild species ancestor of the D genome of hexaploid wheat (Triticum aestivum, AABBDD). It contains many disease resistance genes that have been lost in the long-term evolution of wheat and is an important genetic resource for the mining and utilization of wheat disease resistance genes. In recent years, the genome sequence of Aegilops tauschii has been preliminarily completed, which has laid a good foundation for the further exploration of wheat disease resistance genes in Aegilops tauschii. There are many studies on disease resistance genes in Aegilops tauschii; in order to provide better help for the disease resistance breeding of wheat, this paper analyzes and reviews the relationship between Aegilops tauschii and wheat, the research progress of Aegilops tauschii, the discovery of disease resistance genes from Aegilops tauschii, and the application of disease resistance genes from Aegilops tauschii to modern wheat breeding, providing a reference for the further exploration and utilization of Aegilops tauschii in wheat disease resistance breeding.

A chromosome-scale genome assembly of Dasypyrum villosum provides insights into its application as a broad-spectrum disease resistance resource for wheat improvement

Dasypyrum villosum is one of the most valuable gene resources in wheat improvement, especially for disease resistance. The mining of favorable genes from D. villosum is frustrated by the lack of a whole genome sequence. In this study, we generated a doubled-haploid line, 91C43^DH, using microspore culture and obtained a 4.05-GB high-quality, chromosome-scale genome assembly for D. villosum. The assembly contains39 727 high-confidence genes, and 85.31% of the sequences are repetitive. Two reciprocal translocation events were detected, and 7VS-4VL is a unique translocation in D. villosum. The prolamin seed storage protein-coding genes were found to be duplicated; in particular, the genes encoding low-molecular-weight glutenin at the Glu-V3 locus were significantly expanded. RNA sequencing (RNA-seq) analysis indicated that, after Blumeria graminearum f.sp tritici (Bgt) inoculation, there were more upregulated genes involved in the pattern-triggered immunity and effector-triggered immunity defense pathways in D. villosum than in Triticum urartu. MNase hypersensitive sequencing (MH-seq) identified two Bgt-inducible MH sites (MHSs), one in the promoter and one in the 3' terminal region of the powdery mildew resistance (Pm) gene NLR1-V. Each site had two subpeaks and they were termed MHS1 (MHS1.1/1.2) and MHS2 (MHS2.1/2.2). Bgt-inducible MHS2.2 was uniquely present in D. villosum, and MHS1.1 was more inducible in D. villosum than in wheat, suggesting that MHSs may be critical for regulation of NLR1-V expression and plant defense. In summary, this study provides a valuable genome resource for functional genomics studies and wheat-D. villosum introgression breeding. The identified regulatory mechanisms may also be exploited to develop new strategies for enhancing Pm resistance by optimizing gene expression in wheat.

Genomic Characterization of hox Genes in Senegalese Sole (Solea senegalensis, Kaup 1858): Clues to Evolutionary Path in Pleuronectiformes

The Senegalese sole (Solea senegalensis, Kaup 1858), a marine flatfish, belongs to the Pleuronectiformes order. It is a commercially important species for fisheries and aquaculture. However, in aquaculture, several production bottlenecks have still to be resolved, including skeletal deformities and high mortality during the larval and juvenile phase. The study aims to characterize the hox gene clusters in S. senegalensis to understand better the developmental and metamorphosis process in this species. Using a BAC library, the clones that contain hox genes were isolated, sequenced by NGS and used as BAC-FISH probes. Subsequently the hox clusters were studied by sequence analysis, comparative genomics, and cytogenetic and phylogenetic analysis. Cytogenetic analysis demonstrated the localization of four BAC clones on chromosome pairs 4, 12, 13, and 16 of the Senegalese sole cytogenomic map. Comparative and phylogenetic analysis showed a highly conserved organization in each cluster and different phylogenetic clustering in each hox cluster. Analysis of structural and repetitive sequences revealed accumulations of polymorphisms mediated by repetitive elements in the hoxba cluster, mainly retroelements. Therefore, a possible loss of the hoxb7a gene can be established in the Pleuronectiformes lineage. This work allows the organization and regulation of hox clusters to be understood, and is a good base for further studies of expression patterns.

An established protocol for generating transgenic wheat for wheat functional genomics via particle bombardment

Wheat is one of the most important food crops in the world and is considered one of the top targets in crop biotechnology. With the high-quality reference genomes of wheat and its relative species and the recent burst of genomic resources in Triticeae, demands to perform gene functional studies in wheat and genetic improvement have been rapidly increasing, requiring that production of transgenic wheat should become a routine technique. While established for more than 20 years, the particle bombardment-mediated wheat transformation has not become routine yet, with only a handful of labs being proficient in this technique. This could be due to, at least partly, the low transformation efficiency and the technical difficulties. Here, we describe the current version of this method through adaptation and optimization. We report the detailed protocol of producing transgenic wheat by the particle gun, including several critical steps, from the selection of appropriate explants (i.e., immature scutella), the preparation of DNA-coated gold particles, and several established strategies of tissue culture. More importantly, with over 20 years of experience in wheat transformation in our lab, we share the many technical details and recommendations and emphasize that the particle bombardment-mediated approach has fewer limitations in genotype dependency and vector construction when compared with the Agrobacterium-mediated methods. The particle bombardment-mediated method has been successful for over 30 wheat genotypes, from the tetraploid durum wheat to the hexaploid common wheat, from modern elite varieties to landraces. In conclusion, the particle bombardment-mediated wheat transformation has demonstrated its potential and wide applications, and the full set of protocol, experience, and successful reports in many wheat genotypes described here will further its impacts, making it a routine and robust technique in crop research labs worldwide.

Adaptive radiation in Orinus, an endemic alpine grass of the Qinghai-Tibet Plateau, based on comparative transcriptomic analysis

The species of Orinus (Poaceae) are important alpine plants with a variety of phenotypic traits and potential usages in molecular breeding toward drought-tolerant forage crops. However, the genetic basis of evolutionary adaption and diversification in the genus is still unclear. In the present study, we obtained transcriptomes for the two most divergent species, O. thoroldii and O. kokonoricus, using the Illumina platform and de novo assembly. In total, we generated 23,029 and 24,086 unigenes with N50 values of 1188 and 1203 for O. thoroldii and O. kokonoricus respectively, and identified 19,005 pairs of putative orthologs between the two species of Orinus. For these orthologs, estimations of non-synonymous/synonymous substitution rate ratios indicated that 568 pairs may be under strongly positive selection (Ka/Ks > 1), and Gene Ontogeny (GO) enrichment analysis revealed that significantly enriched pathways were in DNA repair and resistance to abiotic stress. Meanwhile, the divergence times of species between O. thoroldii and O. kokonoricus occurred 3.2 million years ago (Mya), and the recent evolutionary branch is an allotetraploid species, Cleistogenes songorica. We also detected a Ks peak of ∼0.60 for Orinus. Additionally, we identified 188 pairs of differentially expressed genes (DEGs) between the two species of Orinus, which were significantly enrich in stress resistance and lateral root development. Thus, we considered that the species diversification and evolutionary adaption of this genus was initiated by environmental selection, followed by phenotypic differentiation, finally leading to niche separation in the Qinghai-Tibet Plateau.

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.

Establishment and identification of six wheat-Thinopyrum ponticum disomic addition lines derived from partial amphiploid Xiaoyan 7430

Six wheat-Thinopyrum ponticum disomic addition lines derived from partial amphiploid Xiaoyan 7430 were identified using in situ hybridization and SNP microarray, the homoeologous group and stripe rust resistance of each alien chromosome were determined, and Th. ponticum chromosome-specific markers were developed. Xiaoyan 7430 is a significant partial amphiploid, which is used to set up a bridge for transferring valuable genes from Thinopyrum ponticum (Podp.) Barkworth & D.R. Dewey into common wheat. To accelerate the application of these useful genes in enriching the genetic variability of cultivated wheat by chromosome engineering, a complete set of derived addition lines has been created from Xiaoyan 7430. The chromosome composition of each line was characterized by the combination of genomic in situ hybridization and multicolor fluorescence in situ hybridization (mc-FISH), and the homoeology of each alien chromosome was determined by wheat SNP microarray analysis. Addition line WTA55 with alien group-6 chromosome was evaluated resistant to stripe rust isolates at both the seedling and grain-filling stages (Zadoks scale at z.11 and z.73). Diagnostic marker analysis proved that it could carry a novel stripe rust resistance gene derived from Th. ponticum. Furthermore, a FISH probe and 45 molecular markers specific for alien chromosomes were developed based on specific-locus amplified fragment sequencing (SLAF-seq). Of which 27 markers were separately located on single alien chromosome, and some of them could be used to identify the derived translocation lines. This set of addition lines as well as the molecular markers and the FISH probe will promote the introgression of abundant variation from Th. ponticum into wheat in wheat improvement programs.

Seed position in spikelet as a contributing factor to the success of the winter annual invasive grass Aegilops tauschii

Seed position - dependent effects on seed dormancy/germination are well documented at the inflorescence/infructescence level, but less is known about seeds at different positions within a dispersal unit. For the invasive winter annual grass Aegilops tauschii, we quantified morphology, mass and dormancy/germination of seeds from basal (1), middle (2), and distal (3) positions in two spikelet types (Left and Right). We also investigated seedling emergence, survival, plant size and seed production of plants from seeds in different spikelet positions of two spikelet types under different soil nutrient and water conditions. We found that these seed, seedling and plant traits performed as mirror images between the Left and Right spikelet types. The middle seed was significantly the longest and had the maximum mass, while the basal seed was the shortest and had medium mass. Middle seeds had the highest increase in mass during imbibition and the highest germination percentages and rates, while basal seeds had the lowest. Seedling emergence and survival, plant size and seed production for each position of seeds were highest in the added fertilizer combined with regular watering treatment and lowest in the no fertilizer combined with natural moisture, while height of plants derived from the middle and the distal seeds was significantly higher than that of plants derived from the basal seeds under all soil nutrient and water conditions. Seedling survival, number of tillers per plant and seed production per plant from the middle and distal seeds were significantly lower than those from basal seeds under all soil nutrient and water treatments. The considerable variation in seedling emergence and survival, plant size and seed production between seeds in different positions in the spikelet results in much flexibility in all stages of the life cycle, thereby likely contributing to the invasiveness of A. tauschii.

Quick mapping and characterization of a co-located kernel length and thousand-kernel weight-related QTL in wheat

A co-located KL and TKW-related QTL with no negative effect on PH and AD was rapidly identified using BSA and wheat 660 K SNP array. Its effect was validated in a panel of 218 wheat accessions. Kernel length (KL) and thousand-kernel weight (TKW) of wheat (Triticum aestivum L.) contribute significantly to kernel yield. In the present study, a recombinant inbred line (RIL) population derived from the cross between the wheat line S849-8 with larger kernels and more spikelets per spike and the line SY95-71 was developed. Further, of both the bulked segregant analysis (BSA) and the wheat 660 K single nucleotide polymorphism (SNP) array were used to rapidly identify genomic regions for kernel-related traits from this RIL population. Kompetitive Allele Specific PCR markers were further developed in the SNP-enriched region on the 2D chromosome to construct a genetic map. Both QKL.sicau-SSY-2D for KL and QTKW.sicau-SSY-2D for TKW were identified at multiple environments on chromosome arm 2DL. These two QTLs explained 9.68-23.02% and 6.73-18.32% of the phenotypic variation, respectively. The effects of this co-located QTL were successfully verified in a natural population consisting of 218 Sichuan wheat accessions. Interestingly, the major QTL was significantly and positively correlated with spike length, but did not negatively affect spikelet number per spike (SNS), plant height, or anthesis date. These results indicated that it is possible to synchronously improve kernel weight and SNS by using this QTL. Additionally, several genes associated with kernel development and filling rate were predicted and sequenced in the QTL-containing physical intervals of reference genomes of 'Chinese spring' and Aegilops tauschii. Collectively, these results provide a QTL with great breeding potential and its linked markers which should be helpful for fine mapping and molecular breeding.

Development and application of universal ND-FISH probes for detecting P-genome chromosomes based on Agropyron cristatum transposable elements

Fluorescence in situ hybridization (FISH) is a basic tool that is widely used in cytogenetic research. The detection efficiency of conventional FISH is limited due to its time-consuming nature. Oligonucleotide (oligo) probes with fluorescent labels have been applied in non-denaturing FISH (ND-FISH) assays, which greatly streamline experimental processes and save costs and time. Agropyron cristatum, which contains one basic genome, "P," is a vital wild relative for wheat improvement. However, oligo probes for detecting P-genome chromosomes based on ND-FISH assays have not been reported. In this study, according to the distribution of transposable elements (TEs) in Triticeae genomes, 94 oligo probes were designed based on three types of A. cristatum sequences. ND-FISH validation showed that 12 single oligo probes generated a stable and obvious hybridization signal on whole P chromosomes in the wheat background. To improve signal intensity, mixed probes (Oligo-pAc) were prepared by using the 12 successful probes and validated in the diploid accession A. cristatum Z1842, a small segmental translocation line and six allopolyploid wild relatives containing the P genome. The signals of Oligo-pAc covered the entire chromosomes of A. cristatum and were more intense than those of single probes. The results indicate that Oligo-pAc can replace conventional genomic in situ hybridization (GISH) probes to identify P chromosomes or segments in non-P-genome backgrounds. Finally, we provide a rapid and efficient method specifically for detecting P chromosomes in wheat backgrounds by combining the Oligo-pAc probe with the Oligo-pSc119.2-1 and Oligo-pTa535-1 probes, which can replace conventional sequential GISH/FISH assays. Altogether, we developed a set of oligo probes based on the ND-FISH assays to identify P-genome chromosomes, which can promote utilization of A. cristatum in wheat improvement programs.

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.

Epistatic interaction effect between chromosome 1BL (Yr29) and a novel locus on 2AL facilitating resistance to stripe rust in Chinese wheat Changwu 357-9

Four stable QTL for adult plant resistance were identified in wheat line Changwu 357-9, including a new QTL on 2AL showing significant interaction with Yr29 to reduce stripe rust severity. Stripe rust (yellow rust) is a serious disease of bread wheat (Triticum aestivum L.) worldwide. Genetic resistance is considered the most economical, effective and environmentally friendly method to control the disease and to minimize the use of fungicides. The current study focused on characterizing the components of stripe rust resistance and understanding the interactions in Changwu 357-9 (CW357-9)/Avocet S RIL population. A genetic linkage map constructed using a new GenoBaits Wheat 16K Panel and the 660K SNP array had 5104 polymorphic SNP markers spanning 3533.11 cM. Four stable QTL, consistently identified across five environments, were detected on chromosome arms 1BL, 2AL, 3DS, and 6BS in Changwu357-9. The most effective QTL QYrCW357-1BL was Yr29. The 6BS QTL was identified as Yr78, which has been combined with the 1BL QTL in many wheat cultivars and breeding lines. The novel QTL on 2AL with moderate effect showed a stable and significant epistatic interaction with Yr29. The QTL on 3DL should be same as QYrsn.nwafu-3DL and enriches the overall stripe rust resistance gene pool for breeding. Polymorphisms of flanking AQP markers AX-110020417 (for QYrCW357-1BL), AX-110974948 (for QYrCW357-2AL), AX-109466386 (for QYrCW357-3DL), and AX-109995005 (for QYrCW357-6BS) were evaluated in a diversity panel including 225 wheat cultivars and breeding lines. These results suggested that these high-throughput markers could be used to introduce QYrCW357-1BL, QYrCW357-2AL, QYrCW357-3DL, and QYrCW357-6BS into commercial wheat cultivars. Combinations of these genes with other APR QTL should lead to higher levels of stripe rust resistance along with the beneficial effects of multi-disease resistance gene Yr29 on improving resistance to other diseases.

Cytogenetic identification and molecular marker development of a novel wheat-Thinopyrum ponticum translocation line with powdery mildew resistance

A new wheat-Thinopyrum ponticum translocation line with excellent powdery mildew resistance was produced, and alien-specific PCR markers and FISH probes were developed by SLAF-seq. Powdery mildew is one of the most threatening diseases in wheat production. Thinopyrum ponticum (Podp.) Barkworth and D. R. Dewey, as a wild relative, has been used for wheat genetic improvement for the better part of a century. In view of the good powdery mildew resistance of Th. ponticum, we have been working to transfer the resistance genes from Th. ponticum to wheat by creating translocation lines. In this study, a new wheat-Th. ponticum translocation line with excellent resistance and agronomic performance was developed and through seedling disease evaluation, gene postulation and diagnostic marker analysis proved to carry a novel Pm gene derived from Th. ponticum. Cytogenetic analysis revealed that a small alien segment was translocated to the terminal of chromosome 1D to form new translocation TTh-1DS·1DL chromosome. The translocation breakpoint was determined to lie in 21.5 Mb region of chromosome 1D by using Wheat660K SNP array analysis. Based on specific-locus amplified fragment sequencing (SLAF-seq) technology, eight molecular markers and one repetitive sequence probe were developed, which were specific for Th. ponticum. Fortunately, the probe could be used in distinguishing six alien chromosome pairs in partial amphiploid Xiaoyan 7430 by fluorescence in situ hybridization (FISH). Furthermore, a Thinopyrum-specific oligonucleotide probe was designed depending on the sequence information of the FISH probe. The novel translocation line could be used in wheat disease resistance breeding, and these specific markers and probes will enable wheat breeders to rapidly trace the alien genome with the novel Pm gene(s).

BS-Seq reveals major role of differential CHH methylation during leaf rust resistance in wheat (Triticum aestivum L.)

Epigenetic regulation of the activity of defense genes during onset of diseases or resistance against diseases in plants is an active area of research. In the present study, a pair of wheat NILs for leaf rust resistance gene Lr28 (R) in the background of an Indian cultivar HD2329 (S) was used for a study of DNA methylation mediated regulation of gene expression. Leaf samples were collected at 0 h before (S0 and R0) and 96 h after inoculation (S96 and R96). The DNA samples were subjected to BS-Seq and sequencing data were used for identification of differentially methylated/demethylated regions/genes (DMRs and DMGs). Following four pairs of comparisons were used for this purpose: S0 vs S96; S0 vs R0; R0 vs R96; S96 vs R96. Major role of CHH methylation relative to that of CG and CHG methylation was observed. Some important observations include the following: (i) abundance of CHH methylation among DMRs; (ii) predominance of DMRs in intergenic region, relative to other genomic regions (promoters, exons, introns, TSS and TTS); (iii) abundance of transposable elements (TEs) in DMRs with CHH context; (iv) demethylation mediated high expression of genes during susceptible reaction (S0 vs S96) and methylation mediated low expression of genes during resistant reaction (R0 vs R96 and S96 vs R96); (v) major genes under regulation encode proteins, which differ from those encoded by genes regulated during susceptible reaction and (vi) ~ 500 DMGs carried differential binding sites for H3K4/K27me3 marks suggesting joint involvement of DNA and H3 methylation. Thus, CHH methylation either alone or in combination with histone methylation plays a major role in regulating the expression of genes involved in wheat-leaf rust interaction.

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.

Comparative transcriptome analysis of two contrasting resistant and susceptible Aegilops tauschii accessions to wheat leaf rust (Puccinia triticina) using RNA-sequencing

Leaf rust, caused by Puccinia triticina Eriks., is the most common rust disease of wheat (Triticum aestivum L.) worldwide. Owing to the rapid evolution of virulent pathotypes, new and effective leaf rust resistance sources must be found. Aegilops tauschii, an excellent source of resistance genes to a wide range of diseases and pests, may provide novel routes for resistance to this disease. In this study, we aimed to elucidate the transcriptome of leaf rust resistance in two contrasting resistant and susceptible Ae. tauschii accessions using RNA-sequencing. Gene ontology, analysis of pathway enrichment and transcription factors provided an apprehensible review of differentially expressed genes and highlighted biological mechanisms behind the Aegilops–P. triticina interaction. The results showed the resistant accession could uniquely recognize pathogen invasion and respond precisely via reducing galactosyltransferase and overexpressing chromatin remodeling, signaling pathways, cellular homeostasis regulation, alkaloid biosynthesis pathway and alpha-linolenic acid metabolism. However, the suppression of photosynthetic pathway and external stimulus responses were observed upon rust infection in the susceptible genotype. In particular, this first report of comparative transcriptome analysis offers an insight into the strength and weakness of Aegilops against leaf rust and exhibits a pipeline for future wheat breeding programs.

K-mer counting and curated libraries drive efficient annotation of repeats in plant genomes

The annotation of repetitive sequences within plant genomes can help in the interpretation of observed phenotypes. Moreover, repeat masking is required for tasks such as whole-genome alignment, promoter analysis, or pangenome exploration. Although homology-based annotation methods are computationally expensive, k-mer strategies for masking are orders of magnitude faster. Here, we benchmarked a two-step approach, where repeats were first called by k-mer counting and then annotated by comparison to curated libraries. This hybrid protocol was tested on 20 plant genomes from Ensembl, with the k-mer-based Repeat Detector (Red) and two repeat libraries (REdat, last updated in 2013, and nrTEplants, curated for this work). Custom libraries produced by RepeatModeler were also tested. We obtained repeated genome fractions that matched those reported in the literature but with shorter repeated elements than those produced directly by sequence homology. Inspection of the masked regions that overlapped genes revealed no preference for specific protein domains. Most Red-masked sequences could be successfully classified by sequence similarity, with the complete protocol taking less than 2 h on a desktop Linux box. A guide to curating your own repeat libraries and the scripts for masking and annotating plant genomes can be obtained at https://github.com/Ensembl/plant-scripts.

Genome assembly of a Mesoamerican derived variety of lima bean: a foundational cultivar in the Mid-Atlantic USA

Lima bean, Phaseolus lunatus, is closely related to common bean and is high in fiber and protein, with a low glycemic index. Lima bean is widely grown in the state of Delaware, where late summer and early fall weather are conducive to pod production. The same weather conditions also promote diseases such as pod rot and downy mildew, the latter of which has caused previous epidemics. A better understanding of the genes underlying resistance to this and other pathogens is needed to keep this industry thriving in the region. Our current study sought to sequence, assemble, and annotate a commercially available cultivar called Bridgeton, which could then serve as a reference genome, a basis of comparison to other Phaseolus taxa, and a resource for the identification of potential resistance genes. Combined efforts of sequencing, linkage, and comparative analysis resulted in a 623 Mb annotated assembly for lima bean, as well as a better understanding of an evolutionarily dynamic resistance locus in legumes.

Aegilops tauschii genome assembly Aet v5.0 features greater sequence contiguity and improved annotation

Aegilops tauschii is the donor of the D subgenome of hexaploid wheat and an important genetic resource. The reference-quality genome sequence Aet v4.0 for Ae. tauschii acc. AL8/78 was therefore an important milestone for wheat biology and breeding. Further advances in sequencing acc. AL8/78 and release of the Aet v5.0 sequence assembly are reported here. Two new optical maps were constructed and used in the revision of pseudomolecules. Gaps were closed with Pacific Biosciences long-read contigs, decreasing the gap number by 38,899. Transposable elements and protein-coding genes were reannotated. The number of annotated high-confidence genes was reduced from 39,635 in Aet v4.0 to 32,885 in Aet v5.0. A total of 2245 biologically important genes, including those affecting plant phenology, grain quality, and tolerance of abiotic stresses in wheat, was manually annotated and disease-resistance genes were annotated by a dedicated pipeline. Disease-resistance genes encoding nucleotide-binding site domains, receptor-like protein kinases, and receptor-like proteins were preferentially located in distal chromosome regions, whereas those encoding transmembrane coiled-coil proteins were dispersed more evenly along the chromosomes. Discovery, annotation, and expression analyses of microRNA (miRNA) precursors, mature miRNAs, and phasiRNAs are reported, including miRNA target genes. Other small RNAs, such as hc-siRNAs and tRFs, were characterized. These advances enhance the utility of the Ae. tauschii genome sequence for wheat genetics, biotechnology, and breeding.

Long-range assembly of sequences helps to unravel the genome structure and small variation of the wheat-Haynaldia villosa translocated chromosome 6VS.6AL

Genomics studies in wild species of wheat have been limited due to the lack of references; however, new technologies and bioinformatics tools have much potential to promote genomic research. The wheat-Haynaldia villosa translocation line T6VS·6AL has been widely used as a backbone parent of wheat breeding in China. Therefore, revealing the genome structure of translocation chromosome 6VS·6AL will clarify how this chromosome formed and will help to determine how it affects agronomic traits. In this study, chromosome flow sorting, NGS sequencing and Chicago long-range linkage assembly were innovatively used to produce the assembled sequences of 6VS·6AL, and gene prediction and genome structure characterization at the molecular level were effectively performed. The analysis discovered that the short arm of 6VS·6AL was actually composed of a large distal segment of 6VS, a small proximal segment of 6AS and the centromere of 6A, while the collinear region in 6VS corresponding to 230-260 Mb of 6AS-Ta was deleted when the recombination between 6VS and 6AS occurred. In addition to the molecular mechanism of the increased grain weight and enhanced spike length produced by the translocation chromosome, it may be correlated with missing GW2-V and an evolved NRT-V cluster. Moreover, a fine physical bin map of 6VS was constructed by the high-throughput developed 6VS-specific InDel markers and a series of newly identified small fragment translocation lines involving 6VS. This study will provide essential information for mining of new alien genes carried by the 6VS·6AL translocation chromosome.

Present and future prospects for wheat improvement through genome editing and advanced technologies

Wheat (Triticum aestivum, 2n = 6x = 42, AABBDD) is one of the most important staple food crops in the world. Despite the fact that wheat production has significantly increased over the past decades, future wheat production will face unprecedented challenges from global climate change, increasing world population, and water shortages in arid and semi-arid lands. Furthermore, excessive applications of diverse fertilizers and pesticides are exacerbating environmental pollution and ecological deterioration. To ensure global food and ecosystem security, it is essential to enhance the resilience of wheat production while minimizing environmental pollution through the use of cutting-edge technologies. However, the hexaploid genome and gene redundancy complicate advances in genetic research and precision gene modifications for wheat improvement, thus impeding the breeding of elite wheat cultivars. In this review, we first introduce state-of-the-art genome-editing technologies in crop plants, especially wheat, for both functional genomics and genetic improvement. We then outline applications of other technologies, such as GWAS, high-throughput genotyping and phenotyping, speed breeding, and synthetic biology, in wheat. Finally, we discuss existing challenges in wheat genome editing and future prospects for precision gene modifications using advanced genome-editing technologies. We conclude that the combination of genome editing and other molecular breeding strategies will greatly facilitate genetic improvement of wheat for sustainable global production.

A transposon in the vacuolar sorting receptor gene TaVSR1-B promoter region is associated with wheat root depth at booting stage

Root depth, as an important component of root architecture, plays a significant role in growth, grain yield determination and abiotic stress tolerance in crop plants, but its genetic basis remains poorly elucidated. In this study, a panel composed of 323 wheat (Triticum aestivum L.) accessions was assessed for variation in root depth and genotyped with the Wheat 660K SNP Array. GWAS (genome-wide association study) detected significant association between a 125 bp miniature inverted-repeat transposable element (MITE) in the promoter of the TaVSR1-B gene with root depth at the booting stage. We showed that the MITE repressed TaVSR1-B expression by DNA methylation and H3K27 tri-methylation. The roles of TaVSR1-B in root growth were verified by altered expression of the gene in transgenic wheat, rice and a tavsr1 TILLING mutant. Increased TaVSR1-B expression made the root elongation zone shorter and the differentiation zone longer, leading to deeper root. This work provides novel insight into the genetic basis of variation in root depth and a promising target for genetic improvement of root architecture in wheat.

Sequence of the supernumerary B chromosome of maize provides insight into its drive mechanism and evolution

B chromosomes are enigmatic elements in thousands of plant and animal genomes that persist in populations despite being nonessential. They circumvent the laws of Mendelian inheritance but the molecular mechanisms underlying this behavior remain unknown. Here we present the sequence, annotation, and analysis of the maize B chromosome providing insight into its drive mechanism. The sequence assembly reveals detailed locations of the elements involved with the cis and trans functions of its drive mechanism, consisting of nondisjunction at the second pollen mitosis and preferential fertilization of the egg by the B-containing sperm. We identified 758 protein-coding genes in 125.9 Mb of B chromosome sequence, of which at least 88 are expressed. Our results demonstrate that transposable elements in the B chromosome are shared with the standard A chromosome set but multiple lines of evidence fail to detect a syntenic genic region in the A chromosomes, suggesting a distant origin. The current gene content is a result of continuous transfer from the A chromosomal complement over an extended evolutionary time with subsequent degradation but with selection for maintenance of this nonvital chromosome.