Search

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.

Potato is the third most important food crop in the world. Diverse pathogens threaten sustainable crop production but can be controlled, in many cases, through the deployment of disease resistance genes belonging to the family of nucleotide-binding, leucine-rich-repeat (NLR) genes.To identify functional NLRs in established varieties, we have successfully established SMRT-AgRenSeq in tetraploid potatoes and have further enhanced the methodology by including dRenSeq in an approach that we term SMRT-AgRenSeq-d. The inclusion of dRenSeq enables the filtering of candidates after the association analysis by establishing a presence/absence matrix across resistant and susceptible potatoes that is translated into an F1 score. Using a SMRT-RenSeq based sequence representation of the NLRome from the cultivar Innovator, SMRT-AgRenSeq-d analyses reliably identified the late blight resistance benchmark genesR1, R2-like, R3aandR3bin a panel of 117 varieties with variable phenotype penetrations. All benchmark genes were identified with an F1 score of 1 which indicates absolute linkage in the panel.When applied to the elusive nematode disease resistance geneGpa5that controls the Potato Cyst Nematode (PCN) speciesGlobodera pallida(pathotypes Pa2/3), SMRT-AgRenSeq-d identified nine strong candidates. These map to the previously established position on potato chromosome 5 and are potential homologs of the late blight resistance geneR1.Assuming that NLRs are involved in controlling many types of resistances, SMRT-AgRenSeq-d can readily be applied to diverse crops and pathogen systems. In potato, SMRT-AgRenSeq-d lends itself, for example, to further study the elusive PCN resistancesH1orH3for which phenotypic data exist.

Diagnostic markers for Rrs1Rh4 have been identified by testing for associations between SNPs within the Rrs1 interval in 150 barley genotypes and their resistance to Rhynchosporium commune isolates recognised by lines containing Rrs1. Rhynchosporium or barley scald, caused by the destructive fungal pathogen Rhynchosporium commune, is one of the most economically important diseases of barley in the world. Barley landraces from Syria and Jordan demonstrated high resistance to rhynchosporium in the field. Genotyping of a wide range of barley cultivars and landraces, including known sources of different Rrs1 genes/alleles, across the Rrs1 interval, followed by association analysis of this genotypic data with resistance phenotypes to R. commune isolates recognised by Rrs1, allowed the identification of diagnostic markers for Rrs1Rh4. These markers are specific to Rrs1Rh4 and do not detect other Rrs1 genes/alleles. The Rrs1Rh4 diagnostic markers represent a resource that can be exploited by breeders for the sustainable deployment of varietal resistance in new cultivars. Thirteen out of the 55 most resistant Syrian and Jordanian landraces were shown to contain markers specific to Rrs1Rh4. One of these lines came from Jordan, with the remaining 12 lines from different locations in Syria. One of the Syrian landraces containing Rrs1Rh4 was also shown to have Rrs2. The remaining landraces that performed well against rhynchosporium in the field are likely to contain other resistance genes and represent an important novel resource yet to be exploited by European breeders.

Accurate characterization of splice junctions as well as transcription start and end sites in reference transcriptomes allows precise quantification of transcripts from RNA-seq data and enable 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 (RTD) 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 splice junctions and transcript start and end sites. BaRTv2.18 shows increased transcript diversity and completeness compared to an earlier version, BaRTv1.0. The accuracy of transcript level quantification, splice junctions and transcript start and end sites has been validated extensively using parallel technologies and analysis, including high resolution RT PCR 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.

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.

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., Measurement(s) gene expression Technology Type(s) transcription profiling assay Factor Type(s) Genotype • Abiotic Stress • Developmental stage • Tissue • Biotic stress Sample Characteristic - Organism Hordeum vulgare Sample Characteristic - Environment sodium chloride salt • drought • increased temperature • decreased temperature Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.13643387

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* – Expression of RNA) 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

