The new uORFdb: integrating literature, sequence, and variation data in a central hub for uORF research.

Upstream open reading frames (uORFs) are initiated by AUG or near-cognate start codons and have been identified in the transcript leader sequences of the majority of eukaryotic transcripts. Functionally, uORFs are implicated in downstream translational regulation of the main protein coding sequence and may serve as a source of non-canonical peptides. Genetic defects in uORF sequences have been linked to the development of various diseases, including cancer. To simplify uORF-related research, the initial release of uORFdb in 2014 provided a comprehensive and manually curated collection of uORF-related literature. Here, we present an updated sequence-based version of uORFdb, accessible at https://www.bioinformatics.uni-muenster.de/tools/uorfdb. The new uORFdb enables users to directly access sequence information, graphical displays, and genetic variation data for over 2.4 million human uORFs. It also includes sequence data of >4.2 million uORFs in 12 additional species. Multiple uORFs can be displayed in transcript- and reading-frame-specific models to visualize the translational context. A variety of filters, sequence-related information, and links to external resources (UCSC Genome Browser, dbSNP, ClinVar) facilitate immediate in-depth analysis of individual uORFs. The database also contains uORF-related somatic variation data obtained from whole-genome sequencing (WGS) analyses of 677 cancer samples collected by the TCGA consortium.

NewickTreeModifier: A simple web tool to prune and modify Newick trees.

Large-scale selection analyses of protein-coding sequences and phylogenetic tree reconstructions require suitable trees in Newick format. We developed the NewickTreeModifier (NTM), a simple web-based tool to trim and modify Newick trees for such analyses. The users can choose provided master trees or upload a tree to prune it to selected species available in FASTA, NEXUS, or PHYLIP sequence format with an internal converter, a simple species list, or directly determined from a checklist interface of the master trees. Plant, insect, and vertebrate master trees comprise the maximum number of species in an up-to-date phylogenetic order directly transferable to the pruned Newick outfile. NTM is available at https://retrogenomics.uni-muenster.de/tools/ntm.

The multicomparative 2-n-way genome suite.

To effectively analyze the increasing amounts of available genomic data, improved comparative analytical tools that are accessible to and applicable by a broad scientific community are essential. We built the "2-n-way" software suite to provide a fundamental and innovative processing framework for revealing and comparing inserted elements among various genomes. The suite comprises two user-friendly web-based modules. The 2-way module generates pairwise whole-genome alignments of target and query species. The resulting genome coordinates of blocks (matching sequences) and gaps (missing sequences) from multiple 2-ways are then transferred to the n-way module and sorted into projects, in which user-defined coordinates from reference species are projected to the block/gap coordinates of orthologous loci in query species to provide comparative information about presence (blocks) or absence (gaps) patterns of targeted elements over many entire genomes and phylogroups. Thus, the 2-n-way software suite is ideal for performing multidirectional, non-ascertainment-biased screenings to extract all possible presence/absence data of user-relevant elements in orthologous sequences. To highlight its applicability and versatility, we used 2-n-way to expose approximately 100 lost introns in vertebrates, analyzed thousands of potential phylogenetically informative bat and whale retrotransposons, and novel human exons as well as thousands of human polymorphic retrotransposons.

GPAC-genome presence/absence compiler: a web application to comparatively visualize multiple genome-level changes.

