BaRTv1.0: an improved barley reference transcript dataset to determine accurate changes in the barley transcriptome using RNA-seq

Advancing Sustainable Malting Practices: Aquaporins as Potential Breeding Targets for Improved Water Uptake during Controlled Germination of Barley (Hordeum vulgare L.)

The conversion of raw barley (Hordeum vulgare L.) to malt requires a process of controlled germination, where the grain is submerged in water to raise the moisture content to >40%. The transmembrane proteins, aquaporins, influence water uptake during the initial stage of controlled germination, yet little is known of their involvement in malting. With the current focus on sustainability, understanding the mechanisms of water uptake and usage during the initial stages of malting has become vital in improving efficient malting practices. In this study, we used quantitative proteomics analysis of two malting barley genotypes demonstrating differing water-uptake phenotypes in the initial stages of malting. Our study quantified 19 transmembrane proteins from nine families, including seven distinct aquaporin isoforms, including the plasma intrinsic proteins (PIPs) PIP1;1, PIP2;1, and PIP2;4 and the tonoplast intrinsic proteins (TIPs) TIP1;1, TIP2;3, TIP3;1, and TIP3;2. Our findings suggest that the presence of TIP1;1, TIP3;1, and TIP3;2 in the mature barley grain proteome is essential for facilitating water uptake, influencing cell turgor and the formation of large central lytic vacuoles aiding storage reserve hydrolysis and endosperm modification efficiency. This study proposes that TIP3s mediate water uptake in malting barley grain, offering potential breeding targets for improving sustainable malting practices.

The Hessian fly resistance gene HvRHF1 is localized in an NBS-LRR gene cluster in barley

Hessian fly (Mayetiola destructor Say) is a significant pest in cereal crops, causing substantial yield losses worldwide. While host resistance is the most efficient method for pest control, research on genetic characterization of Hessian fly resistance in barley (Hordeum vulgare L.) has been limited, and the underlying resistance mechanism remains largely unknown. In this study, we conducted fine mapping of a crucial Hessian fly resistance locus, known as HvRHF1, using a biparental population. Assisted with genetic markers and robust phenotyping assay, we pinpointed the HvRHF1 gene to an ~ 82 kb region on chromosome 4H. Gene prediction and annotation revealed that the HvRHF1 locus comprises three complete NBS-LRR genes, which are characteristic of disease resistance genes. As a result, our study not only provides valuable resources for resistance in barley and genetic tools for breeding, but also identifies candidate genes that lay the foundation for cloning HvRHF1. This endeavor will significantly contribute to our understanding of the molecular mechanisms underlying cereal resistance to Hessian fly.

Unveiling unique alternative splicing responses to low temperature in Zoysia japonica through ZjRTD1.0, a high-quality reference transcript dataset

Inadequate reference databases in RNA-seq analysis can hinder data utilization and interpretation. In this study, we have successfully constructed a high-quality reference transcript dataset, ZjRTD1.0, for Zoysia japonica, a widely-used turfgrass with exceptional tolerance to various abiotic stress, including low temperatures and salinity. This dataset comprises 113,089 transcripts from 57,143 genes. BUSCO analysis demonstrates exceptional completeness (92.4%) in ZjRTD1.0, with reduced proportions of fragmented (3.3%) and missing (4.3%) orthologs compared to prior datasets. ZjRTD1.0 enables more precise analyses, including transcript quantification and alternative splicing assessments using public datasets, which identified a substantial number of differentially expressed transcripts (DETs) and differential alternative splicing (DAS) events, leading to several novel findings on Z. japonica's responses to abiotic stresses. First, spliceosome gene expression influenced alternative splicing significantly under abiotic stress, with a greater impact observed during low-temperature stress. Then, a significant positive correlation was found between the number of differentially expressed genes (DEGs) encoding protein kinases and the frequency of DAS events, suggesting the role of protein phosphorylation in regulating alternative splicing. Additionally, our results suggest possible involvement of serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs) in generating inclusion/exclusion isoforms under low-temperature stress. Furthermore, our investigation revealed a significantly enhanced overlap between DEGs and differentially alternatively spliced genes (DASGs) in response to low-temperature stress, suggesting a unique co-regulatory mechanism governing transcription and splicing in the context of low-temperature response. In conclusion, we have proven that ZjRTD1.0 will serve as a reliable and useful resource for future transcriptomic analyses in Z. japonica.

