Plant pangenomes for crop improvement, biodiversity and evolution.

Plant genome sequences catalogue genes and the genetic elements that regulate their expression. Such inventories further research aims as diverse as mapping the molecular basis of trait diversity in domesticated plants or inquiries into the origin of evolutionary innovations in flowering plants millions of years ago. The transformative technological progress of DNA sequencing in the past two decades has enabled researchers to sequence ever more genomes with greater ease. Pangenomes - complete sequences of multiple individuals of a species or higher taxonomic unit - have now entered the geneticists' toolkit. The genomes of crop plants and their wild relatives are being studied with translational applications in breeding in mind. But pangenomes are applicable also in ecological and evolutionary studies, as they help classify and monitor biodiversity across the tree of life, deepen our understanding of how plant species diverged and show how plants adapt to changing environments or new selection pressures exerted by human beings.

Evolution of the Grain Dispersal System in Barley.

About 12,000 years ago in the Near East, humans began the transition from hunter-gathering to agriculture-based societies. Barley was a founder crop in this process, and the most important steps in its domestication were mutations in two adjacent, dominant, and complementary genes, through which grains were retained on the inflorescence at maturity, enabling effective harvesting. Independent recessive mutations in each of these genes caused cell wall thickening in a highly specific grain "disarticulation zone," converting the brittle floral axis (the rachis) of the wild-type into a tough, non-brittle form that promoted grain retention. By tracing the evolutionary history of allelic variation in both genes, we conclude that spatially and temporally independent selections of germplasm with a non-brittle rachis were made during the domestication of barley by farmers in the southern and northern regions of the Levant, actions that made a major contribution to the emergence of early agrarian societies.

The giant diploid faba genome unlocks variation in a global protein crop.

Increasing the proportion of locally produced plant protein in currently meat-rich diets could substantially reduce greenhouse gas emissions and loss of biodiversity1. However, plant protein production is hampered by the lack of a cool-season legume equivalent to soybean in agronomic value2. Faba bean (Vicia faba L.) has a high yield potential and is well suited for cultivation in temperate regions, but genomic resources are scarce. Here, we report a high-quality chromosome-scale assembly of the faba bean genome and show that it has expanded to a massive 13 Gb in size through an imbalance between the rates of amplification and elimination of retrotransposons and satellite repeats. Genes and recombination events are evenly dispersed across chromosomes and the gene space is remarkably compact considering the genome size, although with substantial copy number variation driven by tandem duplication. Demonstrating practical application of the genome sequence, we develop a targeted genotyping assay and use high-resolution genome-wide association analysis to dissect the genetic basis of seed size and hilum colour. The resources presented constitute a genomics-based breeding platform for faba bean, enabling breeders and geneticists to accelerate the improvement of sustainable protein production across the Mediterranean, subtropical and northern temperate agroecological zones.

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 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.

Dynamic Phytomeric Growth Contributes to Local Adaptation in Barley.

Vascular plants have segmented body axes with iterative nodes and internodes. Appropriate node initiation and internode elongation are fundamental to plant fitness and crop yield; however, how these events are spatiotemporally coordinated remains elusive. We show that in barley (Hordeum vulgare L.), selections during domestication have extended the apical meristematic phase to promote node initiation, but constrained subsequent internode elongation. In both vegetative and reproductive phases, internode elongation displays a dynamic proximal-distal gradient, and among subpopulations of domesticated barleys worldwide, node initiation and proximal internode elongation are associated with latitudinal and longitudinal gradients, respectively. Genetic and functional analyses suggest that, in addition to their converging roles in node initiation, flowering-time genes have been repurposed to specify the timing and duration of internode elongation. Our study provides an integrated view of barley node initiation and internode elongation and suggests that plant architecture should be recognized as a collection of dynamic phytomeric units in the context of crop adaptive evolution.

Helical coiling of metaphase chromatids.

Chromatids of mitotic chromosomes were suggested to coil into a helix in early cytological studies and this assumption was recently supported by chromosome conformation capture (3C) sequencing. Still, direct differential visualization of a condensed chromatin fibre confirming the helical model was lacking. Here, we combined Hi-C analysis of purified metaphase chromosomes, biopolymer modelling and spatial structured illumination microscopy of large fluorescently labeled chromosome segments to reveal the chromonema - a helically-wound, 400 nm thick chromatin thread forming barley mitotic chromatids. Chromatin from adjacent turns of the helix intermingles due to the stochastic positioning of chromatin loops inside the chromonema. Helical turn size varies along chromosome length, correlating with chromatin density. Constraints on the observable dimensions of sister chromatid exchanges further supports the helical chromonema model.

Root angle is controlled by EGT1 in cereal crops employing an antigravitropic mechanism.

