E‐Clic – easy climate data converter

There are an increasing number of studies that are now focussing on the influence of climate change on species’ distributions. However, access to predictive climatic datasets for future scenarios is difficult due to their specific formats and/or the need to be geographically downscaled. The TYN dataset is freely available to users and provides a synthetic format with several climatic models and IPCC future climate scenarios. Moreover, the CRU historical dataset (1901–2000) is also available which allows users to create baseline models for current climatic variables. E‐clic is a free, user‐friendly software package that offers three different ways to convert these two datasets into a spatially explicit raster format which is compatible with the most common geographic information systems and usable on different platforms.     

Comparison of VHH‐Fc antibody production in Arabidopsis thaliana,Nicotiana benthamiana and Pichia pastoris

VHHs or nanobodies are widely acknowledged as interesting diagnostic and therapeutic tools. However, for some applications, multivalent antibody formats, such as the dimeric VHH‐Fc format, are desired to increase the functional affinity. The scope of this study was to compare transient expression of diagnostic VHH‐Fc antibodies in Nicotiana benthamiana leaves with their stable expression in Arabidopsis thaliana seeds and Pichia pastoris. To this end, VHH‐Fc antibodies targeting green fluorescent protein or the A. thaliana seed storage proteins (albumin and globulin) were produced in the three platforms. Differences were mainly observed in the accumulation levels and glycosylation patterns. Interestingly, although in plants oligomannosidic N‐glycans were expected for KDEL‐tagged VHH‐Fcs, several VHH‐Fcs with an intact KDEL‐tag carried complex‐type N‐glycans, suggesting a dysfunctional retention in the endoplasmic reticulum. All VHH‐Fcs were equally functional across expression platforms and several outperformed their corresponding VHH in terms of sensitivity in ELISA.     

Double‐digest RAD sequencing using Ion Proton semiconductor platform (ddRADseq‐ion) with nonmodel organisms

Research in evolutionary biology involving nonmodel organisms is rapidly shifting from using traditional molecular markers such as mtDNA and microsatellites to higher throughput SNP genotyping methodologies to address questions in population genetics, phylogenetics and genetic mapping. Restriction site associated DNA sequencing (RAD sequencing or RADseq) has become an established method for SNP genotyping on Illumina sequencing platforms. Here, we developed a protocol and adapters for double‐digest RAD sequencing for Ion Torrent (Life Technologies; Ion Proton, Ion PGM) semiconductor sequencing. We sequenced thirteen genomic libraries of three different nonmodel vertebrate species on Ion Proton with PI chips: Arctic charr Salvelinus alpinus, European whitefish Coregonus lavaretus and common lizard Zootoca vivipara. This resulted in ~962 million single‐end reads overall and a mean of ~74 million reads per library. We filtered the genomic data using Stacks, a bioinformatic tool to process RAD sequencing data. On average, we obtained ~11 000 polymorphic loci per library of 6–30 individuals. We validate our new method by technical and biological replication, by reconstructing phylogenetic relationships, and using a hybrid genetic cross to track genomic variants. Finally, we discuss the differences between using the different sequencing platforms in the context of RAD sequencing, assessing possible advantages and disadvantages. We show that our protocol can be used for Ion semiconductor sequencing platforms for the rapid and cost‐effective generation of variable and reproducible genetic markers.     

jackalope: A swift,versatile phylogenomic and high‐throughput sequencing simulator

High‐throughput sequencing (HTS) is central to the study of population genomics and has an increasingly important role in constructing phylogenies. Choices in research design for sequencing projects can include a wide range of factors, such as sequencing platform, depth of coverage and bioinformatic tools. Simulating HTS data better informs these decisions, as users can validate software by comparing output to the known simulation parameters. However, current standalone HTS simulators cannot generate variant haplotypes under even somewhat complex evolutionary scenarios, such as recombination or demographic change. This greatly reduces their usefulness for fields such as population genomics and phylogenomics. Here I present the R package jackalope that simply and efficiently simulates (i) sets of variant haplotypes from a reference genome and (ii) reads from both Illumina and Pacific Biosciences platforms. Haplotypes can be simulated using phylogenies, gene trees, coalescent‐simulation output, population‐genomic summary statistics, and Variant Call Format (VCF) files. jackalope can simulate single, paired‐end or mate‐pair Illumina reads, as well as reads from Pacific Biosciences. These simulations include sequencing errors, mapping qualities, multiplexing and optical/PCR duplicates. It can read reference genomes from fasta files and can simulate new ones, and all outputs can be written to standard file formats. jackalope is available for Mac, Windows and Linux systems.     