Plastid terminal oxidase is required for chloroplast biogenesis in barley

Chloroplast biogenesis is critical for crop biomass and economic yield. However, chloroplast development is a very complicated process coordinated by cross-communication between the nucleus and plastids, and the underlying mechanisms have not been fully revealed. To explore the regulatory machinery for chloroplast biogenesis, we conducted map-based cloning of the Grandpa 1 (Gpa1) gene regulating chloroplast development in barley. The spontaneous mutation gpa1.a caused a variegation phenotype of the leaf, dwarfed growth, reduced grain yield, and increased tiller number. Genetic mapping anchored the Gpa1 gene onto 2H within a gene cluster functionally related to photosynthesis or chloroplast differentiation. One gene (HORVU.MOREX.r3.2HG0213170) in the delimited region encodes a putative plastid terminal oxidase (PTOX) in thylakoid membranes, which is homologous to IMMUTANS (IM) of Arabidopsis. The IM gene is required for chloroplast biogenesis and maintenance of functional thylakoids in Arabidopsis. Using CRISPR technology and gene transformation, we functionally validated that the PTOX-encoding gene, HORVU.MOREX.r3.2HG0213170, is the causal gene of Gpa1. Gene expression and chemical analysis revealed that the carotenoid biosynthesis pathway is suppressed by the gpa1 mutation, rendering mutants vulnerable to photobleaching. Our results showed that the overtillering associated with the gpa1 mutation was caused by the lower accumulation of carotenoid-derived strigolactones (SLs) in the mutant. The cloning of Gpa1 not only improves our understanding of the molecular mechanisms underlying chloroplast biosynthesis but also indicates that the PTOX activity is conserved between monocots and dicots for the establishment of the photosynthesis factory.

Proteomic exploration reveals a metabolic rerouting due to low oxygen during controlled germination of malting barley (Hordeum vulgare L.)

Barley (Hordeum vulgare L.) is used in malt production for brewing applications. Barley malting involves a process of controlled germination that modifies the grain by activating enzymes to solubilize starch and proteins for brewing. Initially, the grain is submerged in water to raise grain moisture, requiring large volumes of water. Achieving grain modification at reduced moisture levels can contribute to the sustainability of malting practices. This study combined proteomics, bioinformatics, and biochemical phenotypic analysis of two malting barley genotypes with observed differences in water uptake and modification efficiency. We sought to reveal the molecular mechanisms at play during controlled germination and explore the roles of protein groups at 24 h intervals across the first 72 h. Overall, 3,485 protein groups were identified with 793 significant differentially abundant (DAP) within and between genotypes, involved in various biological processes, including protein synthesis, carbohydrate metabolism, and hydrolysis. Functional integration into metabolic pathways, such as glycolysis, pyruvate, starch and sucrose metabolism, revealed a metabolic rerouting due to low oxygen enforced by submergence during controlled germination. This SWATH-MS study provides a comprehensive proteome reference, delivering new insights into the molecular mechanisms underlying the impacts of low oxygen during controlled germination. It is concluded that continued efficient modification of malting barley subjected to submergence is largely due to the capacity to reroute energy to maintain vital processes, particularly protein synthesis.

Natural variation in HvAT10 underlies grain cell wall-esterified phenolic acid content in cultivated barley

The phenolic acids, ferulic acid and p-coumaric acid, are components of plant cell walls in grasses, including many of our major food crops. They have important health-promoting properties in grain, and influence the digestibility of biomass for industrial processing and livestock feed. Both phenolic acids are assumed to be critical to cell wall integrity and ferulic acid, at least, is important for cross-linking cell wall components, but the role of p-coumaric acid is unclear. Here we identify alleles of a BAHD p-coumaroyl arabinoxylan transferase, HvAT10, as responsible for the natural variation in cell wall-esterified phenolic acids in whole grain within a cultivated two-row spring barley panel. We show that HvAT10 is rendered non-functional by a premature stop codon mutation in half of the genotypes in our mapping panel. This results in a dramatic reduction in grain cell wall-esterifed p-coumaric acid, a moderate rise in ferulic acid, and a clear increase in the ferulic acid to p-coumaric acid ratio. The mutation is virtually absent in wild and landrace germplasm suggesting an important function for grain arabinoxylan p-coumaroylation pre-domestication that is dispensable in modern agriculture. Intriguingly, we detected detrimental impacts of the mutated locus on grain quality traits where it was associated with smaller grain and poorer malting properties. HvAT10 could be a focus for improving grain quality for malting or phenolic acid content in wholegrain foods.

