GEDI: a user-friendly toolbox for analysis of large-scale gene expression data

Abstract Background Several mathematical and statistical methods have been proposed in the last few years to analyze microarray data. Most of those methods involve complicated formulas, and software implementations that require advanced computer programming skills. Researchers from other areas may experience difficulties when they attempting to use those methods in their research. Here we present an user-friendly toolbox which allows large-scale gene expression analysis to be carried out by biomedical researchers with limited programming skills. Results Here, we introduce an user-friendly toolbox called GEDI (Gene Expression Data Interpreter), an extensible, open-source, and freely-available tool that we believe will be useful to a wide range of laboratories, and to researchers with no background in Mathematics and Computer Science, allowing them to analyze their own data by applying both classical and advanced approaches developed and recently published by Fujita et al. Conclusion GEDI is an integrated user-friendly viewer that combines the state of the art SVR, DVAR and SVAR algorithms, previously developed by us. It facilitates the application of SVR, DVAR and SVAR, further than the mathematical formulas present in the corresponding publications, and allows one to better understand the results by means of available visualizations. Both running the statistical methods and visualizing the results are carried out within the graphical user interface, rendering these algorithms accessible to the broad community of researchers in Molecular Biology.

Genomic analysis of ERVWE2 locus in patients with multiple sclerosis: absence of genetic association but potential role of human endogenous retrovirus type W elements in molecular mimicry with myelin antigen

Human endogenous retroviruses (HERVs) arise from ancient infections of the host germline cells by exogenous retroviruses, constituting 8% of the human genome. Elevated level of envelope transcripts from HERVs-W has been detected in CSF, plasma and brain tissues from patients with Multiple Sclerosis (MS), most of them from Xq22.3, 15q21.3, and 6q21 chromosomes. However, since the locus Xq22.3 (ERVWE2) lack the 5' LTR promoter and the putative protein should be truncated due to a stop codon, we investigated the ERVWE2 genomic loci from 84 individuals, including MS patients with active HERV-W expression detected in PBMC. In addition, an automated search for promoter sequences in 20 kb nearby region of ERVWE2 reference sequence was performed. Several putative binding sites for cellular cofactors and enhancers were found, suggesting that transcription may occur via alternative promoters. However, ERVWE2 DNA sequencing of MS and healthy individuals revealed that all of them harbor a stop codon at site 39, undermining the expression of a full-length protein. Finally, since plaque formation in central nervous system (CNS) of MS patients is attributed to immunological mechanisms triggered by autoimmune attack against myelin, we also investigated the level of similarity between envelope protein and myelin oligodendrocyte glycoprotein (MOG). Comparison of the MOG to the envelope identified five retroviral regions similar to the Ig-like domain of MOG. Interestingly, one of them includes T and B cell epitopes, capable to induce T effector functions and circulating Abs in rats. In sum, although no DNA substitutions that would link ERVWE2 to the MS pathogeny was found, the similarity between the envelope protein to MOG extends the idea that ERVEW2 may be involved on the immunopathogenesis of MS, maybe facilitating the MOG recognizing by the immune system. Although awaiting experimental evidences, the data presented here may expand the scope of the endogenous retroviruses involvement on MS pathogenesis

Proteomic Approaches Identify Members of Cofilin Pathway Involved in Oral Tumorigenesis

The prediction of tumor behavior for patients with oral carcinomas remains a challenge for clinicians. The presence of lymph node metastasis is the most important prognostic factor but it is limited in predicting local relapse or survival. This highlights the need for identifying biomarkers that may effectively contribute to prediction of recurrence and tumor spread. In this study, we used one-and two-dimensional gel electrophoresis, mass spectrometry and immunodetection methods to analyze protein expression in oral squamous cell carcinomas. Using a refinement for classifying oral carcinomas in regard to prognosis, we analyzed small but lymph node metastasis-positive versus large, lymph node metastasis-negative tumors in order to contribute to the molecular characterization of subgroups with risk of dissemination. Specific protein patterns favoring metastasis were observed in the "more-aggressive'' group defined by the present study. This group displayed upregulation of proteins involved in migration, adhesion, angiogenesis, cell cycle regulation, anti-apoptosis and epithelial to mesenchymal transition, whereas the "less-aggressive'' group was engaged in keratinocyte differentiation, epidermis development, inflammation and immune response. Besides the identification of several proteins not yet described as deregulated in oral carcinomas, the present study demonstrated for the first time the role of cofilin-1 in modulating cell invasion in oral carcinomas.

Demographic processes in the montane Atlantic rainforest: Molecular and cytogenetic evidence from the endemic frog Proceratophrys boiei