Root angle in crops represents a key trait for efficient capture of soil resources. Root angle is determined by competing gravitropic versus antigravitropic offset (AGO) mechanisms. Here we report a root angle regulatory gene termed ENHANCED GRAVITROPISM1 (EGT1) that encodes a putative AGO component, whose loss-of-function enhances root gravitropism. Mutations in barley and wheat EGT1 genes confer a striking root phenotype, where every root class adopts a steeper growth angle. EGT1 encodes an F-box and Tubby domain-containing protein that is highly conserved across plant species. Haplotype analysis found that natural allelic variation at the barley EGT1 locus impacts root angle. Gravitropic assays indicated that Hvegt1 roots bend more rapidly than wild-type. Transcript profiling revealed Hvegt1 roots deregulate reactive oxygen species (ROS) homeostasis and cell wall-loosening enzymes and cofactors. ROS imaging shows that Hvegt1 root basal meristem and elongation zone tissues have reduced levels. Atomic force microscopy measurements detected elongating Hvegt1 root cortical cell walls are significantly less stiff than wild-type. In situ analysis identified HvEGT1 is expressed in elongating cortical and stele tissues, which are distinct from known root gravitropic perception and response tissues in the columella and epidermis, respectively. We propose that EGT1 controls root angle by regulating cell wall stiffness in elongating root cortical tissue, counteracting the gravitropic machinery's known ability to bend the root via its outermost tissues. We conclude that root angle is controlled by EGT1 in cereal crops employing an antigravitropic mechanism.

The auxin efflux carrier PIN1a regulates vascular patterning in cereal roots.

Barley (Hordeum vulgare) is an important global cereal crop and a model in genetic studies. Despite advances in characterising barley genomic resources, few mutant studies have identified genes controlling root architecture and anatomy, which plays a critical role in capturing soil resources. Our phenotypic screening of a TILLING mutant collection identified line TM5992 exhibiting a short-root phenotype compared with wild-type (WT) Morex background. Outcrossing TM5992 with barley variety Proctor and subsequent SNP array-based bulk segregant analysis, fine mapped the mutation to a cM scale. Exome sequencing pinpointed a mutation in the candidate gene HvPIN1a, further confirming this by analysing independent mutant alleles. Detailed analysis of root growth and anatomy in Hvpin1a mutant alleles exhibited a slower growth rate, shorter apical meristem and striking vascular patterning defects compared to WT. Expression and mutant analyses of PIN1 members in the closely related cereal brachypodium (Brachypodium distachyon) revealed that BdPIN1a and BdPIN1b were redundantly expressed in root vascular tissues but only Bdpin1a mutant allele displayed root vascular defects similar to Hvpin1a. We conclude that barley PIN1 genes have sub-functionalised in cereals, compared to Arabidopsis (Arabidopsis thaliana), where PIN1a sequences control root vascular patterning.

The crane fly glycosylated triketide δ-lactone cornicinine elicits akinete differentiation of the cyanobiont in aquatic Azolla fern symbioses.

The restriction of plant-symbiont dinitrogen fixation by an insect semiochemical had not been previously described. Here we report on a glycosylated triketide δ-lactone from Nephrotoma cornicina crane flies, cornicinine, that causes chlorosis in the floating-fern symbioses from the genus Azolla. Only the glycosylated trans-A form of chemically synthesized cornicinine was active: 500 nM cornicinine in the growth medium turned all cyanobacterial filaments from Nostoc azollae inside the host leaf-cavities into akinetes typically secreting CTB-bacteriocins. Cornicinine further inhibited akinete germination in Azolla sporelings, precluding re-establishment of the symbiosis during sexual reproduction. It did not impact development of the plant Arabidopsis thaliana or several free-living cyanobacteria from the genera Anabaena or Nostoc but affected the fern host without cyanobiont. Fern-host mRNA sequencing from isolated leaf cavities confirmed high NH4-assimilation and proanthocyanidin biosynthesis in this trichome-rich tissue. After cornicinine treatment, it revealed activation of Cullin-RING ubiquitin-ligase-pathways, known to mediate metabolite signaling and plant elicitation consistent with the chlorosis phenotype, and increased JA-oxidase, sulfate transport and exosome formation. The work begins to uncover molecular mechanisms of cyanobiont differentiation in a seed-free plant symbiosis important for wetland ecology or circular crop-production today, that once caused massive CO2 draw-down during the Eocene geological past.

Hi-C techniques: from genome assemblies to transcription regulation.

