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

Non-cell-autonomous signaling associated with barley ALOG1 specifies spikelet meristem determinacy

Inflorescence architecture and crop productivity are often tightly coupled in our major cereal crops. However, the underlying genetic mechanisms controlling cereal inflorescence development remain poorly understood. Here, we identified recessive alleles of barley (Hordeum vulgare L.) HvALOG1 (Arabidopsis thaliana LSH1 and Oryza G1) that produce non-canonical extra spikelets and fused glumes abaxially to the central spikelet from the upper-mid portion until the tip of the inflorescence. Notably, we found that HvALOG1 exhibits a boundary-specific expression pattern that specifically excludes reproductive meristems, implying the involvement of previously proposed localized signaling centers for branch regulation. Importantly, during early spikelet formation, non-cell-autonomous signals associated with HvALOG1 expression may specify spikelet meristem determinacy, while boundary formation of floret organs appears to be coordinated in a cell-autonomous manner. Moreover, barley ALOG family members synergistically modulate inflorescence morphology, with HvALOG1 predominantly governing meristem maintenance and floral organ development. We further propose that spatiotemporal redundancies of expressed HvALOG members specifically in the basal inflorescence may be accountable for proper patterning of spikelet formation in mutant plants. Our research offers new perspectives on regulatory signaling roles of ALOG transcription factors during the development of reproductive meristems in cereal inflorescences.

Transcriptional changes during crown-root development and emergence in barley (Hordeum vulgare L.)

BACKGROUND

Roots play an important role during plant growth and development, ensuring water and nutrient uptake. Understanding the mechanisms regulating their initiation and development opens doors towards root system architecture engineering.

RESULTS

Here, we investigated by RNA-seq analysis the changes in gene expression in the barley stem base of 1 day-after-germination (DAG) and 10DAG seedlings when crown roots are formed. We identified 2,333 genes whose expression was lower in the stem base of 10DAG seedlings compared to 1DAG seedlings. Those genes were mostly related to basal cellular activity such as cell cycle organization, protein biosynthesis, chromatin organization, cytoskeleton organization or nucleotide metabolism. In opposite, 2,932 genes showed up-regulation in the stem base of 10DAG seedlings compared to 1DAG seedlings, and their function was related to phytohormone action, solute transport, redox homeostasis, protein modification, secondary metabolism. Our results highlighted genes that are likely involved in the different steps of crown root formation from initiation to primordia differentiation and emergence, and revealed the activation of different hormonal pathways during this process.

CONCLUSIONS

This whole transcriptomic study is the first study aiming at understanding the molecular mechanisms controlling crown root development in barley. The results shed light on crown root emergence that is likely associated with a strong cell wall modification, death of the cells covering the crown root primordium, and the production of defense molecules that might prevent pathogen infection at the site of root emergence.

Are cereal grasses a single genetic system?

In 1993, a passionate and provocative call to arms urged cereal researchers to consider the taxon they study as a single genetic system and collaborate with each other. Since then, that group of scientists has seen their discipline blossom. In an attempt to understand what unity of genetic systems means and how the notion was borne out by later research, we survey the progress and prospects of cereal genomics: sequence assemblies, population-scale sequencing, resistance gene cloning and domestication genetics. Gene order may not be as extraordinarily well conserved in the grasses as once thought. Still, several recurring themes have emerged. The same ancestral molecular pathways defining plant architecture have been co-opted in the evolution of different cereal crops. Such genetic convergence as much as cross-fertilization of ideas between cereal geneticists has led to a rich harvest of genes that, it is hoped, will lead to improved varieties.

Disentangling hydroxynitrile glucoside biosynthesis in a barley (Hordeum vulgare) metabolon provides access to elite malting barleys for ethyl carbamate-free whisky production

Barley produces several specialized metabolites, including five α-, β-, and γ-hydroxynitrile glucosides (HNGs). In malting barley, presence of the α-HNG epiheterodendrin gives rise to undesired formation of ethyl carbamate in the beverage production, especially after distilling. Metabolite-GWAS identified QTLs and underlying gene candidates possibly involved in the control of the relative and absolute content of HNGs, including an undescribed MATE transporter. By screening 325 genetically diverse barley accessions, we discovered three H. vulgare ssp. spontaneum (wild barley) lines with drastic changes in the relative ratios of the five HNGs. Knock-out (KO)-lines, isolated from the barley FIND-IT resource and each lacking one of the functional HNG biosynthetic genes (CYP79A12, CYP71C103, CYP71C113, CYP71U5, UGT85F22 and UGT85F23) showed unprecedented changes in HNG ratios enabling assignment of specific and mutually dependent catalytic functions to the biosynthetic enzymes involved. The highly similar relative ratios between the five HNGs found across wild and domesticated barley accessions indicate assembly of the HNG biosynthetic enzymes in a metabolon, the functional output of which was reconfigured in the absence of a single protein component. The absence or altered ratios of the five HNGs in the KO-lines did not change susceptibility to the fungal phytopathogen Pyrenophora teres causing net blotch. The study provides a deeper understanding of the organization of HNG biosynthesis in barley and identifies a novel, single gene HNG-0 line in an elite spring barley background for direct use in breeding of malting barley, eliminating HNGs as a source of ethyl carbamate formation in whisky production.

A chromosome level genome assembly of Pseudoroegneria Libanotica reveals a key Kcs gene involves in the cuticular wax elongation for drought resistance

BACKGROUND

The genus Pseudoroegneria (Nevski) Löve (Triticeae, Poaceae), whose genome symbol was designed as "St", accounts for more than 60% of perennial Triticeae species. The diploid species Psudoroegneria libanotica (2n = 14) contains the most ancient St genome, exhibited strong drought resistance, and was morphologically covered by cuticular wax on the aerial part. Therefore, the St-genome sequencing data could provide fundamental information for studies of genome evolution and reveal its mechanisms of cuticular wax and drought resistance.

RESULTS

In this study, we reported the chromosome-level genome assembly for the St genome of Pse. libanotica, with a total size of 2.99 Gb. 46,369 protein-coding genes annotated and 71.62% was repeat sequences. Comparative analyses revealed that the genus Pseudoroegneria diverged during the middle and late Miocene. During this period, unique genes, gene family expansion, and contraction in Pse. libanotica were enriched in biotic and abiotic stresses, such as fatty acid biosynthesis which may greatly contribute to its drought adaption. Furthermore, we investigated genes associated with the cuticular wax formation and water deficit and found a new Kcs gene evm.TU.CTG175.54. It plays a critical role in the very long chain fatty acid (VLCFA) elongation from C18 to C26 in Pse. libanotica. The function needs more evidence to be verified.

CONCLUSIONS

We sequenced and assembled the St genome in Triticeae and discovered a new KCS gene that plays a role in wax extension to cope with drought. Our study lays a foundation for the genome diversification of Triticeae species and deciphers cuticular wax formation genes involved in drought resistance.

An improved chromosome-level genome assembly of perennial ryegrass (Lolium perenne L.)

