The dynamics of touch-responsive gene expression in cereals.

Wind, rain, herbivores, obstacles, neighbouring plants, etc. provide important mechanical cues to steer plant growth and survival. Mechanostimulation to stimulate yield and stress resistance of crops is of significant research interest, yet a molecular understanding of transcriptional responses to touch is largely absent in cereals. To address this, we performed whole-genome transcriptomics following mechanostimulation of wheat, barley, and the recent genome-sequenced oat. The largest transcriptome changes occurred ±25 min after touching, with most of the genes being upregulated. While most genes returned to basal expression level by 1-2 h in oat, many genes retained high expression even 4 h post-treatment in barley and wheat. Functional categories such as transcription factors, kinases, phytohormones, and Ca2+ regulation were affected. In addition, cell wall-related genes involved in (hemi)cellulose, lignin, suberin, and callose biosynthesis were touch-responsive, providing molecular insight into mechanically induced changes in cell wall composition. Furthermore, several cereal-specific transcriptomic footprints were identified that were not observed in Arabidopsis. In oat and barley, we found evidence for systemic spreading of touch-induced signalling. Finally, we provide evidence that both the jasmonic acid-dependent and the jasmonic acid-independent pathways underlie touch-signalling in cereals, providing a detailed framework and marker genes for further study of (a)biotic stress responses in cereals.

Chromosome-level genome assembly and population genomic resource to accelerate orphan crop lablab breeding.

Under-utilised orphan crops hold the key to diversified and climate-resilient food systems. Here, we report on orphan crop genomics using the case of Lablab purpureus (L.) Sweet (lablab) - a legume native to Africa and cultivated throughout the tropics for food and forage. Our Africa-led plant genome collaboration produces a high-quality chromosome-scale assembly of the lablab genome. Our assembly highlights the genome organisation of the trypsin inhibitor genes - an important anti-nutritional factor in lablab. We also re-sequence cultivated and wild lablab accessions from Africa confirming two domestication events. Finally, we examine the genetic and phenotypic diversity in a comprehensive lablab germplasm collection and identify genomic loci underlying variation of important agronomic traits in lablab. The genomic data generated here provide a valuable resource for lablab improvement. Our inclusive collaborative approach also presents an example that can be explored by other researchers sequencing indigenous crops, particularly from low and middle-income countries (LMIC).

A reference-guided TILLING by amplicon-sequencing platform supports forward and reverse genetics in barley.

Barley is a diploid species with a genome smaller than those of other members of the Triticeae tribe, making it an attractive model for genetic studies in Triticeae crops. The recent development of barley genomics has created a need for a high-throughput platform to identify genetically uniform mutants for gene function investigations. In this study, we report an ethyl methanesulfonate (EMS)-mutagenized population consisting of 8525 M3 lines in the barley landrace "Hatiexi" (HTX), which we complement with a high-quality de novo assembly of a reference genome for this genotype. The mutation rate within the population ranged from 1.51 to 4.09 mutations per megabase, depending on the treatment dosage of EMS and the mutation discrimination platform used for genotype analysis. We implemented a three-dimensional DNA pooling strategy combined with multiplexed amplicon sequencing to create a highly efficient and cost-effective TILLING (targeting induced locus lesion in genomes) platform in barley. Mutations were successfully identified from 72 mixed amplicons within a DNA pool containing 64 individual mutants and from 56 mixed amplicons within a pool containing 144 individuals. We discovered abundant allelic mutants for dozens of genes, including the barley Green Revolution contributor gene Brassinosteroid insensitive 1 (BRI1). As a proof of concept, we rapidly determined the causal gene responsible for a chlorotic mutant by following the MutMap strategy, demonstrating the value of this resource to support forward and reverse genetic studies in barley.

Genome analysis in Avena sativa reveals hidden breeding barriers and opportunities for oat improvement.

