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.

Natural variation in meiotic recombination rate shapes introgression patterns in intraspecific hybrids between wild and domesticated barley.

Meiotic recombination rates vary considerably between species, populations and individuals. The genetic exchange between homologous chromosomes plays a major role in evolution by breaking linkage between advantageous and deleterious alleles in the case of introgressions. Identifying recombination rate modifiers is thus of both fundamental and practical interest to understand and utilize variation in meiotic recombination rates. We investigated recombination rate variation in a large intraspecific hybrid population (named HEB-25) derived from a cross between domesticated barley and 25 wild barley accessions. We observed quantitative variation in total crossover number with a maximum of a 1.4-fold difference between subpopulations and increased recombination rates across pericentromeric regions. The meiosis-specific α-kleisin cohesin subunit REC8 was identified as a candidate gene influencing crossover number and patterning. Furthermore, we quantified wild barley introgression patterns and revealed how local and genome-wide recombination rate variation shapes patterns of introgression. The identification of allelic variation in REC8 in combination with the observed changes in crossover patterning suggest a difference in how chromatin loops are tethered to the chromosome axis, resulting in reduced crossover suppression across pericentromeric regions. Local and genome-wide recombination rate variation is shaping patterns of introgressions and thereby directly influences the consequences of linkage drag.

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

BACKGROUND: The time required to analyse RNA-seq data varies considerably, due to discrete steps for computational assembly, quantification of gene expression and splicing analysis. Recent fast non-alignment tools such as Kallisto and Salmon overcome these problems, but these tools require a high quality, comprehensive reference transcripts dataset (RTD), which are rarely available in plants. RESULTS: A high-quality, non-redundant barley gene RTD and database (Barley Reference Transcripts - BaRTv1.0) has been generated. BaRTv1.0, was constructed from a range of tissues, cultivars and abiotic treatments and transcripts assembled and aligned to the barley cv. Morex reference genome (Mascher et al. Nature; 544: 427-433, 2017). Full-length cDNAs from the barley variety Haruna nijo (Matsumoto et al. Plant Physiol; 156: 20-28, 2011) determined transcript coverage, and high-resolution RT-PCR validated alternatively spliced (AS) transcripts of 86 genes in five different organs and tissue. These methods were used as benchmarks to select an optimal barley RTD. BaRTv1.0-Quantification of Alternatively Spliced Isoforms (QUASI) was also made to overcome inaccurate quantification due to variation in 5' and 3' UTR ends of transcripts. BaRTv1.0-QUASI was used for accurate transcript quantification of RNA-seq data of five barley organs/tissues. This analysis identified 20,972 significant differentially expressed genes, 2791 differentially alternatively spliced genes and 2768 transcripts with differential transcript usage. CONCLUSION: A high confidence barley reference transcript dataset consisting of 60,444 genes with 177,240 transcripts has been generated. Compared to current barley transcripts, BaRTv1.0 transcripts are generally longer, have less fragmentation and improved gene models that are well supported by splice junction reads. Precise transcript quantification using BaRTv1.0 allows routine analysis of gene expression and AS.

Contrasting genetic regulation of plant development in wild barley grown in two European environments revealed by nested association mapping.

Barley is cultivated more widely than the other major world crops because it adapts well to environmental constraints, such as drought, heat, and day length. To better understand the genetic control of local adaptation in barley, we studied development in the nested association mapping population HEB-25, derived from crossing 25 wild barley accessions with the cultivar 'Barke'. HEB-25 was cultivated in replicated field trials in Dundee (Scotland) and Halle (Germany), differing in regard to day length, precipitation, and temperature. Applying a genome-wide association study, we located 60 and 66 quantitative trait locus (QTL) regions regulating eight plant development traits in Dundee and Halle, respectively. A number of QTLs could be explained by known major genes such as PHOTOPERIOD 1 (Ppd-H1) and FLOWERING LOCUS T (HvFT-1) that regulate plant development. In addition, we observed that developmental traits in HEB-25 were partly controlled via genotype × environment and genotype × donor interactions, defined as location-specific and family-specific QTL effects. Our findings indicate that QTL alleles are available in the wild barley gene pool that show contrasting effects on plant development, which may be deployed to improve adaptation of cultivated barley to future environmental changes.

Enhancement of Glen Moy x Latham raspberry linkage map using GbS to further understand control of developmental processes leading to fruit ripening.