Genetic and physical localization of a major susceptibility gene to Pyrenophora teres f. maculata in barley

KEY MESSAGE

Genetic characterization of a major spot form net blotch susceptibility locus to using linkage mapping to identify a candidate gene and user-friendly markers in barley. Spot form net blotch (SFNB), caused by the necrotrophic fungal pathogen Pyrenophora teres f. maculata (Ptm), is an economically important foliar diseases in barley. Although various resistance loci have been identified, breeding for SFNB-resistant varieties has been hampered due to the complex virulence profile of Ptm populations. One resistance locus in the host may be effective against one specific isolate, but it may confer susceptibility to other isolates. A major susceptibility QTL on chromosome 7H, named Sptm1, was consistently identified in many studies. In the present study, we conduct fine mapping to localize Sptm1 with high resolution. A segregating population was developed from selected F₂ progenies of the cross Tradition (S) × PI 67381 (R), in which the disease phenotype was determined by the Sptm1 locus alone. Disease phenotypes of critical recombinants were confirmed in the following two consecutive generations. Genetic mapping anchored the Sptm1 gene to an ⁓400 kb region on chromosome 7H. Gene prediction and annotation identified six protein-coding genes in the delimited Sptm1 region, and the gene encoding a putative cold-responsive protein kinase was selected as a strong candidate. Therefore, providing fine localization and candidate of Sptm1 for functional validation, our study will facilitate the understanding of susceptibility mechanism underlying the barley-Ptm interaction and offers a potential target for gene editing to develop valuable materials with broad-spectrum resistance to SFNB.

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.

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.

High-quality reference transcriptome construction improves RNA-seq quantification in Oryza sativa indica

The Reference Transcriptomic Dataset (RTD) is an accurate and comprehensive collection of transcripts originating from a given organism. It holds the key to precise transcript quantification and downstream analysis of differential expressions and regulations. Currently, transcriptome annotations for most crop plants are far from complete. For example, Oryza sativa indica (O. sativa indica) is reported to have 40,759 transcripts in the Ensembl database without alternative transcript isoforms and alternative splicing (AS) events. To generate a high-quality RTD, we conducted RNA sequencing of rice leaf samples collected at various time points during Rhizoctonia solani infection. The obtained reads were analyzed by adopting the recently developed computational analysis pipeline to assemble the RTD with increased transcript and AS diversity for O. sativa indica (IndicaRTD). After stringent quality filtering, the newly constructed transcriptome annotation was comprised of 122,968 non-redundant transcripts from 53,695 genes. This study identified many novel transcripts compared to Ensembl deposited data that are important for regulating molecular and physiological processes in the plant system. Currently, the assembled IndicaRTD must allow fast quantification of transcript and gene expression with high precision.

The evolutionary patterns of barley pericentromeric chromosome regions, as shaped by linkage disequilibrium and domestication

The distribution of recombination events along large cereal chromosomes is uneven and is generally restricted to gene-rich telomeric ends. To understand how the lack of recombination affects diversity in the large pericentromeric regions, we analysed deep exome capture data from a final panel of 815 Hordeum vulgare (barley) cultivars, landraces and wild barleys, sampled from across their eco-geographical ranges. We defined and compared variant data across the pericentromeric and non-pericentromeric regions, observing a clear partitioning of diversity both within and between chromosomes and germplasm groups. Dramatically reduced diversity was found in the pericentromeres of both cultivars and landraces when compared with wild barley. We observed a mixture of completely and partially differentiated single-nucleotide polymorphisms (SNPs) between domesticated and wild gene pools, suggesting that domesticated gene pools were derived from multiple wild ancestors. Patterns of genome-wide linkage disequilibrium, haplotype block size and number, and variant frequency within blocks showed clear contrasts among individual chromosomes and between cultivars and wild barleys. Although most cultivar chromosomes shared a single major pericentromeric haplotype, chromosome 7H clearly differentiated the two-row and six-row types associated with different geographical origins. Within the pericentromeric regions we identified 22 387 non-synonymous SNPs, 92 of which were fixed for alternative alleles in cultivar versus wild accessions. Surprisingly, only 29 SNPs found exclusively in the cultivars were predicted to be 'highly deleterious'. Overall, our data reveal an unconventional pericentromeric genetic landscape among distinct barley gene pools, with different evolutionary processes driving domestication and diversification.