Oat (Avena sativa L.) is an important and nutritious cereal crop, and there is a growing need to identify genes that contribute to improved oat varieties. Here we utilize a newly sequenced and annotated oat reference genome to locate and characterize quantitative trait loci (QTLs) affecting agronomic and grain-quality traits in five oat populations. We find strong and significant associations between the positions of candidate genes and QTL that affect heading date, as well as those that influence the concentrations of oil and ß-glucan in the grain. We examine genome-wide recombination profiles to confirm the presence of a large, unbalanced translocation from chromosome 1 C to 1 A, and a possible inversion on chromosome 7D. Such chromosome rearrangements appear to be common in oat, where they cause pseudo-linkage and recombination suppression, affecting the segregation, localization, and deployment of QTLs in breeding programs.

The mosaic oat genome gives insights into a uniquely healthy cereal crop.

Cultivated oat (Avena sativa L.) is an allohexaploid (AACCDD, 2n = 6x = 42) thought to have been domesticated more than 3,000 years ago while growing as a weed in wheat, emmer and barley fields in Anatolia1,2. Oat has a low carbon footprint, substantial health benefits and the potential to replace animal-based food products. However, the lack of a fully annotated reference genome has hampered efforts to deconvolute its complex evolutionary history and functional gene dynamics. Here we present a high-quality reference genome of A. sativa and close relatives of its diploid (Avena longiglumis, AA, 2n = 14) and tetraploid (Avena insularis, CCDD, 2n = 4x = 28) progenitors. We reveal the mosaic structure of the oat genome, trace large-scale genomic reorganizations in the polyploidization history of oat and illustrate a breeding barrier associated with the genome architecture of oat. We showcase detailed analyses of gene families implicated in human health and nutrition, which adds to the evidence supporting oat safety in gluten-free diets, and we perform mapping-by-sequencing of an agronomic trait related to water-use efficiency. This resource for the Avena genus will help to leverage knowledge from other cereal genomes, improve understanding of basic oat biology and accelerate genomics-assisted breeding and reanalysis of quantitative trait studies.

The Barley and Wheat Pan-Genomes.

To unlock the genetic potential in crops, multi-genome comparisons are an essential tool. Decreasing costs and improved sequencing technologies have democratized plant genome sequencing and led to a vast increase in the amount of available reference sequences on the one hand and enabled the assembly of even the largest and most complex and repetitive crops genomes such as wheat and barley. These developments have led to the era of pan-genomics in recent years. Pan-genome projects enable the definition of the core and dispensable genome for various crop species as well as the analysis of structural and functional variation and hence offer unprecedented opportunities for exploring and utilizing the genetic basis of natural variation in crops. Comparing, analyzing, and visualizing these multiple reference genomes and their diversity requires powerful and specialized computational strategies and tools.

Chromosome-scale assembly of barley cv. 'Haruna Nijo' as a resource for barley genetics.

Cultivated barley (Hordeum vulgare ssp. vulgare) is used for food, animal feed, and alcoholic beverages and is widely grown in temperate regions. Both barley and its wild progenitor (H. vulgare ssp. spontaneum) have large 5.1-Gb genomes. High-quality chromosome-scale assemblies for several representative barley genotypes, both wild and domesticated, have been constructed recently to populate the nascent barley pan-genome infrastructure. Here, we release a chromosome-scale assembly of the Japanese elite malting barley cultivar 'Haruna Nijo' using a similar methodology as in the barley pan-genome project. The 4.28-Gb assembly had a scaffold N50 size of 18.9 Mb. The assembly showed high collinearity with the barley reference genome 'Morex' cultivar, with some inversions. The pseudomolecule assembly was characterized using transcript evidence of gene projection derived from the reference genome and de novo gene annotation achieved using published full-length cDNA sequences and RNA-Seq data for 'Haruna Nijo'. We found good concordance between our whole-genome assembly and the publicly available BAC clone sequence of 'Haruna Nijo'. Interesting phenotypes have since been identified in Haruna Nijo; its genome sequence assembly will facilitate the identification of the underlying genes.

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.

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.

Chromosome-scale assembly of wild barley accession "OUH602".

