We present here a new version of the Arlequin program available under three different forms: a Windows graphical version (Winarl35), a console version of Arlequin (arlecore), and a specific console version to compute summary statistics (arlsumstat). The command-line versions run under both Linux and Windows. The main innovations of the new version include enhanced outputs in XML format, the possibility to embed graphics displaying computation results directly into output files, and the implementation of a new method to detect loci under selection from genome scans. Command-line versions are designed to handle large series of files, and arlsumstat can be used to generate summary statistics from simulated data sets within an Approximate Bayesian Computation framework.

PGDSpider is freely available under the BSD 3-Clause license on http://cmpg.unibe.ch/software/PGDSpider/.

BackgroundThe purpose of gene set enrichment analysis (GSEA) is to find general trends in the huge lists of genes or proteins generated by many functional genomics techniques and bioinformatics analyses.ResultsHere we present SetRank, an advanced GSEA algorithm which is able to eliminate many false positive hits. The key principle of the algorithm is that it discards gene sets that have initially been flagged as significant, if their significance is only due to the overlap with another gene set. The algorithm is explained in detail and its performance is compared to that of other methods using objective benchmarking criteria. Furthermore, we explore how sample source bias can affect the results of a GSEA analysis.ConclusionsThe benchmarking results show that SetRank is a highly specific tool for GSEA. Furthermore, we show that the reliability of results can be improved by taking sample source bias into account. SetRank is available as an R package and through an online web interface.

BackgroundThe development of next-generation sequencing has made it possible to sequence whole genomes at a relatively low cost. However, de novo genome assemblies remain challenging due to short read length, missing data, repetitive regions, polymorphisms and sequencing errors. As more and more genomes are sequenced, reference-guided assembly approaches can be used to assist the assembly process. However, previous methods mostly focused on the assembly of other genotypes within the same species. We adapted and extended a reference-guided de novo assembly approach, which enables the usage of a related reference sequence to guide the genome assembly. In order to compare and evaluate de novo and our reference-guided de novo assembly approaches, we used a simulated data set of a repetitive and heterozygotic plant genome.ResultsThe extended reference-guided de novo assembly approach almost always outperforms the corresponding de novo assembly program even when a reference of a different species is used. Similar improvements can be observed in high and low coverage situations. In addition, we show that a single evaluation metric, like the widely used N50 length, is not enough to properly rate assemblies as it not always points to the best assembly evaluated with other criteria. Therefore, we used the summed z-scores of 36 different statistics to evaluate the assemblies.ConclusionsThe combination of reference mapping and de novo assembly provides a powerful tool to improve genome reconstruction by integrating information of a related genome. Our extension of the reference-guided de novo assembly approach enables the application of this strategy not only within but also between related species. Finally, the evaluation of genome assemblies is often not straight forward, as the truth is not known. Thus one should always use a combination of evaluation metrics, which not only try to assess the continuity but also the accuracy of an assembly.Electronic supplementary materialThe online version of this article (10.1186/s12859-017-1911-6) contains supplementary material, which is available to authorized users.

Myeloproliferative neoplasms (MPN) show dysregulated JAK2 signaling. JAK2 inhibitors provide clinical benefits, but compensatory activation of MAPK pathway signaling impedes efficacy. We hypothesized that dual targeting of JAK2 and ERK1/2 could enhance clone control and therapeutic efficacy. We employed genetic and pharmacologic targeting of ERK1/2 in Jak2V617F MPN mice, cells and patient clinical isolates. Competitive transplantations of Jak2V617F vs. wild-type bone marrow (BM) showed that ERK1/2 deficiency in hematopoiesis mitigated MPN features and reduced the Jak2V617F clone in blood and hematopoietic progenitor compartments. ERK1/2 ablation combined with JAK2 inhibition suppressed MAPK transcriptional programs, normalized cytoses and promoted clone control suggesting dual JAK2/ERK1/2 targeting as enhanced corrective approach. Combined pharmacologic JAK2/ERK1/2 inhibition with ruxolitinib and ERK inhibitors reduced proliferation of Jak2V617F cells and corrected erythrocytosis and splenomegaly of Jak2V617F MPN mice. Longer-term treatment was able to induce clone reductions. BM fibrosis was significantly decreased in MPLW515L-driven MPN to an extent not seen with JAK2 inhibitor monotherapy. Colony formation from JAK2V617F patients’ CD34+ blood and BM was dose-dependently inhibited by combined JAK2/ERK1/2 inhibition in PV, ET, and MF subsets. Overall, we observed that dual targeting of JAK2 and ERK1/2 was able to enhance therapeutic efficacy suggesting a novel treatment approach for MPN.