This work is an update and extension of the previously published article "Ultralong Oxford Nanopore Reads Enable the Development of a Reference-Grade Perennial Ryegrass Genome Assembly" by Frei et al. The published genome assembly of the doubled haploid perennial ryegrass (Lolium perenne L.) genotype Kyuss (Kyuss v1.0) marked a milestone for forage grass research and breeding. However, order and orientation errors may exist in the pseudo-chromosomes of Kyuss, since barley (Hordeum vulgare L.), which diverged 30 million years ago from perennial ryegrass, was used as the reference to scaffold Kyuss. To correct for structural errors possibly present in the published Kyuss assembly, we de novo assembled the genome again and generated 50-fold coverage high-throughput chromosome conformation capture (Hi-C) data to assist pseudo-chromosome construction. The resulting new chromosome-level assembly Kyuss v2.0 showed improved quality with high contiguity (contig N50 = 120 Mb), high completeness (total BUSCO score = 99%), high base-level accuracy (QV = 50), and correct pseudo-chromosome structure (validated by Hi-C contact map). This new assembly will serve as a better reference genome for Lolium spp. and greatly benefit the forage and turf grass research community.

Genome-wide association mapping for seedling and adult resistance to powdery mildew in barley

KEY MESSAGE

Two new major QTL were identified for powdery mildew resistance. We confirmed that the QTL on 7HS contributed mainly to the adult-plant resistance, while another one on chromosome arm 1HS made a significant contribution to the seedling resistance. Powdery mildew (PM), caused by Blumeria hordei, can occur at all post emergent stages of barley and constantly threatens crop production. To identify more genes for effective resistance to powdery mildew for use in breeding programs, 696 barley accessions collected from different regions of the world were evaluated for PM resistance at seedling and adult growth stages in three different states of Australia. These barley accessions were genotyped using DArTSeq with over 18,000 markers for a genome-wide association study (GWAS). Using the FarmCPU model, 54 markers showed significant associations with PM resistance scored at the seedling and adult-plant stages in different states of Australia. Another 40 markers showed tentative associations (LOD > 4.0) with resistance. These markers are distributed across all seven barley chromosomes. Most of them were grouped into eleven QTL regions, coinciding with the locations of most of the reported resistance genes. Two major MTAs were identified on chromosome arms 3HS and 5HL, with one on 3HS contributing to adult plant resistance and the one on 5HL to both seedling and adult plant resistance. An MTA on 7HS contributed mainly to the adult-plant resistance, while another one on chromosome arm 1HS made a significant contribution to the seedling resistance.

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.

Phylogeny and functional diversity of halophilic microbial communities from a thalasso environment

The El-Rawda solar saltern, located in North Sinai, Egypt, is formed through the process of water evaporation from the Bradawil lagoon. This evaporation leads to the precipitation of gypsum, halite minerals, and salt flats, which subsequently cover the southern and eastern areas of the lagoon. This study employed the shotgun metagenomic approach, the illumine platform, and bioinformatic tools to investigate the taxonomic composition and functional diversity of halophilic microbial communities in solar saltern. The metagenomic reads obtained from the brine sample exhibited a greater count compared to those from the sediment sample. Notably, the brine sample was primarily characterized by an abundance of archaea, while the sediment sample displayed a dominant abundance of bacteria. Both samples exhibited a relatively low abundance of eukaryotes, while viruses were only found in the brine sample. Furthermore, the comparative analysis of functional pathways showed many important processes related to central metabolism and protein processing in brine and sediment samples. In brief, this research makes a valuable contribution to the understanding of very halophilic ecosystems in Egypt, providing insights into their microbial biodiversity and functional processes.

UDP-glucosyltransferase HvUGT13248 confers type II resistance to Fusarium graminearum in barley

Fusarium head blight (FHB) of barley (Hordeum vulgare) causes yield losses and accumulation of trichothecene mycotoxins (e.g. deoxynivalenol [DON]) in grains. Glucosylation of DON to the nontoxic DON-3-O-glucoside (D3G) is catalyzed by UDP-glucosyltransferases (UGTs), such as barley UGT13248. We explored the natural diversity of UGT13248 in 496 barley accessions and showed that all carried potential functional alleles of UGT13248, as no genotypes showed strongly increased seedling sensitivity to DON. From a TILLING population, we identified 2 mutant alleles (T368I and H369Y) that, based on protein modeling, likely affect the UDP-glucose binding of UGT13248. In DON feeding experiments, DON-to-D3G conversion was strongly reduced in spikes of these mutants compared to controls, and plants overexpressing UGT13248 showed increased resistance to DON and increased DON-to-D3G conversion. Moreover, field-grown plants carrying the T368I or H369Y mutations inoculated with Fusarium graminearum showed increased FHB disease severity and reduced D3G production. Barley is generally considered to have type II resistance that limits the spread of F. graminearum from the infected spikelet to adjacent spikelets. Point inoculation experiments with F. graminearum showed increased infection spread in T368I and H369Y across the spike compared to wild type, while overexpression plants showed decreased spread of FHB symptoms. Confocal microscopy revealed that F. graminearum spread to distant rachis nodes in T368I and H369Y mutants but was arrested at the rachis node of the inoculated spikelet in wild-type plants. Taken together, our data reveal that UGT13248 confers type II resistance to FHB in barley via conjugation of DON to D3G.

Barley shows reduced Fusarium head blight under drought and modular expression of differentially expressed genes under combined stress

Plants often face simultaneous abiotic and biotic stress conditions; however, physiological and transcriptional responses under such combined stress conditions are still not fully understood. Spring barley (Hordeum vulgare) is susceptible to Fusarium head blight (FHB), which is strongly affected by weather conditions. We therefore studied the potential influence of drought on FHB severity and plant responses in three varieties of different susceptibility. We found strongly reduced FHB severity in susceptible varieties under drought. The number of differentially expressed genes (DEGs) and strength of transcriptomic regulation reflected the concentrations of physiological stress markers such as abscisic acid or fungal DNA contents. Infection-related gene expression was associated with susceptibility rather than resistance. Weighted gene co-expression network analysis revealed 18 modules of co-expressed genes that reflected the pathogen- or drought-response in the three varieties. A generally infection-related module contained co-expressed genes for defence, programmed cell death, and mycotoxin detoxification, indicating that the diverse genotypes used a similar defence strategy towards FHB, albeit with different degrees of success. Further, DEGs showed co-expression in drought- or genotype-associated modules that correlated with measured phytohormones or the osmolyte proline. The combination of drought stress with infection led to the highest numbers of DEGs and resulted in a modular composition of the single-stress responses rather than a specific transcriptional output.

Cytokinin Translocation to, and Biosynthesis and Metabolism within, Cereal and Legume Seeds: Looking Back to Inform the Future

Early in the history of cytokinins, it was clear that Zea mays seeds contained not just trans-zeatin, but its nucleosides and nucleotides. Subsequently, both pods and seeds of legumes and cereal grains have been shown to contain a complex of cytokinin forms. Relative to the very high quantities of cytokinin detected in developing seeds, only a limited amount appears to have been translocated from the parent plant. Translocation experiments, and the detection of high levels of endogenous cytokinin in the maternal seed coat tissues of legumes, indicates that cytokinin does not readily cross the maternal/filial boundary, indicating that the filial tissues are autonomous for cytokinin biosynthesis. Within the seed, trans-zeatin plays a key role in sink establishment and it may also contribute to sink strength. The roles, if any, of the other biologically active forms of cytokinin (cis-zeatin, dihydrozeatin and isopentenyladenine) remain to be elucidated. The recent identification of genes coding for the enzyme that leads to the biosynthesis of trans-zeatin in rice (OsCYP735A3 and 4), and the identification of a gene coding for an enzyme (CPN1) that converts trans-zeatin riboside to trans-zeatin in the apoplast, further cements the key role played by trans-zeatin in plants.