BACKGROUND: The changing climate is altering timing of key fruit ripening processes and increasing the occurrence of fruit defects. To improve our understanding of the genetic control of raspberry fruit development an enhanced genetic linkage map was developed and used to examine ripening phenotypic data. RESULTS: In this study we developed an enhanced genetic linkage map for the raspberry cvs. Glen Moy x Latham reference mapping population using genotyping by sequencing (GbS). Alignment to a newly sequenced draft reference genome of red raspberry, cultivar (cv.) Glen Moy, identified 8019 single nucleotide polymorphisms (SNPs). After stringent filtering to take account of read coverage over all the progeny individuals, association with a single chromosome, heterozygosity and marker regression mapping, 2348 high confidence SNPs were retained and integrated with an existing raspberry genetic map. The linkage map contained many more SNPs segregating in Latham than in Glen Moy. This caused difficulties in quantitative trait loci (QTL) mapping with standard software and a novel analysis based on a hidden Markov model was used to improve the mapping. QTL mapping using the newly generated dense genetic map not only corroborated previously identified genetic locations but also provided additional genetic elements controlling fruit ripening in raspberry. CONCLUSION: The high-density GbS map located the QTL peaks more precisely than in earlier studies, aligned the QTLs with Glen Moy genome scaffolds, narrowed the range of potential candidate genes to these regions that can be utilised in other populations or in gene expression studies to confirm their role and increased the repertoire of markers available to understand the genetic control of fruit ripening traits.

The Dynamics of Transcript Abundance during Cellularization of Developing Barley Endosperm.