Genotyping by Sequencing Advancements in Barley

Barley is considered an ideal crop to study cereal genetics due to its close relationship with wheat and diploid ancestral genome. It plays a crucial role in reducing risks to global food security posed by climate change. Genetic variations in the traits of interest in crops are vital for their improvement. DNA markers have been widely used to estimate these variations in populations. With the advancements in next-generation sequencing, breeders could access different types of genetic variations within different lines, with single-nucleotide polymorphisms (SNPs) being the most common type. However, genotyping barley with whole genome sequencing (WGS) is challenged by the higher cost and computational demand caused by the large genome size (5.5GB) and a high proportion of repetitive sequences (80%). Genotyping-by-sequencing (GBS) protocols based on restriction enzymes and target enrichment allow a cost-effective SNP discovery by reducing the genome complexity. In general, GBS has opened up new horizons for plant breeding and genetics. Though considered a reliable alternative to WGS, GBS also presents various computational difficulties, but GBS-specific pipelines are designed to overcome these challenges. Moreover, a robust design for GBS can facilitate the imputation to the WGS level of crops with high linkage disequilibrium. The complete exploitation of GBS advancements will pave the way to a better understanding of crop genetics and offer opportunities for the successful improvement of barley and its close relatives.

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.

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.

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.

Reciprocal epigenetic remodeling controls testicular cancer hypersensitivity to hypomethylating agents and chemotherapy

Testicular germ cell tumors (TGCTs) are aggressive but sensitive to cisplatin-based chemotherapy. Alternative therapies are needed for tumors refractory to cisplatin with hypomethylating agents providing one possibility. The mechanisms of cisplatin hypersensitivity and resistance in TGCTs remain poorly understood. Recently, it has been shown that TGCTs, even those resistant to cisplatin, are hypersensitive to very low doses of hypomethylating agents including 5-aza deoxy-cytosine (5-aza) and guadecitabine. We undertook a pharmacogenomic approach in order to better understand mechanisms of TGCT hypomethylating agent hypersensitivity by generating a panel of acquired 5-aza-resistant TGCT cells and contrasting these to previously generated acquired isogenic cisplatin-resistant cells from the same parent. Interestingly, there was a reciprocal relationship between cisplatin and 5-aza sensitivity, with cisplatin resistance associated with increased sensitivity to 5-aza and 5-aza resistance associated with increased sensitivity to cisplatin. Unbiased transcriptome analysis revealed 5-aza-resistant cells strongly downregulated polycomb target gene expression, the exact opposite of the finding for cisplatin-resistant cells, which upregulated polycomb target genes. This was associated with a dramatic increase in H3K27me3 and decrease in DNMT3B levels in 5-aza-resistant cells, the exact opposite changes seen in cisplatin-resistant cells. Evidence is presented that reciprocal regulation of polycomb and DNMT3B may be initiated by changes in DNMT3B levels as DNMT3B knockdown alone in parental cells resulted in increased expression of H3K27me3, EZH2, and BMI1, conferred 5-aza resistance and cisplatin sensitization, and mediated genome-wide repression of polycomb target gene expression. Finally, genome-wide analysis revealed that 5-aza-resistant, cisplatin-resistant, and DNMT3B-knockdown cells alter the expression of a common set of polycomb target genes. This study highlights that reciprocal epigenetic changes mediated by DNMT3B and polycomb may be a key driver of the unique cisplatin and 5-aza hypersensitivity of TGCTs and suggests that distinct epigenetic vulnerabilities may exist for pharmacological targeting of TGCTs.

Transcriptome Analyses of Barley Roots Inoculated with Novel Paenibacillus sp. and Erwinia gerundensis Strains Reveal Beneficial Early-Stage Plant-Bacteria Interactions