A gap-free reference genome reveals structural variations associated with flowering time in rapeseed (Brassica napus)

Allopolyploid oilseed rape (Brassica napus) is an important oil crop and vegetable. However, the latest version of its reference genome, with collapsed duplications, gaps, and other issues, prevents comprehensive genomic analysis. Herein, we report a gap-free assembly of the rapeseed cv. Xiang5A genome using a combination of ONT (Oxford Nanopore Technologies) ultra-long reads, PacBio high-fidelity reads, and Hi-C datasets. It includes gap-free assemblies of all 19 chromosomes and telomere-to-telomere assemblies of eight chromosomes. Compared with previously published genomes of B. napus, our gap-free genome, with a contig N50 length of 50.70 Mb, has complete assemblies of 9 of 19 chromosomes without manual intervention, and greatly improves contiguity and completeness, thereby representing the highest quality genome assembly to date. Our results revealed that B. napus Xiang5A underwent nearly complete triplication and allotetraploidy relative to Arabidopsis thaliana. Using the gap-free assembly, we found that 917 flowering-related genes were affected by structural variation, including BnaA03.VERNALIZATION INSENSITIVE 3 and BnaC04.HIGH EXPRESSION OF OSMOTICALLY RESPONSIVE GENES 1. These genes may play crucial roles in regulating flowering time and facilitating the adaptation of Xiang5A in the Yangtze River Basin of China. This reference genome provides a valuable genetic resource for rapeseed functional genomic studies and breeding.

A pathogen-induced putative NAC transcription factor mediates leaf rust resistance in barley

Leaf rust, caused by Puccinia hordei, is one of the most widespread and damaging foliar diseases affecting barley. The barley leaf rust resistance locus Rph7 has been shown to have unusually high sequence and haplotype divergence. In this study, we isolate the Rph7 gene using a fine mapping and RNA-Seq approach that is confirmed by mutational analysis and transgenic complementation. Rph7 is a pathogen-induced, non-canonical resistance gene encoding a protein that is distinct from other known plant disease resistance proteins in the Triticeae. Structural analysis using an AlphaFold2 protein model suggests that Rph7 encodes a putative NAC transcription factor with a zinc-finger BED domain with structural similarity to the N-terminal DNA-binding domain of the NAC transcription factor (ANAC019) from Arabidopsis. A global gene expression analysis suggests Rph7 mediates the activation and strength of the basal defence response. The isolation of Rph7 highlights the diversification of resistance mechanisms available for engineering disease control in crops.

Genome-wide association and targeted transcriptomic analyses reveal loci and candidate genes regulating preharvest sprouting in barley

KEY MESSAGE

Genome-wide association study of diverse barley genotypes identified loci, single nucleotide polymorphisms and candidate genes that control seed dormancy and therefore enhance resistance to preharvest sprouting. Preharvest sprouting (PHS) causes significant yield and quality loss in barley and it is strongly associated with the level of seed dormancy. This study performed genome-wide association study using a collection of 255 diverse barley genotypes grown over four environments to identify loci controlling dormancy/PHS. Our phenotypic analysis revealed substantial variation in germination index/dormancy levels among the barley genotypes. Marker-trait association and linkage disequilibrium (LD) decay analyses identified 16 single nucleotide polymorphisms (SNPs) and two QTLs associated with dormancy/PHS, respectively, on chromosome 3H and 5H explaining 6.9% to 11.1% of the phenotypic variation. QTL.5H consist of 14 SNPs of which 12 SNPs satisfy the FDR threshold of α = 0.05, and it may represent the SD2 locus. The QTL on 3H consists of one SNP that doesn't satisfy FDR (α = 0.05). Genes harbouring the significant SNPs were analyzed for their expression pattern in the seeds of selected dormant and non-dormant genotypes. Of these genes, HvRCD1, HvPSRP1 and HvF3H exhibited differential expression between the dormant and non-dormant seed samples, suggesting their role in controlling seed dormancy/PHS. Three SNPs located within the differentially expressed genes residing in QTL.5H explained considerable phenotypic variation (≥ 8.6%), suggesting their importance in regulating PHS resistance. Analysis of the SNP marker data in QTL.5H identified a haplotype for PHS resistance. Overall, the study identified loci, SNPs and candidate genes that control dormancy and therefore play important roles in enhancing PHS resistance in barley through marker-assisted breeding.

Investigating microglia-neuron crosstalk by characterizing microglial contamination in human and mouse patch-seq datasets

Microglia are cells with diverse roles, including the regulation of neuronal excitability. We leveraged Patch-seq to assess the presence and effects of microglia in the local microenvironment of recorded neurons. We first quantified the amounts of microglial transcripts in three Patch-seq datasets of human and mouse neocortical neurons, observing extensive contamination. Variation in microglial contamination was explained foremost by donor identity, particularly in human samples, and additionally by neuronal cell type identity in mice. Gene set enrichment analysis suggests that microglial contamination is reflective of activated microglia, and that these transcriptional signatures are distinct from those captured via single-nucleus RNA-seq. Finally, neurons with greater microglial contamination differed markedly in their electrophysiological characteristics, including lowered input resistances and more depolarized action potential thresholds. Our results generalize beyond Patch-seq to suggest that activated microglia may be widely present across brain slice preparations and contribute to neuron- and donor-related electrophysiological variability in vitro.

Triticeae crop genome biology: an endless frontier

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

Genetic resources and precise gene editing for targeted improvement of barley abiotic stress tolerance

Abiotic stresses, predominately drought, heat, salinity, cold, and waterlogging, adversely affect cereal crops. They limit barley production worldwide and cause huge economic losses. In barley, functional genes under various stresses have been identified over the years and genetic improvement to stress tolerance has taken a new turn with the introduction of modern gene-editing platforms. In particular, clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) is a robust and versatile tool for precise mutation creation and trait improvement. In this review, we highlight the stress-affected regions and the corresponding economic losses among the main barley producers. We collate about 150 key genes associated with stress tolerance and combine them into a single physical map for potential breeding practices. We also overview the applications of precise base editing, prime editing, and multiplexing technologies for targeted trait modification, and discuss current challenges including high-throughput mutant genotyping and genotype dependency in genetic transformation to promote commercial breeding. The listed genes counteract key stresses such as drought, salinity, and nutrient deficiency, and the potential application of the respective gene-editing technologies will provide insight into barley improvement for climate resilience.

Genome analyses reveal population structure and a purple stigma color gene candidate in finger millet

Finger millet is a key food security crop widely grown in eastern Africa, India and Nepal. Long considered a 'poor man's crop', finger millet has regained attention over the past decade for its climate resilience and the nutritional qualities of its grain. To bring finger millet breeding into the 21^st century, here we present the assembly and annotation of a chromosome-scale reference genome. We show that this ~1.3 million years old allotetraploid has a high level of homoeologous gene retention and lacks subgenome dominance. Population structure is mainly driven by the differential presence of large wild segments in the pericentromeric regions of several chromosomes. Trait mapping, followed by variant analysis of gene candidates, reveals that loss of purple coloration of anthers and stigma is associated with loss-of-function mutations in the finger millet orthologs of the maize R1/B1 and Arabidopsis GL3/EGL3 anthocyanin regulatory genes. Proanthocyanidin production in seed is not affected by these gene knockouts.

