Aplicação da análise de agrupamento de dados de expressão gênica temporal a dados em painel

O objetivo deste trabalho foi determinar a melhor alternativa, entre os métodos de agrupamento hierárquico (Ward) e de otimização (Tocher), para a formação de grupos homogêneos de séries de expressão gênica, e realizar previsões quanto à expressão gênica dessas séries, a partir de pequeno número de observações temporais. Os dados utilizados referem-se à expressão de genes que atuam sobre o ciclo celular de Saccharomyces cerevisiae e corresponderam a 114 séries de expressão gênica, cada uma com dez valores de "fold-change" (medida da expressão gênica) ao longo do tempo (0, 15, 30, 45, 60, 75, 90, 105, 120 e 135 min). As estimativas dos parâmetros dos modelos autorregressivos AR(p) foram previamente ajustadas a séries individuais (de cada gene) de dados "microarray time series" e utilizadas, como variáveis, no processo de agrupamento. As previsões da expressão gênica foram feitas dentro de cada grupo formado, a partir dos ajustes no modelo AR(p) para dados em painel. O método de Ward foi o mais apropriado para a formação de grupos de genes com séries homogêneas. Uma vez obtidos esses grupos, é possível ajustar o modelo AR(2) para dados em painel e predizer a expressão gênica em um tempo futuro (135 min), a partir de um pequeno número de observações temporais (os outros nove valores de "fold-change").

BadiRate: estimating family turnover rates by likelihood-based methods

Motivation: The comparative analysis of gene gain and loss rates is critical for understanding the role of natural selection and adaptation in shaping gene family sizes. Studying complete genome data from closely related species allows accurate estimation of gene family turnover rates. Current methods and software tools, however, are not well designed for dealing with certain kinds of functional elements, such as microRNAs or transcription factor binding sites. Results: Here, we describe BadiRate, a new software tool to estimate family turnover rates, as well as the number of elements in internal phylogenetic nodes, by likelihood-based methods and parsimony. It implements two stochastic population models, which provide the appropriate statistical framework for testing hypothesis, such as lineage-specific gene family expansions or contractions. We have assessed the accuracy of BadiRate by computer simulations, and have also illustrated its functionality by analyzing a representative empirical dataset.

Contamination of genomic databases by HIV-1 and its possible consequences. A study in Bioinformatics.

Peer-reviewed

Desenho e validação de iniciadores microssatélites SSR para mamoneira

O objetivo deste trabalho foi desenhar, validar e otimizar pares de iniciadores microssatélites SSR para mamoneira. O desenho dos pares de iniciadores foi feito por meio do aplicativo Websat, a partir de sequências depositadas no GenBank do National Center for Biotechnology Information (GenBank/NCBI), e a sua qualidade foi aferida com uso do aplicativo web NetPrimer. Foram utilizadas diferentes concentrações de DNA, cloreto de magnésio, pares de iniciadores, dNTPs e temperatura de anelamento para otimização das condições de PCR. Um total de 30 pares de iniciadores SSR foi desenhado, sintetizado e otimizado. O gel de agarose foi utilizado para detecção dos produtos amplificados, e o gel desnaturante de poliacrilamida, na otimização das condições de PCR e na identificação de polimorfismo. Os pares de iniciadores apresentaram percentagem média de guanina/citosina (GC) igual a 47,29% e produtos amplificados com tamanhos entre 128 e 381 pb. Vinte e nove pares de iniciadores SSR (96,7%) foram validados, dos quais nove foram polimórficos (23,3%). As concentrações otimizadas para amplificação são: DNA, 25 ng; cloreto de magnésio, 1,2 mmol L-1; iniciadores Forward e Reverse, 0,4 mmol L-1; dNTPs, 0,1 mmol L-1; e temperatura de anelamento, 62 a 64ºC. As ferramentas de bioinformática Websat e Net Primer podem ser utilizadas para desenvolver iniciadores microssatélites de qualidade, para a mamoneira, a partir de sequências depositadas no GenBank/NCBI.

Transfer entropy reconstruction and labeling of neuronal connections from simulated calcium imaging

