Noisy: Identification of problematic columns in multiple sequence alignments

<p>Abstract</p> <p>Motivation</p> <p>Sequence-based methods for phylogenetic reconstruction from (nucleic acid) sequence data are notoriously plagued by two effects: homoplasies and alignment errors. Large evolutionary distances imply a large number of homoplastic sites. As most protein-coding genes show dramatic variations in substitution rates that are not uncorrelated across the sequence, this often leads to a patchwork pattern of (i) phylogenetically informative and (ii) effectively randomized regions. In highly variable regions, furthermore, alignment errors accumulate resulting in sometimes misleading signals in phylogenetic reconstruction.</p> <p>Results</p> <p>We present here a method that, based on assessing the distribution of character states along a cyclic ordering of the taxa, allows the identification of phylogenetically uninformative homoplastic sites in a multiple sequence alignment. Removal of these sites appears to improve the performance of phylogenetic reconstruction algorithms as measured by various indices of "tree quality". In particular, we obtain more stable trees due to the exclusion of phylogenetically incompatible sites that most likely represent strongly randomized characters.</p> <p>Software</p> <p>The computer program noisy implements this approach. It can be employed to improving phylogenetic reconstruction capability with quite a considerable success rate whenever (1) the average bootstrap support obtained from the original alignment is low, and (2) there are sufficiently many taxa in the data set – at least, say, 12 to 15 taxa. The software can be obtained under the GNU Public License from <url>http://www.bioinf.uni-leipzig.de/Software/noisy/</url>.</p

Penile metastasis from primary transitional cell carcinoma of the renal pelvis: first manifestation of systemic spread

BACKGROUND: Almost one-third of all penile metastases are detected at the same time as a primary tumor, whereas the remaining two-thirds are detected a mean of 18 months after the discovery of the primary tumor. Cutaneous metastasis of transitional cell carcinoma (TCC) is extremely rare and generally accepted as the late manifestation of a systemic spread. CASE PRESENTATION: We report the first case of simultaneous penile and lung metastases from a primary TCC of the renal pelvis in a 76-year-old man, that occurred 8 years after a left nephroureterectomy. CONCLUSIONS: This case report underscores the importance of physical examinations of the skin of patients who undergo surgical procedures for TCC from bladder as well as from the upper urinary tract, including those seemingly without metastatic disease, because of the possibility of skin and penile metastatic spread

Combining estimates of interest in prognostic modelling studies after multiple imputation: current practice and guidelines

Background: Multiple imputation (MI) provides an effective approach to handle missing covariate data within prognostic modelling studies, as it can properly account for the missing data uncertainty. The multiply imputed datasets are each analysed using standard prognostic modelling techniques to obtain the estimates of interest. The estimates from each imputed dataset are then combined into one overall estimate and variance, incorporating both the within and between imputation variability. Rubin's rules for combining these multiply imputed estimates are based on asymptotic theory. The resulting combined estimates may be more accurate if the posterior distribution of the population parameter of interest is better approximated by the normal distribution. However, the normality assumption may not be appropriate for all the parameters of interest when analysing prognostic modelling studies, such as predicted survival probabilities and model performance measures. Methods: Guidelines for combining the estimates of interest when analysing prognostic modelling studies are provided. A literature review is performed to identify current practice for combining such estimates in prognostic modelling studies. Results: Methods for combining all reported estimates after MI were not well reported in the current literature. Rubin's rules without applying any transformations were the standard approach used, when any method was stated. Conclusion: The proposed simple guidelines for combining estimates after MI may lead to a wider and more appropriate use of MI in future prognostic modelling studies

Molecular evolution of HoxA13 and the multiple origins of limbless morphologies in amphibians and reptiles