During plant growth, sodium (Na+) in the soil is transported via the xylem from the root to the shoot. While excess Na+ is toxic to most plants, non-toxic concentrations have been shown to improve crop yields under certain conditions, such as when soil K+ is low. We quantified grain Na+ across a barley genome-wide association study panel grown under non-saline conditions and identified variants of a Class 1 HIGH-AFFINITY-POTASSIUM-TRANSPORTER (HvHKT1;5)-encoding gene responsible for Na+ content variation under these conditions. A leucine to proline substitution at position 189 (L189P) in HvHKT1;5 disturbs its characteristic plasma membrane localisation and disrupts Na+ transport. Under low and moderate soil Na+, genotypes containing HvHKT1:5P189 accumulate high concentrations of Na+ but exhibit no evidence of toxicity. As the frequency of HvHKT1:5P189 increases significantly in cultivated European germplasm, we cautiously speculate that this non-functional variant may enhance yield potential in non-saline environments, possibly by offsetting limitations of low available K+.

We quantified grain sodium (Na+) content across a barley GWAS panel grown under optimal conditions. We identified a strong association with a region containing two low and one high Na+ accumulating haplotypes of a Class 1 HIGH-AFFINITY POTASSIUM TRANSPORTER (HKT1;5) known to be involved in regulating plant Na+ homeostasis. The haplotypes exhibited an average 1.8-fold difference in grain Na+ content. We show that an L189P substitution disrupts Na+ transport in the high Na+ lines, disturbs the plasma membrane localisation typical of HKT1;5 and induces a conformational change in the protein predicted to compromise function. Under NaCl stress, lines containing P189 accumulate high levels of Na+, but show no significant difference in biomass. P189 increases in frequency from wild-species to elite cultivars leading us to speculate that the compromised haplotype is undergoing directional selection possibly due to the value of Na+ as a functional nutrient in non-saline environments.

During plant growth, sodium (Na+) in the soil is transported via the xylem from the root to the shoot. While excess Na+ is toxic to most plants, non-toxic concentrations have been shown to improve crop yields under certain conditions, such as when soil K+ is low. We quantified grain Na+ across a barley genome-wide association study panel grown under non-saline conditions and identified variants of a Class 1 HIGH-AFFINITY-POTASSIUM-TRANSPORTER (HvHKT1;5)-encoding gene responsible for Na+ content variation under these conditions. A leucine to proline substitution at position 189 (L189P) in HvHKT1;5 disturbs its characteristic plasma membrane localisation and disrupts Na+ transport. Under low and moderate soil Na+, genotypes containing HvHKT1:5P189 accumulate high concentrations of Na+ but exhibit no evidence of toxicity. As the frequency of HvHKT1:5P189 increases significantly in cultivated European germplasm, we cautiously speculate that this non-functional variant may enhance yield potential in non-saline environments, possibly by offsetting limitations of low available K+., Kelly Houston et al. report a genome-wide association study for sodium content in barley to find genetic variants that may improve yield under low soil K + levels. They identify variants of the Na+ transporter-encoding gene HvHKT1;5 as important for sodium content variation in non-saline conditions.

Key message A broad-spectrum late blight disease-resistance gene from Solanum verrucosum has been mapped to potato chromosome 9. The gene is distinct from previously identified-resistance genes. Abstract We have identified and characterised a broad-spectrum resistance to Phytophthora infestans from the wild Mexican species Solanum verrucosum. Diagnostic resistance gene enrichment (dRenSeq) revealed that the resistance is not conferred by previously identified nucleotide-binding, leucine-rich repeat genes. Utilising the sequenced potato genome as a reference, two complementary enrichment strategies that target resistance genes (RenSeq) and single/low-copy number genes (Generic-mapping enrichment Sequencing; GenSeq), respectively, were deployed for the rapid, SNP-based mapping of the resistance through bulked-segregant analysis. Both approaches independently positioned the resistance, referred to as Rpi-ver1, to the distal end of potato chromosome 9. Stringent post-enrichment read filtering identified a total of 64 informative SNPs that corresponded to the expected ratio for significant polymorphisms in the parents as well as the bulks. Of these, 61 SNPs are located on potato chromosome 9 and reside within 27 individual genes, which in the sequenced potato clone DM locate to positions 45.9 to 60.9 Mb. RenSeq- and GenSeq-derived SNPs within the target region were converted into allele-specific PCR-based KASP markers and further defined the position of the resistance to a 4.3 Mb interval at the bottom end of chromosome 9 between positions 52.62–56.98 Mb. Electronic supplementary material The online version of this article (10.1007/s00122-018-3078-6) contains supplementary material, which is available to authorized users.