Our understanding of genome-wide and comparative sequence information has been broadened considerably by the databases available from the University of California Santa Cruz (UCSC) Genome Bioinformatics Department. In particular, the identification and visualization of genomic sequences, present in some species but absent in others, led to fundamental insights into gene and genome evolution. However, the UCSC tools currently enable one to visualize orthologous genomic loci for a range of species in only a single locus. For large-scale comparative analyses of such presence/absence patterns a multilocus view would be more desirable. Such a tool would enable us to compare thousands of relevant loci simultaneously and to resolve many different questions about, for example, phylogeny, specific aspects of genome and gene evolution, such as the gain or loss of exons and introns, the emergence of novel transposed elements, nonprotein-coding RNAs, and viral genomic particles. Here, we present the first tool to facilitate the parallel analysis of thousands of genomic loci for cross-species presence/absence patterns based on multiway genome alignments. This genome presence/absence compiler uses annotated or other compilations of coordinates of genomic locations and compiles all presence/absence patterns in a flexible, color-coded table linked to the individual UCSC Genome Browser alignments. We provide examples of the versatile information content of such a screening system especially for 7SL-derived transposed elements, nuclear mitochondrial DNA, DNA transposons, and miRNAs in primates (http://www.bioinformatics.uni-muenster.de/tools/gpac, last accessed October 1, 2014).

Origin of the 1918 pandemic H1N1 influenza A virus as studied by codon usage patterns and phylogenetic analysis.

The pandemic of 1918 was caused by an H1N1 influenza A virus, which is a negative strand RNA virus; however, little is known about the nature of its direct ancestral strains. Here we applied a broad genetic and phylogenetic analysis of a wide range of influenza virus genes, in particular the PB1 gene, to gain information about the phylogenetic relatedness of the 1918 H1N1 virus. We compared the RNA genome of the 1918 strain to many other influenza strains of different origin by several means, including relative synonymous codon usage (RSCU), effective number of codons (ENC), and phylogenetic relationship. We found that the PB1 gene of the 1918 pandemic virus had ENC values similar to the H1N1 classical swine and human viruses, but different ENC values from avian as well as H2N2 and H3N2 human viruses. Also, according to the RSCU of the PB1 gene, the 1918 virus grouped with all human isolates and "classical" swine H1N1 viruses. The phylogenetic studies of all eight RNA gene segments of influenza A viruses may indicate that the 1918 pandemic strain originated from a H1N1 swine virus, which itself might be derived from a H1N1 avian precursor, which was separated from the bulk of other avian viruses in toto a long time ago. The high stability of the RSCU pattern of the PB1 gene indicated that the integrity of RNA structure is more important for influenza virus evolution than previously thought.

Co-expression and clinical utility of AR-FL and AR splice variants AR-V3, AR-V7 and AR-V9 in prostate cancer.

BACKGROUND: Androgen receptor (AR) splice variants (AR-Vs) have been discussed as a biomarker in prostate cancer (PC). However, some reports question the predictive property of AR-Vs. From a mechanistic perspective, the connection between AR full length (AR-FL) and AR-Vs is not fully understood. Here, we aimed to investigate the dependence of AR-FL and AR-V expression levels on AR gene activity. Additionally, we intended to comprehensively analyze presence of AR-FL and three clinically relevant AR-Vs (AR-V3, AR-V7 and AR-V9) in different stages of disease, especially with respect to clinical utility in PC patients undergoing AR targeted agent (ARTA) treatment. METHODS: AR-FL and AR-V levels were analyzed in PC and non-PC cell lines upon artificial increase of AR pre-mRNA using either drug treatment or AR gene activation. Furthermore, expression of AR-FL and AR-Vs was determined in PC specimen at distinct stages of disease (primary (n = 10) and metastatic tissues (n = 20), liquid biopsy samples (n = 422), mCRPC liquid biopsy samples of n = 96 patients starting novel treatment). Finally, baseline AR-FL and AR-V status was correlated with clinical outcome in a defined cohort of n = 65 mCRPC patients undergoing ARTA treatment. RESULTS: We revealed rising levels of AR-FL accompanied with appearance and increase of AR-Vs in dependence of elevated AR pre-mRNA levels. We also noticed increase in AR-FL and AR-V levels throughout disease progression. AR-V expression was always associated with high AR-FL levels without any sample being solely AR-V positive. In patients undergoing ARTA treatment, AR-FL did show prognostic, yet not predictive validity. Additionally, we observed a substantial clinical response to ARTA treatment even in AR-V positive patients. Accordingly, multivariate analysis did not demonstrate independent significance of AR-Vs in neither predictive nor prognostic clinical utility. CONCLUSION: We demonstrate a correlation between AR-FL and AR-V expression during PC progression; with AR-V expression being a side-effect of elevated AR pre-mRNA levels. Clinically, AR-V positivity relies on high levels of AR-FL, making cells still vulnerable to ARTA treatment, as demonstrated by AR-FL and AR-V positive patients responding to ARTA treatment. Thus, AR-FL and AR-V might be considered as a prognostic, yet not predictive biomarker in mCRPC patients.

paPAML: An Improved Computational Tool to Explore Selection Pressure on Protein-Coding Sequences.

Evolution is change over time. Although neutral changes promoted by drift effects are most reliable for phylogenetic reconstructions, selection-relevant changes are of only limited use to reconstruct phylogenies. On the other hand, comparative analyses of neutral and selected changes of protein-coding DNA sequences (CDS) retrospectively tell us about episodic constrained, relaxed, and adaptive incidences. The ratio of sites with nonsynonymous (amino acid altering) versus synonymous (not altering) mutations directly measures selection pressure and can be analysed by using the Phylogenetic Analysis by Maximum Likelihood (PAML) software package. We developed a CDS extractor for compiling protein-coding sequences (CDS-extractor) and parallel PAML (paPAML) to simplify, amplify, and accelerate selection analyses via parallel processing, including detection of negatively selected sites. paPAML compiles results of site, branch-site, and branch models and detects site-specific negative selection with the output of a codon list labelling significance values. The tool simplifies selection analyses for casual and inexperienced users and accelerates computing speeds up to the number of allocated computer threads. We then applied paPAML to examine the evolutionary impact on a new GINS Complex Subunit 3 exon, and neutrophil-associated as well as lysin and apolipoprotein genes. Compared with codeml (PAML version 4.9j) and HyPhy (HyPhy FEL version 2.5.26), all paPAML test runs performed with 10 computing threads led to identical selection pressure results, whereas the total selection analysis via paPAML, including all model comparisons, was about 3 to 5 times faster than the longest running codeml model and about 7 to 15 times faster than the entire processing time of these codeml runs.

Comparison of circulating tumor cells and AR-V7 as clinical biomarker in metastatic castration-resistant prostate cancer patients.

Biomarker in metastatic castration resistant prostate cancer (mCRPC) treatment are rare. We aimed to compare the clinical value of circulating tumor cells (CTCs) and androgen receptor splice variant 7 (AR-V7) as biomarker in mCRPC patients undergoing androgen receptor-targeted agent (ARTA) treatment. Overall cohort (65 patients) was stratified regarding either CTC or AR-V7 status followed by further sub-stratification of the respective other marker. Subsequently, prostate specific antigen (PSA) response, progression free survival (PFS) and overall survival (OS)) of subgroups was compared. CTCs and AR-V7 were detected in 54 (83%) and 33 (61%) patients, respectively. All AR-V7 + were CTC +. We detected PSA response in all subgroups. For PFS and OS, biomarker stratification revealed differences between all subgroups. Interestingly, no significant differences of AR-V7 transcript copy numbers were detected between responding and non-responding patients. Additionally, multivariable analysis revealed no independent prognostic value of AR-V7 positivity. Both biomarkers show clinical value in prognosticating clinical outcome. Nonetheless, AR-V7 stratification underestimates the heterogenous subgroup of CTC - and CTC + patient, the latter requiring more intense clinical surveillance. Additionally, AR-V7 level does not correlate with clinical response. Thus, the value of AR-V7 as a clinical biomarker must be considered skeptically.

