BlaSTorage: a fast package to parse, manage and store BLAST results

Background

Large-scale sequence studies requiring BLAST-based analysis produce huge amounts of data to be parsed. BLAST parsers are available, but they are often missing some important features, such as keeping all information from the raw BLAST output, allowing direct access to single results, and performing logical operations over them.

Findings

We implemented BlaSTorage, a Python package that parses multi BLAST results and returns them in a purpose-built object-database format. Unlike other BLAST parsers, BlaSTorage retains and stores all parts of BLAST results, including alignments, without loss of information; a complete API allows access to all the data components.

Conclusions

BlaSTorage shows comparable speed of more basic parser written in compiled languages as C++ and can be easily integrated into web applications or software pipelines.     

RootScan: Software for high-throughput analysis of root anatomical traits

Background and aims

RootScan is a program for semi-automated image analysis of anatomical traits in root cross-sections.

Methods

RootScan uses pixel thresholds to separate the cross-section from its background and to divide it into tissue regions. Area measurements and object counts are performed within various regions of interest. A graphical user interface permits the user to see which regions are selected, to edit those selections, and to rate and comment on the data. The structure of the program allows for organized workflow and increased data collection efficiency.

Results

The program collects data on more than 20 variables per image including areas of the cross-section, stele, cortex, aerenchyma lacunae, xylem vessels, and counts of cortical cells and cell files. An increased rate of data collection allows collection of four times more variables in less time than is possible with current methods. Correlation analysis shows that RootScan data is equal or greater in accuracy than data collected with Photoshop.

Conclusions

Compared with currently available tools, this software offers considerable improvements in the amount and quality of data, ease of use, and time needed for data collection. RootScan permits phenotypic scoring of physiologically and agronomically important traits on a large number of genotypes.     

Towards better digital pathology workflows: programming libraries for high-speed sharpness assessment of Whole Slide Images

Background

Since microscopic slides can now be automatically digitized and integrated in the clinical workflow, quality assessment of Whole Slide Images (WSI) has become a crucial issue. We present a no-reference quality assessment method that has been thoroughly tested since 2010 and is under implementation in multiple sites, both public university-hospitals and private entities. It is part of the FlexMIm R&D project which aims to improve the global workflow of digital pathology. For these uses, we have developed two programming libraries, in Java and Python, which can be integrated in various types of WSI acquisition systems, viewers and image analysis tools.

Methods

Development and testing have been carried out on a MacBook Pro i7 and on a bi-Xeon 2.7GHz server. Libraries implementing the blur assessment method have been developed in Java, Python, PHP5 and MySQL5. For web applications, JavaScript, Ajax, JSON and Sockets were also used, as well as the Google Maps API. Aperio SVS files were converted into the Google Maps format using VIPS and Openslide libraries.

Results

We designed the Java library as a Service Provider Interface (SPI), extendable by third parties. Analysis is computed in real-time (3 billion pixels per minute). Tests were made on 5000 single images, 200 NDPI WSI, 100 Aperio SVS WSI converted to the Google Maps format.

Conclusions

Applications based on our method and libraries can be used upstream, as calibration and quality control tool for the WSI acquisition systems, or as tools to reacquire tiles while the WSI is being scanned. They can also be used downstream to reacquire the complete slides that are below the quality threshold for surgical pathology analysis. WSI may also be displayed in a smarter way by sending and displaying the regions of highest quality before other regions. Such quality assessment scores could be integrated as WSI's metadata shared in clinical, research or teaching contexts, for a more efficient medical informatics workflow.

     

Simple tools for assembling and searching high-density picolitre pyrophosphate sequence data

Background

The advent of pyrophosphate sequencing makes large volumes of sequencing data available at a lower cost than previously possible. However, the short read lengths are difficult to assemble and the large dataset is difficult to handle. During the sequencing of a virus from the tsetse fly, Glossina pallidipes, we found the need for tools to search quickly a set of reads for near exact text matches.

Methods

A set of tools is provided to search a large data set of pyrophosphate sequence reads under a "live" CD version of Linux on a standard PC that can be used by anyone without prior knowledge of Linux and without having to install a Linux setup on the computer. The tools permit short lengths of de novo assembly, checking of existing assembled sequences, selection and display of reads from the data set and gathering counts of sequences in the reads.

Results