Neuronal dynamics are fundamentally constrained by the underlying structural network architecture, yet much of the details of this synaptic connectivity are still unknown even in neuronal cultures in vitro. Here we extend a previous approach based on information theory, the Generalized Transfer Entropy, to the reconstruction of connectivity of simulated neuronal networks of both excitatory and inhibitory neurons. We show that, due to the model-free nature of the developed measure, both kinds of connections can be reliably inferred if the average firing rate between synchronous burst events exceeds a small minimum frequency. Furthermore, we suggest, based on systematic simulations, that even lower spontaneous inter-burst rates could be raised to meet the requirements of our reconstruction algorithm by applying a weak spatially homogeneous stimulation to the entire network. By combining multiple recordings of the same in silico network before and after pharmacologically blocking inhibitory synaptic transmission, we show then how it becomes possible to infer with high confidence the excitatory or inhibitory nature of each individual neuron.

DnaSP Version 5. DNA Sequence Polymorphism

DnaSP, DNA Sequence Polymorphism, is a software package for the analysis of nucleotide polymorphism from aligned DNA sequence data. DnaSP can estimate several measures of DNA sequence variation within and between populations (in noncoding, synonymous or nonsynonymous sites, or in various sorts of codon positions), as well as linkage disequilibrium, recombination, gene flow and gene conversion parameters. DnaSP can also carry out several tests of neutrality: Hudson, Kreitman and Aguadé (1987), Tajima (1989), McDonald and Kreitman (1991), Fu and Li (1993), and Fu (1997) tests. Additionally, DnaSP can estimate the confidence intervals of some test-statistics by the coalescent. The results of the analyses are displayed on tabular and graphic form.

VariScan Software

VariScan is a software package for the analysis of DNA sequence polymorphisms at the whole genome scale. Among other features, the software:(1) can conduct many population genetic analyses; (2) incorporates a multiresolution wavelet transform-based method that allows capturing relevant information from DNA polymorphism data; and (3) it facilitates the visualization of the results in the most commonly used genome browsers.

miR-143 Interferes with ERK5 Signaling, and Abrogates Prostate Cancer Progression in Mice

Abstract Background: Micro RNAs are small, non-coding, single-stranded RNAs that negatively regulate gene expression at the post-transcriptional level. Since miR-143 was found to be down-regulated in prostate cancer cells, we wanted to analyze its expression in human prostate cancer, and test the ability of miR-43 to arrest prostate cancer cell growth in vitro and in vivo. Results: Expression of miR-143 was analyzed in human prostate cancers by quantitative PCR, and by in situ hybridization. miR-143 was introduced in cancer cells in vivo by electroporation. Bioinformatics analysis and luciferase-based assays were used to determine miR-143 targets. We show in this study that miR-143 levels are inversely correlated with advanced stages of prostate cancer. Rescue of miR-143 expression in cancer cells results in the arrest of cell proliferation and the abrogation of tumor growth in mice. Furthermore, we show that the effects of miR-143 are mediated, at least in part by the inhibition of extracellular signal-regulated kinase-5 (ERK5) activity. We show here that ERK5 is a miR-143 target in prostate cancer. Conclusions: miR-143 is as a new target for prostate cancer treatment.

Kernel-PCA data integration with enhanced interpretability

Background Nowadays, combining the different sources of information to improve the biological knowledge available is a challenge in bioinformatics. One of the most powerful methods for integrating heterogeneous data types are kernel-based methods. Kernel-based data integration approaches consist of two basic steps: firstly the right kernel is chosen for each data set; secondly the kernels from the different data sources are combined to give a complete representation of the available data for a given statistical task. Results We analyze the integration of data from several sources of information using kernel PCA, from the point of view of reducing dimensionality. Moreover, we improve the interpretability of kernel PCA by adding to the plot the representation of the input variables that belong to any dataset. In particular, for each input variable or linear combination of input variables, we can represent the direction of maximum growth locally, which allows us to identify those samples with higher/lower values of the variables analyzed. Conclusions The integration of different datasets and the simultaneous representation of samples and variables together give us a better understanding of biological knowledge.

Genetic Variation in the TP53 Pathway and Bladder Cancer Risk. A Comprehensive Analysis

