Tablet: Visualizing Next-Generation Sequence Assemblies and Mappings.

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.

tropiTree: an NGS-based EST-SSR resource for 24 tropical tree species.

The development of genetic tools for non-model organisms has been hampered by cost, but advances in next-generation sequencing (NGS) have created new opportunities. In ecological research, this raises the prospect for developing molecular markers to simultaneously study important genetic processes such as gene flow in multiple non-model plant species within complex natural and anthropogenic landscapes. Here, we report the use of bar-coded multiplexed paired-end Illumina NGS for the de novo development of expressed sequence tag-derived simple sequence repeat (EST-SSR) markers at low cost for a range of 24 tree species. Each chosen tree species is important in complex tropical agroforestry systems where little is currently known about many genetic processes. An average of more than 5,000 EST-SSRs was identified for each of the 24 sequenced species, whereas prior to analysis 20 of the species had fewer than 100 nucleotide sequence citations. To make results available to potential users in a suitable format, we have developed an open-access, interactive online database, tropiTree (http://bioinf.hutton.ac.uk/tropiTree), which has a range of visualisation and search facilities, and which is a model for the efficient presentation and application of NGS data.

Resistance gene enrichment sequencing (RenSeq) enables reannotation of the NB-LRR gene family from sequenced plant genomes and rapid mapping of resistance loci in segregating populations.

RenSeq is a NB-LRR (nucleotide binding-site leucine-rich repeat) gene-targeted, Resistance gene enrichment and sequencing method that enables discovery and annotation of pathogen resistance gene family members in plant genome sequences. We successfully applied RenSeq to the sequenced potato Solanum tuberosum clone DM, and increased the number of identified NB-LRRs from 438 to 755. The majority of these identified R gene loci reside in poorly or previously unannotated regions of the genome. Sequence and positional details on the 12 chromosomes have been established for 704 NB-LRRs and can be accessed through a genome browser that we provide. We compared these NB-LRR genes and the corresponding oligonucleotide baits with the highest sequence similarity and demonstrated that ~80% sequence identity is sufficient for enrichment. Analysis of the sequenced tomato S. lycopersicum 'Heinz 1706' extended the NB-LRR complement to 394 loci. We further describe a methodology that applies RenSeq to rapidly identify molecular markers that co-segregate with a pathogen resistance trait of interest. In two independent segregating populations involving the wild Solanum species S. berthaultii (Rpi-ber2) and S. ruiz-ceballosii (Rpi-rzc1), we were able to apply RenSeq successfully to identify markers that co-segregate with resistance towards the late blight pathogen Phytophthora infestans. These SNP identification workflows were designed as easy-to-adapt Galaxy pipelines.

Using Tablet for visual exploration of second-generation sequencing 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 Tablet-a high-performance graphical viewer for visualization of 2GS assemblies and read mappings-plays 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.

Alternative splicing and nonsense-mediated decay modulate expression of important regulatory genes in Arabidopsis.

Alternative splicing (AS) coupled to nonsense-mediated decay (NMD) is a post-transcriptional mechanism for regulating gene expression. We have used a high-resolution AS RT-PCR panel to identify endogenous AS isoforms which increase in abundance when NMD is impaired in the Arabidopsis NMD factor mutants, upf1-5 and upf3-1. Of 270 AS genes (950 transcripts) on the panel, 102 transcripts from 97 genes (32%) were identified as NMD targets. Extrapolating from these data around 13% of intron-containing genes in the Arabidopsis genome are potentially regulated by AS/NMD. This cohort of naturally occurring NMD-sensitive AS transcripts also allowed the analysis of the signals for NMD in plants. We show the importance of AS in introns in 5' or 3'UTRs in modulating NMD-sensitivity of mRNA transcripts. In particular, we identified upstream open reading frames overlapping the main start codon as a new trigger for NMD in plants and determined that NMD is induced if 3'-UTRs were >350 nt. Unexpectedly, although many intron retention transcripts possess NMD features, they are not sensitive to NMD. Finally, we have shown that AS/NMD regulates the abundance of transcripts of many genes important for plant development and adaptation including transcription factors, RNA processing factors and stress response genes.

Identification, utilisation and mapping of novel transcriptome-based markers from blackcurrant (Ribes nigrum).

BACKGROUND: Deep-level second generation sequencing (2GS) technologies are now being applied to non-model species as a viable and favourable alternative to Sanger sequencing. Large-scale SNP discovery was undertaken in blackcurrant (Ribes nigrum L.) using transcriptome-based 2GS 454 sequencing on the parental genotypes of a reference mapping population, to generate large numbers of novel markers for the construction of a high-density linkage map. RESULTS: Over 700,000 reads were produced, from which a total of 7,000 SNPs were found. A subset of polymorphic SNPs was selected to develop a 384-SNP OPA assay using the Illumina BeadXpress platform. Additionally, the data enabled identification of 3,000 novel EST-SSRs. The selected SNPs and SSRs were validated across diverse Ribes germplasm, including mapping populations and other selected Ribes species.SNP-based maps were developed from two blackcurrant mapping populations, incorporating 48% and 27% of assayed SNPs respectively. A relatively high proportion of visually monomorphic SNPs were investigated further by quantitative trait mapping of theta score outputs from BeadStudio analysis, and this enabled additional SNPs to be placed on the two maps. CONCLUSIONS: The use of 2GS technology for the development of markers is superior to previously described methods, in both numbers of markers and biological informativeness of those markers. Whilst the numbers of reads and assembled contigs were comparable to similar sized studies of other non-model species, here a high proportion of novel genes were discovered across a wide range of putative function and localisation. The potential utility of markers developed using the 2GS approach in downstream breeding applications is discussed.

Highly parallel gene-to-BAC addressing using microarrays.

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.

Comparative visualization of genetic and physical maps with Strudel.

UNLABELLED: Data visualization can play a key role in comparative genomics, for example, underpinning the investigation of conserved synteny patterns. Strudel is a desktop application that allows users to easily compare both genetic and physical maps interactively and efficiently. It can handle large datasets from several genomes simultaneously, and allows all-by-all comparisons between these. AVAILABILITY AND IMPLEMENTATION: Installers for Strudel are available for Windows, Linux, Solaris and Mac OS X at http://bioinf.scri.ac.uk/strudel/.

Flapjack--graphical genotype visualization.

SUMMARY: New software tools for graphical genotyping are required that can routinely handle the large data volumes generated by the high-throughput single-nucleotide polymorphism (SNP) platforms, genotyping-by-sequencing and other comparable genotyping technologies. Flapjack has been developed to facilitate analysis of these data, providing real time rendering with rapid navigation and comparisons between lines, markers and chromosomes, with visualization, sorting and querying based on associated data, such as phenotypes, quantitative trait loci or other mappable features. AVAILABILITY: Flapjack is freely available for Microsoft Windows, Mac OS X, Linux and Solaris, and can be downloaded from http://bioinf.scri.ac.uk/flapjack .

Candidate genes associated with bud dormancy release in blackcurrant (Ribes nigrum L.).

BACKGROUND: The detrimental effects of mild winter temperatures on the consistency of cropping of blackcurrant (Ribes nigrum L.) in parts of Europe have led to increasing interest in the genetic control of dormancy release in this species. This study examined patterns of gene expression in leaf buds of blackcurrant to identify key differential changes in these profiles around the time of budbreak. RESULTS: Using leaf bud tissue of blackcurrant, a cDNA library was generated as a source of blackcurrant ESTs for construction of a custom microarray, which was used to identify differential gene expression during dormancy release. Gene activity was lowest in early stages of dormancy, increasing to reach a maximum around the time of budbreak. Genes with significantly changing expression profiles were clustered and evidence is provided for the transient activity of genes previously associated with dormancy processes in other species. Expression profiling identified candidate genes which were mapped onto a blackcurrant genetic linkage map containing budbreak-related QTL. Three genes, which putatively encode calmodulin-binding protein, beta tubulin and acetyl CoA carboxylase respectively, were found to co-localise with budbreak QTL. CONCLUSIONS: This study provides insight into the genetic control of dormancy transition in blackcurrant, identifying key changes in gene expression around budbreak. Genetic mapping of ESTs enabled the identification of genes which co-localise with previously-characterised blackcurrant QTL, and it is concluded that these genes have probable roles in release of dormancy and can therefore provide a basis for the development of genetic markers for future breeding deployment.

Exploiting induced variation to dissect quantitative traits in barley.

The identification of genes underlying complex quantitative traits such as grain yield by means of conventional genetic analysis (positional cloning) requires the development of several large mapping populations. However, it is possible that phenotypically related, but more extreme, allelic variants generated by mutational studies could provide a means for more efficient cloning of QTLs (quantitative trait loci). In barley (Hordeum vulgare), with the development of high-throughput genome analysis tools, efficient genome-wide identification of genetic loci harbouring mutant alleles has recently become possible. Genotypic data from NILs (near-isogenic lines) that carry induced or natural variants of genes that control aspects of plant development can be compared with the location of QTLs to potentially identify candidate genes for development--related traits such as grain yield. As yield itself can be divided into a number of allometric component traits such as tillers per plant, kernels per spike and kernel size, mutant alleles that both affect these traits and are located within the confidence intervals for major yield QTLs may represent extreme variants of the underlying genes. In addition, the development of detailed comparative genomic models based on the alignment of a high-density barley gene map with the rice and sorghum physical maps, has enabled an informed prioritization of 'known function' genes as candidates for both QTLs and induced mutant genes.

An eQTL analysis of partial resistance to Puccinia hordei in barley.

BACKGROUND: Genetic resistance to barley leaf rust caused by Puccinia hordei involves both R genes and quantitative trait loci. The R genes provide higher but less durable resistance than the quantitative trait loci. Consequently, exploring quantitative or partial resistance has become a favorable alternative for controlling disease. Four quantitative trait loci for partial resistance to leaf rust have been identified in the doubled haploid Steptoe (St)/Morex (Mx) mapping population. Further investigations are required to study the molecular mechanisms underpinning partial resistance and ultimately identify the causal genes. METHODOLOGY/PRINCIPAL FINDINGS: We explored partial resistance to barley leaf rust using a genetical genomics approach. We recorded RNA transcript abundance corresponding to each probe on a 15K Agilent custom barley microarray in seedlings from St and Mx and 144 doubled haploid lines of the St/Mx population. A total of 1154 and 1037 genes were, respectively, identified as being P. hordei-responsive among the St and Mx and differentially expressed between P. hordei-infected St and Mx. Normalized ratios from 72 distant-pair hybridisations were used to map the genetic determinants of variation in transcript abundance by expression quantitative trait locus (eQTL) mapping generating 15685 eQTL from 9557 genes. Correlation analysis identified 128 genes that were correlated with resistance, of which 89 had eQTL co-locating with the phenotypic quantitative trait loci (pQTL). Transcript abundance in the parents and conservation of synteny with rice allowed us to prioritise six genes as candidates for Rphq11, the pQTL of largest effect, and highlight one, a phospholipid hydroperoxide glutathione peroxidase (HvPHGPx) for detailed analysis. CONCLUSIONS/SIGNIFICANCE: The eQTL approach yielded information that led to the identification of strong candidate genes underlying pQTL for resistance to leaf rust in barley and on the general pathogen response pathway. The dataset will facilitate a systems appraisal of this host-pathogen interaction and, potentially, for other traits measured in this population.

Tablet--next generation sequence assembly visualization.

SUMMARY: Tablet is a lightweight, high-performance graphical viewer for next-generation sequence assemblies and alignments. Supporting a range of input assembly formats, Tablet provides high-quality visualizations showing data in packed or stacked views, allowing instant access and navigation to any region of interest, and whole contig overviews and data summaries. Tablet is both multi-core aware and memory efficient, allowing it to handle assemblies containing millions of reads, even on a 32-bit desktop machine. AVAILABILITY: Tablet is freely available for Microsoft Windows, Apple Mac OS X, Linux and Solaris. Fully bundled installers can be downloaded from http://bioinf.scri.ac.uk/tablet in 32- and 64-bit versions.

Use of genome sequence data in the design and testing of SSR markers for Phytophthora species.

BACKGROUND: Microsatellites or single sequence repeats (SSRs) are a powerful choice of marker in the study of Phytophthora population biology, epidemiology, ecology, genetics and evolution. A strategy was tested in which the publicly available unigene datasets extracted from genome sequences of P. infestans, P. sojae and P. ramorum were mined for candidate SSR markers that could be applied to a wide range of Phytophthora species. RESULTS: A first approach, aimed at the identification of polymorphic SSR loci common to many Phytophthora species, yielded 171 reliable sequences containing 211 SSRs. Microsatellites were identified from 16 target species representing the breadth of diversity across the genus. Repeat number ranged from 3 to 16 with most having seven repeats or less and four being the most commonly found. Trinucleotide repeats such as (AAG)n, (AGG)n and (AGC)n were the most common followed by pentanucleotide, tetranucleotide and dinucleotide repeats. A second approach was aimed at the identification of useful loci common to a restricted number of species more closely related to P. sojae (P. alni, P. cambivora, P. europaea and P. fragariae). This analysis yielded 10 trinucleotide and 2 tetranucleotide SSRs which were repeated 4, 5 or 6 times. CONCLUSION: Key studies on inter- and intra-specific variation of selected microsatellites remain. Despite the screening of conserved gene coding regions, the sequence diversity between species was high and the identification of useful SSR loci applicable to anything other than the most closely related pairs of Phytophthora species was challenging. That said, many novel SSR loci for species other than the three 'source species' (P. infestans, P. sojae and P. ramorum) are reported, offering great potential for the investigation of Phytophthora populations. In addition to the presence of microsatellites, many of the amplified regions may represent useful molecular marker regions for other studies as they are highly variable and easily amplifiable from different Phytophthora species.

Isolation, analysis and marker utility of novel miniature inverted repeat transposable elements from the barley genome.

Four previously undescribed families of miniature inverted repeat transposable elements (MITEs) were isolated by searching barley genomic DNA using structure-based criteria. Putative MITEs were confirmed by PCR to determine their insertional polymorphism in a panel of diverse barley germplasm. Copy numbers for all these familes are somewhat low (less than 1,000 copies per family per haploid genome). In contrast to previous studies, a higher proportion of insertions of the new MITEs are found within known transposable elements (27%) than are associated with genes (15%). Preliminary studies were conducted on two of the new MITE families to test their utility as molecular markers. Insertional polymorphism levels for both the families are high and diversity trees produced by both the families are similar and congruent with known relationships among the germplasm studied, suggesting that both the MITE families are useful markers of barley genetic diversity.

Co-ordinated gene expression during phases of dormancy release in raspberry (Rubus idaeus L.) buds.

Bud break in raspberry (Rubus idaeus L.) is often poor and uneven, with many of the subapical buds remaining in a dormant state. In order to determine the dormancy status of raspberry buds, an empirical measure of bud burst in a growth-permissive environment following exposure to chilling (4 degrees C cold storage) was developed. For cv. Glen Ample, percentage bud burst in intact canes and isolated nodes was recorded after 14 d. Isolated nodes (a measure of endodormancy) achieved 100% bud burst after approximately 1500 h chilling whereas buds on intact plants (combined endo- and paradormancy) required an additional 1000 h chilling. A microarray approach was used to follow changes in gene expression that occurred during dormancy transition. The probes for the microarrays were obtained from endodormant and paradormant raspberry bud cDNA libraries. The expression profiles of 5300 clones from these libraries were subjected to principal component analysis to determine the most significant expression patterns. Sequence analysis of these clones, in many cases, enabled their functional categorization and the development of hypotheses concerning the mechanisms of bud dormancy release. Thus a set of novel candidates for key dormancy-related genes from raspberry buds have been identified. Bud dormancy is fundamental to the study of plant developmental processes and, in addition, its regulation is of significant economic importance to fruit and horticultural industries.

Recent history of artificial outcrossing facilitates whole-genome association mapping in elite inbred crop varieties.

Genomewide association studies depend on the extent of linkage disequilibrium (LD), the number and distribution of markers, and the underlying structure in populations under study. Outbreeding species generally exhibit limited LD, and consequently, a very large number of markers are required for effective whole-genome association genetic scans. In contrast, several of the world's major food crops are self-fertilizing inbreeding species with narrow genetic bases and theoretically extensive LD. Together these are predicted to result in a combination of low resolution and a high frequency of spurious associations in LD-based studies. However, inbred elite plant varieties represent a unique human-induced pseudo-outbreeding population that has been subjected to strong selection for advantageous alleles. By assaying 1,524 genomewide SNPs we demonstrate that, after accounting for population substructure, the level of LD exhibited in elite northwest European barley, a typical inbred cereal crop, can be effectively exploited to map traits by using whole-genome association scans with several hundred to thousands of biallelic SNPs.

Genome-wide SNP discovery and linkage analysis in barley based on genes responsive to abiotic stress.

More than 2,000 genome-wide barley single nucleotide polymorphisms (SNPs) were developed by resequencing unigene fragments from eight diverse accessions. The average genome-wide SNP frequency observed in 877 unigenes was 1 SNP per 200 bp. However, SNP frequency was highly variable with the least number of SNP and SNP haplotypes observed within European cultivated germplasm reflecting effects of breeding history on genetic diversity. More than 300 SNP loci were mapped genetically in three experimental mapping populations which allowed the construction of an integrated SNP map incorporating a large number of RFLP, AFLP and SSR markers (1,237 loci in total). The genes used for SNP discovery were selected based on their transcriptional response to a variety of abiotic stresses. A set of known barley abiotic stress QTL was positioned on the linkage map, while the available sequence and gene expression information facilitated the identification of genes potentially associated with these traits. Comparison of the sequenced SNP loci to the rice genome sequence identified several regions of highly conserved gene order providing a framework for marker saturation in barley genomic regions of interest. The integration of genome-wide SNP and expression data with available genetic and phenotypic information will facilitate the identification of gene function in barley and other non-model organisms.

Single-feature polymorphism discovery in the barley transcriptome.

A probe-level model for analysis of GeneChip gene-expression data is presented which identified more than 10,000 single-feature polymorphisms (SFP) between two barley genotypes. The method has good sensitivity, as 67% of known single-nucleotide polymorphisms (SNP) were called as SFPs. This method is applicable to all oligonucleotide microarray data, accounts for SNP effects in gene-expression data and represents an efficient and versatile approach for highly parallel marker identification in large genomes.

Characterisation and functional analysis of two barley caleosins expressed during barley caryopsis development.

Two full-length cDNA sequences homologous to caleosin, a seed-storage oil-body protein from sesame, were identified from a series of barley grain development cDNA libraries and further characterised. The cDNAs, subsequently termed HvClo1 and HvClo2, encode proteins of 34 kDa and 28 kDa, respectively. Real-time RT-PCR indicated that HvClo1 is expressed abundantly during the later stages of embryogenesis and is seed-specific, accumulating in the scutellum of mature embryos. HvClo2 is expressed mainly in the endosperm tissues of the developing grain. We show that HvClo1 and HvClo2 are paralogs that co-segregate on barley chromosome 2HL. Transient expression of HvClo1 in lipid storage and non-storage cells of barley using biolistic particle bombardment indicates that caleosins have different subcellular locations from the structural oil-body protein oleosin, and by inference participate in different sorting pathways. We observe that caleosin sorts via small vesicles, suggesting a likely association with lipid trafficking, membrane expansion and oil-body biogenesis.