Genetic Dissection of Spike Productivity Traits in the Siberian Collection of Spring Barley

Barley (Hordeum vulgare L.) is one of the most commonly cultivated cereals worldwide. Its local varieties can represent a valuable source of unique genetic variants useful for crop improvement. The aim of this study was to reveal loci contributing to spike productivity traits in Siberian spring barley and to develop diagnostic DNA markers for marker-assisted breeding programs. For this purpose we conducted a genome-wide association study using a panel of 94 barley varieties. In total, 64 SNPs significantly associated with productivity traits were revealed. Twenty-three SNP markers were validated by genotyping in an independent sample set using competitive allele-specific PCR (KASP). Finally, fourteen markers associated with spike productivity traits on chromosomes 2H, 4H and 5H can be suggested for use in breeding programs.

SWEET11b transports both sugar and cytokinin in developing barley grains

Even though Sugars Will Eventually be Exported Transporters (SWEETs) have been found in every sequenced plant genome, a comprehensive understanding of their functionality is lacking. In this study, we focused on the SWEET family of barley (Hordeum vulgare). A radiotracer assay revealed that expressing HvSWEET11b in African clawed frog (Xenopus laevis) oocytes facilitated the bidirectional transfer of not only just sucrose and glucose, but also cytokinin. Barley plants harboring a loss-of-function mutation of HvSWEET11b could not set viable grains, while the distribution of sucrose and cytokinin was altered in developing grains of plants in which the gene was knocked down. Sucrose allocation within transgenic grains was disrupted, which is consistent with the changes to the cytokinin gradient across grains, as visualized by magnetic resonance imaging and Fourier transform infrared spectroscopy microimaging. Decreasing HvSWEET11b expression in developing grains reduced overall grain size, sink strength, the number of endopolyploid endosperm cells, and the contents of starch and protein. The control exerted by HvSWEET11b over sugars and cytokinins likely predetermines their synergy, resulting in adjustments to the grain's biochemistry and transcriptome.

Barley's inability to germinate after submergence depends on hypoxia-induced secondary dormancy

Global climate change has dramatically increased flooding events, which have a strong impact on crop production. Barley is one of the most important cereals and its cultivation includes a broad range of different environments. We tested the capacity to germinate of a large barley panel after a short period of submergence followed by a recovery phase. We demonstrated that sensitive barley varieties activate underwater secondary dormancy because of a lower permeability to oxygen dissolved in water. In sensitive barley accessions, secondary dormancy is removed by nitric oxide donors. Our genome wide association study results uncovered a laccase gene located in a region of significant marker-trait association that is differently regulated during grain development and plays a key role in this process. We believe that our findings will help to improve the genetics of barley thereby increasing the capacity of seeds to germinate after a short period of flooding.

Novel exotic alleles of EARLY FLOWERING 3 determine plant development in barley

EARLY FLOWERING 3 (ELF3) is an important regulator of various physiological and developmental processes and hence may serve to improve plant adaptation which will be substantial for future plant breeding. To expand the limited knowledge on barley ELF3 in determining agronomic traits, we conducted field studies with heterogeneous inbred families (HIFs) derived from selected lines of the wild barley nested association mapping population HEB-25. During two growing seasons, phenotypes of nearly isogenic HIF sister lines, segregating for exotic and cultivated alleles at the ELF3 locus, were compared for ten developmental and yield-related traits. We determine novel exotic ELF3 alleles and show that HIF lines, carrying the exotic ELF3 allele, accelerated plant development compared to the cultivated ELF3 allele, depending on the genetic background. Remarkably, the most extreme effects on phenology could be attributed to one exotic ELF3 allele differing from the cultivated Barke ELF3 allele in only one SNP. This SNP causes an amino acid substitution (W669G), which predictively has an impact on the protein structure of ELF3, thereby possibly affecting phase separation behaviour and nano-compartment formation of ELF3 and, potentially, also affecting its local cellular interactions causing significant trait differences between HIF sister lines.

BarleyExpDB: an integrative gene expression database for barley

BACKGROUND

RNA-sequencing (RNA-seq) has been widely used to study the dynamic expression patterns of transcribed genes, which can lead to new biological insights. However, processing and analyzing these huge amounts of histological data remains a great challenge for wet labs and field researchers who lack bioinformatics experience and computational resources.

RESULTS

We present BarleyExpDB, an easy-to-operate, free, and web-accessible database that integrates transcriptional profiles of barley at different growth and developmental stages, tissues, and stress conditions, as well as differential expression of mutants and populations to build a platform for barley expression and visualization. The expression of a gene of interest can be easily queried by searching by known gene ID or sequence similarity. Expression data can be displayed as a heat map, along with functional descriptions as well as Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Proteins Families Database, and Simple Modular Architecture Research Tool annotations.

CONCLUSIONS

BarleyExpDB will serve as a valuable resource for the barley research community to leverage the vast publicly available RNA-seq datasets for functional genomics research and crop molecular breeding.

Genetic Localization and Homologous Genes Mining for Barley Grain Size

Grain size is an important agronomic trait determining barley yield and quality. An increasing number of QTLs (quantitative trait loci) for grain size have been reported due to the improvement in genome sequencing and mapping. Elucidating the molecular mechanisms underpinning barley grain size is vital for producing elite cultivars and accelerating breeding processes. In this review, we summarize the achievements in the molecular mapping of barley grain size over the past two decades, highlighting the results of QTL linkage analysis and genome-wide association studies. We discuss the QTL hotspots and predict candidate genes in detail. Moreover, reported homologs that determine the seed size clustered into several signaling pathways in model plants are also listed, providing the theoretical basis for mining genetic resources and regulatory networks of barley grain size.

A molecular framework for grain number determination in barley

Flowering plants with indeterminate inflorescences often produce more floral structures than they require. We found that floral primordia initiations in barley (Hordeum vulgare L.) are molecularly decoupled from their maturation into grains. While initiation is dominated by flowering-time genes, floral growth is specified by light signaling, chloroplast, and vascular developmental programs orchestrated by barley CCT MOTIF FAMILY 4 (HvCMF4), which is expressed in the inflorescence vasculature. Consequently, mutations in HvCMF4 increase primordia death and pollination failure, mainly through reducing rachis greening and limiting plastidial energy supply to developing heterotrophic floral tissues. We propose that HvCMF4 is a sensory factor for light that acts in connection with the vascular-localized circadian clock to coordinate floral initiation and survival. Notably, stacking beneficial alleles for both primordia number and survival provides positive implications on grain production. Our findings provide insights into the molecular underpinnings of grain number determination in cereal crops.

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.

Characterization of Leaf Rust Resistance in International Barley Germplasm Using Genome-Wide Association Studies