The Quantification of Representative Sequences pipeline for amplicon sequencing: case study on within‐population ITS1 sequence variation in a microparasite infecting Daphnia

Next generation sequencing (NGS) platforms are replacing traditional molecular biology protocols like cloning and Sanger sequencing. However, accuracy of NGS platforms has rarely been measured when quantifying relative frequencies of genotypes or taxa within populations. Here we developed a new bioinformatic pipeline (QRS) that pools similar sequence variants and estimates their frequencies in NGS data sets from populations or communities. We tested whether the estimated frequency of representative sequences, generated by 454 amplicon sequencing, differs significantly from that obtained by Sanger sequencing of cloned PCR products. This was performed by analysing sequence variation of the highly variable first internal transcribed spacer (ITS1) of the ichthyosporean Caullerya mesnili, a microparasite of cladocerans of the genus Daphnia. This analysis also serves as a case example of the usage of this pipeline to study within‐population variation. Additionally, a public Illumina data set was used to validate the pipeline on community‐level data. Overall, there was a good correspondence in absolute frequencies of C. mesnili ITS1 sequences obtained from Sanger and 454 platforms. Furthermore, analyses of molecular variance (amova ) revealed that population structure of C. mesnili differs across lakes and years independently of the sequencing platform. Our results support not only the usefulness of amplicon sequencing data for studies of within‐population structure but also the successful application of the QRS pipeline on Illumina‐generated data. The QRS pipeline is freely available together with its documentation under GNU Public Licence version 3 at http://code.google.com/p/quantification-representative-sequences .     

PipeCraft: Flexible open‐source toolkit for bioinformatics analysis of custom high‐throughput amplicon sequencing data

High‐throughput sequencing methods have become a routine analysis tool in environmental sciences as well as in public and private sector. These methods provide vast amount of data, which need to be analysed in several steps. Although the bioinformatics may be applied using several public tools, many analytical pipelines allow too few options for the optimal analysis for more complicated or customized designs. Here, we introduce PipeCraft, a flexible and handy bioinformatics pipeline with a user‐friendly graphical interface that links several public tools for analysing amplicon sequencing data. Users are able to customize the pipeline by selecting the most suitable tools and options to process raw sequences from Illumina, Pacific Biosciences, Ion Torrent and Roche 454 sequencing platforms. We described the design and options of PipeCraft and evaluated its performance by analysing the data sets from three different sequencing platforms. We demonstrated that PipeCraft is able to process large data sets within 24 hr. The graphical user interface and the automated links between various bioinformatics tools enable easy customization of the workflow. All analytical steps and options are recorded in log files and are easily traceable.     

SimRAD: an R package for simulation‐based prediction of the number of loci expected in RADseq and similar genotyping by sequencing approaches