Introduction: Germline variants in TP63 have been consistently associated with several tumors, including bladder cancer, indicating the importance of TP53 pathway in cancer genetic susceptibility. However, variants in other related genes, including TP53 rs1042522 (Arg72Pro), still present controversial results. We carried out an in depth assessment of associations between common germline variants in the TP53 pathway and bladder cancer risk. Material and Methods: We investigated 184 tagSNPs from 18 genes in 1,058 cases and 1,138 controls from the Spanish Bladder Cancer/EPICURO Study. Cases were newly-diagnosed bladder cancer patients during 1998–2001. Hospital controls were age-gender, and area matched to cases. SNPs were genotyped in blood DNA using Illumina Golden Gate and TaqMan assays. Cases were subphenotyped according to stage/grade and tumor p53 expression. We applied classical tests to assess individual SNP associations and the Least Absolute Shrinkage and Selection Operator (LASSO)-penalized logistic regression analysis to assess multiple SNPs simultaneously. Results: Based on classical analyses, SNPs in BAK1 (1), IGF1R (5), P53AIP1 (1), PMAIP1 (2), SERINPB5 (3), TP63 (3), and TP73 (1) showed significant associations at p-value#0.05. However, no evidence of association, either with overall risk or with specific disease subtypes, was observed after correction for multiple testing (p-value$0.8). LASSO selected the SNP rs6567355 in SERPINB5 with 83% of reproducibility. This SNP provided an OR = 1.21, 95%CI 1.05–1.38, p-value = 0.006, and a corrected p-value = 0.5 when controlling for over-estimation. Discussion: We found no strong evidence that common variants in the TP53 pathway are associated with bladder cancer susceptibility. Our study suggests that it is unlikely that TP53 Arg72Pro is implicated in the UCB in white Europeans. SERPINB5 and TP63 variation deserve further exploration in extended studies.

Methylglyoxal Produced by Amyloid- Peptide-Induced Nitrotyrosination of Triosephosphate Isomerase Triggers Neuronal Death in Alzheimer’s Disease

Amyloid-β peptide (Aβ) aggregates induce nitro-oxidative stress, contributing to the characteristic neurodegeneration found in Alzheimer's disease (AD). One of the most strongly nitrotyrosinated proteins in AD is the triosephosphate isomerase (TPI) enzyme which regulates glycolytic flow, and its efficiency decreased when it is nitrotyrosinated. The main aims of this study were to analyze the impact of TPI nitrotyrosination on cell viability and to identify the mechanism behind this effect. In human neuroblastoma cells (SH-SY5Y), we evaluated the effects of Aβ42 oligomers on TPI nitrotyrosination. We found an increased production of methylglyoxal (MG), a toxic byproduct of the inefficient nitro-TPI function. The proapoptotic effects of Aβ42 oligomers, such as decreasing the protective Bcl2 and increasing the proapoptotic caspase-3 and Bax, were prevented with a MG chelator. Moreover, we used a double mutant TPI (Y165F and Y209F) to mimic nitrosative modifications due to Aβ action. Neuroblastoma cells transfected with the double mutant TPI consistently triggered MG production and a decrease in cell viability due to apoptotic mechanisms. Our data show for the first time that MG is playing a key role in the neuronal death induced by Aβ oligomers. This occurs because of TPI nitrotyrosination, which affects both tyrosines associated with the catalytic center.

Risk Prediction Scores for Recurrence and Progression of Non-Muscle Invasive Bladder Cancer: An International Validation in Primary Tumours

Abstract Objective: We aimed to determine the validity of two risk scores for patients with non-muscle invasive bladder cancer in different European settings, in patients with primary tumours. Methods: We included 1,892 patients with primary stage Ta or T1 non-muscle invasive bladder cancer who underwent a transurethral resection in Spain (n = 973), the Netherlands (n = 639), or Denmark (n = 280). We evaluated recurrence-free survival and progression-free survival according to the European Organisation for Research and Treatment of Cancer (EORTC) and the Spanish Urological Club for Oncological Treatment (CUETO) risk scores for each patient and used the concordance index (c-index) to indicate discriminative ability. Results: The 3 cohorts were comparable according to age and sex, but patients from Denmark had a larger proportion of patients with the high stage and grade at diagnosis (p,0.01). At least one recurrence occurred in 839 (44%) patients and 258 (14%) patients had a progression during a median follow-up of 74 months. Patients from Denmark had the highest 10- year recurrence and progression rates (75% and 24%, respectively), whereas patients from Spain had the lowest rates (34% and 10%, respectively). The EORTC and CUETO risk scores both predicted progression better than recurrence with c-indices ranging from 0.72 to 0.82 while for recurrence, those ranged from 0.55 to 0.61. Conclusion: The EORTC and CUETO risk scores can reasonably predict progression, while prediction of recurrence is more difficult. New prognostic markers are needed to better predict recurrence of tumours in primary non-muscle invasive bladder cancer patients.

