High-throughput genotyping arrays continue to be an attractive, cost-effective alternative to sequencing based approaches. We have developed a new 50k Illumina Infinium iSelect genotyping array for barley, a cereal crop species of major international importance. The majority of SNPs on the array have been extracted from variants called in exome capture data of a wide range of European barley germplasm. We used the recently published barley pseudomolecule assembly to map the exome capture data, which allowed us to generate markers with accurate physical positions and detailed gene annotation. Markers from an existing and widely used barley 9k Infinium iSelect array were carried over onto the 50k chip for backward compatibility. The array design featured 49,267 SNP markers that converted into 44,040 working assays, of which 43,461 were scorable in GenomeStudio. Of the working assays, 6,251 are from the 9k iSelect platform. We validated the SNPs by comparing the genotype calls from the new array to legacy datasets. Rates of agreement averaged 98.1 and 93.9% respectively for the legacy 9k iSelect SNP set (Comadran et al., 2012) and the exome capture SNPs. To test the utility of the 50k chip for genetic mapping, we genotyped a segregating population derived from a Golden Promise × Morex cross (Liu et al., 2014) and mapped over 14,000 SNPs to genetic positions which showed a near exact correspondence to their known physical positions. Manual adjustment of the cluster files used by the interpreting software for genotype scoring improved results substantially, but migration of cluster files between sites led to a deterioration of results, suggesting that local adjustment of cluster files is required on a site-per-site basis. Information relating to the markers on the chip is available online at https://ics.hutton.ac.uk/50k.

The barley inflorescence (spike) comprises a multi-noded central stalk (rachis) with tri-partite clusters of uni-floretted spikelets attached alternately along its length. Relative fertility of lateral spikelets within each cluster leads to spikes with two or six rows of grain, or an intermediate morphology. Understanding the mechanisms controlling this key developmental step could provide novel solutions to enhanced grain yield. Classical genetic studies identified five major SIX-ROWED SPIKE (VRS) genes, with four now known to encode transcription factors. Here we identify and characterise the remaining major VRS gene, VRS3, as encoding a putative Jumonji C-type H3K9me2/me3 demethylase, a regulator of chromatin state. Exploring the expression network modulated by VRS3 reveals specific interactions, both with other VRS genes and genes involved in stress, hormone and sugar metabolism. We show that combining a vrs3 mutant allele with natural six-rowed alleles of VRS1 and VRS5 leads to increased lateral grain size and greater grain uniformity.

BackgroundThe 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.ResultsA 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.ConclusionA 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.

A high-quality, non-redundant barley gene reference transcript dataset and database (Barley Reference Transcripts -BaRTv1.0) has been generated. BaRTv1.0, was constructed from 11 RNA-seq experimental datasets consisting of 808 samples from a range of tissues, cultivars and abiotic treatments. Transcripts were assembled and aligned to the barley cv. Morex reference genome (Mascher et al., 2017). Full-length cDNAs from the barley variety Haruna nijo (Matsumoto et al., 2011) were used to determine transcript coverage, and high-resolutionRT-PCR to validate alternatively spliced (AS) transcripts of 86 genes in five different organs and tissues. These methods were used as benchmarks to select optimal parameters for transcript assembly with StringTie. Further filtering to remove fragmented, poorly supported and lowly expressed transcripts generated BaRTv1.0, which consists of 60,444 genes with 177,240 transcripts. Overall, BaRTv1.0 contains higher confidence transcripts when compared to the current barley (HORVU) transcriptome; they are generally longer, have less fragmentation and improved gene models that are well supported by splice junction reads.We also generated a version of the RTD specifically for quantification of alternative spliced isoforms (BaRTv1.0-Quantification of Alternatively Spliced Isoforms -QUASI) for accurate expression and AS analysis. BaRTv1.0-QUASI overcomes inaccurate quantification due to variation in 5' and 3' UTR ends of transcripts and was used for accurate transcript quantification of RNA-seq data of five barley organs and tissues. Precise transcript quantification allows routine analysis of AS and we identified 20,972 significant differentially expressed genes, 2,791 differentially alternatively spliced genes and 2,768 transcripts with differential transcript usage.