Within the cereal grain, the endosperm and its nutrient reserves are critical for successful germination and in the context of grain utilization. The identification of molecular determinants of early endosperm development, particularly regulators of cell division and cell wall deposition, would help predict end-use properties such as yield, quality, and nutritional value. Custom microarray data have been generated using RNA isolated from developing barley grain endosperm 3 d to 8 d after pollination (DAP). Comparisons of transcript abundance over time revealed 47 gene expression modules that can be clustered into 10 broad groups. Superimposing these modules upon cytological data allowed patterns of transcript abundance to be linked with key stages of early grain development. Here, attention was focused on how the datasets could be mined to explore and define the processes of cell wall biosynthesis, remodeling, and degradation. Using a combination of spatial molecular network and gene ontology enrichment analyses, it is shown that genes involved in cell wall metabolism are found in multiple modules, but cluster into two main groups that exhibit peak expression at 3 DAP to 4 DAP and 5 DAP to 8 DAP. The presence of transcription factor genes in these modules allowed candidate genes for the control of wall metabolism during early barley grain development to be identified. The data are publicly available through a dedicated web interface (https://ics.hutton.ac.uk/barseed/), where they can be used to interrogate co- and differential expression for any other genes, groups of genes, or transcription factors expressed during early endosperm development.

Differential expression of microRNAs and potential targets under drought stress in barley.

Drought is a crucial environmental constraint limiting crop production in many parts of the world. microRNA (miRNA) based gene regulation has been shown to act in several pathways, including crop response to drought stress. Sequence based profiling and computational analysis have revealed hundreds of miRNAs and their potential targets in different plant species under various stress conditions, but few have been biologically verified. In this study, 11 candidate miRNAs were tested for their expression profiles in barley. Differences in accumulation of only four miRNAs (Ath-miR169b, Osa-miR1432, Hv-miRx5 and Hv-miR166b/c) were observed between drought-treated and well-watered barley in four genotypes. miRNA targets were predicted using degradome analysis of two, different genotypes, and genotype-specific target cleavage was observed. Inverse correlation of mature miRNA accumulation with miRNA target transcripts was also genotype dependent under drought treatment. Drought-responsive miRNAs accumulated predominantly in mesophyll tissues. Our results demonstrate genotype-specific miRNA regulation under drought stress and evidence for their role in mediating expression of target genes for abiotic stress response in barley.

Chromatin state analysis of the barley epigenome reveals a higher-order structure defined by H3K27me1 and H3K27me3 abundance.

Combinations of histones carrying different covalent modifications are a major component of epigenetic variation. We have mapped nine modified histones in the barley seedling epigenome by chromatin immunoprecipitation next-generation sequencing (ChIP-seq). The chromosomal distributions of the modifications group them into four different classes, and members of a given class also tend to coincide at the local DNA level, suggesting that global distribution patterns reflect local epigenetic environments. We used this peak sharing to define 10 chromatin states representing local epigenetic environments in the barley genome. Five states map mainly to genes and five to intergenic regions. Two genic states involving H3K36me3 are preferentially associated with constitutive gene expression, while an H3K27me3-containing genic state is associated with differentially expressed genes. The 10 states display striking distribution patterns that divide barley chromosomes into three distinct global environments. First, telomere-proximal regions contain high densities of H3K27me3 covering both genes and intergenic DNA, together with very low levels of the repressive H3K27me1 modification. Flanking these are gene-rich interior regions that are rich in active chromatin states and have greatly decreased levels of H3K27me3 and increasing amounts of H3K27me1 and H3K9me2. Lastly, H3K27me3-depleted pericentromeric regions contain gene islands with active chromatin states separated by extensive retrotransposon-rich regions that are associated with abundant H3K27me1 and H3K9me2 modifications. We propose an epigenomic framework for barley whereby intergenic H3K27me3 specifies facultative heterochromatin in the telomere-proximal regions and H3K27me1 is diagnostic for constitutive heterochromatin elsewhere in the barley genome.

An ultra-high density genetic linkage map of perennial ryegrass (Lolium perenne) using genotyping by sequencing (GBS) based on a reference shotgun genome assembly.

BACKGROUND AND AIMS: High density genetic linkage maps that are extensively anchored to assembled genome sequences of the organism in question are extremely useful in gene discovery. To facilitate this process in perennial ryegrass (Lolium perenne L.), a high density single nucleotide polymorphism (SNP)- and presence/absence variant (PAV)-based genetic linkage map has been developed in an F2 mapping population that has been used as a reference population in numerous studies. To provide a reference sequence to which to align genotyping by sequencing (GBS) reads, a shotgun assembly of one of the grandparents of the population, a tenth-generation inbred line, was created using Illumina-based sequencing. METHODS: The assembly was based on paired-end Illumina reads, scaffolded by mate pair and long jumping distance reads in the range of 3-40 kb, with >200-fold initial genome coverage. A total of 169 individuals from an F2 mapping population were used to construct PstI-based GBS libraries tagged with unique 4-9 nucleotide barcodes, resulting in 284 million reads, with approx. 1·6 million reads per individual. A bioinformatics pipeline was employed to identify both SNPs and PAVs. A core genetic map was generated using high confidence SNPs, to which lower confidence SNPs and PAVs were subsequently fitted in a straightforward binning approach. KEY RESULTS: The assembly comprises 424 750 scaffolds, covering 1·11 Gbp of the 2·5 Gbp perennial ryegrass genome, with a scaffold N50 of 25 212 bp and a contig N50 of 3790 bp. It is available for download, and access to a genome browser has been provided. Comparison of the assembly with available transcript and gene model data sets for perennial ryegrass indicates that approx. 570 Mbp of the gene-rich portion of the genome has been captured. An ultra-high density genetic linkage map with 3092 SNPs and 7260 PAVs was developed, anchoring just over 200 Mb of the reference assembly. CONCLUSIONS: The combined genetic map and assembly, combined with another recently released genome assembly, represent a significant resource for the perennial ryegrass genetics community.

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

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis.

BACKGROUND: Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis. RESULTS: We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts-twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage. CONCLUSIONS: AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

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.

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.

High-Resolution RT-PCR Analysis of Alternative Barley Transcripts.

Assembly of the barley genome and extensive use of RNA-seq has resulted in an abundance of gene expression data and the recognition of wide-scale production of alternatively spliced transcripts. Here, we describe in detail a high-resolution reverse transcription-PCR based panel (HR RT-PCR) that confirms the accuracy of alternatively spliced transcripts from RNA-seq and allows quantification of changes in the proportion of splice isoforms between different experimental conditions, time points, tissues, genotypes, ecotypes, and treatments. By validating a selection of barley genes, use of the panel gives confidence or otherwise to the genome-wide global changes in alternatively spliced transcripts reported by RNA-seq. This simple assay can readily be applied to perform detailed transcript isoform analysis for any gene in any species.

A Transcript and Metabolite Atlas of Blackcurrant Fruit Development Highlights Hormonal Regulation and Reveals the Role of Key Transcription Factors.

Blackcurrant fruit collected at six stages of development were assessed for changes in gene expression using custom whole transcriptome microarrays and for variation in metabolite content using a combination of liquid chromatography-mass spectrometry and gas chromatography-mass spectrometry. Principal components analysis demonstrated that fruit development could be clearly defined according to their transcript or metabolite profiles. During early developmental stages, metabolite profiles were dominated by amino acids and tannins, whilst transcript profiles were enriched in functions associated with cell division, anatomical structure morphogenesis and cell wall metabolism. During mid fruit development, fatty acids accumulated and transcript profiles were consistent with seed and embryo development. At the later stages, sugars and anthocyanins accumulated consistent with transcript profiles that were associated with secondary metabolism. Transcript data also indicated active signaling during later stages of fruit development. A targeted analysis of signaling networks revealed a dynamic activation and repression of almost 60 different transcripts encoding transcription factors across the course of fruit development, many of which have been demonstrated as pivotal to controlling such processes in other species. Transcripts associated with cytokinin and gibberellin were highly abundant at early fruit development, whilst those associated with ABA and ethylene tended to be more abundant at later stages. The data presented here provides an insight into fruit development in blackcurrant and provides a foundation for further work in the elucidation of the genetic basis of fruit quality.

A transcriptional reference map of defence hormone responses in potato.

Phytohormones are involved in diverse aspects of plant life including the regulation of plant growth, development and reproduction, as well as governing biotic and abiotic stress responses. We have generated a comprehensive transcriptional reference map of the early potato responses to exogenous application of the defence hormones abscisic acid, brassinolides (applied as epibrassinolide), ethylene (applied as the ethylene precursor aminocyclopropanecarboxylic acid), salicylic acid and jasmonic acid (applied as methyl jasmonate). Of the 39000 predicted genes on the microarray, a total of 2677 and 2473 genes were significantly differentially expressed at 1 h and 6 h after hormone treatment, respectively. Specific marker genes newly identified for the early hormone responses in potato include: a homeodomain 20 transcription factor (DMG400000248) for abscisic acid; a SAUR gene (DMG400016561) induced in epibrassinolide treated plants; an osmotin gene (DMG400003057) specifically enhanced by aminocyclopropanecarboxylic acid; a gene weakly similar to AtWRKY40 (DMG402007388) that was induced by salicylic acid; and a jasmonate ZIM-domain protein 1 (DMG400002930) which was specifically activated by methyl jasmonate. An online database has been set up to query the expression patterns of potato genes represented on the microarray that can also incorporate future microarray or RNAseq-based expression studies.

Construction of reference chromosome-scale pseudomolecules for potato: integrating the potato genome with genetic and physical maps.

The genome of potato, a major global food crop, was recently sequenced. The work presented here details the integration of the potato reference genome (DM) with a new sequence-tagged site marker-based linkage map and other physical and genetic maps of potato and the closely related species tomato. Primary anchoring of the DM genome assembly was accomplished by the use of a diploid segregating population, which was genotyped with several types of molecular genetic markers to construct a new ~936 cM linkage map comprising 2469 marker loci. In silico anchoring approaches used genetic and physical maps from the diploid potato genotype RH89-039-16 (RH) and tomato. This combined approach has allowed 951 superscaffolds to be ordered into pseudomolecules corresponding to the 12 potato chromosomes. These pseudomolecules represent 674 Mb (~93%) of the 723 Mb genome assembly and 37,482 (~96%) of the 39,031 predicted genes. The superscaffold order and orientation within the pseudomolecules are closely collinear with independently constructed high density linkage maps. Comparisons between marker distribution and physical location reveal regions of greater and lesser recombination, as well as regions exhibiting significant segregation distortion. The work presented here has led to a greatly improved ordering of the potato reference genome superscaffolds into chromosomal "pseudomolecules".

Hope: a construct central to living with chronic obstructive pulmonary disease.

Background. Hope plays an integral role in health and illness and may assist individuals to cope in difficult and adverse circumstances, for instance when living with an illness such as chronic obstructive pulmonary disease (COPD), which can demand continuous adaptation. Aim. This paper reports the meaning of hope in people living with COPD as described by seven participants involved in a home-based pulmonary maintenance program. Methods. Using an interpretive phenomenological approach a purposive sample of seven participants were interviewed to understand participants' experiences of the phenomena of hope. An interpretative description is provided. Results. Thematic analysis revealed a number of themes, including that hope persists despite chronic illness and the unpredictable dilemmas of living with COPD. Many benefits were found to be gained from involvement in a home-based pulmonary maintenance program, including increasing exercise capacity, hope and wellbeing. Conclusion. Despite the limitations imposed by living with COPD participants revealed a determination to live as normal a life as possible. The pulmonary maintenance program was pivotal in assisting participants to improve exercise capacity, hope and wellbeing. Nurses may have a role to play in helping people with COPD maintain or regain hope.