Demonstrations are given of the use of the tools to help with checking an assembly against the fragment data set; investigating homopolymer lengths, repeat regions and polymorphisms; and resolving inserted bases caused by incomplete chain extension.

Conclusion

The additional information contained in a pyrophosphate sequencing data set beyond a basic assembly is difficult to access due to a lack of tools. The set of simple tools presented here would allow anyone with basic computer skills and a standard PC to access this information.     

PedGenie: meta genetic association testing in mixed family and case-control designs

Background

With the increased availability of high throughput data, such as DNA microarray data, researchers are capable of producing large amounts of biological data. During the analysis of such data often there is the need to further explore the similarity of genes not only with respect to their expression, but also with respect to their functional annotation which can be obtained from Gene Ontology (GO).

Results

We present the freely available software package GOSim, which allows to calculate the functional similarity of genes based on various information theoretic similarity concepts for GO terms. GOSim extends existing tools by providing additional lately developed functional similarity measures for genes. These can e.g. be used to cluster genes according to their biological function. Vice versa, they can also be used to evaluate the homogeneity of a given grouping of genes with respect to their GO annotation. GOSim hence provides the researcher with a flexible and powerful tool to combine knowledge stored in GO with experimental data. It can be seen as complementary to other tools that, for instance, search for significantly overrepresented GO terms within a given group of genes.

Conclusion

GOSim is implemented as a package for the statistical computing environment R and is distributed under GPL within the CRAN project.     

Generation of selectable marker-free transgenic eggplant resistant to Alternaria solani using the R/RS site-specific recombination system

Key message

Marker-free transgenic eggplants, exhibiting enhanced resistance to Alternaria solani , can be generated on plant growth regulators (PGRs)- and antibiotic-free MS medium employing the multi-auto-transformation (MAT) vector, pMAT21 - wasabi defensin , wherein isopentenyl transferase ( ipt ) gene is used as a positive selection marker.

Abstract

Use of the selection marker genes conferring antibiotic or herbicide resistance in transgenic plants has been considered a serious problem for environment and the public. Multi-auto-transformation (MAT) vector system has been one of the tools to excise the selection marker gene and produce marker-free transgenic plants. Ipt gene was used as a selection marker gene. Wasabi defensin gene, isolated from Wasabia japonica (a Japanese horseradish which has been a potential source of antimicrobial proteins), was used as a gene of interest. Wasabi defensin gene was cloned from the binary vector, pEKH-WD, to an ipt-type MAT vector, pMAT21, by gateway cloning technology and transferred to Agrobacterium tumefaciens strain EHA105. Infected cotyledon explants of eggplant were cultured on PGRs- and antibiotic-free MS medium. Extreme shooty phenotype/ipt shoots were produced by the explants infected with the pMAT21-wasabi defensin (WD). The same PGRs- and antibiotic-free MS medium was used in subcultures of the ipt shoots. Subsequently, morphologically normal shoots emerged from the Ipt shoots. Molecular analyses of genomic DNA from transgenic plants confirmed the integration of the WD gene and excision of the selection marker (ipt gene). Expression of the WD gene was confirmed by RT-PCR and Northern blot analyses. In vitro whole plant and detached leaf assay of the marker-free transgenic plants exhibited enhanced resistance against Alternaria solani.     

Transparent mediation-based access to multiple yeast data sources using an ontology driven interface

Background

Over the past decade the workflow system paradigm has evolved as an efficient and user-friendly approach for developing complex bioinformatics applications. Two popular workflow systems that have gained acceptance by the bioinformatics community are Taverna and Galaxy. Each system has a large user-base and supports an ever-growing repository of application workflows. However, workflows developed for one system cannot be imported and executed easily on the other. The lack of interoperability is due to differences in the models of computation, workflow languages, and architectures of both systems. This lack of interoperability limits sharing of workflows between the user communities and leads to duplication of development efforts.

Results

In this paper, we present Tavaxy, a stand-alone system for creating and executing workflows based on using an extensible set of re-usable workflow patterns. Tavaxy offers a set of new features that simplify and enhance the development of sequence analysis applications: It allows the integration of existing Taverna and Galaxy workflows in a single environment, and supports the use of cloud computing capabilities. The integration of existing Taverna and Galaxy workflows is supported seamlessly at both run-time and design-time levels, based on the concepts of hierarchical workflows and workflow patterns. The use of cloud computing in Tavaxy is flexible, where the users can either instantiate the whole system on the cloud, or delegate the execution of certain sub-workflows to the cloud infrastructure.