Barley (Hordeum vulgare) was domesticated from its wild ancestral form ca. 10,000 years ago in the Fertile Crescent and is widely cultivated throughout the world, except for in tropical areas. The genome size of both cultivated barley and its conspecific wild ancestor is approximately 5 Gb. High-quality chromosome-level assemblies of 19 cultivated and one wild barley genotype were recently established by pan-genome analysis. Here, we release another equivalent short-read assembly of the wild barley accession "OUH602." A series of genetic and genomic resources were developed for this genotype in prior studies. Our assembly contains more than 4.4 Gb of sequence, with a scaffold N50 value of over 10 Mb. The haplotype shows high collinearity with the most recently updated barley reference genome, "Morex" V3, with some inversions. Gene projections based on "Morex" gene models revealed 46,807 protein-coding sequences and 43,375 protein-coding genes. Alignments to publicly available sequences of bacterial artificial chromosome (BAC) clones of "OUH602" confirm the high accuracy of the assembly. Since more loci of interest have been identified in "OUH602," the release of this assembly, with detailed genomic information, should accelerate gene identification and the utilization of this key wild barley accession.

Chromosome-scale genome assembly of the transformation-amenable common wheat cultivar 'Fielder'.

We have established a high-quality, chromosome-level genome assembly for the hexaploid common wheat cultivar 'Fielder', an American, soft, white, pastry-type wheat released in 1974 and known for its amenability to Agrobacterium tumefaciens-mediated transformation and genome editing. Accurate, long-read sequences were obtained using PacBio circular consensus sequencing with the HiFi approach. Sequence reads from 16 SMRT cells assembled using the hifiasm assembler produced assemblies with N50 greater than 20 Mb. We used the Omni-C chromosome conformation capture technique to order contigs into chromosome-level assemblies, resulting in 21 pseudomolecules with a cumulative size of 14.7 and 0.3 Gb of unanchored contigs. Mapping of published short reads from a transgenic wheat plant with an edited seed-dormancy gene, TaQsd1, identified four positions of transgene insertion into wheat chromosomes. Detection of guide RNA sequences in pseudomolecules provided candidates for off-target mutation induction. These results demonstrate the efficiency of chromosome-scale assembly using PacBio HiFi reads and their application in wheat genome-editing studies.

Long-read sequence assembly: a technical evaluation in barley.

Sequence assembly of large and repeat-rich plant genomes has been challenging, requiring substantial computational resources and often several complementary sequence assembly and genome mapping approaches. The recent development of fast and accurate long-read sequencing by circular consensus sequencing (CCS) on the PacBio platform may greatly increase the scope of plant pan-genome projects. Here, we compare current long-read sequencing platforms regarding their ability to rapidly generate contiguous sequence assemblies in pan-genome studies of barley (Hordeum vulgare). Most long-read assemblies are clearly superior to the current barley reference sequence based on short-reads. Assemblies derived from accurate long reads excel in most metrics, but the CCS approach was the most cost-effective strategy for assembling tens of barley genomes. A downsampling analysis indicated that 20-fold CCS coverage can yield very good sequence assemblies, while even five-fold CCS data may capture the complete sequence of most genes. We present an updated reference genome assembly for barley with near-complete representation of the repeat-rich intergenic space. Long-read assembly can underpin the construction of accurate and complete sequences of multiple genomes of a species to build pan-genome infrastructures in Triticeae crops and their wild relatives.

De Novo Genome Assembly of the Japanese Wheat Cultivar Norin 61 Highlights Functional Variation in Flowering Time and Fusarium-Resistant Genes in East Asian Genotypes.

Bread wheat is a major crop that has long been the focus of basic and breeding research. Assembly of its genome has been difficult because of its large size and allohexaploid nature (AABBDD genome). Following the first reported assembly of the genome of the experimental strain Chinese Spring (CS), the 10+ Wheat Genomes Project was launched to produce multiple assemblies of worldwide modern cultivars. The only Asian cultivar in the project is Norin 61, a representative Japanese cultivar adapted to grow across a broad latitudinal range, mostly characterized by a wet climate and a short growing season. Here, we characterize the key aspects of its chromosome-scale genome assembly spanning 15 Gb with a raw scaffold N50 of 22 Mb. Analysis of the repetitive elements identified chromosomal regions unique to Norin 61 that encompass a tandem array of the pathogenesis-related 13 family. We report novel copy-number variations in the B homeolog of the florigen gene FT1/VRN3, pseudogenization of its D homeolog and the association of its A homeologous alleles with the spring/winter growth habit. Furthermore, the Norin 61 genome carries typical East Asian functional variants different from CS, ranging from a single nucleotide to multi-Mb scale. Examples of such variation are the Fhb1 locus, which confers Fusarium head-blight resistance, Ppd-D1a, which confers early flowering, Glu-D1f for Asian noodle quality and Rht-D1b, which introduced semi-dwarfism during the green revolution. The adoption of Norin 61 as a reference assembly for functional and evolutionary studies will enable comprehensive characterization of the underexploited Asian bread wheat diversity.