Plant growth-promoting bacteria can improve host plant traits including nutrient uptake and metabolism and tolerance to biotic and abiotic stresses. Understanding the molecular basis of plant–bacteria interactions using dual RNA-seq analyses provides key knowledge of both host and bacteria simultaneously, leading to future enhancements of beneficial interactions. In this study, dual RNA-seq analyses were performed to provide insights into the early-stage interactions between barley seedlings and three novel bacterial strains (two Paenibacillus sp. strains and one Erwinia gerundensis strain) isolated from the perennial ryegrass seed microbiome. Differentially expressed bacterial and barley genes/transcripts involved in plant–bacteria interactions were identified, with varying species- and strain-specific responses. Overall, transcriptome profiles suggested that all three strains improved stress response, signal transduction, and nutrient uptake and metabolism of barley seedlings. Results also suggested potential improvements in seedling root growth via repressing ethylene biosynthesis in roots. Bacterial secondary metabolite gene clusters producing compounds that are potentially associated with interactions with the barley endophytic microbiome and associated with stress tolerance of plants under nutrient limiting conditions were also identified. The results of this study provided the molecular basis of plant growth-promoting activities of three novel bacterial strains in barley, laid a solid foundation for the future development of these three bacterial strains as biofertilisers, and identified key differences between bacterial strains of the same species in their responses to plants.

Proteome Analysis of Condensed Barley Mitotic Chromosomes

Proteins play a major role in the three-dimensional organization of nuclear genome and its function. While histones arrange DNA into a nucleosome fiber, other proteins contribute to higher-order chromatin structures in interphase nuclei, and mitotic/meiotic chromosomes. Despite the key role of proteins in maintaining genome integrity and transferring hereditary information to daughter cells and progenies, the knowledge about their function remains fragmentary. This is particularly true for the proteins of condensed chromosomes and, in particular, chromosomes of plants. Here, we purified barley mitotic metaphase chromosomes by a flow cytometric sorting and characterized their proteins. Peptides from tryptic protein digests were fractionated either on a cation exchanger or reversed-phase microgradient system before liquid chromatography coupled to tandem mass spectrometry. Chromosomal proteins comprising almost 900 identifications were classified based on a combination of software prediction, available database localization information, sequence homology, and domain representation. A biological context evaluation indicated the presence of several groups of abundant proteins including histones, topoisomerase 2, POLYMERASE 2, condensin subunits, and many proteins with chromatin-related functions. Proteins involved in processes related to DNA replication, transcription, and repair as well as nucleolar proteins were found. We have experimentally validated the presence of FIBRILLARIN 1, one of the nucleolar proteins, on metaphase chromosomes, suggesting that plant chromosomes are coated with proteins during mitosis, similar to those of human and animals. These results improve significantly the knowledge of plant chromosomal proteins and provide a basis for their functional characterization and comparative phylogenetic analyses.

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.

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.

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.

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.

FINDER: an automated software package to annotate eukaryotic genes from RNA-Seq data and associated protein sequences

BACKGROUND

Gene annotation in eukaryotes is a non-trivial task that requires meticulous analysis of accumulated transcript data. Challenges include transcriptionally active regions of the genome that contain overlapping genes, genes that produce numerous transcripts, transposable elements and numerous diverse sequence repeats. Currently available gene annotation software applications depend on pre-constructed full-length gene sequence assemblies which are not guaranteed to be error-free. The origins of these sequences are often uncertain, making it difficult to identify and rectify errors in them. This hinders the creation of an accurate and holistic representation of the transcriptomic landscape across multiple tissue types and experimental conditions. Therefore, to gauge the extent of diversity in gene structures, a comprehensive analysis of genome-wide expression data is imperative.

RESULTS

We present FINDER, a fully automated computational tool that optimizes the entire process of annotating genes and transcript structures. Unlike current state-of-the-art pipelines, FINDER automates the RNA-Seq pre-processing step by working directly with raw sequence reads and optimizes gene prediction from BRAKER2 by supplementing these reads with associated proteins. The FINDER pipeline (1) reports transcripts and recognizes genes that are expressed under specific conditions, (2) generates all possible alternatively spliced transcripts from expressed RNA-Seq data, (3) analyzes read coverage patterns to modify existing transcript models and create new ones, and (4) scores genes as high- or low-confidence based on the available evidence across multiple datasets. We demonstrate the ability of FINDER to automatically annotate a diverse pool of genomes from eight species.

CONCLUSIONS

FINDER takes a completely automated approach to annotate genes directly from raw expression data. It is capable of processing eukaryotic genomes of all sizes and requires no manual supervision-ideal for bench researchers with limited experience in handling computational tools.

EORNA, a barley gene and transcript abundance database