The shoot apical and axillary meristems control shoot development, effectively influencing lateral branch and leaf formation. The barley (Hordeum vulgare) uniculm2 (cul2) mutation blocks axillary meristem development, and mutant plants lack lateral branches (tillers) that normally develop from the crown. A genetic screen for cul2 suppressors recovered two recessive alleles of ELIGULUM-A (ELI-A) that partially rescued the cul2 tillering phenotype. Mutations in ELI-A produce shorter plants with fewer tillers and disrupt the leaf blade-sheath boundary, producing liguleless leaves and reduced secondary cell wall development in stems and leaves. ELI-A is predicted to encode an unannotated protein containing an RNaseH-like domain that is conserved in land plants. ELI-A transcripts accumulate at the preligule boundary, the developing ligule, leaf margins, cells destined to develop secondary cell walls, and cells surrounding leaf vascular bundles. Recent studies have identified regulatory similarities between boundary development in leaves and lateral organs. Interestingly, we observed ELI-A transcripts at the preligule boundary, suggesting that ELI-A contributes to boundary formation between the blade and sheath. However, we did not observe ELI-A transcripts at the axillary meristem boundary in leaf axils, suggesting that ELI-A is not involved in boundary development for axillary meristem development. Our results show that ELI-A contributes to leaf and lateral branch development by acting as a boundary gene during ligule development but not during lateral branch development.

Robbie Waugh, Nils Stein, Gary Muehlbauer and colleagues report the exome sequencing of 267 landraces and wild accessions of barley from diverse regions to study adaptations to different agricultural environments. They observe correlations of days to heading and height with environment and find that variation in flowering-associated genes has strong geographical structuring. After domestication, during a process of widespread range extension, barley adapted to a broad spectrum of agricultural environments. To explore how the barley genome responded to the environmental challenges it encountered, we sequenced the exomes of a collection of 267 georeferenced landraces and wild accessions. A combination of genome-wide analyses showed that patterns of variation have been strongly shaped by geography and that variant-by-environment associations for individual genes are prominent in our data set. We observed significant correlations of days to heading (flowering) and height with seasonal temperature and dryness variables in common garden experiments, suggesting that these traits were major drivers of environmental adaptation in the sampled germplasm. A detailed analysis of known flowering-associated genes showed that many contain extensive sequence variation and that patterns of single- and multiple-gene haplotypes exhibit strong geographical structuring. This variation appears to have substantially contributed to range-wide ecogeographical adaptation, but many factors key to regional success remain unidentified.

Exome capture is a reduced representation approach that selectively captures sequence from only the gene-bearing regions of a genome. It is based on probes targeted at these regions and, compared with whole genome shotgun sequencing, leads to a significant reduction in cost and data processing effort while still providing insights into the most relevant part of a genome. An exome capture array for barley was released in 2013 and this has opened the door to numerous studies that have put this technology to good use. In this chapter we detail the laboratory protocols required for enrichment and sequencing, and provide detailed step-by-step instructions for the bioinformatics analysis of the resulting data.

Exome capture is a reduced representation approach that selectively captures sequence from only the gene-bearing regions of a genome. It is based on probes targeted at these regions and, compared with whole genome shotgun sequencing, leads to a significant reduction in cost and data processing effort while still providing insights into the most relevant part of a genome. An exome capture array for barley was released in 2013 and this has opened the door to numerous studies that have put this technology to good use. In this chapter we detail the laboratory protocols required for enrichment and sequencing, and provide detailed step-by-step instructions for the bioinformatics analysis of the resulting data.

Major resistance gene to rhynchosporium, Rrs18, maps close to the telomere on the short arm of chromosome 6H in barley. Rhynchosporium or barley scald caused by a fungal pathogen Rhynchosporium commune is one of the most destructive and economically important diseases of barley in the world. Testing of Steptoe × Morex and CIho 3515 × Alexis doubled haploid populations has revealed a large effect QTL for resistance to R. commune close to the telomere on the short arm of chromosome 6H, present in both populations. Mapping markers flanking the QTL from both populations onto the 2017 Morex genome assembly revealed a rhynchosporium resistance locus independent of Rrs13 that we named Rrs18. The causal gene was fine mapped to an interval of 660 Kb using Steptoe × Morex backcross 1 S2 and S3 lines with molecular markers developed from Steptoe exome capture variant calling. Sequencing RNA from CIho 3515 and Alexis revealed that only 4 genes within the Rrs18 interval were transcribed in leaf tissue with a serine/threonine protein kinase being the most likely candidate for Rrs18.