A panel of 114 genetically diverse barley lines were assessed in the greenhouse and field for resistance to the pathogen Puccinia hordei, the causal agent of barley leaf rust. Multi-pathotype tests revealed that 16.6% of the lines carried the all-stage resistance (ASR) gene Rph3, followed by Rph2 (4.4%), Rph1 (1.7%), Rph12 (1.7%) or Rph19 (1.7%). Five lines (4.4%) were postulated to carry the gene combinations Rph2+9.am, Rph2+19 and Rph8+19. Three lines (2.6%) were postulated to carry Rph15 based on seedling rust tests and genotyping with a marker linked closely to this gene. Based on greenhouse seedling tests and adult-plant field tests, 84 genotypes (73.7%) were identified as carrying APR, and genotyping with molecular markers linked closely to three known APR genes (Rph20, Rph23 and Rph24) revealed that 48 of the 84 genotypes (57.1%) likely carry novel (uncharacterized) sources of APR. Seven lines were found to carry known APR gene combinations (Rph20+Rph23, Rph23+Rph24 and Rph20+Rph24), and these lines had higher levels of field resistance compared to those carrying each of these three APR genes singly. GWAS identified 12 putative QTLs; strongly associated markers located on chromosomes 1H, 2H, 3H, 5H and 7H. Of these, the QTL on chromosome 7H had the largest effect on resistance response to P. hordei. Overall, these studies detected several potentially novel genomic regions associated with resistance. The findings provide useful information for breeders to support the utilization of these sources of resistance to diversify resistance to leaf rust in barley and increase resistance durability.

Transcriptome profiling disclosed the effect of single and combined drought and heat stress on reprogramming of genes expression in barley flag leaf

Despite numerous studies aimed at unraveling the genetic background of barley's response to abiotic stress, the modulation of the transcriptome induced by combinatorial drought and increased temperature remains largely unrecognized. Very limited studies were done, especially on the flag leaf, which plays an important role in grain filling in cereals. In the present study, transcriptome profiles, along with chlorophyll fluorescence parameters and yield components, were compared between barley genotypes with different flag leaf sizes under single and combined drought and heat stress. High-throughput mRNA sequencing revealed 2,457 differentially expressed genes, which were functionally interpreted using Gene Ontology term enrichment analysis. The transcriptomic signature under double stress was more similar to effects caused by drought than by elevated temperature; it was also manifested at phenotypic and chlorophyll fluorescence levels. Both common and stress-specific changes in transcript abundance were identified. Genes regulated commonly across stress treatments, determining universal stress responses, were associated, among others, with responses to drought, heat, and oxidative stress. In addition, changes specific to the size of the flag leaf blade were found. Our study allowed us to identify sets of genes assigned to various processes underlying the response to drought and heat, including photosynthesis, the abscisic acid pathway, and lipid transport. Genes encoding LEA proteins, including dehydrins and heat shock proteins, were especially induced by stress treatments. Some association between genetic composition and flag leaf size was confirmed. However, there was no general coincidence between SNP polymorphism of genotypes and differential expression of genes induced by stress factors. This research provided novel insight into the molecular mechanisms of barley flag leaf that determine drought and heat response, as well as their co-occurrence.

Sequencing and Assembly of Polyploid Genomes

Polyploidy has been observed throughout major eukaryotic clades and has played a vital role in the evolution of angiosperms. Recent polyploidizations often result in highly complex genome structures, posing challenges to genome assembly and phasing. Recent advances in sequencing technologies and genome assembly algorithms have enabled high-quality, near-complete chromosome-level assemblies of polyploid genomes. Advances in novel sequencing technologies include highly accurate single-molecule sequencing with HiFi reads, chromosome conformation capture with Hi-C technique, and linked reads sequencing. Additionally, new computational approaches have also significantly improved the precision and reliability of polyploid genome assembly and phasing, such as HiCanu, hifiasm, ALLHiC, and PolyGembler. Herein, we review recently published polyploid genomes and compare the various sequencing, assembly, and phasing approaches that are utilized in these genome studies. Finally, we anticipate that accurate and telomere-to-telomere chromosome-level assembly of polyploid genomes could ultimately become a routine procedure in the near future.

Super-resolution microscopy reveals the number and distribution of topoisomerase IIα and CENH3 molecules within barley metaphase chromosomes

Topoisomerase IIα (Topo IIα) and the centromere-specific histone H3 variant CENH3 are key proteins involved in chromatin condensation and centromere determination, respectively. Consequently, they are required for proper chromosome segregation during cell divisions. We combined two super-resolution techniques, structured illumination microscopy (SIM) to co-localize Topo IIα and CENH3, and photoactivated localization microscopy (PALM) to determine their molecule numbers in barley metaphase chromosomes. We detected a dispersed Topo IIα distribution along chromosome arms but an accumulation at centromeres, telomeres, and nucleolus-organizing regions. With a precision of 10-50 nm, we counted ~ 20,000-40,000 Topo IIα molecules per chromosome, 28% of them within the (peri)centromere. With similar precision, we identified ~13,500 CENH3 molecules per centromere where Topo IIα proteins and CENH3-containing chromatin intermingle. In short, we demonstrate PALM as a useful method to count and localize single molecules with high precision within chromosomes. The ultrastructural distribution and the detected amount of Topo IIα and CENH3 are instrumental for a better understanding of their functions during chromatin condensation and centromere determination.

Genome-Wide Analysis of Snf2 Gene Family Reveals Potential Role in Regulation of Spike Development in Barley

Sucrose nonfermenting 2 (Snf2) family proteins, as the catalytic core of ATP-dependent chromatin remodeling complexes, play important roles in nuclear processes as diverse as DNA replication, transcriptional regulation, and DNA repair and recombination. The Snf2 gene family has been characterized in several plant species; some of its members regulate flower development in Arabidopsis. However, little is known about the members of the family in barley (Hordeum vulgare). Here, 38 Snf2 genes unevenly distributed among seven chromosomes were identified from the barley (cv. Morex) genome. Phylogenetic analysis categorized them into 18 subfamilies. They contained combinations of 21 domains and consisted of 3 to 34 exons. Evolution analysis revealed that segmental duplication contributed predominantly to the expansion of the family in barley, and the duplicated gene pairs have undergone purifying selection. About eight hundred Snf2 family genes were identified from 20 barley accessions, ranging from 38 to 41 genes in each. Most of these genes were subjected to purification selection during barley domestication. Most were expressed abundantly during spike development. This study provides a comprehensive characterization of barley Snf2 family members, which should help to improve our understanding of their potential regulatory roles in barley spike development.

Defining Composition and Function of the Rhizosphere Microbiota of Barley Genotypes Exposed to Growth-Limiting Nitrogen Supplies