The SIB Swiss Institute of Bioinformatics (www.isb-sib.ch) provides world-class bioinformatics databases, software tools, services and training to the international life science community in academia and industry. These solutions allow life scientists to turn the exponentially growing amount of data into knowledge. Here, we provide an overview of SIB's resources and competence areas, with a strong focus on curated databases and SIB's most popular and widely used resources. In particular, SIB's Bioinformatics resource portal ExPASy features over 150 resources, including UniProtKB/Swiss-Prot, ENZYME, PROSITE, neXtProt, STRING, UniCarbKB, SugarBindDB, SwissRegulon, EPD, arrayMap, Bgee, SWISS-MODEL Repository, OMA, OrthoDB and other databases, which are briefly described in this article.

SummaryIn domestic goats, the polled intersex syndrome (PIS) refers to XX female‐to‐male sex reversal associated with the absence of horn growth (polled). The causal variant was previously reported as a 11.7 kb deletion at approximately 129 Mb on chromosome 1 that affects the transcription of both FOXL2 and several long non‐coding RNAs. In the meantime the presence of different versions of the PIS deletion was postulated and trials to establish genetic testing with the existing molecular genetic information failed. Therefore, we revisited this variant by long‐read whole‐genome sequencing of two genetically female (XX) goats, a PIS‐affected and a horned control. This revealed the presence of a more complex structural variant consisting of a deletion with a total length of 10 159 bp and an inversely inserted approximately 480 kb‐sized duplicated segment of a region located approximately 21 Mb further downstream on chromosome 1 containing two genes, KCNJ15 and ERG. Publicly available short‐read whole‐genome sequencing data, Sanger sequencing of the breakpoints and FISH using BAC clones corresponding to both involved genome regions confirmed this structural variant. A diagnostic PCR was developed for simultaneous genotyping of carriers for this variant and determination of their genetic sex. We showed that the variant allele was present in all 334 genotyped polled goats of diverse breeds and that all analyzed 15 PIS‐affected XX goats were homozygous. Our findings enable for the first time a precise genetic diagnosis for polledness and PIS in goats and add a further genomic feature to the complexity of the PIS phenomenon.

Dendritic cells (DC) and monocytes are vital for the initiation of innate and adaptive immune responses. Recently, we identified bona fide DC subsets in blood of cattle, revealing subset‐ and species‐specific transcription of toll‐like receptors (TLR). In the present study, we analyzed phenotypic and transcriptional responses of bovine DC subsets and monocytes to in vitro stimulation with four to six different TLR ligands. Bovine DC subsets, especially plasmacytoid DC (pDC), showed a clear increase of CCR7, CD25, CD40, CD80, CD86, and MHC‐II expression both on mRNA and protein level. Flow cytometric detection of p38 MAPK phosphorylation 15 min after stimulation confirmed activation of DC subsets and monocytes in accordance with TLR gene expression. Whole‐transcriptome sequencing of sorted and TLR‐stimulated subsets revealed potential ligand‐ and subset‐specific regulation of genes associated with inflammation, T‐cell co‐stimulation, migration, metabolic reprogramming, and antiviral activity. Gardiquimod was found to evoke strong responses both in DC subsets and monocytes, while Poly(I:C) and CpG preferentially triggered responses in cDC1 and pDC, respectively. This in‐depth analysis of ligand responsiveness is essential for the rational design of vaccine adjuvants in cattle, and provides a solid basis for comparative studies on DC and monocyte biology across species.