In this chapter, we refer to the expressed portion of the barley genome as the relatively small fraction of the total cellular DNA that either contains the genes that ultimately produce proteins, or that directly/indirectly controls the level, location and/or timing of when these genes are expressed and proteins are produced. We start by describing the dynamics of tissue and time-dependent gene expression and how common patterns across multiple samples can provide clues about gene networks involved in common biological processes. We then describe some of the complexities of how a single mRNA template can be differentially processed by alternative splicing to generate multiple different proteins or provide a mechanism to regulate the amount of functional gene product in a cell at a given point in time. We extend our analysis, using a number of biological examples, to address how diverse families of small non-coding microRNAs specifically regulate gene expression, and complete our appraisal by looking at the physical/molecular environment around genes that can result in either the promotion or repression of gene expression. We conclude by assessing some of the issues that remain around our ability to fully exploit the depth and power of current approaches for analysing gene expression and propose improvements that could be made using new but available sequencing and bioinformatics technologies.

Gene prediction in large and highly repetitive grass genomes like barley is complicated by large numbers of transposable elements (TEs), pseudogenes and often incomplete or un-/miss-oriented genomic sequence. In this chapter, we describe the automated gene prediction and annotation pipeline used for the latest barley reference genome sequence, as well as the genomic evidence used to predict gene models. Additional topics cover the (automated) functional annotation, the evaluation of the gene models, and a comprehensive discussion about shortcomings of the current annotation and ways to improve it further.

Summary 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., Significance Statement Triticeae cereals are among the most important crops worldwide yet their epigenomics are largely unstudied. Here we have used ChIP‐Seq to explore the epigenome of barley, identifying and mapping chromatin states of nine modified histones at both gene and genome levels.

Background /Methodology: Taking advantage of a newly available barley (Hordeum vulgare) genome sequence assembly and increasing number of miRNAs identified in other plant species, we carried out computational analyses to identify and characterize miRNAs conserved in barley. We investigate the locations of miRNAs on the barley genome assembly and provide annotation of the functions of their predicted target genes. We compare our results to previous miRNA studies and publicly available barley small RNA libraries. Results: 116 mature miRNA sequences from 60 miRNA families have been found in the barley genome assembly by our miRNA identification pipeline. Closely related cereal crops contain most of the miRNA families that we found in the barley genome assembly. Most miRNA genes were located in intergenic regions or introns. Among the 116 mature miRNAs predicted, 80 have been reported in previous barley miRNA studies. Eight mature miRNA sequences have never reported in the previous barley miRNA studies. Conclusions: This in silico study has provided updated information in characterizing plant miRNAs in barley. The identified miRNA and precursor sequences, their genomic locations as well as predicted target transcripts will serve as valuable resources for future studies.

A genotyping by sequencing (GbS) approach is reported in blackcurrant (Ribes nigrum L.) using a de novo read assembly method developed because of the current absence of a reference genome sequence for this species. A new approach to single nucleotide polymorphism (SNP) genotype calling is described, where individual genotypes for a large number of SNPs were characterised from the GbS counts using a novel method based on a functional regression of major and minor allele read counts. The high-quality GbS SNPs were combined with SNPs and simple sequence repeats generated from other technologies to develop a linkage map with increased marker density and improved genome coverage, containing up to 204 SNPs on each linkage group. SNPs of lower quality were then located on the map using quantitative trait locus (QTL) interval mapping of the proportion of the major allele. Two QTL each for 100-berry weight and Brix scores, measured over three years, were identified using the map. The use of this approach to identify and map a significant number of novel SNPs in a woody species with hitherto limited genomic resources may have generic application to other under-resourced and minor species in the development of cost-effective and efficient high-density genetic maps.