Developmental processes and their results, morphological characters, are inherited through transmission of genes regulating development. While there is ample evidence that cis-regulatory elements tend to be modular, with sequence segments dedicated to different roles, the situation for proteins is less clear, being particularly complex for transcription factors with multiple functions. Some motifs mediating protein-protein interactions may be exclusive to particular developmental roles, but it is also possible that motifs are mostly shared among different processes. Here we focus on HoxA13, a protein essential for limb development. We asked whether the HoxA13 amino acid sequence evolved similarly in three limbless clades: Gymnophiona, Amphisbaenia and Serpentes. We explored variation in ω (dN/dS) using a maximum-likelihood framework and HoxA13sequences from 47 species. Comparisons of evolutionary models provided low ω global values and no evidence that HoxA13 experienced relaxed selection in limbless clades. Branch-site models failed to detect evidence for positive selection acting on any site along branches of Amphisbaena and Gymnophiona, while three sites were identified in Serpentes. Examination of alignments did not reveal consistent sequence differences between limbed and limbless species. We conclude that HoxA13 has no modules exclusive to limb development, which may be explained by its involvement in multiple developmental processes

Proteinortho: Detection of (Co-)orthologs in large-scale analysis

<p>Abstract</p> <p>Background</p> <p>Orthology analysis is an important part of data analysis in many areas of bioinformatics such as comparative genomics and molecular phylogenetics. The ever-increasing flood of sequence data, and hence the rapidly increasing number of genomes that can be compared simultaneously, calls for efficient software tools as brute-force approaches with quadratic memory requirements become infeasible in practise. The rapid pace at which new data become available, furthermore, makes it desirable to compute genome-wide orthology relations for a given dataset rather than relying on relations listed in databases.</p> <p>Results</p> <p>The program <monospace>Proteinortho</monospace> described here is a stand-alone tool that is geared towards large datasets and makes use of distributed computing techniques when run on multi-core hardware. It implements an extended version of the reciprocal best alignment heuristic. We apply <monospace>Proteinortho</monospace> to compute orthologous proteins in the complete set of all 717 eubacterial genomes available at NCBI at the beginning of 2009. We identified thirty proteins present in 99% of all bacterial proteomes.</p> <p>Conclusions</p> <p><monospace>Proteinortho</monospace> significantly reduces the required amount of memory for orthology analysis compared to existing tools, allowing such computations to be performed on off-the-shelf hardware.</p

“Genes”

In order to describe a cell at molecular level, a notion of a “gene” is neither necessary nor helpful. It is sufficient to consider the molecules (i.e., chromosomes, transcripts, proteins) and their interactions to describe cellular processes. The downside of the resulting high resolution is that it becomes very tedious to address features on the organismal and phenotypic levels with a language based on molecular terms. Looking for the missing link between biological disciplines dealing with different levels of biological organization, we suggest to return to the original intent behind the term “gene”. To this end, we propose to investigate whether a useful notion of “gene” can be constructed based on an underlying notion of function, and whether this can serve as the necessary link and embed the various distinct gene concepts of biological (sub)disciplines in a coherent theoretical framework. In reply to the Genon Theory recently put forward by Klaus Scherrer and Jürgen Jost in this journal, we shall discuss a general approach to assess a gene definition that should then be tested for its expressiveness and potential cross-disciplinary relevance

Tumour suppressor microRNA-584 directly targets oncogene Rock-1 and decreases invasion ability in human clear cell renal cell carcinoma

BackgroundThe purpose of this study was to identify new tumour suppressor microRNAs (miRs) in clear cell renal cell carcinoma (ccRCC), carry out functional analysis of their suppressive role and identify their specific target genes.MethodsTo explore suppressor miRs in RCC, miR microarray and real-time PCR were performed using HK-2 and A-498 cells. Cell viability, invasion and wound healing assays were carried out for functional analysis after miR transfection. To determine target genes of miR, we used messenger RNA (mRNA) microarray and target scan algorithms to identify target oncogenes. A 3'UTR luciferase assay was also performed. Protein expression of target genes in ccRCC tissues was confirmed by immunohistochemistry and was compared with miR-584 expression in ccRCC tissues.ResultsExpression of miR-584 in RCC (A-498 and 769-P) cells was downregulated compared with HK-2 cells. Transfection of miR-584 dramatically decreased cell motility. The ROCK-1 mRNA was inhibited by miR-584 and predicted to be target gene. The miR-584 decreased 3'UTR luciferase activity of ROCK-1 and ROCK-1 protein expression. Low expression of miR-584 in ccRCC tissues was correlated with high expression of ROCK-1 protein. The knockdown of ROCK-1 by siRNA inhibited cell motility.ConclusionmiR-584 is a new tumour suppressor miR in ccRCC and inhibits cell motility through downregulation of ROCK-1