A high-quality, barley gene reference transcript dataset (BaRTv1.0), was used to quantify gene and transcript abundances from 22 RNA-seq experiments, covering 843 separate samples. Using the abundance data we developed a Barley Expression Database (EORNA*) to underpin a visualisation tool that displays comparative gene and transcript abundance data on demand as transcripts per million (TPM) across all samples and all the genes. EORNA provides gene and transcript models for all of the transcripts contained in BaRTV1.0, and these can be conveniently identified through either BaRT or HORVU gene names, or by direct BLAST of query sequences. Browsing the quantification data reveals cultivar, tissue and condition specific gene expression and shows changes in the proportions of individual transcripts that have arisen via alternative splicing. TPM values can be easily extracted to allow users to determine the statistical significance of observed transcript abundance variation among samples or perform meta analyses on multiple RNA-seq experiments. * Eòrna is the Scottish Gaelic word for Barley.

EORNA, a barley gene and transcript abundance database

A high-quality, barley gene reference transcript dataset (BaRTv1.0), was used to quantify gene and transcript abundances from 22 RNA-seq experiments, covering 843 separate samples. Using the abundance data we developed a Barley Expression Database (EORNA*) to underpin a visualisation tool that displays comparative gene and transcript abundance data on demand as transcripts per million (TPM) across all samples and all the genes. EORNA provides gene and transcript models for all of the transcripts contained in BaRTV1.0, and these can be conveniently identified through either BaRT or HORVU gene names, or by direct BLAST of query sequences. Browsing the quantification data reveals cultivar, tissue and condition specific gene expression and shows changes in the proportions of individual transcripts that have arisen via alternative splicing. TPM values can be easily extracted to allow users to determine the statistical significance of observed transcript abundance variation among samples or perform meta analyses on multiple RNA-seq experiments. * Eòrna is the Scottish Gaelic word for Barley.

Barley Anther and Meiocyte Transcriptome Dynamics in Meiotic Prophase I

In flowering plants, successful germinal cell development and meiotic recombination depend upon a combination of environmental and genetic factors. To gain insights into this specialized reproductive development program we used short- and long-read RNA-sequencing (RNA-seq) to study the temporal dynamics of transcript abundance in immuno-cytologically staged barley (Hordeum vulgare) anthers and meiocytes. We show that the most significant transcriptional changes in anthers occur at the transition from pre-meiosis to leptotene-zygotene, which is followed by increasingly stable transcript abundance throughout prophase I into metaphase I-tetrad. Our analysis reveals that the pre-meiotic anthers are enriched in long non-coding RNAs (lncRNAs) and that entry to meiosis is characterized by their robust and significant down regulation. Intriguingly, only 24% of a collection of putative meiotic gene orthologs showed differential transcript abundance in at least one stage or tissue comparison. Argonautes, E3 ubiquitin ligases, and lys48 specific de-ubiquitinating enzymes were enriched in prophase I meiocyte samples. These developmental, time-resolved transcriptomes demonstrate remarkable stability in transcript abundance in meiocytes throughout prophase I after the initial and substantial reprogramming at meiosis entry and the complexity of the regulatory networks involved in early meiotic processes.

Re-examination of the APETALA2/Ethylene-Responsive Factor Gene Family in Barley (Hordeum vulgare L.) Indicates a Role in the Regulation of Starch Synthesis

The APETALA2/Ethylene-Responsive factor (AP2/ERF) gene family is a large plant-specific transcription factor family, which plays important roles in regulating plant growth and development. A role in starch synthesis is among the multiple functions of this family of transcription factors. Barley (Hordeum vulgare L.) is one of the most important cereals for starch production. However, there are limited data on the contribution of AP2 transcription factors in barley. In this study, we used the recently published barley genome database (Morex) to identify 185 genes of the HvAP2/ERF family. Compared with previous work, we identified 64 new genes in the HvAP2/ERF gene family and corrected some previously misannotated and duplicated genes. After phylogenetic analysis, HvAP2/ERF genes were classified into four subfamilies and 18 subgroups. Expression profiling showed different patterns of spatial and temporal expression for HvAP2/ERF genes. Most of the 12 HvAP2/ERF genes analyzed using quantitative reverse transcription-polymerase chain reaction had similar expression patterns when compared with those of starch synthase genes in barley, except for HvAP2-18 and HvERF-73. HvAP2-18 is homologous to OsRSR1, which negatively regulates the synthesis of rice starch. Luciferase reporter gene, and yeast one-hybrid assays showed that HvAP2-18 bound the promoter of AGP-S and SBE1 in vitro. Thus, HvAP2-18 might be an interesting candidate gene to further explore the mechanisms involved in the regulation of starch synthesis in barley.