Conclusions

Tavaxy reduces the workflow development cycle by introducing the use of workflow patterns to simplify workflow creation. It enables the re-use and integration of existing (sub-) workflows from Taverna and Galaxy, and allows the creation of hybrid workflows. Its additional features exploit recent advances in high performance cloud computing to cope with the increasing data size and complexity of analysis. The system can be accessed either through a cloud-enabled web-interface or downloaded and installed to run within the user's local environment. All resources related to Tavaxy are available at http://www.tavaxy.org.     

Managing potato wart: a review of present research status and future perspective

Key message

Identification of resistance genes to potato wart disease caused by Synchytrium endobioticum is the key for developing diagnostic markers for breeding resistant cultivars. We present an overview on the current knowledge of this host-pathogen system and molecular advances while highlighting future research focus.

Abstract

Potato wart is a quarantined disease of cultivated potato (Solanum tuberosum L.) caused by the obligate biotrophic, soil-borne fungus Synchytrium endobioticum (Schilb.) Perc. Since its discovery by Schilberszky in 1896, the management of wart disease was enabled by research efforts focusing on understanding and classifying the causative agent, its mode of infection, pathogenesis, geographical distribution, detection and chemical control, on developing screening methods for host resistance and on genetic analyses, which led to the development of resistant cultivars. These early successes are currently challenged by new S. endobioticum pathotypes evolving and the increased risk of dissemination by potato tuber trade. New research efforts are therefore required to ensure continuation of effective and sustainable management of the potato wart disease. Advances in molecular biology and genomic tools offer potential for innovations. This review presents an overview on what we know about this complex host-pathogen interaction, highlights recent molecular work and embarks on an outlook towards future research directions.     

GelJ – a tool for analyzing DNA fingerprint gel images

Background

DNA fingerprinting is a technique for comparing DNA patterns that has applications in a wide variety of contexts. Several commercial and freely-available tools can be used to analyze DNA fingerprint gel images; however, commercial tools are expensive and usually difficult to use; and, free tools support the basic functionality for DNA fingerprint analysis, but lack some instrumental features to obtain accurate results.

Results

In this paper, we present GelJ, a feather-weight, user-friendly, platform-independent, open-source and free tool for analyzing DNA fingerprint gel images. Some of the outstanding features of GelJ are mechanisms for accurate lane- and band-detection, several options for computing migration models, a number of band- and curve-based similarity methods, different techniques for generating dendrograms, comparison of banding patterns from different experiments, and database support.

Conclusions

GelJ is an easy to use tool for analyzing DNA fingerprint gel images. It combines the best characteristics of both free and commercial tools: GelJ is light and simple to use (as free programs), but it also includes the necessary features to obtain precise results (as commercial programs). In addition, GelJ incorporates new functionality that is not supported by any other tool.

     

Quantifying the effect of soil moisture content on segmenting root system architecture in X-ray computed tomography images

Aims

A commonly accepted challenge when visualising plant roots in X-ray micro Computed Tomography (μCT) images is the similar X-ray attenuation of plant roots and soil phases. Soil moisture content remains a recognised, yet currently uncharacterised source of segmentation error. This work sought to quantify the effect of soil moisture content on the ability to segment roots from soil in μCT images.

Methods

Rice (Oryza sativa) plants grown in contrasting soils (loamy sand and clay loam) were μCT scanned daily for nine days whilst drying from saturation. Root volumes were segmented from μCT images and compared with volumes derived by root washing.

Results

At saturation the overlapping attenuation values of root material, water-filled soil pores and soil organic matter significantly hindered segmentation. However, in dry soil (ca. six days of drying post-saturation) the air-filled pores increased image noise adjacent to roots and impeded accurate visualisation of root material. The root volume was most accurately segmented at field capacity.

Conclusions

Root volumes can be accurately segmented from μCT images of undisturbed soil without compromising the growth requirements of the plant providing soil moisture content is kept at field capacity. We propose all future studies in this area should consider the error associated with scanning at different soil moisture contents.     

Differential regulation of cytokinin oxidase genes and cytokinin-induced auxin biosynthesis by cellular cytokinin level in transgenic poplars

Key message