Historical climatic refugia predict genetic diversity in lowland endemics of the Brazilian Atlantic rainforest. Yet, available data reveal distinct biological responses to the Last Glacial Maximum (LGM) conditions across species of different altitudinal ranges. We show that species occupying Brazil's montane forests were significantly less affected by LGM conditions relative to lowland specialists, but that pre-Pleistocene tectonics greatly influenced their geographic variation. Our conclusions are based on palaeoclimatic distribution models, molecular sequences of the cytochrome b, 16S, and RAG-1 genes, and karyotype data for the endemic frog Proceratophrys boiei. DNA and chromosomal data identify in P. boiei at least two broadly divergent phylogroups, which have not been distinguished morphologically. Cytogenetic results also indicate an area of hybridization in southern Sao Paulo. The location of the phylogeographic break broadly matches the location of a NW-SE fault, which underwent reactivation in the Neogene and led to remarkable landscape changes in southeastern Brazil. Our results point to different mechanisms underpinning diversity patterns in lowland versus montane tropical taxa, and help us to understand the processes responsible for the large number of narrow endemics currently observed in montane areas of the southern Atlantic forest hotspot. (C) 2011 Elsevier Inc. All rights reserved.

Structural role of the active-site metal in the conformation of Trypanosoma brucei phosphoglycerate mutase

Phosphoglycerate mutases (PGAMs) participate in both the glycolytic and the gluconeogenic pathways in reversible isomerization of 3-phosphoglycerate and 2-phosphoglycerate. PGAMs are members of two distinct protein families: enzymes that are dependent on or independent of the 2,3-bisphosphoglycerate cofactor. We determined the X-ray structure of the monomeric Trypanosoma brucei independent PGAM (TbiPGAM) in its apoenzyme form, and confirmed this observation by small angle X-ray scattering data. Comparing the TbiPGAM structure with the Leishmania mexicana independent PGAM structure, previously reported with a phosphoglycerate molecule bound to the active site, revealed the domain movement resulting from active site occupation. The structure reported here shows the interaction between Asp319 and the metal bound to the active site, and its contribution to the domain movement. Substitution of the metal-binding residue Asp319 by Ala resulted in complete loss of independent PGAM activity, and showed for the first time its involvement in the enzymes function. As TbiPGAM is an attractive molecular target for drug development, the apoenzyme conformation described here provides opportunities for its use in structure-based drug design approaches.

A molecular linkage map for Drosophila mediopunctata confirms synteny with Drosophila melanogaster and suggests a region that controls the variation in the number of abdominal spots

The classic approach to gene discovery relies on the construction of linkage maps. We report the first molecular-based linkage map for Drosophila mediopunctata, a neotropical species of the tripunctata group. Eight hundred F2 individuals were genotyped at 49 microsatellite loci, resulting in a map that is approximate to 450 centimorgans long. Five linkage groups were detected, and the species' chromosomes were identified through cross-references to BLASTn searches and Muller elements. Strong synteny was observed when compared with the Drosophila melanogaster chromosome arms, but little conservation in the gene order was seen. The incorporation of morphological data corresponding to the number of central abdominal spots on the map was consistent with the expected location of a genomic region responsible for the phenotype on the second chromosome.

Molecular characterization of a mariner-like element in the Atta sexdens rubropilosa genome

Mobile elements are widely present in eukaryotic genomes. They are repeated DNA segments that are able to move from one locus to another within the genome. They are divided into two main categories, depending on their mechanism of transposition, involving RNA (class I) or DNA (class II) molecules. The mariner-like elements are class II transposons. They encode their own transposase, which is necessary and sufficient for transposition in the absence of host factors. They are flanked by a short inverted terminal repeat and a TA dinucleotide target site, which is duplicated upon insertion. The transposase consists of two domains, an N-terminal inverted terminal repeat binding domain and a C-terminal catalytic domain. We identified a transposable element with molecular characteristics of a mariner-like element in Atta sexdens rubropilosa genome. Identification started from a PCR with degenerate primers and queen genomic DNA templates, with which it was possible to amplify a fragment with mariner transposable-element homology. Phylogenetic analysis demonstrated that this element belongs to the mauritiana subfamily of mariner-like elements and it was named Asmar1. We found that Asmar1 is homologous to a transposon described from another ant, Messor bouvieri. The predicted transposase sequence demonstrated that Asmar1 has a truncated transposase ORF. This study is part of a molecular characterization of mobile elements in the Atta spp genome. Our finding of mariner-like elements in all castes of this ant could be useful to help understand the dynamics of mariner-like element distribution in the Hymenoptera.

Involvement of TSSA (trypomastigote small surface antigen) in Trypanosoma cruzi invasion of mammalian cells