Application of high‐throughput sequencing platforms in the field of ecology and evolutionary biology is developing quickly with the introduction of efficient methods to reduce genome complexity. Numerous approaches for genome complexity reduction have been developed using different combinations of restriction enzymes, library construction strategies and fragment size selection. As a result, the choice of which techniques to use may become cumbersome, because it is difficult to anticipate the number of loci resulting from each method. We developed SimRAD, an R package that performs in silico restriction enzyme digests and fragment size selection as implemented in most restriction site associated DNA polymorphism and genotyping by sequencing methods. In silico digestion is performed on a reference genome or on a randomly generated DNA sequence when no reference genome sequence is available. SimRAD accurately predicts the number of loci under alternative protocols when a reference genome sequence is available for the targeted species (or a close relative) but may be unreliable when no reference genome is available. SimRAD is also useful for fine‐tuning a given protocol to adjust the number of targeted loci. Here, we outline the functionality of SimRAD and provide an illustrative example of the use of the package (available on the CRAN at http://cran.r-project.org/web/packages/SimRAD ).     

multiplierz v2.0: A Python‐based ecosystem for shared access and analysis of native mass spectrometry data

The continued evolution of modern mass spectrometry instrumentation and associated methods represents a critical component in efforts to decipher the molecular mechanisms which underlie normal physiology and understand how dysregulation of biological pathways contributes to human disease. The increasing scale of these experiments combined with the technological diversity of mass spectrometers presents several challenges for community‐wide data access, analysis, and distribution. Here we detail a redesigned version of multiplierz, our Python software library which leverages our common application programming interface (mzAPI) for analysis and distribution of proteomic data. New features include support for a wider range of native mass spectrometry file types, interfaces to additional database search engines, compatibility with new reporting formats, and high‐level tools to perform post‐search proteomic analyses. A GUI desktop environment, mzDesktop, provides access to multiplierz functionality through a user friendly interface. multiplierz is available for download from: https://github.com/BlaisProteomics/multiplierz ; and mzDesktop is available for download from: https://sourceforge.net/projects/multiplierz/     

FERONIA receptor kinase interacts with S‐adenosylmethionine synthetase and suppresses S‐adenosylmethionine production and ethylene biosynthesis in Arabidopsis

Environmental inputs such as stress can modulate plant cell metabolism, but the detailed mechanism remains unclear. We report here that FERONIA (FER), a plasma membrane receptor‐like kinase, may negatively regulate the S‐adenosylmethionine (SAM) synthesis by interacting with two S‐adenosylmethionine synthases (SAM1 and SAM2). SAM participates in ethylene, nicotianamine and polyamine biosynthetic pathways and provides the methyl group for protein and DNA methylation reactions. The Arabidopsis fer mutants contained a higher level of SAM and ethylene in plant tissues and displayed a dwarf phenotype. Such phenotype in the fer mutants was mimicked by over‐expressing the S‐adenosylmethionine synthetase in transgenic plants, whereas sam1/2 double mutant showed an opposite phenotype. We propose that FER receptor kinase, in response to environmental stress and plant hormones such as auxin and BR, interacts with SAM synthases and down‐regulates ethylene biosynthesis.     

Getting started in mapping-by-sequencing

Next-generation sequencing (NGS) technologies al ow the cost-effective sequencing of whole genomes and have expanded the scope of genomics to novel applications, such as the genome-wide characterizatio...     

Increased T‐cell stimulating activity by mutated SEC2 correlates with its improved antitumour potency

Aims: To investigate the improved antitumour activity of SAM‐3 compared with recombinant staphylococcal enterotoxins C2 (rSEC2). Methods and Results: Methylthiazol tetrazolium and flow cytometry assays showed that the antitumour activity of SAM‐3 in vivo was improved because of enhanced T‐cell stimulating potency, resulting in massive activation of T cells, particularly CD4 ⁺ and CD8 ⁺ T cells, and subsequent cytokine release. Quantitative real‐time PCR assay showed that despite similar Vβ specificities induced by rSEC2 and SAM‐3, the quantities of activated T cells bearing specific Vβin vitro were different. Conclusions: The results strongly suggested that the increased SAM‐3–T‐cell receptor (TCR) binding affinity contributed to massive T‐cell activation and cytokine release, substantially amplifying antitumour immune response in vivo. Significance and Impact of the Study: This study provided evidence for the mechanism of SAM‐3 antitumour activity improvement compared with rSEC2. Results indicated that SAM‐3 could be used as a potent powerful candidate agent for tumour treatment in clinics.     

gmconvert: file conversion for genemapper output files

gmconvert is a platform‐independent program provided in GUI (for Apple OS X and Windows XP) and command‐line versions (for other platforms). gmconvert allows rapid reformatting of microsatellite data from output files produced by Applied Biosystems genemapper software (version 3.x). The program will re‐array data into three formats commonly used in downstream analysis: genepop , cervus , and gerud . gmconvert will greatly increase the speed of data preparation prior to analysis and aid in reducing transpositional errors associated with manual re‐arraying and reformatting steps. gmconvert is available from http://gallus.forestry.uga.edu/software/ .     

Conversion of array‐based single nucleotide polymorphic markers for use in targeted genotyping by sequencing in hexaploid wheat (Triticum aestivum)

Wheat breeders and academics alike use single nucleotide polymorphisms (SNP s) as molecular markers to characterize regions of interest within the hexaploid wheat genome. A number of SNP ‐based genotyping platforms are available, and their utility depends upon factors such as the available technologies, number of data points required, budgets and the technical expertise required. Unfortunately, markers can rarely be exchanged between existing and newly developed platforms, meaning that previously generated data cannot be compared, or combined, with more recently generated data sets. We predict that genotyping by sequencing will become the predominant genotyping technology within the next 5–10 years. With this in mind, to ensure that data generated from current genotyping platforms continues to be of use, we have designed and utilized SNP ‐based capture probes from several thousand existing and publicly available probes from Axiom® and KASP ? genotyping platforms. We have validated our capture probes in a targeted genotyping by sequencing protocol using 31 previously genotyped UK elite hexaploid wheat accessions. Data comparisons between targeted genotyping by sequencing, Axiom® array genotyping and KASP ? genotyping assays, identified a set of 3256 probes which reliably bring together targeted genotyping by sequencing data with the previously available marker data set. As such, these probes are likely to be of considerable value to the wheat community. The probe details, full probe sequences and a custom built analysis pipeline may be freely downloaded from the CerealsDB website (http://www.cerealsdb.uk.net/cerealgenomics/CerealsDB /sequence_capture.php).     

Assessing the importance of artificial nest‐sites in the population dynamics of endangered Northern Aplomado Falcons Falco femoralis septentrionalis in South Texas using stochastic simulation models

Habitat availability might be the most important determinant of success for a species reintroduction programme, making investigation of the quality and quantity of habitat needed to produce self‐sustaining populations a research priority for reintroduction ecologists. We used a stochastic model of population dynamics to predict whether attempts to improve existing breeding territories using artificial nest platforms improved the population growth rate and persistence of a reintroduced population of Northern Aplomado Falcons Falco femoralis septentrionalis in South Texas. We further assessed whether the creation of new territories, i.e. conversion of entire areas to suitable habitat and not simply the erection of nest platforms, would lead to a subsequent increase in the nesting population. Our model was able to reproduce several characteristics of the wild population and predicted the number of breeding pairs per year strikingly well (R² = 0.97). Simulations revealed that the addition of nest platforms improved productivity such that the population would decline to extinction without them but is stable since their installation. Moreover, the model predicted that the increase in productivity due to nest platforms would cause the population to saturate available breeding territories, at which point the population would contain a moderate proportion of non‐territorial birds that could occupy territories if new ones become available. Population size would therefore be proportional to the increase in available territories. Our study demonstrates that artificial nest‐sites can be an effective tool for the management of reintroduced species.     

An aminoacylase activity from Streptomyces ambofaciens catalyzes the acylation of lysine on α‐position and peptides on N‐terminal position

The presence of aminoacylase activities was investigated in a crude extract of Streptomyces ambofaciens ATCC23877. First activities catalyzing the hydrolysis of N‐α or ε‐acetyl‐L‐lysine were identified. Furthermore, the acylation of lysine and different peptides was studied and compared with results obtained with lipase B of Candida antarctica (CALB). Different regioselectivities were demonstrated for the two classes of enzymes. CALB was able to catalyze acylation only on the ε‐position whereas the crude extract from S. ambofaciens possessed the rare ability to catalyze the N‐acylation on the α‐position of the lysine or of the amino‐acid in N‐terminal position of peptides. Two genes, SAM23877_1485 and SAM23877_1734, were identified in the genome of Streptomyces ambofaciens ATCC23877 whose products show similarities with the previously identified aminoacylases from Streptomyces mobaraensis. The proteins encoded by these two genes were responsible for the major aminoacylase hydrolytic activities. Furthermore, we show that the hydrolysis of N‐α‐acetyl‐L‐lysine could be attributed to the product of SAM23877_1734 gene.     

Evaluation of gene expression measurements from commercial microarray platforms

Multiple commercial microarrays for measuring genome-wide gene expression levels are currently available, including oligonucleotide and cDNA, single- and two-channel formats. This study reports on the results of gene expression measurements generated from identical RNA preparations that were obtained using three commercially available microarray platforms. RNA was collected from PANC-1 cells grown in serum-rich medium and at 24 h following the removal of serum. Three biological replicates were prepared for each condition, and three experimental replicates were produced for the first biological replicate. RNA was labeled and hybridized to microarrays from three major suppliers according to manufacturers’ protocols, and gene expression measurements were obtained using each platform’s standard software. For each platform, gene targets from a subset of 2009 common genes were compared. Correlations in gene expression levels and comparisons for significant gene expression changes in this subset were calculated, and showed considerable divergence across the different platforms, suggesting the need for establishing industrial manufacturing standards, and further independent and thorough validation of the technology.     

Reconstructing A/B compartments as revealed by Hi-C using long-range correlations in epigenetic data

Analysis of Hi-C data has shown that the genome can be divided into two compartments called A/B compartments. These compartments are cell-type specific and are associated with open and closed chromatin. We show that A/B compartments can reliably be estimated using epigenetic data from several different platforms: the Illumina 450 k DNA methylation microarray, DNase hypersensitivity sequencing, single-cell ATAC sequencing and single-cell whole-genome bisulfite sequencing. We do this by exploiting that the structure of long-range correlations differs between open and closed compartments. This work makes A/B compartment assignment readily available in a wide variety of cell types, including many human cancers.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0741-y) contains supplementary material, which is available to authorized users.