CRISPR/Cas9-Induced Loss-of-Function Mutation in the Barley Mitogen-Activated Protein Kinase 6 Gene Causes Abnormal Embryo Development Leading to Severely Reduced Grain Germination and Seedling Shootless Phenotype

The diverse roles of mitogen-activated protein kinases (MAPKs, MPKs) in plant development could be efficiently revealed by reverse genetic studies. In Arabidopsis, mpk6 knockout mutants complete the life cycle; however, ~40% of their embryos show defects in the development leading to abnormal phenotypes of seeds and seedlings' roots. Contrary to the Arabidopsis MPK6, the rice MPK6 (OsMPK6) is an essential gene as transfer DNA (T-DNA) insertion and CRISPR/Cas9 induced loss-of-function mutations in the OsMPK6 cause early embryo arrest. In this study, we successfully developed a viable transgenic barley line with the CRISPR/Cas9-induced heterozygous single base pair cytosine-guanine (CG) deletion [wild type (WT)/-1C] in the third exon of the HvMPK6 gene, a barley ortholog of the Arabidopsis and rice MPK6. There were no obvious macroscopic phenotype differences between the WT/-1C plants and WT plants. All the grains collected from the WT/-1C plants were of similar size and appearance. However, seedling emergence percentage (SEP) from these grains was substantially decreased in the soil in the T2 and T3 generation. The mutation analysis of the 248 emerged T2 and T3 generation plants showed that none of them was a biallelic mutant in the HvMPK6 gene, suggesting lethality of the -1C/-1C homozygous knockout mutation. In the soil, the majority of the -1C/-1C grains did not germinate and the minority of them developed into abnormal seedlings with a shootless phenotype and a reduced root system. Some of the -1C/-1C seedlings also developed one or more small chlorotic leaf blade-like structure/structures. The -1C/-1C grains contained the late-stage developed abnormal embryos with the morphologically obvious scutellum and root part of the embryonic axis but with the missing or substantially reduced shoot part of the embryonic axis. The observed embryonic abnormalities correlated well with the shootless phenotype of the seedlings and suggested that the later-stage defect is predetermined already during the embryo development. In conclusion, our results indicate that barley MPK6 is essential for the embryologically predetermined shoot formation, but not for the most aspects of the embryo and early seedling development.

TB1: from domestication gene to tool for many trades

This article comments on: Dixon LE, Pasquariello M, Boden SA. 2020. TEOSINTE BRANCHED1 regulates height and stem internode length in bread wheat. Journal of Experimental Botany 71, 4742–4750.

Barley's Second Spring as A Model Organism for Chloroplast Research

Barley (Hordeum vulgare) has been widely used as a model crop for studying molecular and physiological processes such as chloroplast development and photosynthesis. During the second half of the 20th century, mutants such as albostrians led to the discovery of the nuclear-encoded, plastid-localized RNA polymerase and the retrograde (chloroplast-to-nucleus) signalling communication pathway, while chlorina-f2 and xantha mutants helped to shed light on the chlorophyll biosynthetic pathway, on the light-harvesting proteins and on the organization of the photosynthetic apparatus. However, during the last 30 years, a large fraction of chloroplast research has switched to the more “user-friendly” model species Arabidopsis thaliana, the first plant species whose genome was sequenced and published at the end of 2000. Despite its many advantages, Arabidopsis has some important limitations compared to barley, including the lack of a real canopy and the absence of the proplastid-to-chloroplast developmental gradient across the leaf blade. These features, together with the availability of large collections of natural genetic diversity and mutant populations for barley, a complete genome assembly and protocols for genetic transformation and gene editing, have relaunched barley as an ideal model species for chloroplast research. In this review, we provide an update on the genomics tools now available for barley, and review the biotechnological strategies reported to increase photosynthesis efficiency in model species, which deserve to be validated in barley.

A Genome Assembly of the Barley 'Transformation Reference' Cultivar Golden Promise