Cereal grasses of the Triticeae tribe have been the major food source in temperate regions since the dawn of agriculture. Their large genomes are characterized by a high content of repetitive elements and large pericentromeric regions that are virtually devoid of meiotic recombination. Here we present a high-quality reference genome assembly for barley (Hordeum vulgare L.). We use chromosome conformation capture mapping to derive the linear order of sequences across the pericentromeric space and to investigate the spatial organization of chromatin in the nucleus at megabase resolution. The composition of genes and repetitive elements differs between distal and proximal regions. Gene family analyses reveal lineage-specific duplications of genes involved in the transport of nutrients to developing seeds and the mobilization of carbohydrates in grains. We demonstrate the importance of the barley reference sequence for breeding by inspecting the genomic partitioning of sequence variation in modern elite germplasm, highlighting regions vulnerable to genetic erosion.

This chapter is designed to be a practical guide to using Tablet for the visualization of next/second-generation (NGS) sequencing data. NGS data is being produced more frequently and in greater data volumes every year. As such, it is increasingly important to have tools which enable biologists and bioinformaticians to understand and gain key insights into their data. Visualization can play a key role in the exploration of such data as well as aid in the visual validation of sequence assemblies and features such as single nucleotide polymorphisms (SNPs). We aim to show several use cases which demonstrate Tablet's ability to visually highlight various situations of interest which can arise in NGS data.

The advent of second-generation sequencing (2GS) has provided a range of significant new challenges for the visualization of sequence assemblies. These include the large volume of data being generated, short-read lengths and different data types and data formats associated with the diversity of new sequencing technologies. This article illustrates how Tabletca high-performance graphical viewer for visualization of 2GS assemblies and read mappingscplays an important role in the analysis of these data. We present Tablet, and through a selection of use cases, demonstrate its value in quality assurance and scientific discovery, through features such as whole-reference coverage overviews, variant highlighting, paired-end read mark-up, GFF3-based feature tracks and protein translations. We discuss the computing and visualization techniques utilized to provide a rich and responsive graphical environment that enables users to view a range of file formats with ease. Tablet installers can be freely downloaded from http://bioinf.hutton.ac.uk/tablet in 32 or 64 -bit versions for Windows, OS X, Linux or Solaris. For further details on the Tablet, contact tablet@hutton.ac.uk.

Second-generation sequencing now provides the potential for low-cost generation of whole-genome sequences. However, for large-genome organisms with high repetitive DNA content, genome-wide short read sequence assembly is currently impossible, with accurate ordering and localization of genes still relying heavily on integration with physical and genetic maps. To facilitate this process, we have used Agilent microarrays to simultaneously address thousands of gene sequences to individual BAC clones and contiguous sequences that form part of an emerging physical map of the large and currently unsequenced 5.3-Gb barley genome. The approach represents a cost-effective, highly parallel alternative to traditional addressing methods. By coupling the gene-to-BAC address data with gene-based molecular markers, thousands of BACs can be anchored directly to the genetic map, thereby generating a framework for orientating and ordering genes, and providing direct links to phenotypic traits.

Species- and genus-level taxonomy of desmids depends largely on shape and detail of the cell wall and chloroplast morphology. The depth of most desmid semicells, relative to the focal depth of conventional light microscopes, means that morphological characteristics are usually illustrated by drawings, made from material that is mounted in water to allow reorientation of specimens to different aspects of shape and pattern. Though a productive approach for two centuries, this has the disadvantages that features not initially detected or thought irrelevant are not recorded, drawing quality is variable, and individual specimens are rarely retained for further study. We describe methods for making permanent preparations of desmid cell walls and using these to produce extended depth of focus summary images and three-dimensional (3-D) reconstructions. Together with World-Wide Web dissemination of image stacks, these advances make it practical to make a desirable change from typification via drawings to typification via single or multiple preserved specimens. They will also facilitate standardization of taxon concepts and identification.