The invention of chromosome-conformation capture (3C) techniques, in particular the key method Hi-C providing genome-wide information about chromatin contacts, revolutionized the way we study the three-dimensional (3D) organization of the nuclear genome and how it impacts transcription, replication and DNA repair. Since the frequency of chromatin contacts between pairs of genomic segments predictably relates to the distance in the linear genome, the Hi-C information has also proved useful for scaffolding genomic sequences. Here, we review recent enhancements in experimental procedures of Hi-C and its various derivatives such as Micro-C, HiChIP, and Capture Hi-C. We assess advantages and limitations of the techniques, and present examples of their use in recent plant studies. We also report on progress in computational tools used in assembling genome sequences.

Capitalizing genebank core collections for rare and novel disease resistance loci to enhance barley resilience.

In the realm of agricultural sustainability, the utilization of plant genetic resources (PGRs) for enhanced disease resistance is paramount. Preservation efforts in genebanks are justified by their potential contributions to future crop improvement. To capitalize on the potential of PGRs, we focused on a barley core collection from the German ex situ genebank, and contrasted it with a European elite collection. The phenotypic assessment included 812 PGRs and 298 elites with a particular emphasis on four disease traits (Puccinia hordei, Blumeria graminis hordei, Ramularia collo-cygni, and Rhynchosporium commune). An integrated genome-wide association study, employing both Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK) and a linear mixed model, was performed to unravel the genetic underpinnings of disease resistance. A total of 932 marker-trait associations were identified and assigned to 49 quantitative trait loci. The accumulation of novel and rare resistance alleles significantly bolstered the overall resistance level in PGRs. Three PGR donors with high counts of novel/rare alleles and exhibited exceptional resistance to leaf rust and powdery mildew were identified, offering promise for targeted pre-breeding goals and enhanced resilience in forthcoming varieties. Our findings underscore the critical contribution of PGRs to strengthening crop resilience and advancing sustainable agricultural practices.

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.

ENHANCED GRAVITROPISM 2 encodes a STERILE ALPHA MOTIF-containing protein that controls root growth angle in barley and wheat.

The root growth angle defines how roots grow toward the gravity vector and is among the most important determinants of root system architecture. It controls water uptake capacity, nutrient use efficiency, stress resilience, and, as a consequence, yield of crop plants. We demonstrated that the egt2 (enhanced gravitropism 2) mutant of barley exhibits steeper root growth of seminal and lateral roots and an auxin-independent higher responsiveness to gravity compared to wild-type plants. We cloned the EGT2 gene by a combination of bulked-segregant analysis and whole genome sequencing. Subsequent validation experiments by an independent CRISPR/Cas9 mutant allele demonstrated that egt2 encodes a STERILE ALPHA MOTIF domain-containing protein. In situ hybridization experiments illustrated that EGT2 is expressed from the root cap to the elongation zone. We demonstrated the evolutionary conserved role of EGT2 in root growth angle control between barley and wheat by knocking out the EGT2 orthologs in the A and B genomes of tetraploid durum wheat. By combining laser capture microdissection with RNA sequencing, we observed that seven expansin genes were transcriptionally down-regulated in the elongation zone. This is consistent with a role of EGT2 in this region of the root where the effect of gravity sensing is executed by differential cell elongation. Our findings suggest that EGT2 is an evolutionary conserved regulator of root growth angle in barley and wheat that could be a valuable target for root-based crop improvement strategies in cereals.

Six-rowed wild-growing barleys are hybrids of diverse origins.

Crop-wild gene flow is common when domesticated plants and their wild relatives grow close to each other. The resultant hybrid forms appear as semi-domesticates and were sometimes considered as missing links between crops and their wild progenitors. Wild-growing barleys in Central and Eastern Asia, named Hordeum agriocrithon, show hallmark characters of both wild and domesticated forms. Their spikes disintegrate at maturity to disperse without human intervention, but bear lateral grains, which were favored by early farmers and are absent from other wild barleys. As an intermediate form, H. agriocrithon has been proposed several times as a progenitor of domesticated barley. Here, we used genome-wide marker data and whole-genome resequencing to show that all H. agriocrithon accessions of a major germplasm collection are hybrid forms that arose multiple times by admixture of diverse domesticated and wild populations. Although H. agriocrithon barleys have not played a special role in barley domestication, future analysis of the adaptative potential of bi-directional crop-wild gene flow in extant barleys may prove a fertile research field.

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.

Barley Ror1 encodes a class XI myosin required for mlo-based broad-spectrum resistance to the fungal powdery mildew pathogen.