Barley (Hordeum vulgare) is one of the most important crops worldwide and is also considered a research model for the large-genome small grain temperate cereals. Despite genomic resources improving all the time, they are limited for the cv Golden Promise, the most efficient genotype for genetic transformation. We have developed a barley cv Golden Promise reference assembly integrating Illumina paired-end reads, long mate-pair reads, Dovetail Chicago in vitro proximity ligation libraries and chromosome conformation capture sequencing (Hi-C) libraries into a contiguous reference assembly. The assembled genome of 7 chromosomes and 4.13Gb in size, has a super-scaffold N50 after Chicago libraries of 4.14Mb and contains only 2.2% gaps. Using BUSCO (benchmarking universal single copy orthologous genes) as evaluation the genome assembly contains 95.2% of complete and single copy genes from the plant database. A high-quality Golden Promise reference assembly will be useful and utilized by the whole barley research community but will prove particularly useful for CRISPR-Cas9 experiments.

A Genome Assembly of the Barley 'Transformation Reference' Cultivar Golden Promise

Barley (Hordeum vulgare) is one of the most important crops worldwide and is also considered a research model for the large-genome small grain temperate cereals. Despite genomic resources improving all the time, they are limited for the cv Golden Promise, the most efficient genotype for genetic transformation. We have developed a barley cv Golden Promise reference assembly integrating Illumina paired-end reads, long mate-pair reads, Dovetail Chicago in vitro proximity ligation libraries and chromosome conformation capture sequencing (Hi-C) libraries into a contiguous reference assembly. The assembled genome of 7 chromosomes and 4.13Gb in size, has a super-scaffold N50 after Chicago libraries of 4.14Mb and contains only 2.2% gaps. Using BUSCO (benchmarking universal single copy orthologous genes) as evaluation the genome assembly contains 95.2% of complete and single copy genes from the plant database. A high-quality Golden Promise reference assembly will be useful and utilized by the whole barley research community but will prove particularly useful for CRISPR-Cas9 experiments.

The Effect of Drought on Transcriptome and Hormonal Profiles in Barley Genotypes With Contrasting Drought Tolerance

Like many cereal crops, barley is also negatively affected by drought stress. However, due to its simple genome as well as enhanced stress resilient nature compared to rice and wheat, barley has been considered as a model to decipher drought tolerance in cereals. In the present study, transcriptomic and hormonal profiles along with several biochemical features were compared between drought-tolerant (Otis) and drought-sensitive (Baronesse) barley genotypes subjected to drought to identify molecular and biochemical differences between the genotypes. The drought-induced decrease in the leaf relative water content, net photosynthesis, and biomass accumulation was relatively low in Otis compared to Baronesse. The hormonal profiles did not reveal significant differences for majority of the compounds other than the GA20 and the cis-zeatin-o-glucoside (c-ZOG), whose levels were greatly increased in Otis compared to Baronesse under drought. The major differences that emerged from the transcriptome analysis are; (1), the overall number of differentially expressed genes was relatively low in drought-tolerant Otis compared to drought-sensitive Baronesse; (2), a wax biosynthesis gene (CER1), and NAC transcription factors were specifically induced in Otis but not in Baronesse; (3), the degree of upregulation of betaine aldehyde dehydrogenase and a homeobox transcription factor (genes with proven roles in imparting drought tolerance), was greater in Otis compared to Baronesse; (4) the extent of downregulation of gene expression profiles for proteins of the reaction center photosystem II (PSII) (D1 and D2) was low in Otis compared to Baronesse; and, (5), alternative splicing (AS) was also found to differ between the genotypes under drought. Taken together, the overall transcriptional responses were low in drought-tolerant Otis but the genes that could confer drought tolerance were either specifically induced or greatly upregulated in the tolerant genotype and these differences could be important for drought tolerance in barley.

Alternative Splicing and DNA Damage Response in Plants

Plants are exposed to a variety of abiotic and biotic stresses that may result in DNA damage. Endogenous processes - such as DNA replication, DNA recombination, respiration, or photosynthesis - are also a threat to DNA integrity. It is therefore essential to understand the strategies plants have developed for DNA damage detection, signaling, and repair. Alternative splicing (AS) is a key post-transcriptional process with a role in regulation of gene expression. Recent studies demonstrate that the majority of intron-containing genes in plants are alternatively spliced, highlighting the importance of AS in plant development and stress response. Not only does AS ensure a versatile proteome and influence the abundance and availability of proteins greatly, it has also emerged as an important player in the DNA damage response (DDR) in animals. Despite extensive studies of DDR carried out in plants, its regulation at the level of AS has not been comprehensively addressed. Here, we provide some insights into the interplay between AS and DDR in plants.