TSSA (trypomastigote small surface antigen) is a polymorphic mucin-like molecule displayed on the surface of Trypanosoma cruzi trypomastigote forms. To evaluate its functional properties, we undertook comparative biochemical and genetic approaches on isoforms present in parasite stocks from extant evolutionary lineages (CL Brener and Sylvio X-10). We show that CL Brener TSSA, but not the Sylvio X-10 counterpart, exhibits dose-dependent and saturable binding towards non-macrophagic cell lines. This binding triggers Ca2+-based signalling responses in the target cell while providing an anchor for the invading parasite. Accordingly, exogenous addition of either TSSA-derived peptides or specific antibodies significantly inhibits invasion of CL Brener, but not Sylvio X-10, trypomastigotes. Non-infective epimastigote forms, which do not express detectable levels of TSSA, were stably transfected with TSSA cDNA from either parasite stock. Although both transfectants produced a surface-associated mucin-like TSSA product, epimastigotes expressing CL Brener TSSA showed a similar to 2-fold increase in their attachment to mammalian cells. Overall, these findings indicate that CL Brener TSSA functions as a parasite adhesin, engaging surface receptor(s) and inducing signalling pathways on the host cell as a prerequisite for parasite internalization. More importantly, the contrasting functional features of TSSA isoforms provide one appealing mechanism underlying the differential infectivity of T. cruzi stocks.

Prion potency in stem cells biology

Prion protein (PrP) can be considered a pivotal molecule because it interacts with several partners to perform a diverse range of critical biological functions that might differ in embryonic and adult cells. In recent years, there have been major advances in elucidating the putative role of PrP in the basic biology of stem cells in many different systems. Here, we review the evidence indicating that PrP is a key molecule involved in driving different aspects of the potency of embryonic and tissue-specific stem cells in self-perpetuation and differentiation in many cell types. It has been shown that PrP is involved in stem cell self-renewal, controlling pluripotency gene expression, proliferation and neural and cardiomyocyte differentiation. PrP also has essential roles in distinct processes that regulate tissue-specific stem cell biology in nervous and hematopoietic systems and during muscle regeneration. Results from our own investigations have shown that PrP is able to modulate self-renewal and proliferation in neural stem cells, processes that are enhanced by PrP interactions with stress inducible protein 1 (STI1). Thus, the available data reveal the influence of PrP in acting upon the maintenance of pluripotent status or the differentiation of stem cells from the early embryogenesis through adulthood.

Effect of the thermostat in the molecular dynamics simulation on the folding of the model protein chignolin

Molecular dynamics simulations of the model protein chignolin with explicit solvent were carried out, in order to analyze the influence of the Berendsen thermostat on the evolution and folding of the peptide. The dependence of the peptide behavior on temperature was tested with the commonly employed thermostat scheme consisting of one thermostat for the protein and another for the solvent. The thermostat coupling time of the protein was increased to infinity, when the protein is not in direct contact with the thermal bath, a situation known as minimally invasive thermostat. In agreement with other works, it was observed that only in the last situation the instantaneous temperature of the model protein obeys a canonical distribution. As for the folding studies, it was shown that, in the applications of the commonly utilized thermostat schemes, the systems are trapped in local minima regions from which it has difficulty escaping. With the minimally invasive thermostat the time that the protein needs to fold was reduced by two to three times. These results show that the obstacles to the evolution of the extended peptide to the folded structure can be overcome when the temperature of the peptide is not directly controlled.

Molecular Dynamics Investigations of PRODAN in a DLPC Bilayer

Molecular dynamics computer simulations have been performed to identify preferred positions of the fluorescent probe PRODAN in a fully hydrated DLPC bilayer in the fluid phase. In addition to the intramolecular charge-transfer first vertical excited state, we considered different charge distributions for the electronic ground state of the PRODAN molecule by distinct atomic charge models corresponding to the probe molecule in vacuum as well as polarized in a weak and a strong dielectric solvent (cyclohexane and water). Independent on the charge distribution model of PRODAN, we observed a preferential orientation of this molecule in the bilayer with the dimethylamino group pointing toward the membrane's center and the carbonyl oxygen toward the membrane's interface. However, changing the charge distribution model of PRODAN, independent of its initial position in the equilibrated DLPC membrane, we observed different preferential positions. For the ground state representation without polarization and the in-cyclohexane polarization, the probe maintains its position close to the membrane's center. Considering the in-water polarization model, the probe approaches more of the polar headgroup region of the bilayer, with a strong structural correlation with the choline group, exposing its oxygen atom to water molecules. PRODAN's representation of the first vertical excited state with the in-water polarization also approaches the polar region of the membrane with the oxygen atom exposed to the bilayer's hydration shell. However, this model presents a stronger structural correlation with the phosphate groups than the ground state. Therefore, we conclude that the orientation of the PRODAN molecule inside the DLPC membrane is well-defined, but its position is very sensitive to the effect of the medium polarization included here by different models for the atomic charge distribution of the probe.