The present work with transgenic poplar lines producing varying levels of trans -zeatin suggests the existence of a switching threshold for triggering ckx gene expression or suppressing cytokinin-induced auxin.

Abstract

Cytokinins have an important role in growth and developmental processes of plants. Transgenic plants with varying levels of cellular cytokinin are convenient tools for studying its role in morphogenetic as well as molecular responses. In this work, the transgenic lines producing either high level of cellular trans-zeatin (HX lines) or moderate level (MX lines) were compared with regard to their cytokinin oxidase activities and cellular auxin content. The HX lines showed typical cytokinin phenotypes including leafy shoots and spontaneous shoot formation on hormone free medium. In contrast, the MX lines did not show any striking phenotypes. However, in leaf disk culture on hormone free medium, they regenerated roots and subsequently formed shoots from the roots. Determination of cellular IAA content revealed a significant increase in the level in MX lines but not in HX lines. Of nine cytokinin oxidase genes (ckx) examined by qPCR, five were activated in HX lines but not in MX lines. Among them, ckx4 appeared to play a key role in maintaining cellular cytokinin level since it showed more than 1,000-fold increase in HX lines and in the leaf disks of untransformed control exposed to exogenous cytokinins. Although low level of cellular cytokinin did not induce the expression of ckx genes, it appeared to trigger cellular IAA biosynthesis.     

Fine mapping and candidate gene analysis of the novel thermo-sensitive genic male sterility tms9-1 gene in rice

Key message

Fine mapping of the novel thermo-sensitive genic male sterility locus tms9 - 1 in the traditional TGMS line HengnongS-1 revealed that the MALE STERILITY1 homolog OsMS1 is the candidate gene.

Abstract

Photoperiod-thermo-sensitive genic male sterility (P/TGMS) has been widely used in the two-line hybrid rice breeding system. HengnongS-1 is one of the oldest TGMS lines and is often used in indica two-line breeding programs in China. In this study, our genetic analysis showed that the TGMS gene in HengnongS-1 was controlled by a single recessive gene that was non-allelic with the other TGMS loci identified, including C815S, Zhu1S and Y58S. Using SSR markers and bulked segregant analysis, we located the TGMS locus on chromosome 9 and named the gene tms9-1. Fine mapping further narrowed the tms9-1 loci to a 162 kb interval between two dCAPS markers. Sequence analysis revealed that a T to C substitution results in an amino acid change in the tms9-1 candidate gene (Os09g27620) in HengnongS-1 as compared to Minghui63. Sequencing of other rice accessions, including six P/TGMS lines, seven indica varieties and nine japonica varieties, showed that this SNP was exclusive to HengnongS-1. With multiple sequence alignment and expression pattern analyses, the rice MALE STERILITY1 homolog OsMS1 gene was identified as the candidate gene for tms9-1. Therefore, our study identified a novel TGMS locus and will facilitate the functional identification of the tms9-1 gene. Moreover, the markers linked to the tms9-1 gene will provide useful tools for the development of new TGMS lines by marker-assisted selection in two-line hybrid rice breeding programs.     

MulSatDB: a first online database for mulberry microsatellites

Key message

Simple sequence repeat motifs were mined from the genome and EST sequences of Morus notabilis and archived in MulSatDB. Bioinformatics tools were integrated with the database for the analysis of genomic datasets.

Abstract

Mulberry is a crop of economic importance in sericulture, which shapes the lives of millions of rural people among different Eurasian and Latin American countries. Limited availability of genomic resources has constrained the molecular breeding efforts in mulberry, a poorly studied crop. Microsatellite or simple sequence repeat (SSR) has revolutionized the plant breeding and is used in linkage mapping, association studies, diversity, and parentage analysis, etc. Recent availability of mulberry whole genome assembly provided an opportunity for the development of mulberry-specific DNA markers. In this study, we mined a total of 217,312 microsatellites from whole genome and 961 microsatellites from EST sequences of Morus notabilis. Mono-repeats were predominant among both whole genome and EST sequences. The SSR containing EST sequences were functionally annotated, and SSRs mined from whole genome were mapped on chromosomes of the phylogenetically related genus—Fragaria vesca, to aid the selection of markers based on the function and location. All the mined markers were archived in the mulberry microsatellite database (MulSatDB), and the markers can be retrieved based on different criteria like marker location, repeat kind, motif type and size. Primer3plus and CMap tools are integrated with the database to design primers for PCR amplification and to visualize markers on F. vesca chromosomes, respectively. A blast tool is also integrated to collate new markers with the database. MulSatDB is the first and complete destination for mulberry researchers to browse SSR markers, design primers, and locate markers on strawberry chromosomes. MulSatDB is freely accessible at http://btismysore.in/mulsatdb.     

