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

High-Resolution Mapping of Barley mild mosaic virus Resistance Gene rym15

Barley yellow mosaic virus (BaYMV) and Barley mild mosaic virus (BaMMV), which are transmitted by the soil-borne plasmodiophorid Polymyxa graminis, cause high yield losses in barley. In previous studies, the recessive BaMMV resistance gene rym15, derived from the Japanese landrace Chikurin Ibaraki 1, was mapped on chromosome 6HS of Hordeum vulgare. In this study, 423 F₄ segmental recombinant inbred lines (RILs) were developed from crosses of Chikurin Ibaraki 1 with two BaMMV-susceptible cultivars, Igri (139 RILs) and Uschi (284 RILs). A set of 32 competitive allele-specific PCR (KASP) assays, designed using single nucleotide polymorphisms (SNPs) from the barley 50 K Illumina Infinium iSelect SNP chip, genotyping by sequencing (GBS) and whole-genome sequencing (WGS), was used as a backbone for construction of two high-resolution maps. Using this approach, the target locus was narrowed down to 0.161 cM and 0.036 cM in the Igri × Chikurin Ibaraki 1 (I × C) and Chikurin Ibaraki 1 × Uschi (C × U) populations, respectively. Corresponding physical intervals of 11.3 Mbp and 0.281 Mbp were calculated for I × C and C × U, respectively, according to the Morex v3 genome sequence. In the 0.281 Mbp target region, six high confidence (HC) and two low confidence (LC) genes were identified. Genome assemblies of BaMMV-susceptible cultivars Igri and Golden Promise from the barley pan-genome, and a HiFi assembly of Chikurin Ibaraki 1 together with re-sequencing data for the six HC and two LC genes in susceptible parental cultivar Uschi revealed functional SNPs between resistant and susceptible genotypes only in two of the HC genes. These SNPs are the most promising candidates for the development of functional markers and the two genes represent promising candidates for functional analysis.

Genomic insights into positive selection during barley domestication

BACKGROUND

Cultivated barley (Hordeum vulgare) is widely used in animal feed, beverages, and foods and has become a model crop for molecular evolutionary studies. Few studies have examined the evolutionary fates of different types of genes in barley during the domestication process.

RESULTS

The rates of nonsynonymous substitution (Ka) to synonymous substitution (Ks) were calculated by comparing orthologous genes in different barley groups (wild vs. landrace and landrace vs. improved cultivar). The rates of evolution, properties, expression patterns, and diversity of positively selected genes (PSGs) and negatively selected genes (NSGs) were compared. PSGs evolved more rapidly, possessed fewer exons, and had lower GC content than NSGs; they were also shorter and had shorter intron, exon, and first exon lengths. Expression levels were lower, the tissue specificity of expression was higher, and codon usage bias was weaker for PSGs than for NSGs. Nucleotide diversity analysis revealed that PSGs have undergone a more severe genetic bottleneck than NSGs. Several candidate PSGs were involved in plant growth and development, which might make them as excellent targets for the molecular breeding of barley.

CONCLUSIONS

Our comprehensive analysis of the evolutionary, structural, and functional divergence between PSGs and NSGs in barley provides new insight into the evolutionary trajectory of barley during domestication. Our findings also aid future functional studies of PSGs in barley.

An interolog-based barley interactome as an integration framework for immune signaling

The barley MLA nucleotide-binding leucine-rich-repeat (NLR) receptor and its orthologs confer recognition specificity to many fungal diseases, including powdery mildew, stem-, and stripe rust. We used interolog inference to construct a barley protein interactome (Hordeum vulgare predicted interactome, HvInt) comprising 66,133 edges and 7,181 nodes, as a foundation to explore signaling networks associated with MLA. HvInt was compared with the experimentally validated Arabidopsis interactome of 11,253 proteins and 73,960 interactions, verifying that the 2 networks share scale-free properties, including a power-law distribution and small-world network. Then, by successive layering of defense-specific "omics" datasets, HvInt was customized to model cellular response to powdery mildew infection. Integration of HvInt with expression quantitative trait loci (eQTL) enabled us to infer disease modules and responses associated with fungal penetration and haustorial development. Next, using HvInt and infection-time-course RNA sequencing of immune signaling mutants, we assembled resistant and susceptible subnetworks. The resulting differentially coexpressed (resistant - susceptible) interactome is essential to barley immunity, facilitates the flow of signaling pathways and is linked to mildew resistance locus a (Mla) through trans eQTL associations. Lastly, we anchored HvInt with new and previously identified interactors of the MLA coiled coli + nucleotide-binding domains and extended these to additional MLA alleles, orthologs, and NLR outgroups to predict receptor localization and conservation of signaling response. These results link genomic, transcriptomic, and physical interactions during MLA-specified immunity.

GrainGenes: a data-rich repository for small grains genetics and genomics