The microbiota populating the rhizosphere, the interface between roots and soil, can modulate plant growth, development, and health. These microbial communities are not stochastically assembled from the surrounding soil, but their composition and putative function are controlled, at least partially, by the host plant. Here, we use the staple cereal barley as a model to gain novel insights into the impact of differential applications of nitrogen, a rate-limiting step for global crop production, on the host genetic control of the rhizosphere microbiota. Using a high-throughput amplicon sequencing survey, we determined that nitrogen availability for plant uptake is a factor promoting the selective enrichment of individual taxa in the rhizosphere of wild and domesticated barley genotypes. Shotgun sequencing and metagenome-assembled genomes revealed that this taxonomic diversification is mirrored by a functional specialization, manifested by the differential enrichment of multiple Gene Ontology terms, of the microbiota of plants exposed to nitrogen conditions limiting barley growth. Finally, a plant soil feedback experiment revealed that host control of the barley microbiota underpins the assembly of a phylogenetically diverse group of bacteria putatively required to sustain plant performance under nitrogen-limiting supplies. Taken together, our observations indicate that under nitrogen conditions limiting plant growth, host-microbe and microbe-microbe interactions fine-tune the host genetic selection of the barley microbiota at both taxonomic and functional levels. The disruption of these recruitment cues negatively impacts plant growth. IMPORTANCE The microbiota inhabiting the rhizosphere, the thin layer of soil surrounding plant roots, can promote the growth, development, and health of their host plants. Previous research indicated that differences in the genetic composition of the host plant coincide with variations in the composition of the rhizosphere microbiota. This is particularly evident when looking at the microbiota associated with input-demanding modern cultivated varieties and their wild relatives, which have evolved under marginal conditions. However, the functional significance of these differences remains to be fully elucidated. We investigated the rhizosphere microbiota of wild and cultivated genotypes of the global crop barley and determined that nutrient conditions limiting plant growth amplify the host control on microbes at the root-soil interface. This is reflected in a plant- and genotype-dependent functional specialization of the rhizosphere microbiota, which appears to be required for optimal plant growth. These findings provide novel insights into the significance of the rhizosphere microbiota for plant growth and sustainable agriculture.

A technical guide to TRITEX, a computational pipeline for chromosome-scale sequence assembly of plant genomes

BACKGROUND

As complete and accurate genome sequences are becoming easier to obtain, more researchers wish to get one or more of them to support their research endeavors. Reliable and well-documented sequence assembly workflows find use in reference or pangenome projects.

RESULTS

We describe modifications to the TRITEX genome assembly workflow motivated by the rise of fast and easy long-read contig assembly of inbred plant genomes and the routine deployment of the toolchains in pangenome projects. New features include the use as surrogates of or complements to dense genetic maps and the introduction of user-editable tables to make the curation of contig placements easier and more intuitive.

CONCLUSION

Even maximally contiguous sequence assemblies of the telomere-to-telomere sort, and to a yet greater extent, the fragmented kind require validation, correction, and comparison to reference standards. As pangenomics is burgeoning, these tasks are bound to become more widespread and TRITEX is one tool to get them done. This technical guide is supported by a step-by-step computational tutorial accessible under https://tritexassembly.bitbucket.io/ . The TRITEX source code is hosted under this URL: https://bitbucket.org/tritexassembly .

Discovery of DNA polymorphisms via genome-resequencing and development of molecular markers between two barley cultivars

Genome resequencing uncovers genome-wide DNA polymorphisms that are useful for the development of high-density InDel markers between two barley cultivars. Discovering genomic variations and developing genetic markers are crucial for genetics studies and molecular breeding in cereal crops. Although InDels (insertions and deletions) have become popular because of their abundance and ease of detection, discovery of genome-wide DNA polymorphisms and development of InDel markers in barley have lagged behind other cereal crops such as rice, maize and wheat. In this study, we re-sequenced two barley cultivars, Golden Promise (GP, a classic British spring barley variety) and Hua30 (a Chinese spring barley variety), and mapped clean reads to the reference Morex genome, and identified in total 13,933,145 single nucleotide polymorphisms (SNPs) and 1,240,456 InDels for GP with Morex, 11,297,100 SNPs and 781,687 InDels for Hua30 with Morex, and 13,742,399 SNPs and 1,191,597 InDels for GP with Hua30. We further characterized distinct types, chromosomal distribution patterns, genome location, functional effect, and other features of these DNA polymorphisms. Additionally, we revealed the functional relevance of these identified SNPs/InDels regarding different flowering times between Hua30 and GP within 17 flowering time genes. Furthermore, we developed a series of InDel markers and validated them experimentally in 43 barley core accessions, respectively. Finally, we rebuilt population structure and phylogenetic tree of these 43 barley core accessions. Collectively, all of these genetic resources will facilitate not only the basic research but also applied research in barley.

Evolutionary Conservation and Transcriptome Analyses Attribute Perenniality and Flowering to Day-Length Responsive Genes in Bulbous Barley (Hordeum bulbosum)

Rapid population growth and dramatic climatic turnovers are challenging global crop production. These challenges are spurring plant breeders to enhance adaptation and sustainability of major crops. One intriguing approach is to turn annual systems into perennial ones, yet long-term classical breeding efforts to induce perenniality have achieved limited success. Here, we report the results of our investigation of the genetic basis of bulb formation in the nonmodel organism Hordeum bulbosum, a perennial species closely related to barley. To identify candidate genes that regulate bulb formation in H. bulbosum, we applied two complementary approaches. First, we explored the evolutionary conservation of expressed genes among annual Poaceae species. Next, we assembled a reference transcriptome for H. bulbosum and conducted a differential expression (DE) analysis before and after stimulating bulb initiation. Low conservation was identified in genes related to perenniality in H. bulbosum compared with other species, including bulb development and sugar accumulation genes. We also inspected these genes using a DE analysis, which enabled identification of additional genes responsible for bulb initiation and flowering regulation. We propose a molecular model for the regulation of bulb formation involving storage organ development and starch biosynthesis genes. The high conservation observed along a major part of the pathway between H. bulbosum and barley suggests a potential for the application of biotechnological techniques to accelerate breeding toward perenniality in barley.

The barley DIR gene family: An expanded gene family that is involved in stress responses

Gene family expansion plays a central role in adaptive divergence and, ultimately, speciation is influenced by phenotypic diversity in different environments. Barley (Hordeum vulgare) is the fourth most important cereal crop in the world and is used for brewing purposes, animal feed, and human food. Systematic characterization of expanded gene families is instrumental in the research of the evolutionary history of barley and understanding of the molecular function of their gene products. A total of 31,750 conserved orthologous groups (OGs) were identified using eight genomes/subgenomes, of which 1,113 and 6,739 were rapidly expanded and contracted OGs in barley, respectively. Five expanded OGs containing 20 barley dirigent genes (HvDIRs) were identified. HvDIRs from the same OG were phylogenetically clustered with similar gene structure and domain organization. In particular, 7 and 5 HvDIRs from OG0000960 and OG0001516, respectively, contributed greatly to the expansion of the DIR-c subfamily. Tandem duplication was the driving force for the expansion of the barley DIR gene family. Nucleotide diversity and haplotype network analysis revealed that the expanded HvDIRs experienced severe bottleneck events during barley domestication, and can thus be considered as potential domestication-related candidate genes. The expression profile and co-expression network analysis revealed the critical roles of the expanded HvDIRs in various biological processes, especially in stress responses. HvDIR18, HvDIR19, and HvDIR63 could serve as excellent candidates for further functional genomics studies to improve the production of barley products. Our study revealed that the HvDIR family was significantly expanded in barley and might be involved in different developmental processes and stress responses. Thus, besides providing a framework for future functional genomics and metabolomics studies, this study also identified HvDIRs as candidates for use in improving barley crop resistance to biotic and abiotic stresses.

Conserved signalling components coordinate epidermal patterning and cuticle deposition in barley