The Aegilops ventricosa 2NvS segment in bread wheat: cytology, genomics and breeding.

KEY MESSAGE: The first cytological characterization of the 2NvS segment in hexaploid wheat; complete de novo assembly and annotation of 2NvS segment; 2NvS frequency is increasing 2NvS and is associated with higher yield. The Aegilops ventricosa 2NvS translocation segment has been utilized in breeding disease-resistant wheat crops since the early 1990s. This segment is known to possess several important resistance genes against multiple wheat diseases including root knot nematode, stripe rust, leaf rust and stem rust. More recently, this segment has been associated with resistance to wheat blast, an emerging and devastating wheat disease in South America and Asia. To date, full characterization of the segment including its size, gene content and its association with grain yield is lacking. Here, we present a complete cytological and physical characterization of this agronomically important translocation in bread wheat. We de novo assembled the 2NvS segment in two wheat varieties, 'Jagger' and 'CDC Stanley,' and delineated the segment to be approximately 33 Mb. A total of 535 high-confidence genes were annotated within the 2NvS region, with > 10% belonging to the nucleotide-binding leucine-rich repeat (NLR) gene families. Identification of groups of NLR genes that are potentially N genome-specific and expressed in specific tissues can fast-track testing of candidate genes playing roles in various disease resistances. We also show the increasing frequency of 2NvS among spring and winter wheat breeding programs over two and a half decades, and the positive impact of 2NvS on wheat grain yield based on historical datasets. The significance of the 2NvS segment in wheat breeding due to resistance to multiple diseases and a positive impact on yield highlights the importance of understanding and characterizing the wheat pan-genome for better insights into molecular breeding for wheat improvement.

Comparative Genomics and Functional Studies of Wheat BED-NLR Loci.

Nucleotide-binding leucine-rich-repeat (LRR) receptors (NLRs) with non-canonical integrated domains (NLR-IDs) are widespread in plant genomes. Zinc-finger BED (named after the Drosophila proteins Boundary Element-Associated Factor and DNA Replication-related Element binding Factor, named BED hereafter) are among the most frequently found IDs. Five BED-NLRs conferring resistance against bacterial and fungal pathogens have been characterized. However, it is unknown whether BED-NLRs function in a manner similar to other NLR-IDs. Here, we used chromosome-level assemblies of wheat to explore the Yr7 and Yr5a genomic regions and show that, unlike known NLR-ID loci, there is no evidence for a NLR-partner in their vicinity. Using neighbor-network analyses, we observed that BED domains from BED-NLRs share more similarities with BED domains from single-BED proteins and from BED-containing proteins harboring domains that are conserved in transposases. We identified a nuclear localization signal (NLS) in Yr7, Yr5, and the other characterized BED-NLRs. We thus propose that this is a feature of BED-NLRs that confer resistance to plant pathogens. We show that the NLS was functional in truncated versions of the Yr7 protein when expressed in N. benthamiana. We did not observe cell-death upon the overexpression of Yr7 full-length, truncated, and 'MHD' variants in N. benthamiana. This suggests that either this system is not suitable to study BED-NLR signaling or that BED-NLRs require additional components to trigger cell death. These results define novel future directions to further understand the role of BED domains in BED-NLR mediated resistance.

A haplotype-led approach to increase the precision of wheat breeding.