As one of the US Department of Agriculture-Agricultural Research Service flagship databases, GrainGenes (https://wheat.pw.usda.gov) serves the data and community needs of globally distributed small grains researchers for the genetic improvement of the Triticeae family and Avena species that include wheat, barley, rye and oat. GrainGenes accomplishes its mission by continually enriching its cross-linked data content following the findable, accessible, interoperable and reusable principles, enhancing and maintaining an intuitive web interface, creating tools to enable easy data access and establishing data connections within and between GrainGenes and other biological databases to facilitate knowledge discovery. GrainGenes operates within the biological database community, collaborates with curators and genome sequencing groups and contributes to the AgBioData Consortium and the International Wheat Initiative through the Wheat Information System (WheatIS). Interactive and linked content is paramount for successful biological databases and GrainGenes now has 2917 manually curated gene records, including 289 genes and 254 alleles from the Wheat Gene Catalogue (WGC). There are >4.8 million gene models in 51 genome browser assemblies, 6273 quantitative trait loci and >1.4 million genetic loci on 4756 genetic and physical maps contained within 443 mapping sets, complete with standardized metadata. Most notably, 50 new genome browsers that include outputs from the Wheat and Barley PanGenome projects have been created. We provide an example of an expression quantitative trait loci track on the International Wheat Genome Sequencing Consortium Chinese Spring wheat browser to demonstrate how genome browser tracks can be adapted for different data types. To help users benefit more from its data, GrainGenes created four tutorials available on YouTube. GrainGenes is executing its vision of service by continuously responding to the needs of the global small grains community by creating a centralized, long-term, interconnected data repository. Database URL:https://wheat.pw.usda.gov.

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 Anatolia^1,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.

Assembly of the hexaploid oat genome and its diploid and tetraploid relatives clarifies the evolutionary history of oat and allows mapping of genes for agronomic traits.

The barley leaf rust resistance gene Rph3 encodes a predicted membrane protein and is induced upon infection by avirulent pathotypes of Puccinia hordei

Leaf rust, caused by Puccinia hordei, is an economically significant disease of barley, but only a few major resistance genes to P. hordei (Rph) have been cloned. In this study, gene Rph3 was isolated by positional cloning and confirmed by mutational analysis and transgenic complementation. The Rph3 gene, which originated from wild barley and was first introgressed into cultivated Egyptian germplasm, encodes a unique predicted transmembrane resistance protein that differs from all known plant disease resistance proteins at the amino acid sequence level. Genetic profiles of diverse accessions indicated limited genetic diversity in Rph3 in domesticated germplasm, and higher diversity in wild barley from the Eastern Mediterranean region. The Rph3 gene was expressed only in interactions with Rph3-avirulent P. hordei isolates, a phenomenon also observed for transcription activator-like effector-dependent genes known as executors conferring resistance to Xanthomonas spp. Like known transmembrane executors such as Bs3 and Xa7, heterologous expression of Rph3 in N. benthamiana induced a cell death response. The isolation of Rph3 highlights convergent evolutionary processes in diverse plant-pathogen interaction systems, where similar defence mechanisms evolved independently in monocots and dicots.

Leaf rust is an economically significant disease of barley. Here the authors describe cloning of the barley Rph3 leaf rust resistance gene and reveal it encodes a predicted transmembrane protein that is expressed upon infection by Rph3-avirulent Puccinia hordei isolates.

The value of genotype-specific reference for transcriptome analyses in barley

We demonstrate in this study that using a common reference genome may lead to loss of genotype-specific information in the assembled Reference Transcript Dataset (RTD) and the generation of erroneous, incomplete, or misleading transcriptomics analysis results in barley.

It is increasingly apparent that although different genotypes within a species share “core” genes, they also contain variable numbers of “specific” genes and different structures of “core” genes that are only present in a subset of individuals. Using a common reference genome may thus lead to a loss of genotype-specific information in the assembled Reference Transcript Dataset (RTD) and the generation of erroneous, incomplete or misleading transcriptomics analysis results. In this study, we assembled genotype-specific RTD (sRTD) and common reference–based RTD (cRTD) from RNA-seq data of cultivated Barke and Morex barley, respectively. Our quantitative evaluation showed that the sRTD has a significantly higher diversity of transcripts and alternative splicing events, whereas the cRTD missed 40% of transcripts present in the sRTD and it only has ∼70% accurate transcript assemblies. We found that the sRTD is more accurate for transcript quantification as well as differential expression analysis. However, gene-level quantification is less affected, which may be a reasonable compromise when a high-quality genotype-specific reference is not available.

Triticeae genome sequences reveal huge expansions of gene families implicated in fertility restoration

Breakthroughs in assembly of whole-genome sequencing and targeted sequence capture data have accelerated comparative genomics analyses in cereals with big and complex genomes such as wheat. This newly acquired information has revealed unexpected expansions in two large gene families linked to restoration of fertility in species that exhibit cytoplasmic male sterility. Extreme levels of copy-number and structural variation detected within and between species illustrate the genetic diversity among the family members and reveal the evolutionary mechanisms at work. This new knowledge will greatly facilitate the development of hybrid production strategies in wheat and related species.

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

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

Genomic prediction of the recombination rate variation in barley - A route to highly recombinogenic genotypes

Meiotic recombination is not only fundamental to the adaptation of sexually reproducing eukaryotes in nature but increased recombination rates facilitate the combination of favourable alleles into a single haplotype in breeding programmes. The main objectives of this study were to (i) assess the extent and distribution of the recombination rate variation in cultivated barley (Hordeum vulgare L.), (ii) quantify the importance of the general and specific recombination effects, and (iii) evaluate a genomic selection approach's ability to predict the recombination rate variation. Genetic maps were created for the 45 segregating populations that were derived from crosses among 23 spring barley inbreds with origins across the world. The genome-wide recombination rate among populations ranged from 0.31 to 0.73 cM/Mbp. The crossing design used in this study allowed to separate the general recombination effects (GRE) of individual parental inbreds from the specific recombination effects (SRE) caused by the combinations of parental inbreds. The variance of the genome-wide GRE was found to be about eight times the variance of the SRE. This finding indicated that parental inbreds differ in the efficiency of their recombination machinery. The ability to predict the chromosome or genome-wide recombination rate of an inbred ranged from 0.80 to 0.85. These results suggest that a reliable screening of large genetic materials for their potential to cause a high extent of genetic recombination in their progeny is possible, allowing to systematically manipulate the recombination rate using natural variation.

Aegilops sharonensis genome-assisted identification of stem rust resistance gene Sr62

The wild relatives and progenitors of wheat have been widely used as sources of disease resistance (R) genes. Molecular identification and characterization of these R genes facilitates their manipulation and tracking in breeding programmes. Here, we develop a reference-quality genome assembly of the wild diploid wheat relative Aegilops sharonensis and use positional mapping, mutagenesis, RNA-Seq and transgenesis to identify the stem rust resistance gene Sr62, which has also been transferred to common wheat. This gene encodes a tandem kinase, homologues of which exist across multiple taxa in the plant kingdom. Stable Sr62 transgenic wheat lines show high levels of resistance against diverse isolates of the stem rust pathogen, highlighting the utility of Sr62 for deployment as part of a polygenic stack to maximize the durability of stem rust resistance.

Aegilops sharonensis is a wild diploid relative of wheat. Here, the authors assemble the genome of Ae. sharonensis and use the assembly as an aid to clone the Ae. sharonensis-derived stem rust resistance gene Sr62 in the allohexaploid genome of wheat.

Screening of stable resistant accessions and identification of resistance loci to Barley yellow mosaic virus disease

Background

The disease caused by Barley yellow mosaic virus (BaYMV) infection is a serious threat to autumn-sown barley (Hordeum vulgare L.) production in Europe, East Asia and Iran. Due to the rapid diversification of BaYMV strains, it is urgent to discover novel germplasm and genes to assist breeding new varieties with resistance to different BaYMV strains, thus minimizing the effect of BaYMV disease on barley cropping.

Methods

A natural population consisting of 181 barley accessions with different levels of resistance to BaYMV disease was selected for field resistance identification in two separate locations (Yangzhou and Yancheng, Jiangsu Province, China). Additive main effects and multiplicative interaction (AMMI) analysis was used to identify accessions with stable resistance. Genome-wide association study (GWAS) of BaYMV disease resistance was broadly performed by combining both single nucleotide polymorphisms (SNPs) and specific molecular markers associated with the reported BaYMV disease resistance genes. Furthermore, the viral protein genome linked (VPg) sequences of the virus were amplified and analyzed to assess the differences between the BaYMV strains sourced from the different experimental sites.

Results

Seven barley accessions with lower standardized Area Under the Disease Progress Steps (sAUDPS) index in every environment were identified and shown to have stable resistance to BaYMV disease in each assessed location. Apart from the reported BaYMV disease resistance genes rym4 and rym5, one novel resistance locus explaining 24.21% of the phenotypic variation was identified at the Yangzhou testing site, while two other novel resistance loci that contributed 19.23% and 19.79% of the phenotypic variation were identified at the Yancheng testing site, respectively. Further analysis regarding the difference in the VPg sequence of the predominant strain of BaYMV collected from these two testing sites may explain the difference of resistance loci differentially identified under geographically distinct regions. Our research provides novel genetic resources and resistance loci for breeding barley varieties for BaMYV disease resistance.

Auxin boosts energy generation pathways to fuel pollen maturation in barley

Pollen grains become increasingly independent of the mother plant as they reach maturity through poorly understood developmental programs. We report that the hormone auxin is essential during barley pollen maturation to boost the expression of genes encoding almost every step of heterotrophic energy production pathways. Accordingly, auxin is necessary for the flux of sucrose and hexoses into glycolysis and to increase the levels of pyruvate and two tricarboxylic (TCA) cycle metabolites (citrate and succinate). Moreover, bioactive auxin is synthesized by the pollen-localized enzyme HvYUCCA4, supporting that pollen grains autonomously produce auxin to stimulate a specific cellular output, energy generation, that fuels maturation processes such as starch accumulation. Our results demonstrate that auxin can shift central carbon metabolism to drive plant cell development, which suggests a direct mechanism for auxin's ability to promote growth and differentiation.

Genetic analysis toward more nutritious barley grains for a food secure world

BACKGROUND

Understanding the relationships between nutrition, human health and plant food source is among the highest priorities for public health. Therefore, enhancing the minerals content such as iron (Fe), zinc (Zn) and selenium (Se) in barley (Hordeum vulgare L.) grains is an urgent need to improve the nutritive value of barley grains in overcoming malnutrition and its potential consequencing. This study aimed to expedite biofortification of barley grains by elucidating the genetic basis of Zn, Fe, and Se accumulation in the grains, which will contribute to improved barley nutritional quality.

RESULTS

A genome-wide association study (GWAS) was conducted to detect the genetic architecture for grain Zn, Fe, and Se accumulations in 216 spring barley accessions across two years. All the accessions were genotyped by single nucleotide polymorphisms (SNPs) molecular markers. Mineral heritability values ranging from moderate to high were revealed in both environments. Remarkably, there was a high natural phenotypic variation for all micronutrient accumulation in the used population. High-LD SNP markers (222 SNPs) were detected to be associated with all micronutrients in barley grains across the two environments plus BLUEs. Three genomic regions were detected based on LD, which were identified for the most effective markers that had associations with more than one trait. The strongest SNP-trait associations were found to be physically located within genes that may be involved in grain Zn and Fe homeostasis. Two putative candidate genes were annotated as Basic helix loop helix (BHLH) family transcription factor and Squamosa promoter binding-like protein, respectively, and have been suggested as candidates for increased grain Zn, Fe, and Se accumulation.

CONCLUSIONS

These findings shed a light on the genetic basis of Zn, Fe, and Se accumulation in barley grains and have the potential to assist plant breeders in selecting accessions with high micronutrient concentrations to enhance grain quality and, ultimately human health.

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.

Immunity-associated volatile emissions of β-ionone and nonanal propagate defence responses in neighbouring barley plants

Plants activate biochemical responses to combat stress. (Hemi-)biotrophic pathogens are fended off by systemic acquired resistance (SAR), a primed state allowing plants to respond faster and more strongly upon subsequent infection. Here, we show that SAR-like defences in barley (Hordeum vulgare) are propagated between neighbouring plants, which respond with enhanced resistance to the volatile cues from infected senders. The emissions of the sender plants contained 15 volatile organic compounds (VOCs) associated with infection. Two of these, β-ionone and nonanal, elicited resistance upon plant exposure. Whole-genome transcriptomics analysis confirmed that interplant propagation of defence in barley is established as a form of priming. Although gene expression changes were more pronounced after challenge infection of the receiver plants with Blumeria graminis f. sp. hordei, differential gene expression in response to the volatile cues of the sender plants included an induction of HISTONE DEACETYLASE 2 (HvHDA2) and priming of TETRATRICOPEPTIDE REPEAT-LIKE superfamily protein (HvTPL). Because HvHDA2 and HvTPL transcript accumulation was also enhanced by exposure of barley to β-ionone and nonanal, our data identify both genes as possible defence/priming markers in barley. Our results suggest that VOCs and plant-plant interactions are relevant for possible crop protection strategies priming defence responses in barley.

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

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

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.

Calmodulin and calmodulin-like gene family in barley: Identification, characterization and expression analyses

Calmodulin (CaM) and calmodulin-like (CML) proteins are Ca²⁺ relays and play diverse and multiple roles in plant growth, development and stress responses. However, CaM/CML gene family has not been identified in barley (Hordeum vulgare). In the present study, 5 HvCaMs and 80 HvCMLs were identified through a genome-wide analysis. All HvCaM proteins possessed 4 EF-hand motifs, whereas HvCMLs contained 1 to 4 EF-hand motifs. HvCaM2, HvCaM3 and HvCaM5 coded the same polypeptide although they differed in nucleotide sequence, which was identical to the polypeptides coded by OsCaM1-1, OsCaM1-2 and OsCaM1-3. HvCaMs/CMLs were unevenly distributed over barley 7 chromosomes, and could be phylogenetically classified into 8 groups. HvCaMs/CMLs differed in gene structure, cis-acting elements and tissue expression patterns. Segmental and tandem duplication were observed among HvCaMs/CMLs during evolution. HvCML16, HvCML18, HvCML50 and HvCML78 were dispensable genes and the others were core genes in barley pan-genome. In addition, 14 HvCaM/CML genes were selected to examine their responses to salt, osmotic and low potassium stresses by qRT-PCR, and their expression were stress-and time-dependent. These results facilitate our understanding and further functional identification of HvCaMs/CMLs.

Novel Design of Imputation-Enabled SNP Arrays for Breeding and Research Applications Supporting Multi-Species Hybridization

Array-based single nucleotide polymorphism (SNP) genotyping platforms have low genotype error and missing data rates compared to genotyping-by-sequencing technologies. However, design decisions used to create array-based SNP genotyping assays for both research and breeding applications are critical to their success. We describe a novel approach applicable to any animal or plant species for the design of cost-effective imputation-enabled SNP genotyping arrays with broad utility and demonstrate its application through the development of the Illumina Infinium Wheat Barley 40K SNP array Version 1.0. We show that the approach delivers high quality and high resolution data for wheat and barley, including when samples are jointly hybridised. The new array aims to maximally capture haplotypic diversity in globally diverse wheat and barley germplasm while minimizing ascertainment bias. Comprising mostly biallelic markers that were designed to be species-specific and single-copy, the array permits highly accurate imputation in diverse germplasm to improve the statistical power of genome-wide association studies (GWAS) and genomic selection. The SNP content captures tetraploid wheat (A- and B-genome) and Aegilops tauschii Coss. (D-genome) diversity and delineates synthetic and tetraploid wheat from other wheat, as well as tetraploid species and subgroups. The content includes SNP tagging key trait loci in wheat and barley, as well as direct connections to other genotyping platforms and legacy datasets. The utility of the array is enhanced through the web-based tool, Pretzel (https://plantinformatics.io/) which enables the content of the array to be visualized and interrogated interactively in the context of numerous genetic and genomic resources to be connected more seamlessly to research and breeding. The array is available for use by the international wheat and barley community.

Mutations in a barley cytochrome P450 gene enhances pathogen induced programmed cell death and cutin layer instability

Disease lesion mimic mutants (DLMMs) are characterized by the spontaneous development of necrotic spots with various phenotypes designated as necrotic (nec) mutants in barley. The nec mutants were traditionally considered to have aberrant regulation of programmed cell death (PCD) pathways, which have roles in plant immunity and development. Most barley nec3 mutants express cream to orange necrotic lesions contrasting them from typical spontaneous DLMMs that develop dark pigmented lesions indicative of serotonin/phenolics deposition. Barley nec3 mutants grown under sterile conditions did not exhibit necrotic phenotypes until inoculated with adapted pathogens, suggesting that they are not typical DLMMs. The F₂ progeny of a cross between nec3-γ1 and variety Quest segregated as a single recessive susceptibility gene post-inoculation with Bipolaris sorokiniana, the causal agent of the disease spot blotch. Nec3 was genetically delimited to 0.14 cM representing 16.5 megabases of physical sequence containing 149 annotated high confidence genes. RNAseq and comparative analysis of the wild type and five independent nec3 mutants identified a single candidate cytochrome P450 gene (HORVU.MOREX.r2.6HG0460850) that was validated as nec3 by independent mutations that result in predicted nonfunctional proteins. Histology studies determined that nec3 mutants had an unstable cutin layer that disrupted normal Bipolaris sorokiniana germ tube development.

At the site of pathogen infection, plant defense mechanisms rely on controlled programmed cell death (PCD) to sequester biotrophic pathogens that require living cells to extract nutrients from the host. However, these defense mechanisms are hijacked by necrotrophic plant pathogens that purposefully induce PCD to feed on the dead cells, thus facilitating further disease development. Thus, understanding PCD responses is important for resistance to both classes of pathogens. We characterized five independent disease lesion mimic mutants of barley designated necrotic 3 (nec3) that show aberrant regulation of PCD responses upon pathogen challenge. A cytochrome P450 gene was identified as Nec3 encoding a Tryptamine 5-Hydroxylase that functions as a terminal serotonin biosynthetic enzyme in the Tryptophan pathway of plants. We posit that nec3 mutants have disrupted serotonin biosynthesis resulting in expanded PCD, necrotrophic pathogen susceptibility and cutin layer instability. The nec3 mutants show expanded PCD and disease susceptibility of pathogen-induced necrotic lesions, suggesting a role of serotonin to sequester PCD and suppress pathogen colonization. The identification of Nec3 will facilitate functional analysis to elucidate the role that serotonin plays in the elicitation or suppression of PCD immunity responses to diverse pathogens and the effects it has on cutin layer biosynthesis.

Delineating the elusive BaMMV resistance gene rym15 in barley by medium-resolution mapping

Barley mild mosaic virus (BaMMV), transmitted by the soil-borne protist Polymyxa graminis, has a serious impact on winter barley production. Previously, the BaMMV resistance gene rym15 was mapped on chromosome 6HS, but the order of flanking markers was non-collinear between different maps. To resolve the position of the flanking markers and to enable map-based cloning of rym15, two medium-resolution mapping populations Igri (susceptible) × Chikurin Ibaraki 1 (resistant) (I × C) and Chikurin Ibaraki 1 × Uschi (susceptible) (C × U), consisting of 342 and 180 F2 plants, respectively, were developed. Efficiency of the mechanical inoculation of susceptible standards varied from 87.5 to 100% and in F2 populations from 90.56 to 93.23%. Phenotyping of F2 plants and corresponding F3 families revealed segregation ratios of 250 s:92r (I × C, χ2 = 0.659) and 140 s:40r (C × U, χ2 = 0.741), suggesting the presence of a single recessive resistance gene. After screening the parents with the 50 K Infinium chip and anchoring corresponding SNPs to the barley reference genome, 8 KASP assays were developed and used to remap the gene. Newly constructed maps revealed a collinear order of markers, thereby allowing the identification of high throughput flanking markers. This study demonstrates how construction of medium-resolution mapping populations in combination with robust phenotyping can efficiently resolve conflicting marker ordering and reduce the size of the target interval. In the reference genome era and genome-wide genotyping era, medium-resolution mapping will help accelerate candidate gene identification for traits where phenotyping is difficult.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01270-9.

Gain from genomic selection for a selection index in two-row spring barley

New breeding programs are faced with many challenges including evaluation of unknown germplasm, initiation of breeding populations that will satisfy short- and long-term breeding goals, and implementation of efficient phenotyping strategies for multiple traits. Genomic selection (GS) is a potentially valuable tool for recently established breeding programs to quickly accelerate genetic gain. Genomic selection on selection index (SI) values may increase gain over phenotypic selection but empirical studies remain limited. We compared gain in overall SI value for height, heading date, preharvest sprouting (PHS) resistance, and spot blotch resistance and component traits in two cycles of GS with one round of phenotypic selection (PS) in two-row spring malting barley (Hordeum vulgare L.). Higher realized gain for SI value, height, and PHS was observed with GS compared with PS but GS did not result in significant gain for heading date and spot blotch. Genetic variances for height and heading date, which had small index weights, were not reduced with GS but variances were substantially reduced for heavily weighted PHS and correlated seed germination traits. Inbreeding was increased by GS compared with PS but restricted mating of high breeding value individuals limited potential inbreeding. Our results indicate GS is a useful method to improve selection on index values with different weights.

Proteins, Small Peptides and Other Signaling Molecules Identified as Inconspicuous but Possibly Important Players in Microspores Reprogramming Toward Embryogenesis

In this review, we describe and integrate the latest knowledge on the signaling role of proteins and peptides in the stress-induced microspore embryogenesis (ME) in some crop plants with agricultural importance (i.e., oilseed rape, tobacco, barley, wheat, rice, triticale, rye). Based on the results received from the most advanced omix analyses, we have selected some inconspicuous but possibly important players in microspores reprogramming toward embryogenic development. We provide an overview of the roles and downstream effect of stress-related proteins (e.g., β-1,3-glucanases, chitinases) and small signaling peptides, especially cysteine—(e.g., glutathione, γ-thionins, rapid alkalinization factor, lipid transfer, phytosulfokine) and glycine-rich peptides and other proteins (e.g., fasciclin-like arabinogalactan protein) on acclimation ability of microspores and the cell wall reconstruction in a context of ME induction and haploids/doubled haploids (DHs) production. Application of these molecules, stimulating the induction and proper development of embryo-like structures and green plant regeneration, brings significant improvement of the effectiveness of DHs procedures and could result in its wider incorporation on a commercial scale. Recent advances in the design and construction of synthetic peptides–mainly cysteine-rich peptides and their derivatives–have accelerated the development of new DNA-free genome-editing techniques. These new systems are evolving incredibly fast and soon will find application in many areas of plant science and breeding.

Genome-wide association mapping of Pyrenophora teres f. maculata and Pyrenophora teres f. teres resistance loci utilizing natural Turkish wild and landrace barley populations

Unimproved landraces and wild relatives of crops are sources of genetic diversity that were lost post domestication in modern breeding programs. To tap into this rich resource, genome-wide association studies in large plant genomes have enabled the rapid genetic characterization of desired traits from natural landrace and wild populations. Wild barley (Hordeum spontaneum), the progenitor of domesticated barley (Hordeum vulgare), is dispersed across Asia and North Africa, and has co-evolved with the ascomycetous fungal pathogens Pyrenophora teres f. teres and P. teres f. maculata, the causal agents of the diseases net form of net blotch and spot form of net blotch, respectively. Thus, these wild and local adapted barley landraces from the region of origin of both the host and pathogen represent a diverse gene pool to identify new sources of resistance, due to millions of years of co-evolution. The barley—P. teres pathosystem is governed by complex genetic interactions with dominant, recessive, and incomplete resistances and susceptibilities, with many isolate-specific interactions. Here, we provide the first genome-wide association study of wild and landrace barley from the Fertile Crescent for resistance to both forms of P. teres. A total of 14 loci, four against P. teres f. maculata and 10 against P. teres f. teres, were identified in both wild and landrace populations, showing that both are genetic reservoirs for novel sources of resistance. We also highlight the importance of using multiple algorithms to both identify and validate additional loci.

Genetic Loci Underlying Awn Morphology in Barley

Barley awns are highly active in photosynthesis and account for 30–50% of grain weight in barley. They are diverse in length, ranging from long to awnless, and in shape from straight to hooded or crooked. Their diversity and importance have intrigued geneticists for several decades. A large collection of awnness mutants are available—over a dozen of them have been mapped on chromosomes and a few recently cloned. Different awnness genes interact with each other to produce diverse awn phenotypes. With the availability of the sequenced barley genome and application of new mapping and gene cloning strategies, it will now be possible to identify and clone more awnness genes. A better understanding of the genetic basis of awn diversity will greatly facilitate development of new barley cultivars with improved yield, adaptability and sustainability.

Chorus2: design of genome-scale oligonucleotide-based probes for fluorescence in situ hybridization

Oligonucleotide (oligo)-fluorescence in situ hybridization (FISH) has rapidly becoming the new generation of FISH technique in plant molecular cytogenetics research. Genome-scale identification of single-copy oligos is the foundation of successful oligo-FISH experiments. Here, we introduce Chorus2, a software that is developed specifically for oligo selection. We demonstrate that Chorus2 is highly effective to remove all repetitive elements in selection of single-copy oligos, which is critical for the development of successful FISH probes. Chorus2 is more effective than Chorus, the original version of the pipeline, and OligoMiner for repeat removal. Chorus2 allows to select oligos that are conserved among related species, which extends the usage of oligo-FISH probes among phylogenetically related plant species. We also implemented a new function in Chorus2 that allows development of FISH probes from plant species without an assembled genome. We anticipate that Chorus2 can be used in plants as well as in mammalian and other non-plant species. Chorus2 will broadly facilitate the design of FISH probes for various types of application in molecular cytogenetics research.

Genetic and molecular characterization of determinant of six-rowed spike of barley carrying vrs1.a4

Decisive role of reduced vrs1 transcript abundance in six-rowed spike of barley carrying vrs1.a4 was genetically proved and its potential causes were preliminarily analyzed. Six-rowed spike 1 (vrs1) is the major determinant of the six-rowed spike phenotype of barley (Hordeum vulgare L.). Alleles of Vrs1 have been extensively investigated. Allele vrs1.a4 in six-rowed barley is unique in that it has the same coding sequence as Vrs1.b4 in two-rowed barley. The determinant of row-type in vrs1.a4 carriers has not been experimentally identified. Here, we identified Vrs1.b4 in two-rowed accessions and vrs1.a4 in six-rowed accessions from the Qinghai-Tibet Plateau at high frequency. Genetic analyses revealed a single nuclear gene accounting for row-type alteration in these accessions. Physical mapping identified a 0.08-cM (~ 554-kb) target interval on chromosome 2H, wherein Vrs1 was the most likely candidate gene. Further analysis of Vrs1 expression in offspring of the mapping populations or different Vrs1.b4 and vrs1.a4 lines confirmed that downregulated expression of vrs1.a4 causes six-rowed spike. Regulatory sequence analysis found a single 'TA' dinucleotide deletion in vrs1.a4 carriers within a 'TA' tandem-repeat-enriched region ~ 1 kb upstream of the coding region. DNA methylation levels did not correspond to the expression difference and therefore did not affect Vrs1 expression. More evidence is needed to verify the causal link between the 'TA' deletion and the downregulated Vrs1 expression and hence the six-rowed spike phenotype.

Downregulation of Barley Regulator of Telomere Elongation Helicase 1 Alters the Distribution of Meiotic Crossovers

Programmed meiotic DNA double-strand breaks (DSBs), necessary for proper chromosomal segregation and viable gamete formation, are repaired by homologous recombination (HR) as crossovers (COs) or non-crossovers (NCOs). The mechanisms regulating the number and distribution of COs are still poorly understood. The regulator of telomere elongation helicase 1 (RTEL1) DNA helicase was previously shown to enforce the number of meiotic COs in Caenorhabditis elegans but its function in plants has been studied only in the vegetative phase. Here, we characterised barley RTEL1 gene structure and expression using RNA-seq data previously obtained from vegetative and reproductive organs and tissues. Using RNAi, we downregulated RTEL1 expression specifically in reproductive tissues and analysed its impact on recombination using a barley 50k iSelect SNP Array. Unlike in C. elegans, in a population segregating for RTEL1 downregulated by RNAi, high resolution genome-wide genetic analysis revealed a significant increase of COs at distal chromosomal regions of barley without a change in their total number. Our data reveal the important role of RTEL1 helicase in plant meiosis and control of recombination.

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.

Transcript Profiling of MIKCc MADS-Box Genes Reveals Conserved and Novel Roles in Barley Inflorescence Development

MADS-box genes have a wide range of functions in plant reproductive development and grain production. The ABCDE model of floral organ development shows that MADS-box genes are central players in these events in dicotyledonous plants but the applicability of this model remains largely unknown in many grass crops. Here, we show that transcript analysis of all MIKCc MADS-box genes through barley (Hordeum vulgare L.) inflorescence development reveals co-expression groups that can be linked to developmental events. Thirty-four MIKCc MADS-box genes were identified in the barley genome and single-nucleotide polymorphism (SNP) scanning of 22,626 barley varieties revealed that the natural variation in the coding regions of these genes is low and the sequences have been extremely conserved during barley domestication. More detailed transcript analysis showed that MADS-box genes are generally expressed at key inflorescence developmental phases and across various floral organs in barley, as predicted by the ABCDE model. However, expression patterns of some MADS genes, for example HvMADS58 (AGAMOUS subfamily) and HvMADS34 (SEPALLATA subfamily), clearly deviate from predicted patterns. This places them outside the scope of the classical ABCDE model of floral development and demonstrates that the central tenet of antagonism between A- and C-class gene expression in the ABC model of other plants does not occur in barley. Co-expression across three correlation sets showed that specifically grouped members of the barley MIKCc MADS-box genes are likely to be involved in developmental events driving inflorescence meristem initiation, floral meristem identity and floral organ determination. Based on these observations, we propose a potential floral ABCDE working model in barley, where the classic model is generally upheld, but that also provides new insights into the role of MIKCc MADS-box genes in the developing barley inflorescence.

Barley Viridis-k links an evolutionarily conserved C-type ferredoxin to chlorophyll biosynthesis

Ferredoxins are single-electron carrier proteins involved in various cellular reactions. In chloroplasts, the most abundant ferredoxin accepts electrons from photosystem I and shuttles electrons via ferredoxin NADP+ oxidoreductase to generate NADPH or directly to ferredoxin dependent enzymes. In addition, plants contain other isoforms of ferredoxins. Two of these, named FdC1 and FdC2 in Arabidopsis thaliana, have C-terminal extensions and functions that are poorly understood. Here we identified disruption of the orthologous FdC2 gene in barley (Hordeum vulgare L.) mutants at the Viridis-k locus; these mutants are deficient in the aerobic cyclase reaction of chlorophyll biosynthesis. The magnesium-protoporphyrin IX monomethyl ester cyclase is one of the least characterized enzymes of the chlorophyll biosynthetic pathway and its electron donor has long been sought. Agroinfiltrations showed that the viridis-k phenotype could be complemented in vivo by Viridis-k but not by canonical ferredoxin. VirK could drive the cyclase reaction in vitro and analysis of cyclase mutants showed that in vivo accumulation of VirK is dependent on cyclase enzyme levels. The chlorophyll deficient phenotype of viridis-k mutants suggests that VirK plays an essential role in chlorophyll biosynthesis that cannot be replaced by other ferredoxins, thus assigning a specific function to this isoform of C-type ferredoxins.

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

To date, the potential of utilizing root traits in plant breeding remains largely untapped. In this study, we cloned and characterized the ENHANCED GRAVITROPISM2 (EGT2) gene of barley that encodes a STERILE ALPHA MOTIF domain–containing protein. We demonstrated that EGT2 is a key gene of root growth angle regulation in response to gravity, which is conserved in barley and wheat and could be a promising target for crop improvement in cereals.

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.

RNA-Seq-based DNA marker analysis of the genetics and molecular evolution of Triticeae species

The release of high-quality chromosome-level genome sequences of members of the Triticeae tribe has greatly facilitated genetic and genomic analyses of important crops such as wheat (Triticum aestivum) and barley (Hordeum vulgare). Due to the large diploid genome size of Triticeae plants (ca. 5 Gbp), transcript analysis is an important method for identifying genetic and genomic differences among Triticeae species. In this review, we summarize our results of RNA-Seq analyses of diploid wheat accessions belonging to the genera Aegilops and Triticum. We also describe studies of the molecular relationships among these accessions and provide insight into the evolution of common hexaploid wheat. DNA markers based on polymorphisms within species can be used to map loci of interest. Even though the genome sequence of diploid Aegilops tauschii, the D-genome donor of common wheat, has been released, the diploid barley genome continues to provide key information about the physical structures of diploid wheat genomes. We describe how a series of RNA-Seq analyses of wheat relatives has helped uncover the structural and evolutionary features of genomic and genetic systems in wild and cultivated Triticeae species.

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.

Rapid Generation and Analysis of a Barley Doubled Haploid Line with Higher Nitrogen Use Efficiency Than Parental Lines by F1 Microspore Embryogenesis

Creating varieties with high nitrogen use efficiency (NUE) is crucial for sustainable agriculture development. In this study, a superior barley doubled haploid line (named DH45) with improved NUE was produced via F₁ microspore embryogenesis with three rounds of screening in different nitrogen levels by hydroponic and field experiments. The molecular mechanisms responsible for the NUE of DH45 surpassing that of its parents were investigated by RNA-seq analysis. A total of 1027 differentially expressed genes (DEGs) were identified that were up- or down-regulated in DH45 under low nitrogen conditions but showed no significant differences in the parents. GO analysis indicated that genes involved in nitrogen compound metabolic processes were significantly enriched in DH45 compared with the parents. KEGG analysis showed the MAPK signaling pathway plant to be highly enriched in DH45 relative to its parents, as well as genes involved in alanine, aspartate and glutamate metabolism, and arginine biosynthesis. In conclusion, our study revealed the potential to fix trait superiority in a line by combining crossing with F₁ microspore culture technologies in future crop breeding and also identified several candidate genes that are expressed in shoots and may enable barley to cope with low-nitrogen stress.

Identification and candidate gene mining of HvSS1, a novel qualitative locus on chromosome 6H, regulating the uppermost internode elongation in barley (Hordeum vulgare L.)

A novel qualitative locus regulating the uppermost internode elongation of barley was identified and mapped on 6H, and the candidate gene mining was performed by employing various barley genomic resources. The stem of grass crops, such as barley and wheat, is composed of several interconnected internodes. The extent of elongation of these internodes determines stem height, and hence lodging, canopy architecture, and grain yield. The uppermost internode (UI) is the last internode to elongate. Its elongation contributes largely to stem height and facilitates spike exsertion, which is crucial for final grain yield. Despite the molecular mechanism underlying regulation of UI elongation was extensively investigated in rice, little is known in barley. In this study, we characterized a barley spontaneous mutant, Sheathed Spike 1 (SS1), showing significantly shortened UI and sheathed spike (SS). The extension of UI parenchyma cell in SS1 was significantly suppressed. Exogenous hormone treatments and RNA-seq analysis indicated that the suppression of UI elongation is possibly related to insufficient content of endogenous bioactive gibberellin. Genetic analysis showed that SS1 is possibly controlled by a qualitative dominant nuclear factor. Bulked segregant analysis and further molecular marker mapping identified a novel major locus, HvSS1, in a recombination cold spot expanding 173.44-396.33 Mb on chromosome 6H. The candidate gene mining was further conducted by analyzing sequence differences, spatiotemporal expression patterns, and variant distributions of genes in the candidate interval by employing various barley genomic resources of worldwide collections of barley accessions. This study made insight into genetic control of UI elongation in barley and laid a solid foundation for further gene cloning and functional characterization. The results obtained here also provided valuable information for similar research in wheat.

Genome editing in cereal crops: an overview

Genome-editing technologies offer unprecedented opportunities for crop improvement with superior precision and speed. This review presents an analysis of the current state of genome editing in the major cereal crops- rice, maize, wheat and barley. Genome editing has been used to achieve important agronomic and quality traits in cereals. These include adaptive traits to mitigate the effects of climate change, tolerance to biotic stresses, higher yields, more optimal plant architecture, improved grain quality and nutritional content, and safer products. Not all traits can be achieved through genome editing, and several technical and regulatory challenges need to be overcome for the technology to realize its full potential. Genome editing, however, has already revolutionized cereal crop improvement and is poised to shape future agricultural practices in conjunction with other breeding innovations.

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.

Genome Assembly of the Canadian two-row Malting Barley cultivar AAC Synergy

Barley (Hordeum vulgare L.) is one of the most important global crops. The six-row barley cultivar Morex reference genome has been used by the barley research community worldwide. However, this reference genome can have limitations when used for genomic and genetic diversity analysis studies, gene discovery, and marker development when working in two-row germplasm that is more common to Canadian barley. Here we assembled, for the first time, the genome sequence of a Canadian two-row malting barley, cultivar AAC Synergy. We applied deep Illumina paired-end reads, long mate-pair reads, PacBio sequences, 10X chromium linked read libraries, and chromosome conformation capture sequencing (Hi-C) to generate a contiguous assembly. The genome assembled from super-scaffolds had a size of 4.85 Gb, N50 of 2.32 Mb, and an estimated 93.9% of complete genes from a plant database (BUSCO, benchmarking universal single-copy orthologous genes). After removal of small scaffolds (< 300 Kb), the assembly was arranged into pseudomolecules of 4.14 Gb in size with seven chromosomes plus unanchored scaffolds. The completeness and annotation of the assembly were assessed by comparing it with the updated version of six-row Morex and recently released two-row Golden Promise genome assemblies.

Building pan-genome infrastructures for crop plants and their use in association genetics

Pan-genomic studies aim at representing the entire sequence diversity within a species to provide useful resources for evolutionary studies, functional genomics and breeding of cultivated plants. Cost reductions in high-throughput sequencing and advances in sequence assembly algorithms have made it possible to create multiple reference genomes along with a catalogue of all forms of genetic variations in plant species with large and complex or polyploid genomes. In this review, we summarize the current approaches to building pan-genomes as an in silico representation of plant sequence diversity and outline relevant methods for their effective utilization in linking structural with phenotypic variation. We propose as future research avenues (i) transcriptomic and epigenomic studies across multiple reference genomes and (ii) the development of user-friendly and feature-rich pan-genome browsers.

Improved genome assembly and pan-genome provide key insights into eggplant domestication and breeding

Eggplant (Solanum melongena L.) is an important horticultural crop and one of the most widely grown vegetables from the Solanaceae family. It was domesticated from a wild, prickly progenitor carrying small, round, non-anthocyanic fruits. We obtained a novel, highly contiguous genome assembly of the eggplant '67/3' reference line, by Hi-C retrofitting of a previously released short read- and optical mapping-based assembly. The sizes of the 12 chromosomes and the fraction of anchored genes in the improved assembly were comparable to those of a chromosome-level assembly. We resequenced 23 accessions of S. melongena representative of the worldwide phenotypic, geographic, and genetic diversity of the species, and one each from the closely related species Solanum insanum and Solanum incanum. The eggplant pan-genome contained approximately 51.5 additional megabases and 816 additional genes compared with the reference genome, while the pan-plastome showed little genetic variation. We identified 53 selective sweeps related to fruit color, prickliness, and fruit shape in the nuclear genome, highlighting selection leading to the emergence of present-day S. melongena cultivars from its wild ancestors. Candidate genes underlying the selective sweeps included a MYBL1 repressor and CHALCONE ISOMERASE (for fruit color), homologs of Arabidopsis GLABRA1 and GLABROUS INFLORESCENCE STEMS2 (for prickliness), and orthologs of tomato FW2.2, OVATE, LOCULE NUMBER/WUSCHEL, SUPPRESSOR OF OVATE, and CELL SIZE REGULATOR (for fruit size/shape), further suggesting that selection for the latter trait relied on a common set of orthologous genes in tomato and eggplant.

Identification and spatio-temporal expression analysis of barley genes that encode putative modular xylanolytic enzymes

Arabinoxylans are cell wall polysaccharides whose re-modelling and degradation during plant development are mediated by several classes of xylanolytic enzymes. Here, we present the identification and new annotation of twelve putative (1,4)-β-xylanase and six β-xylosidase genes, and their spatio-temporal expression patterns during vegetative and reproductive growth of barley (Hordeum vulgare cv. Navigator). The encoded xylanase proteins are all predicted to contain a conserved carbohydrate-binding module (CBM) and a catalytic glycoside hydrolase (GH) 10 domain. Additional domains in some xylanases define three discrete phylogenetic clades: one clade contains proteins with an additional N-terminal signal sequence, while another clade contains proteins with multiple CBMs. Homology modelling revealed that all fifteen xylanases likely contain a third domain, a β-sandwich folded from two non-contiguous sequence segments that bracket the catalytic GH domain, which may explain why the full length protein is required for correct folding of the active enzyme. Similarly, predicted xylosidase proteins share a highly conserved domain structure, each with an N-terminal signal peptide, a split GH 3 domain, and a C-terminal fibronectin-like domain. Several genes appear to be ubiquitously expressed during barley growth and development, while four newly annotated xylanase and xylosidase genes are expressed at extremely high levels, which may be of broader interest for industrial applications where cell wall degradation is necessary.

Optical maps refine the bread wheat Triticum aestivum cv. Chinese Spring genome assembly

Until recently, achieving a reference-quality genome sequence for bread wheat was long thought beyond the limits of genome sequencing and assembly technology, primarily due to the large genome size and > 80% repetitive sequence content. The release of the chromosome scale 14.5-Gb IWGSC RefSeq v1.0 genome sequence of bread wheat cv. Chinese Spring (CS) was, therefore, a milestone. Here, we used a direct label and stain (DLS) optical map of the CS genome together with a prior nick, label, repair and stain (NLRS) optical map, and sequence contigs assembled with Pacific Biosciences long reads, to refine the v1.0 assembly. Inconsistencies between the sequence and maps were reconciled and gaps were closed. Gap filling and anchoring of 279 unplaced scaffolds increased the total length of pseudomolecules by 168 Mb (excluding Ns). Positions and orientations were corrected for 233 and 354 scaffolds, respectively, representing 10% of the genome sequence. The accuracy of the remaining 90% of the assembly was validated. As a result of the increased contiguity, the numbers of transposable elements (TEs) and intact TEs have increased in IWGSC RefSeq v2.1 compared with v1.0. In total, 98% of the gene models identified in v1.0 were mapped onto this new assembly through development of a dedicated approach implemented in the MAGAAT pipeline. The numbers of high-confidence genes on pseudomolecules have increased from 105 319 to 105 534. The reconciled assembly enhances the utility of the sequence for genetic mapping, comparative genomics, gene annotation and isolation, and more general studies on the biology of wheat.

Analysis of barley mutants ert-c.1 and ert-d.7 reveals two loci with additive effect on plant architecture

Both mutant ert-c.1 and ert-d.7 carry T2-T3 translocations in the Ert-c gene. Principal coordinate analyses revealed the translocation types and translocation breakpoints. Mutant ert-d.7 is an Ert-c Ert-d double mutant. Mutations in the Ert-c and Ert-d loci are among the most common barley mutations affecting plant architecture. The mutants have various degrees of erect and compact spikes, often accompanied with short and stiff culms. In the current study, complementation tests, linkage mapping, principal coordinate analyses and fine mapping were conducted. We conclude that the original ert-d.7 mutant does not only carry an ert-d mutation but also an ert-c mutation. Combined, mutations in Ert-c and Ert-d cause a pyramid-dense spike phenotype, whereas mutations in only Ert-c or Ert-d give a pyramid and dense phenotype, respectively. Associations between the Ert-c gene and T2-T3 translocations were detected in both mutant ert-c.1 and ert-d.7. Different genetic association patterns indicate different translocation breakpoints in these two mutants. Principal coordinate analysis based on genetic distance and screening of recombinants from all four ends of polymorphic regions was an efficient way to narrow down the region of interest in translocation-involved populations. The Ert-c gene was mapped to the marker interval of 2_0801to1_0224 on 3HL near the centromere. The results illuminate a complex connection between two single genes having additive effects on barley spike architecture and will facilitate the identification of the Ert-c and Ert-d genes.

Improvements in the sequencing and assembly of plant genomes

Advances in DNA sequencing have made it easier to sequence and assemble plant genomes. Here, we extend an earlier study, and compare recent methods for long read sequencing and assembly. Updated Oxford Nanopore Technology software improved assemblies. Using more accurate sequences produced by repeated sequencing of the same molecule (Pacific Biosciences HiFi) resulted in less fragmented assembly of sequencing reads. Using data for increased genome coverage resulted in longer contigs, but reduced total assembly length and improved genome completeness. The original model species, Macadamia jansenii, was also compared with three other Macadamia species, as well as avocado (Persea americana) and jojoba (Simmondsia chinensis). In these angiosperms, increasing sequence data volumes caused a linear increase in contig size, decreased assembly length and further improved already high completeness. Differences in genome size and sequence complexity influenced the success of assembly. Advances in long read sequencing technology continue to improve plant genome sequencing and assembly. However, results were improved by greater genome coverage, with the amount needed to achieve a particular level of assembly being species dependent.

Genome-wide identification of bZIP transcription factor genes related to starch synthesis in barley (Hordeum vulgare L.)

The basic leucine zipper (bZIP) family of genes encode transcription factors that play key roles in plant growth and development. In this study, a total of 92 HvbZIP genes were identified and compared with previous studies using recently released barley genome data. Two novel genes were characterized in this study, and some misannotated and duplicated genes from previous studies have been corrected. Phylogenetic analysis results showed that 92 HvbZIP genes were classified into 10 groups and three unknown groups. The gene structure and motif distribution of the three unknown groups implied that the genes of the three groups may be functionally different. Expression profiling indicated that the HvbZIP genes exhibited different patterns of spatial and temporal expression. Using qRT-PCR, more than 10 HvbZIP genes were identified with expression patterns similar to those of starch synthase genes in barley. Yeast one-hybrid analysis revealed that two of the HvbZIP genes exhibited in vitro binding activity to the promoter of HvAGP-S. The two HvbZIP genes may be candidate genes for further study to explore the mechanism by which they regulate the synthesis of barley starch.

Genomic Prediction of Grain Yield in a Barley MAGIC Population Modeling Genotype per Environment Interaction

Multi-parent Advanced Generation Inter-crosses (MAGIC) lines have mosaic genomes that are generated shuffling the genetic material of the founder parents following pre-defined crossing schemes. In cereal crops, these experimental populations have been extensively used to investigate the genetic bases of several traits and dissect the genetic bases of epistasis. In plants, genomic prediction models are usually fitted using either diverse panels of mostly unrelated accessions or individuals of biparental families and several empirical analyses have been conducted to evaluate the predictive ability of models fitted to these populations using different traits. In this paper, we constructed, genotyped and evaluated a barley MAGIC population of 352 individuals developed with a diverse set of eight founder parents showing contrasting phenotypes for grain yield. We combined phenotypic and genotypic information of this MAGIC population to fit several genomic prediction models which were cross-validated to conduct empirical analyses aimed at examining the predictive ability of these models varying the sizes of training populations. Moreover, several methods to optimize the composition of the training population were also applied to this MAGIC population and cross-validated to estimate the resulting predictive ability. Finally, extensive phenotypic data generated in field trials organized across an ample range of water regimes and climatic conditions in the Mediterranean were used to fit and cross-validate multi-environment genomic prediction models including G×E interaction, using both genomic best linear unbiased prediction and reproducing kernel Hilbert space along with a non-linear Gaussian Kernel. Overall, our empirical analyses showed that genomic prediction models trained with a limited number of MAGIC lines can be used to predict grain yield with values of predictive ability that vary from 0.25 to 0.60 and that beyond QTL mapping and analysis of epistatic effects, MAGIC population might be used to successfully fit genomic prediction models. We concluded that for grain yield, the single-environment genomic prediction models examined in this study are equivalent in terms of predictive ability while, in general, multi-environment models that explicitly split marker effects in main and environmental-specific effects outperform simpler multi-environment models.

Genome-Wide Characterization of WRKY Transcription Factors Revealed Gene Duplication and Diversification in Populations of Wild to Domesticated Barley

The WRKY transcription factors (WRKYs) are known for their crucial roles in biotic and abiotic stress responses, and developmental and physiological processes. In barley, early studies revealed their importance, whereas their diversity at the population scale remains hardly estimated. In this study, 98 HsWRKYs and 103 HvWRKYs have been identified from the reference genome of wild and cultivated barley, respectively. The tandem duplication and segmental duplication events from the cultivated barley were observed. By taking advantage of early released exome-captured sequencing datasets in 90 wild barley accessions and 137 landraces, the diversity analysis uncovered synonymous and non-synonymous variants instead of loss-of-function mutations that had occurred at all WRKYs. For majority of WRKYs, the haplotype and nucleotide diversity both decreased in cultivated barley relative to the wild population. Five WRKYs were detected to have undergone selection, among which haplotypes of WRKY9 were enriched, correlating with the geographic collection sites. Collectively, profiting from the state-of-the-art barley genomic resources, this work represented the characterization and diversity of barley WRKY transcription factors, shedding light on future deciphering of their roles in barley domestication and adaptation.