Faced with terrestrial threats, land plants seal their aerial surfaces with a lipid-rich cuticle. To breathe, plants interrupt their cuticles with adjustable epidermal pores, called stomata, that regulate gas exchange, and develop other specialised epidermal cells such as defensive hairs. Mechanisms coordinating epidermal features remain poorly understood. Addressing this, we studied two loci whose allelic variation causes both cuticular wax-deficiency and misarranged stomata in barley, identifying the underlying genes, Cer-g/ HvYDA1, encoding a YODA-like (YDA) MAPKKK, and Cer-s/ HvBRX-Solo, encoding a single BREVIS-RADIX (BRX) domain protein. Both genes control cuticular integrity, the spacing and identity of epidermal cells, and barley's distinctive epicuticular wax blooms, as well as stomatal patterning in elevated CO₂ conditions. Genetic analyses revealed epistatic and modifying relationships between HvYDA1 and HvBRX-Solo, intimating that their products participate in interacting pathway(s) linking epidermal patterning with cuticular properties in barley. This may represent a mechanism for coordinating multiple adaptive features of the land plant epidermis in a cultivated cereal.

Structural variants in the barley gene pool: precision and sensitivity to detect them using short-read sequencing and their association with gene expression and phenotypic variation

Structural variants (SV) of 23 barley inbreds, detected by the best combination of SV callers based on short-read sequencing, were associated with genome-wide and gene-specific gene expression and, thus, were evaluated to predict agronomic traits. In human genetics, several studies have shown that phenotypic variation is more likely to be caused by structural variants (SV) than by single nucleotide variants. However, accurate while cost-efficient discovery of SV in complex genomes remains challenging. The objectives of our study were to (i) facilitate SV discovery studies by benchmarking SV callers and their combinations with respect to their sensitivity and precision to detect SV in the barley genome, (ii) characterize the occurrence and distribution of SV clusters in the genomes of 23 barley inbreds that are the parents of a unique resource for mapping quantitative traits, the double round robin population, (iii) quantify the association of SV clusters with transcript abundance, and (iv) evaluate the use of SV clusters for the prediction of phenotypic traits. In our computer simulations based on a sequencing coverage of 25x, a sensitivity > 70% and precision > 95% was observed for all combinations of SV types and SV length categories if the best combination of SV callers was used. We observed a significant (P < 0.05) association of gene-associated SV clusters with global gene-specific gene expression. Furthermore, about 9% of all SV clusters that were within 5 kb of a gene were significantly (P < 0.05) associated with the gene expression of the corresponding gene. The prediction ability of SV clusters was higher compared to that of single-nucleotide polymorphisms from an array across the seven studied phenotypic traits. These findings suggest the usefulness of exploiting SV information when fine mapping and cloning the causal genes underlying quantitative traits as well as the high potential of using SV clusters for the prediction of phenotypes in diverse germplasm sets.

Multi-environment genome -wide association mapping of culm morphology traits in barley

In cereals with hollow internodes, lodging resistance is influenced by morphological characteristics such as internode diameter and culm wall thickness. Despite their relevance, knowledge of the genetic control of these traits and their relationship with lodging is lacking in temperate cereals such as barley. To fill this gap, we developed an image analysis-based protocol to accurately phenotype culm diameters and culm wall thickness across 261 barley accessions. Analysis of culm trait data collected from field trials in seven different environments revealed high heritability values (>50%) for most traits except thickness and stiffness, as well as genotype-by-environment interactions. The collection was structured mainly according to row-type, which had a confounding effect on culm traits as evidenced by phenotypic correlations. Within both row-type subsets, outer diameter and section modulus showed significant negative correlations with lodging (<-0.52 and <-0.45, respectively), but no correlation with plant height, indicating the possibility of improving lodging resistance independent of plant height. Using 50k iSelect SNP genotyping data, we conducted multi-environment genome-wide association studies using mixed model approach across the whole panel and row-type subsets: we identified a total of 192 quantitative trait loci (QTLs) for the studied traits, including subpopulation-specific QTLs and 21 main effect loci for culm diameter and/or section modulus showing effects on lodging without impacting plant height. Providing insights into the genetic architecture of culm morphology in barley and the possible role of candidate genes involved in hormone and cell wall-related pathways, this work supports the potential of loci underpinning culm features to improve lodging resistance and increase barley yield stability under changing environments.

BGFD: an integrated multi-omics database of barley gene families

BACKGROUND

A gene family comprises a group of genes with similar functional domains that play various roles in plant growth, development, and responses to environmental stimuli. Barley (Hordeum vulgare L.) is the fourth most cultivated cereal crop worldwide, and it is an important model species for genetic studies. Systematic identification and annotation of gene families are key for studies of molecular function and evolutionary history.

RESULTS