ReCount: A multi-experiment resource of analysis-ready RNA-seq gene count datasets

Background

With next-generation sequencing technologies, experiments that were considered prohibitive only a few years ago are now possible. However, while these technologies have the ability to produce enormous volumes of data, the sequence reads are prone to error. This poses fundamental hurdles when genetic diversity is investigated.

Results

We developed ShoRAH, a computational method for quantifying genetic diversity in a mixed sample and for identifying the individual clones in the population, while accounting for sequencing errors. The software was run on simulated data and on real data obtained in wet lab experiments to assess its reliability.

Conclusions

ShoRAH is implemented in C++, Python, and Perl and has been tested under Linux and Mac OS X. Source code is available under the GNU General Public License at http://www.cbg.ethz.ch/software/shorah.     

IJ_Rhizo: an open-source software to measure scanned images of root samples

Background and aims

This paper provides an overview of the measuring capabilities of IJ_Rhizo, an ImageJ macro that measures scanned images of washed root samples. IJ_Rhizo is open-source, platform-independent and offers a simple graphic user interface (GUI) for a main audience of non-programmer scientists. Being open-source based, it is also fully modifiable to accommodate the specific needs of the more computer-literate users. A comparison of IJ_Rhizo’s performance with that of the widely used commercial package WinRHIZO? is discussed.

Methods

We compared IJ_Rhizo’s performance with that of the commercial package WinRHIZO? using two sets of images, one comprising test-line images, the second consisting of images of root samples collected in the field. IJ_Rhizo and WinRHIZO? estimates were compared by means of correlation and regression analysis.

Results

IJ_Rhizo “Kimura” and WinRHIZO? “Tennant” were the length estimates that were best linearly correlated with each other. Correlation between average root diameter estimates was weaker, due to the sensitivity of this parameter to thresholding and filtering of image background noise.

Conclusions

Overall, IJ_Rhizo offers new opportunities for researchers who cannot afford the cost of commercial software packages to carry out automated measurement of scanned images of root samples, without sacrificing accuracy.     

The MERG Suite: Tools for discovering competencies and associated learning resources

Background

As the demands for competency-based education grow, the need for standards-based tools to allow for publishing and discovery of competency-based learning content is more pressing. This project focused on developing federated discovery services for competency-based medical e-learning content.

Methods

We built a tool suite for authoring and discovery of medical e-learning metadata. The end-user usability of the tool suite was evaluated through a web-based survey.

Results

The suite, implemented as an open-source system, was evaluated to identify areas for improvement.

Conclusion

The MERG suite is a starting point for organizations implementing competency-based e-learning resources.     

Flexible network reconstruction from relational databases with Cytoscape and CytoSQL

Background

Bisulfite sequencing using next generation sequencers yields genome-wide measurements of DNA methylation at single nucleotide resolution. Traditional aligners are not designed for mapping bisulfite-treated reads, where the unmethylated Cs are converted to Ts. We have developed BS Seeker, an approach that converts the genome to a three-letter alphabet and uses Bowtie to align bisulfite-treated reads to a reference genome. It uses sequence tags to reduce mapping ambiguity. Post-processing of the alignments removes non-unique and low-quality mappings.

Results

We tested our aligner on synthetic data, a bisulfite-converted Arabidopsis library, and human libraries generated from two different experimental protocols. We evaluated the performance of our approach and compared it to other bisulfite aligners. The results demonstrate that among the aligners tested, BS Seeker is more versatile and faster. When mapping to the human genome, BS Seeker generates alignments significantly faster than RMAP and BSMAP. Furthermore, BS Seeker is the only alignment tool that can explicitly account for tags which are generated by certain library construction protocols.

Conclusions

BS Seeker provides fast and accurate mapping of bisulfite-converted reads. It can work with BS reads generated from the two different experimental protocols, and is able to efficiently map reads to large mammalian genomes. The Python program is freely available at http://pellegrini.mcdb.ucla.edu/BS_Seeker/BS_Seeker.html.