Loss-of-function alleles of plant MLO genes confer broad-spectrum resistance to powdery mildews in many eudicot and monocot species. Although barley (Hordeum vulgare) mlo mutants have been used in agriculture for more than 40 years, understanding of the molecular principles underlying this type of disease resistance remains fragmentary. Forward genetic screens in barley have revealed mutations in two Required for mlo resistance (Ror) genes that partially impair immunity conferred by mlo mutants. While Ror2 encodes a soluble N-ethylmaleimide-sensitive factor-attached protein receptor (SNARE), the identity of Ror1, located at the pericentromeric region of barley chromosome 1H, remained elusive. We report the identification of Ror1 based on combined barley genomic sequence information and transcriptomic data from ror1 mutant plants. Ror1 encodes the barley class XI myosin Myo11A (HORVU.MOREX.r3.1HG0046420). Single amino acid substitutions of this myosin, deduced from non-functional ror1 mutant alleles, map to the nucleotide-binding region and the interface between the relay-helix and the converter domain of the motor protein. Ror1 myosin accumulates transiently in the course of powdery mildew infection. Functional fluorophore-labeled Ror1 variants associate with mobile intracellular compartments that partially colocalize with peroxisomes. Single-cell expression of the Ror1 tail region causes a dominant-negative effect that phenocopies ror1 loss-of-function mutants. We define a myosin motor for the establishment of mlo-mediated resistance, suggesting that motor protein-driven intracellular transport processes are critical for extracellular immunity, possibly through the targeted transfer of antifungal and/or cell wall cargoes to pathogen contact sites.

Recombination landscape divergence between populations is marked by larger low-recombining regions in domesticated rye.

The genomic landscape of recombination plays an essential role in evolution. Patterns of recombination are highly variable along chromosomes, between sexes, individuals, populations, and species. In many eukaryotes, recombination rates are elevated in sub-telomeric regions and drastically reduced near centromeres, resulting in large low-recombining (LR) regions. The processes of recombination are influenced by genetic factors, such as different alleles of genes involved in meiosis and chromatin structure, as well as external environmental stimuli like temperature and overall stress. In this work, we focused on the genomic landscapes of recombination in a collection of 916 rye (Secale cereale) individuals. By analysing population structure among individuals of different domestication status and geographic origin, we detected high levels of admixture, reflecting the reproductive biology of a self-incompatible, wind-pollinating grass species. We then analysed patterns of recombination in overlapping subpopulations, which revealed substantial variation in the physical size of LR regions, with a tendency for larger LR regions in domesticated subpopulations. Genome-wide association scans (GWAS) for LR region size revealed a major quantitative-trait-locus (QTL) at which, among 18 annotated genes, an ortholog of histone H4 acetyltransferase ESA1 was located. Rye individuals belonging to domesticated subpopulations showed increased synaptonemal complex length, but no difference in crossover frequency, indicating that only the recombination landscape is different. Furthermore, the genomic region harbouring rye ScESA1 showed moderate patterns of selection in domesticated subpopulations, suggesting that larger LR regions were indirectly selected for during domestication to achieve more homogeneous populations for agricultural use.

Genetic mapping of the powdery mildew resistance gene Pm7 on oat chromosome 5D.

KEY MESSAGE: Three independent experiments with different genetic backgrounds mapped the resistance gene Pm7 in the oat genome to the distal part of the long arm of chromosome 5D. Resistance of oat to Blumeria graminis DC. f. sp. avenae is an important breeding goal in Central and Western Europe. In this study, the position of the effective and widely used resistance gene Pm7 in the oat genome was determined based on three independent experiments with different genetic backgrounds: genome-wide association mapping in a diverse set of inbred oat lines and binary phenotype mapping in two bi-parental populations. Powdery mildew resistance was assessed in the field as well as by detached leaf tests in the laboratory. Genotyping-by-sequencing was conducted to establish comprehensive genetic fingerprints for subsequent genetic mapping experiments. All three mapping approaches located the gene to the distal part of the long arm of chromosome 5D in the hexaploid oat genome sequences of OT3098 and 'Sang.' Markers from this region were homologous to a region of chromosome 2Ce of the C-genome species, Avena eriantha, the donor of Pm7, which appears to be the ancestral source of a translocated region on the hexaploid chromosome 5D.

A chromosome conformation capture ordered sequence of the barley genome.

Cereal grasses of the Triticeae tribe have been the major food source in temperate regions since the dawn of agriculture. Their large genomes are characterized by a high content of repetitive elements and large pericentromeric regions that are virtually devoid of meiotic recombination. Here we present a high-quality reference genome assembly for barley (Hordeum vulgare L.). We use chromosome conformation capture mapping to derive the linear order of sequences across the pericentromeric space and to investigate the spatial organization of chromatin in the nucleus at megabase resolution. The composition of genes and repetitive elements differs between distal and proximal regions. Gene family analyses reveal lineage-specific duplications of genes involved in the transport of nutrients to developing seeds and the mobilization of carbohydrates in grains. We demonstrate the importance of the barley reference sequence for breeding by inspecting the genomic partitioning of sequence variation in modern elite germplasm, highlighting regions vulnerable to genetic erosion.