Crop productivity must increase at unprecedented rates to meet the needs of the growing worldwide population. Exploiting natural variation for the genetic improvement of crops plays a central role in increasing productivity. Although current genomic technologies can be used for high-throughput identification of genetic variation, methods for efficiently exploiting this genetic potential in a targeted, systematic manner are lacking. Here, we developed a haplotype-based approach to identify genetic diversity for crop improvement using genome assemblies from 15 bread wheat (Triticum aestivum) cultivars. We used stringent criteria to identify identical-by-state haplotypes and distinguish these from near-identical sequences (~99.95% identity). We showed that each cultivar shares ~59 % of its genome with other sequenced cultivars and we detected the presence of extended haplotype blocks containing hundreds to thousands of genes across all wheat chromosomes. We found that genic sequence alone was insufficient to fully differentiate between haplotypes, as were commonly used array-based genotyping chips due to their gene centric design. We successfully used this approach for focused discovery of novel haplotypes from a landrace collection and documented their potential for trait improvement in modern bread wheat. This study provides a framework for defining and exploiting haplotypes to increase the efficiency and precision of wheat breeding towards optimising the agronomic performance of this crucial crop.

The barley pan-genome reveals the hidden legacy of mutation breeding.

Genetic diversity is key to crop improvement. Owing to pervasive genomic structural variation, a single reference genome assembly cannot capture the full complement of sequence diversity of a crop species (known as the 'pan-genome'1). Multiple high-quality sequence assemblies are an indispensable component of a pan-genome infrastructure. Barley (Hordeum vulgare L.) is an important cereal crop with a long history of cultivation that is adapted to a wide range of agro-climatic conditions2. Here we report the construction of chromosome-scale sequence assemblies for the genotypes of 20 varieties of barley-comprising landraces, cultivars and a wild barley-that were selected as representatives of global barley diversity. We catalogued genomic presence/absence variants and explored the use of structural variants for quantitative genetic analysis through whole-genome shotgun sequencing of 300 gene bank accessions. We discovered abundant large inversion polymorphisms and analysed in detail two inversions that are frequently found in current elite barley germplasm; one is probably the product of mutation breeding and the other is tightly linked to a locus that is involved in the expansion of geographical range. This first-generation barley pan-genome makes previously hidden genetic variation accessible to genetic studies and breeding.

European maize genomes highlight intraspecies variation in repeat and gene content.

The diversity of maize (Zea mays) is the backbone of modern heterotic patterns and hybrid breeding. Historically, US farmers exploited this variability to establish today's highly productive Corn Belt inbred lines from blends of dent and flint germplasm pools. Here, we report de novo genome sequences of four European flint lines assembled to pseudomolecules with scaffold N50 ranging from 6.1 to 10.4 Mb. Comparative analyses with two US Corn Belt lines explains the pronounced differences between both germplasms. While overall syntenic order and consolidated gene annotations reveal only moderate pangenomic differences, whole-genome alignments delineating the core and dispensable genome, and the analysis of heterochromatic knobs and orthologous long terminal repeat retrotransposons unveil the dynamics of the maize genome. The high-quality genome sequences of the flint pool complement the maize pangenome and provide an important tool to study maize improvement at a genome scale and to enhance modern hybrid breeding.

TRITEX: chromosome-scale sequence assembly of Triticeae genomes with open-source tools.

Chromosome-scale genome sequence assemblies underpin pan-genomic studies. Recent genome assembly efforts in the large-genome Triticeae crops wheat and barley have relied on the commercial closed-source assembly algorithm DeNovoMagic. We present TRITEX, an open-source computational workflow that combines paired-end, mate-pair, 10X Genomics linked-read with chromosome conformation capture sequencing data to construct sequence scaffolds with megabase-scale contiguity ordered into chromosomal pseudomolecules. We evaluate the performance of TRITEX on publicly available sequence data of tetraploid wild emmer and hexaploid bread wheat, and construct an improved annotated reference genome sequence assembly of the barley cultivar Morex as a community resource.

Tracing the ancestry of modern bread wheats.

For more than 10,000 years, the selection of plant and animal traits that are better tailored for human use has shaped the development of civilizations. During this period, bread wheat (Triticum aestivum) emerged as one of the world's most important crops. We use exome sequencing of a worldwide panel of almost 500 genotypes selected from across the geographical range of the wheat species complex to explore how 10,000 years of hybridization, selection, adaptation and plant breeding has shaped the genetic makeup of modern bread wheats. We observe considerable genetic variation at the genic, chromosomal and subgenomic levels, and use this information to decipher the likely origins of modern day wheats, the consequences of range expansion and the allelic variants selected since its domestication. Our data support a reconciled model of wheat evolution and provide novel avenues for future breeding improvement.