BioSWR - Semantic Web services Registry for Bioinformatics

Article About the Authors Metrics Comments Related Content Abstract Introduction Functionality Implementation Discussion Acknowledgments Author Contributions References Reader Comments (0) Figures Abstract Despite of the variety of available Web services registries specially aimed at Life Sciences, their scope is usually restricted to a limited set of well-defined types of services. While dedicated registries are generally tied to a particular format, general-purpose ones are more adherent to standards and usually rely on Web Service Definition Language (WSDL). Although WSDL is quite flexible to support common Web services types, its lack of semantic expressiveness led to various initiatives to describe Web services via ontology languages. Nevertheless, WSDL 2.0 descriptions gained a standard representation based on Web Ontology Language (OWL). BioSWR is a novel Web services registry that provides standard Resource Description Framework (RDF) based Web services descriptions along with the traditional WSDL based ones. The registry provides Web-based interface for Web services registration, querying and annotation, and is also accessible programmatically via Representational State Transfer (REST) API or using a SPARQL Protocol and RDF Query Language. BioSWR server is located at http://inb.bsc.es/BioSWR/and its code is available at https://sourceforge.net/projects/bioswr/​under the LGPL license.

p21 as a Transcriptional Co-Repressor of S-Phase and Mitotic Control Genes

It has been previously described that p21 functions not only as a CDK inhibitor but also as a transcriptional co-repressor in some systems. To investigate the roles of p21 in transcriptional control, we studied the gene expression changes in two human cell systems. Using a human leukemia cell line (K562) with inducible p21 expression and human primary keratinocytes with adenoviral-mediated p21 expression, we carried out microarray-based gene expression profiling. We found that p21 rapidly and strongly repressed the mRNA levels of a number of genes involved in cell cycle and mitosis. One of the most strongly down-regulated genes was CCNE2 (cyclin E2 gene). Mutational analysis in K562 cells showed that the N-terminal region of p21 is required for repression of gene expression of CCNE2 and other genes. Chromatin immunoprecipitation assays indicated that p21 was bound to human CCNE2 and other p21-repressed genes gene in the vicinity of the transcription start site. Moreover, p21 repressed human CCNE2 promoter-luciferase constructs in K562 cells. Bioinformatic analysis revealed that the CDE motif is present in most of the promoters of the p21-regulated genes. Altogether, the results suggest that p21 exerts a repressive effect on a relevant number of genes controlling S phase and mitosis. Thus, p21 activity as inhibitor of cell cycle progression would be mediated not only by the inhibition of CDKs but also by the transcriptional down-regulation of key genes.

An integrative computational analysis provides evidence for FBN1-associated network deregulation in trisomy 21

Although approximately 50% of Down Syndrome (DS) patients have heart abnormalities, they exhibit an overprotection against cardiac abnormalities related with the connective tissue, for example a lower risk of coronary artery disease. A recent study reported a case of a person affected by DS who carried mutations in FBN1, the gene causative for a connective tissue disorder called Marfan Syndrome (MFS). The fact that the person did not have any cardiac alterations suggested compensation effects due to DS. This observation is supported by a previous DS meta-analysis at the molecular level where we have found an overall upregulation of FBN1 (which is usually downregulated in MFS). Additionally, that result was cross-validated with independent expression data from DS heart tissue. The aim of this work is to elucidate the role of FBN1 in DS and to establish a molecular link to MFS and MFS-related syndromes using a computational approach. To reach that, we conducted different analytical approaches over two DS studies (our previous meta-analysis and independent expression data from DS heart tissue) and revealed expression alterations in the FBN1 interaction network, in FBN1 co-expressed genes and FBN1-related pathways. After merging the significant results from different datasets with a Bayesian approach, we prioritized 85 genes that were able to distinguish control from DS cases. We further found evidence for several of these genes (47%), such as FBN1, DCN, and COL1A2, being dysregulated in MFS and MFS-related diseases. Consequently, we further encourage the scientific community to take into account FBN1 and its related network for the study of DS cardiovascular characteristics.