A novel web-based TinT application and the chronology of the Primate Alu retroposon activity.

BACKGROUND: DNA sequences afford access to the evolutionary pathways of life. Particularly mobile elements that constantly co-evolve in genomes encrypt recent and ancient information of their host's history. In mammals there is an extraordinarily abundant activity of mobile elements that occurs in a dynamic succession of active families, subfamilies, types, and subtypes of retroposed elements. The high frequency of retroposons in mammals implies that, by chance, such elements also insert into each other. While inactive elements are no longer able to retropose, active elements retropose by chance into other active and inactive elements. Thousands of such directional, element-in-element insertions are found in present-day genomes. To help analyze these events, we developed a computational algorithm (Transpositions in Transpositions, or TinT) that examines the different frequencies of nested transpositions and reconstructs the chronological order of retroposon activities. RESULTS: By examining the different frequencies of such nested transpositions, the TinT application reconstructs the chronological order of retroposon activities. We use such activity patterns as a comparative tool to (1) delineate the historical rise and fall of retroposons and their relations to each other, (2) understand the retroposon-induced complexity of recent genomes, and (3) find selective informative homoplasy-free markers of phylogeny. The efficiency of the new application is demonstrated by applying it to dimeric Alu Short INterspersed Elements (SINE) to derive a complete chronology of such elements in primates. CONCLUSION: The user-friendly, web-based TinT interface presented here affords an easy, automated screening for nested transpositions from genome assemblies or trace data, assembles them in a frequency-matrix, and schematically displays their chronological activity history.

TEclass--a tool for automated classification of unknown eukaryotic transposable elements.

MOTIVATION: The large number of sequenced genomes required the development of software that reconstructs the consensus sequences of transposons and other repetitive elements. However, the available tools usually focus on the accurate identification of raw repeats and provide no information about the taxonomic position of the reconstructed consensi. TEclass is a tool to classify unknown transposable elements into their four main functional categories, which reflect their mode of transposition: DNA transposons, long terminal repeats (LTRs), long interspersed nuclear elements (LINEs) and short interspersed nuclear elements (SINEs). TEclass uses machine learning support vector machine (SVM) for classification based on oligomer frequencies. It achieves 90-97% accuracy in the classification of novel DNA and LTR repeats, and 75% for LINEs and SINEs. AVAILABILITY: http://www.compgen.uni-muenster.de/teclass, stand alone program upon request.

NanoPipe-a web server for nanopore MinION sequencing data analysis.

BACKGROUND: The fast-moving progress of the third-generation long-read sequencing technologies will soon bring the biological and medical sciences to a new era of research. Altogether, the technique and experimental procedures are becoming more straightforward and available to biologists from diverse fields, even without any profound experience in DNA sequencing. Thus, the introduction of the MinION device by Oxford Nanopore Technologies promises to "bring sequencing technology to the masses" and also allows quick and operative analysis in field studies. However, the convenience of this sequencing technology dramatically contrasts with the available analysis tools, which may significantly reduce enthusiasm of a "regular" user. To really bring the sequencing technology to every biologist, we need a set of user-friendly tools that can perform a powerful analysis in an automatic manner. FINDINGS: NanoPipe was developed in consideration of the specifics of the MinION sequencing technologies, providing accordingly adjusted alignment parameters. The range of the target species/sequences for the alignment is not limited, and the descriptive usage page of NanoPipe helps a user to succeed with NanoPipe analysis. The results contain alignment statistics, consensus sequence, polymorphisms data, and visualization of the alignment. Several test cases are used to demonstrate the efficiency of the tool. CONCLUSIONS: Freely available NanoPipe software allows effortless and reliable analysis of MinION sequencing data for experienced bioinformaticians, as well for wet-lab biologists with minimum bioinformatics knowledge. Moreover, for the latter group, we describe the basic algorithm necessary for MinION sequencing analysis from the first to last step.