Background Single Nucleotide Polymorphisms (SNPs) are widely used molecular markers, and their use has increased massively since the inception of Next Generation Sequencing (NGS) technologies, which allow detection of large numbers of SNPs at low cost. However, both NGS data and their analysis are error-prone, which can lead to the generation of false positive (FP) SNPs. We explored the relationship between FP SNPs and seven factors involved in mapping-based variant calling — quality of the reference sequence, read length, choice of mapper and variant caller, mapping stringency and filtering of SNPs by read mapping quality and read depth. This resulted in 576 possible factor level combinations. We used error- and variant-free simulated reads to ensure that every SNP found was indeed a false positive. Results The variation in the number of FP SNPs generated ranged from 0 to 36,621 for the 120 million base pairs (Mbp) genome. All of the experimental factors tested had statistically significant effects on the number of FP SNPs generated and there was a considerable amount of interaction between the different factors. Using a fragmented reference sequence led to a dramatic increase in the number of FP SNPs generated, as did relaxed read mapping and a lack of SNP filtering. The choice of reference assembler, mapper and variant caller also significantly affected the outcome. The effect of read length was more complex and suggests a possible interaction between mapping specificity and the potential for contributing more false positives as read length increases. Conclusions The choice of tools and parameters involved in variant calling can have a dramatic effect on the number of FP SNPs produced, with particularly poor combinations of software and/or parameter settings yielding tens of thousands in this experiment. Between-factor interactions make simple recommendations difficult for a SNP discovery pipeline but the quality of the reference sequence is clearly of paramount importance. Our findings are also a stark reminder that it can be unwise to use the relaxed mismatch settings provided as defaults by some read mappers when reads are being mapped to a relatively unfinished reference sequence from e.g. a non-model organism in its early stages of genomic exploration. Electronic supplementary material The online version of this article (doi:10.1186/s12859-015-0801-z) contains supplementary material, which is available to authorized users.

This chapter is designed to be a practical guide to using Tablet for the visualization of next/second-generation (NGS) sequencing data. NGS data is being produced more frequently and in greater data volumes every year. As such, it is increasingly important to have tools which enable biologists and bioinformaticians to understand and gain key insights into their data. Visualization can play a key role in the exploration of such data as well as aid in the visual validation of sequence assemblies and features such as single nucleotide polymorphisms (SNPs). We aim to show several use cases which demonstrate Tablet's ability to visually highlight various situations of interest which can arise in NGS data.

WHIRLY1 is largely targeted to plastids, where it is a major constituent of the nucleoids. To explore WHIRLY1 functions in barley (Hordeum vulgare), RNA interference-knockdown lines (W1-1, W1-7, and W1-9) that have very low levels of HvWHIRLY1 transcripts were characterized in plants grown under optimal and stress conditions. The WHIRLY1-1 (W1-1), W1-7, and W1-9 plants were phenotypically similar to the wild type but produced fewer tillers and seeds. Photosynthesis rates were similar in all lines, but W1-1, W1-7, and W1-9 leaves had significantly more chlorophyll and less sucrose than the wild type. Transcripts encoding specific subsets of chloroplast-localized proteins, such as ribosomal proteins, subunits of the RNA polymerase, and thylakoid nicotinamide adenine dinucleotide (reduced) and cytochrome b6/f complexes, were much more abundant in the W1-7 leaves than the wild type. Although susceptibility of aphid (Myzus persicae) infestation was similar in all lines, the WHIRLY1-deficient plants showed altered responses to nitrogen deficiency, maintaining higher photosynthetic CO2 assimilation rates than the wild type under limiting nitrogen. Although all lines showed globally similar low nitrogen-dependent changes in transcripts and metabolites, the increased abundance of FAR-RED IMPAIRED RESPONSE1-like transcripts in nitrogen-deficient W1-7 leaves infers that WHIRLY1 has a role in communication between plastid and nuclear genes encoding photosynthetic proteins during abiotic stress.

A method for automatic identification of diatoms (single-celled algae with silica shells) based on extraction of features on the contour of the cells by multi-scale mathematical morphology is presented. After extracting the contour of the cell, it is smoothed adaptively, encoded using Freeman chain code, and converted into a curvature representation which is invariant under translation and scale change. A curvature scale space is built from these data, and the most important features are extracted from it by unsupervised cluster analysis. The resulting pattern vectors, which are also rotation-invariant, provide the input for automatic identification of diatoms by decision trees and k-nearest neighbor classifiers. The method is tested on two large sets of diatom images. The techniques used are applicable to other shapes besides diatoms.