Interleukin-6, tumour necrosis factor α and interleukin-1β in patients with renal cell carcinoma

As regulators of malignant cell behaviour and communication with stroma, cytokines have proved useful in understanding cancer biology and developing novel therapies. In renal cell carcinoma, patients with inflammatory reactions are known to have poor prognosis. In order to elucidate the relation between renal cell carcinoma and the host, serum levels of inflammatory cytokines, interleukin-6, tumour necrosis factor α, interleukin-1β, were measured. One hundred and twenty-two patients with renal cell carcinoma and 21 healthy control subjects were studied, and serum cytokine levels were measured using a highly sensitive ELISA kit. As a result, in the control group, interleukin-6, tumour necrosis factor α and interleukin-1β levels were 1.79±2.03, 2.74±0.94 and 0.16±0.17 pg ml−1, respectively. In the renal cell carcinoma patients, they were 8.91±13.12, 8.44±4.15 and 0.53±0.57 pg ml−1, respectively, and significantly higher. In the comparison of stage, interleukin-6 level was significantly higher in the stage IV group compared to the other stage groups including the control group, while tumour necrosis factor α level was significantly higher in each stage group compared to the control group. As for grade, interleukin-6 level was significantly higher in the grade 3 group compared to the control, grade 1 and grade 2 groups, while tumour necrosis factor α level was significantly higher in each grade group compared to the control group. All cytokines had a positive correlation with tumour size. In regard to the correlation with CRP, all cytokines had a positive correlation with CRP, while interleukin-6 had a particularly strong correlation. In conclusion, interleukin-6 may be one of the factors for the poor prognosis of patients with renal cell carcinoma. In addition, tumour necrosis factor α may be useful in the early diagnosis of renal cell carcinoma and post-operative follow-up

Interleukin-2/interferon-α2a/13-retinoic acid-based chemoimmunotherapy in advanced renal cell carcinoma: results of a prospectively randomised trial of the German Cooperative Renal Carcinoma Chemoimmunotherapy Group (DGCIN)

We performed a prospectively randomised clinical trial to compare the efficacy of four subcutaneous interleukin-2-(sc-IL-2) and sc interferon-α2a (sc-IFN-α2a)-based outpatient regimens in 379 patients with progressive metastatic renal cell carcinoma. Patients with lung metastases, an erythrocyte sedimentation rate ⩽70 mm h−1 and neutrophil counts ⩽6000 μl−1 (group I) were randomised to arm A: sc-IL-2, sc-IFN-α2a, peroral 13-cis-retinoic acid (po-13cRA) (n=78), or arm B: arm A plus inhaled-IL-2 (n=65). All others (group II) were randomised to arm C: arm A plus intravenous 5-fluorouracil (iv-5-FU) (n=116), or arm D: arm A plus po-Capecitabine (n=120). Median overall survival (OS) was 22 months (arm A; 3-year OS: 29.7%) and 18 months (arm B; 3-year OS: 29.2%) in group I, and 18 months (arm C; 3-year OS: 25.7%) and 16 months (arm D; 3-year OS: 32.6%) in group II. There were no statistically significant differences in OS, progression-free survival, and objective response between arms A and B, and between arms C and D, respectively. Given the known therapeutic efficacy of sc-IL-2/sc-INF-α2a/po-13cRA-based outpatient chemoimmunotherapies, our results did not establish survival advantages in favour of po-Capecitabine vs iv-5-FU, and in favour of short-term inhaled-IL-2 in patients with advanced renal cell carcinoma receiving systemic cytokines