We constructed a multi-omics database containing 5593 genes of 77 gene families called the Barley Gene Family Database (BGFD: http://barleygfdb.com ). BGFD is a free, user-friendly, and web-accessible platform that provides data on barley family genes. BGFD provides intuitive visual displays to facilitate studies of the physicochemical properties, gene structure, phylogenetic relationships, and motif organization of genes. Massive multi-omics datasets have been acquired and processed to generate an atlas of expression pattern profiles and genetic variation in BGFD. The platform offers several practical toolkits to conduct searches, browse, and employ BLAST functions, and the data are downloadable.

CONCLUSIONS

BGFD will aid research on the domestication and adaptive evolution of barley; it will also facilitate the screening of candidate genes and exploration of important agronomic traits in barley.

The TCP transcription factor HvTB2 heterodimerizes with VRS5 and controls spike architecture in barley

Understanding the molecular network, including protein-protein interactions, of VRS5 provide new routes towards the identification of other key regulators of plant architecture in barley. The TCP transcriptional regulator TEOSINTE BRANCHED 1 (TB1) is a key regulator of plant architecture. In barley, an important cereal crop, HvTB1 (also referred to as VULGARE SIX-ROWED spike (VRS) 5), inhibits the outgrowth of side shoots, or tillers, and grains. Despite its key role in barley development, there is limited knowledge on the molecular network that is utilized by VRS5. In this work, we performed protein-protein interaction studies of VRS5. Our analysis shows that VRS5 potentially interacts with a diverse set of proteins, including other class II TCP's, NF-Y TF, but also chromatin remodelers. Zooming in on the interaction capacity of VRS5 with other TCP TFs shows that VRS5 preferably interacts with other class II TCP TFs in the TB1 clade. Induced mutagenesis through CRISPR-Cas of one of the putative VRS5 interactors, HvTB2 (also referred to as COMPOSITUM 1 and BRANCHED AND INDETERMINATE SPIKELET 1), resulted in plants that have lost their characteristic unbranched spike architecture. More specifically, hvtb2 mutants exhibited branches arising at the main spike, suggesting that HvTB2 acts as inhibitor of branching. Our protein-protein interaction studies of VRS5 resulted in the identification of HvTB2 as putative interactor of VRS5, another key regulator of spike architecture in barley. The study presented here provides a first step to underpin the protein-protein interactome of VRS5 and to identify other, yet unknown, key regulators of barley plant architecture.

A catalogue of resistance gene homologs and a chromosome-scale reference sequence support resistance gene mapping in winter wheat

A resistance gene atlas is an integral component of the breeder's arsenal in the fight against evolving pathogens. Thanks to high-throughput sequencing, catalogues of resistance genes can be assembled even in crop species with large and polyploid genomes. Here, we report on capture sequencing and assembly of resistance gene homologs in a diversity panel of 907 winter wheat genotypes comprising ex situ genebank accessions and current elite cultivars. In addition, we use accurate long-read sequencing and chromosome conformation capture sequencing to construct a chromosome-scale genome sequence assembly of cv. Attraktion, an elite variety representative of European winter wheat. We illustrate the value of our resource for breeders and geneticists by (i) comparing the resistance gene complements in plant genetic resources and elite varieties and (ii) conducting genome-wide associations scans (GWAS) for the fungal diseases yellow rust and leaf rust using reference-based and reference-free GWAS approaches. The gene content under GWAS peaks was scrutinized in the assembly of cv. Attraktion.

FIND-IT: Accelerated trait development for a green evolution

Improved agricultural and industrial production organisms are required to meet the future global food demands and minimize the effects of climate change. A new resource for crop and microbe improvement, designated FIND-IT (Fast Identification of Nucleotide variants by droplet DigITal PCR), provides ultrafast identification and isolation of predetermined, targeted genetic variants in a screening cycle of less than 10 days. Using large-scale sample pooling in combination with droplet digital PCR (ddPCR) greatly increases the size of low-mutation density and screenable variant libraries and the probability of identifying the variant of interest. The method is validated by screening variant libraries totaling 500,000 barley (Hordeum vulgare) individuals and isolating more than 125 targeted barley gene knockout lines and miRNA or promoter variants enabling functional gene analysis. FIND-IT variants are directly applicable to elite breeding pipelines and minimize time-consuming technical steps to accelerate the evolution of germplasm.

BaRTv2: a highly resolved barley reference transcriptome for accurate transcript-specific RNA-seq quantification

Accurate characterisation of splice junctions (SJs) as well as transcription start and end sites in reference transcriptomes allows precise quantification of transcripts from RNA-seq data, and enables detailed investigations of transcriptional and post-transcriptional regulation. Using novel computational methods and a combination of PacBio Iso-seq and Illumina short-read sequences from 20 diverse tissues and conditions, we generated a comprehensive and highly resolved barley reference transcript dataset from the European 2-row spring barley cultivar Barke (BaRTv2.18). Stringent and thorough filtering was carried out to maintain the quality and accuracy of the SJs and transcript start and end sites. BaRTv2.18 shows increased transcript diversity and completeness compared with an earlier version, BaRTv1.0. The accuracy of transcript level quantification, SJs and transcript start and end sites have been validated extensively using parallel technologies and analysis, including high-resolution reverse transcriptase-polymerase chain reaction and 5'-RACE. BaRTv2.18 contains 39 434 genes and 148 260 transcripts, representing the most comprehensive and resolved reference transcriptome in barley to date. It provides an important and high-quality resource for advanced transcriptomic analyses, including both transcriptional and post-transcriptional regulation, with exceptional resolution and precision.

Chromosome-scale assembly and annotation of the perennial ryegrass genome

BACKGROUND

The availability of chromosome-scale genome assemblies is fundamentally important to advance genetics and breeding in crops, as well as for evolutionary and comparative genomics. The improvement of long-read sequencing technologies and the advent of optical mapping and chromosome conformation capture technologies in the last few years, significantly promoted the development of chromosome-scale genome assemblies of model plants and crop species. In grasses, chromosome-scale genome assemblies recently became available for cultivated and wild species of the Triticeae subfamily. Development of state-of-the-art genomic resources in species of the Poeae subfamily, which includes important crops like fescues and ryegrasses, is lagging behind the progress in the cereal species.

RESULTS

Here, we report a new chromosome-scale genome sequence assembly for perennial ryegrass, obtained by combining PacBio long-read sequencing, Illumina short-read polishing, BioNano optical mapping and Hi-C scaffolding. More than 90% of the total genome size of perennial ryegrass (approximately 2.55 Gb) is covered by seven pseudo-chromosomes that show high levels of collinearity to the orthologous chromosomes of Triticeae species. The transposon fraction of perennial ryegrass was found to be relatively low, approximately 35% of the total genome content, which is less than half of the genome repeat content of cultivated cereal species. We predicted 54,629 high-confidence gene models, 10,287 long non-coding RNAs and a total of 8,393 short non-coding RNAs in the perennial ryegrass genome.

CONCLUSIONS

The new reference genome sequence and annotation presented here are valuable resources for comparative genomic studies in grasses, as well as for breeding applications and will expedite the development of productive varieties in perennial ryegrass and related species.

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

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

Prospects of telomere-to-telomere assembly in barley: Analysis of sequence gaps in the MorexV3 reference genome

The first gapless, telomere-to-telomere (T2T) sequence assemblies of plant chromosomes were reported recently. However, sequence assemblies of most plant genomes remain fragmented. Only recent breakthroughs in accurate long-read sequencing have made it possible to achieve highly contiguous sequence assemblies with a few tens of contigs per chromosome, that is a number small enough to allow for a systematic inquiry into the causes of the remaining sequence gaps and the approaches and resources needed to close them. Here, we analyse sequence gaps in the current reference genome sequence of barley cv. Morex (MorexV3). Optical map and sequence raw data, complemented by ChIP-seq data for centromeric histone variant CENH3, were used to estimate the abundance of centromeric, ribosomal DNA, and subtelomeric repeats in the barley genome. These estimates were compared with copy numbers in the MorexV3 pseudomolecule sequence. We found that almost all centromeric sequences and 45S ribosomal DNA repeat arrays were absent from the MorexV3 pseudomolecules and that the majority of sequence gaps can be attributed to assembly breakdown in long stretches of satellite repeats. However, missing sequences cannot fully account for the difference between assembly size and flow cytometric genome size estimates. We discuss the prospects of gap closure with ultra-long sequence reads.

Accurate recombination estimation from pooled genotyping and sequencing: a case study on barley

Sexual reproduction involves meiotic recombination and the creation of crossing over between homologous chromosomes, which leads to new allele combinations. We present a new approach that uses the allele frequency differences and the physical distance of neighboring polymorphisms to estimate the recombination rate from pool genotyping or sequencing. This allows a considerable cost reduction compared to conventional mapping based on genotyping or sequencing data of single individuals. We evaluated the approach based on computer simulations at various genotyping depths and population sizes as well as applied it to experimental data of 45 barley populations, comprising 4182 RIL. High correlations between the recombination rates from this new pool genetic mapping approach and conventional mapping in simulated and experimental barley populations were observed. The proposed method therefore provides a reliable genetic map position and recombination rate estimation in defined genomic windows.

G4Boost: a machine learning-based tool for quadruplex identification and stability prediction

Background

G-quadruplexes (G4s), formed within guanine-rich nucleic acids, are secondary structures involved in important biological processes. Although every G4 motif has the potential to form a stable G4 structure, not every G4 motif would, and accurate energy-based methods are needed to assess their structural stability. Here, we present a decision tree-based prediction tool, G4Boost, to identify G4 motifs and predict their secondary structure folding probability and thermodynamic stability based on their sequences, nucleotide compositions, and estimated structural topologies.

Results

G4Boost predicted the quadruplex folding state with an accuracy greater then 93% and an F1-score of 0.96, and the folding energy with an RMSE of 4.28 and R² of 0.95 only by the means of sequence intrinsic feature. G4Boost was successfully applied and validated to predict the stability of experimentally-determined G4 structures, including for plants and humans.

Conclusion

G4Boost outperformed the three machine-learning based prediction tools, DeepG4, Quadron, and G4RNA Screener, in terms of both accuracy and F1-score, and can be highly useful for G4 prediction to understand gene regulation across species including plants and humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04782-z.