A reliable measure of similarity based on dependency for short time series: an application to gene expression networks

Background: Microarray techniques have become an important tool to the investigation of genetic relationships and the assignment of different phenotypes. Since microarrays are still very expensive, most of the experiments are performed with small samples. This paper introduces a method to quantify dependency between data series composed of few sample points. The method is used to construct gene co-expression subnetworks of highly significant edges. Results: The results shown here are for an adapted subset of a Saccharomyces cerevisiae gene expression data set with low temporal resolution and poor statistics. The method reveals common transcription factors with a high confidence level and allows the construction of subnetworks with high biological relevance that reveals characteristic features of the processes driving the organism adaptations to specific environmental conditions. Conclusion: Our method allows a reliable and sophisticated analysis of microarray data even under severe constraints. The utilization of systems biology improves the biologists ability to elucidate the mechanisms underlying celular processes and to formulate new hypotheses.

Genomics and proteomics approaches to the study of cancer-stroma interactions

Background: The development and progression of cancer depend on its genetic characteristics as well as on the interactions with its microenvironment. Understanding these interactions may contribute to diagnostic and prognostic evaluations and to the development of new cancer therapies. Aiming to investigate potential mechanisms by which the tumor microenvironment might contribute to a cancer phenotype, we evaluated soluble paracrine factors produced by stromal and neoplastic cells which may influence proliferation and gene and protein expression. Methods: The study was carried out on the epithelial cancer cell line (Hep-2) and fibroblasts isolated from a primary oral cancer. We combined a conditioned-medium technique with subtraction hybridization approach, quantitative PCR and proteomics, in order to evaluate gene and protein expression influenced by soluble paracrine factors produced by stromal and neoplastic cells. Results: We observed that conditioned medium from fibroblast cultures (FCM) inhibited proliferation and induced apoptosis in Hep-2 cells. In neoplastic cells, 41 genes and 5 proteins exhibited changes in expression levels in response to FCM and, in fibroblasts, 17 genes and 2 proteins showed down-regulation in response to conditioned medium from Hep-2 cells (HCM). Nine genes were selected and the expression results of 6 down-regulated genes (ARID4A, CALR, GNB2L1, RNF10, SQSTM1, USP9X) were validated by real time PCR. Conclusions: A significant and common denominator in the results was the potential induction of signaling changes associated with immune or inflammatory response in the absence of a specific protein.

Mucin pattern reflects the origin of the adenocarcinoma in Barrett's esophagus: a retrospective clinical and laboratorial study

Background: Mucin immunoexpression in adenocarcinoma arising in Barrett's esophagus (BE) may indicate the carcinogenesis pathway. The aim of this study was to evaluate resected specimens of adenocarcinoma in BE for the pattern of mucins and to correlate to the histologic classification. Methods: Specimens were retrospectively collected from thirteen patients who underwent esophageal resection due to adenocarcinoma in BE. Sections were scored for the grade of intestinal metaplasia. The tissues were examined by immunohistochemistry for MUC2 and MUC5AC antibodies. Results: Eleven patients were men. The mean age was 61 years old (varied from 40 to 75 years old). The tumor size had a mean of 4.7 +/- 2.3 cm, and the extension of BE had a mean of 7.7 +/- 1.5 cm. Specialized epithelium with intestinal metaplasia was present in all adjacent mucosas. Immunohistochemistry for MUC2 showed immunoreactivity in goblet cells, while MUC5AC was extensively expressed in the columnar gastric cells, localizing to the surface epithelium and extending to a variable degree into the glandular structures in BE. Tumors were classified according to the mucins in gastric type in 7/13 (MUC5AC positive) and intestinal type in 4/13 (MUC2 positive). Two tumors did not express MUC2 or MUC5AC proteins. The pattern of mucin predominantly expressed in the adjacent epithelium was associated to the mucin expression profile in the tumors, p = 0.047. Conclusion: Barrett's esophagus adenocarcinoma shows either gastric or intestinal type pattern of mucin expression. The two types of tumors developed in Barrett's esophagus may reflect the original cell type involved in the malignant transformation.

Cytogenetic findings in an epithelioid sarcoma with angiomatoid features. A case report

Epithelioid sarcoma is a rare, aggressive soft tissue tumor of unknown histogenesis showing predominantly epithelioid cytomorphology. We conducted a conventional and molecular cytogenetic study of a 27-year-old male with epithelioid sarcoma with angiomatoid features. Cytogenetic analysis of epithelioid sarcoma metaphase spreads by GTG-banding revealed a diploid chromosome complement with structural and numerical aberrations. Comparative genomic hybridization analysis demonstrated the amplification of 3p24 similar to pter, 4p15.2-p16 and 18q23, while chromosome losses involved 3p13-p14, 3q24-q26.1, 9q21, and 11q21. Fluorescence in situ hybridization assessment showed normal hybridization patterns for the C-MYC and CCND1 loci; CCND1 RNA overexpression was detected by real-time polymerase chain reaction analysis. Genetic evaluation of this rare condition may be useful in determining if epithelioid sarcoma is associated with a distinct genetic background.

Expression and cellular localization of microRNA-29b and RAX, an activator of the RNA-dependent protein kinase (PKR), in the retina of streptozotocin-induced diabetic rats

Purpose: The apoptosis of retinal neurons plays a critical role in the pathogenesis of diabetic retinopathy (DR), but the molecular mechanisms underlying this phenomenon remain unclear. The purpose of this study was to investigate the cellular localization and the expression of microRNA-29b (miR-29b) and its potential target PKR associated protein X (RAX), an activator of the pro-apoptotic RNA-dependent protein kinase (PKR) signaling pathway, in the retina of normal and diabetic rats. Methods: Retinas were obtained from normal and diabetic rats within 35 days after streptozotocin (STZ) injection. In silico analysis indicated that RAX is a potential target of miR-29b. The cellular localization of miR-29b and RAX was assessed by in situ hybridization and immunofluorescence, respectively. The expression levels of miR-29b and RAX mRNA were evaluated by quantitative reverse transcription PCR (qRT-PCR), and the expression of RAX protein was evaluated by western blot. A luciferase reporter assay and inhibition of endogenous RAX were performed to confirm whether RAX is a direct target of miR-29b as predicted by the in silico analysis. Results: We found that miR-29b and RAX are localized in the retinal ganglion cells (RGCs) and the cells of the inner nuclear layer (INL) of the retinas from normal and diabetic rats. Thus, the expression of miR-29b and RAX, as assessed in the retina by quantitative RT-PCR, reflects their expression in the RGCs and the cells of the INL. We also revealed that RAX protein is upregulated (more than twofold) at 3, 6, 16, and 22 days and downregulated (70%) at 35 days, whereas miR-29b is upregulated (more than threefold) at 28 and 35 days after STZ injection. We did not confirm the computational prediction that RAX is a direct target of miR-29b. Conclusions: Our results suggest that RAX expression may be indirectly regulated by miR-29b, and the upregulation of this miRNA at the early stage of STZ-induced diabetes may have a protective effect against the apoptosis of RGCs and cells of the INL by the pro-apoptotic RNA-dependent protein kinase (PKR) signaling pathway.

Novel Primate-Specific Genes, RMEL 1, 2 and 3, with Highly Restricted Expression in Melanoma, Assessed by New Data Mining Tool

Melanoma is a highly aggressive and therapy resistant tumor for which the identification of specific markers and therapeutic targets is highly desirable. We describe here the development and use of a bioinformatic pipeline tool, made publicly available under the name of EST2TSE, for the in silico detection of candidate genes with tissue-specific expression. Using this tool we mined the human EST (Expressed Sequence Tag) database for sequences derived exclusively from melanoma. We found 29 UniGene clusters of multiple ESTs with the potential to predict novel genes with melanoma-specific expression. Using a diverse panel of human tissues and cell lines, we validated the expression of a subset of three previously uncharacterized genes (clusters Hs.295012, Hs.518391, and Hs.559350) to be highly restricted to melanoma/melanocytes and named them RMEL1, 2 and 3, respectively. Expression analysis in nevi, primary melanomas, and metastatic melanomas revealed RMEL1 as a novel melanocytic lineage-specific gene up-regulated during melanoma development. RMEL2 expression was restricted to melanoma tissues and glioblastoma. RMEL3 showed strong up-regulation in nevi and was lost in metastatic tumors. Interestingly, we found correlations of RMEL2 and RMEL3 expression with improved patient outcome, suggesting tumor and/or metastasis suppressor functions for these genes. The three genes are composed of multiple exons and map to 2q12.2, 1q25.3, and 5q11.2, respectively. They are well conserved throughout primates, but not other genomes, and were predicted as having no coding potential, although primate-conserved and human-specific short ORFs could be found. Hairpin RNA secondary structures were also predicted. Concluding, this work offers new melanoma-specific genes for future validation as prognostic markers or as targets for the development of therapeutic strategies to treat melanoma.

Epigenetic Silencing of CRABP2 and MX1 in Head and Neck Tumors

Head and neck squamous cell carcinoma (HNSCC) is a heterogeneous disease affecting the epithelium of the oral cavity, pharynx and larynx. Conditions of most patients are diagnosed at late stages of the disease, and no sensitive and specific predictors of aggressive behavior have been identified yet. Therefore, early detection and prognostic biomarkers are highly desirable for a more rational management of the disease. Hypermethylation of CpG islands is one of the most important epigenetic mechanisms that leads to gene silencing in tumors and has been extensively used for the identification of biomarkers. In this study, we combined rapid subtractive hybridization and microarray analysis in a hierarchical manner to select genes that are putatively reactivated by the demethylating agent 5-aza-2'-deoxycytidine (5Aza-dC) in HNSCC cell lines (FaDu, UM-SCC-14A, UM-SCC-17A, UM-SCC-38A). This combined analysis identified 78 genes, 35 of which were reactivated in at least 2 cell lines and harbored a CpG island at their 5' region. Reactivation of 3 of these 35 genes (CRABP2, MX1, and SLC15A3) was confirmed by quantitative real-time polymerase chain reaction (PCR; fold change, >= 3). Bisulfite sequencing of their CpG islands revealed that they are indeed differentially methylated in the HNSCC cell lines. Using methylation-specific PCR, we detected a higher frequency of CRABP2 (58.1% for region 1) and MX1 (46.3%) hypermethylation in primary HNSCC when compared with lymphocytes from healthy individuals. Finally, absence of the CRABP2 protein was associated with decreased disease-free survival rates, supporting a potential use of CRABP2 expression as a prognostic biomarker for HNSCC patients.

Characterization of new IS elements and studies of their dispersion in two subspecies of Leifsonia xyli

Background: Leifsonia xyli is a xylem-inhabiting bacterial species comprised of two subspecies: L. xyli subsp. xyli (Lxx) and L. xyli subsp. cynodontis (Lxc). Lxx is the causal agent of ratoon stunting disease in sugarcane commercial fields and Lxc colonizes the xylem of several grasses causing either mild or no symptoms of disease. The completely sequenced genome of Lxx provided insights into its biology and pathogenicity. Since IS elements are largely reported as an important source of bacterial genome diversification and nothing is known about their role in chromosome architecture of L. xyli, a comparative analysis of Lxc and Lxx elements was performed. Results: Sample sequencing of Lxc genome and comparative analysis with Lxx complete DNA sequence revealed a variable number of IS transposable elements acting upon genomic diversity. A detailed characterization of Lxc IS elements and a comparative review with IS elements of Lxx are presented. Each genome showed a unique set of elements although related to same IS families when considering features such as similarity among transposases, inverted and direct repeats, and element size. Most of the Lxc and Lxx IS families assigned were reported to maintain transposition at low levels using translation regulatory mechanisms, consistent with our in silico analysis. Some of the IS elements were found associated with rearrangements and specific regions of each genome. Differences were also found in the effect of IS elements upon insertion, although none of the elements were preferentially associated with gene disruption. A survey of transposases among genomes of Actinobacteria showed no correlation between phylogenetic relatedness and distribution of IS families. By using Southern hybridization, we suggested that diversification of Lxc isolates is also mediated by insertion sequences in probably recent events. Conclusion: Collectively our data indicate that transposable elements are involved in genome diversification of Lxc and Lxx. The IS elements were probably acquired after the divergence of the two subspecies and are associated with genome organization and gene contents. In addition to enhancing understanding of IS element dynamics in general, these data will contribute to our ongoing comparative analyses aimed at understanding the biological differences of the Lxc and Lxx.

Differential expression of genes in salivary glands of male Rhipicephalus (Boophilus)microplus in response to infection with Anaplasma marginale

Background: Bovine anaplasmosis, caused by the rickettsial tick-borne pathogen Anaplasma marginale (Rickettsiales: Anaplasmataceae), is vectored by Rhipicephalus (Boophilus) microplus in many tropical and subtropical regions of the world. A. marginale undergoes a complex developmental cycle in ticks which results in infection of salivary glands from where the pathogen is transmitted to cattle. In previous studies, we reported modification of gene expression in Dermacentor variabilis and cultured Ixodes scapularis tick cells in response to infection with A. marginale. In these studies, we extended these findings by use of a functional genomics approach to identify genes differentially expressed in R. microplus male salivary glands in response to A. marginale infection. Additionally, a R. microplus-derived cell line, BME26, was used for the first time to also study tick cell gene expression in response to A. marginale infection. Results: Suppression subtractive hybridization libraries were constructed from infected and uninfected ticks and used to identify genes differentially expressed in male R. microplus salivary glands infected with A. marginale. A total of 279 ESTs were identified as candidate differentially expressed genes. Of these, five genes encoding for putative histamine-binding protein (22Hbp), von Willebrand factor (94Will), flagelliform silk protein (100Silk), Kunitz-like protease inhibitor precursor (108Kunz) and proline-rich protein BstNI subfamily 3 precursor (7BstNI3) were confirmed by real-time RT-PCR to be down-regulated in tick salivary glands infected with A. marginale. The impact of selected tick genes on A. marginale infections in tick salivary glands and BME26 cells was characterized by RNA interference. Silencing of the gene encoding for putative flagelliform silk protein (100Silk) resulted in reduced A. marginale infection in both tick salivary glands and cultured BME26 cells, while silencing of the gene encoding for subolesin (4D8) significantly reduced infection only in cultured BME26 cells. The knockdown of the gene encoding for putative metallothionein (93 Meth), significantly up-regulated in infected cultured BME26 cells, resulted in higher A. marginale infection levels in tick cells. Conclusions: Characterization of differential gene expression in salivary glands of R. microplus in response to A. marginale infection expands our understanding of the molecular mechanisms at the tick-pathogen interface. Functional studies suggested that differentially expressed genes encoding for subolesin, putative von Willebrand factor and flagelliform silk protein could play a role in A. marginale infection and multiplication in ticks. These tick genes found to be functionally relevant for tick-pathogen interactions will likely be candidates for development of vaccines designed for control of both ticks and tick-borne pathogens.

Nonmonotonic thickness dependence of spin wave energy in ultrathin Fe films: Experiment and theory

High wave-vector spin waves in ultrathin Fe/W(110) films up to 20 monolayers (MLs) thick have been studied using spin-polarized electron energy-loss spectroscopy. An unusual nonmonotonous dependence of the spin wave energies on the film thickness is observed, featuring a pronounced maximum at 2 ML coverage. First-principles theoretical study reveals the origin of this behavior to be in the localization of the spin waves at the surface of the film, as well as in the properties of the interlayer exchange coupling influenced by the hybridization of the electron states of the film and substrate and by the strain.

Anomalous Rashba spin-orbit interaction in InAs/GaSb quantum wells

We theoretically investigate the Rashba spin-orbit interaction in InAs/GaSb quantum wells (QWs). We find that the Rashba spin-splitting (RSS) sensitively depends on the thickness of the InAs layer. The RSS exhibits nonlinear behavior for narrow InAs/GaSb QWs and the oscillating feature for wide InAs/GaSb QWs. The nonlinear and oscillating behaviors arise from the weakened and enhanced interband coupling. The RSS also show asymmetric features respect to the direction of the external electric field. (C) 2008 American Institute of Physics.

Density functional theory and ab initio molecular dynamics study of NO adsorption on Pd(111) and Pt(111) surfaces

The origin of the unique geometry for nitric oxide (NO) adsorption on Pd(111) and Pt(111) surfaces as well as the effect of temperature were studied by density functional theory calculations and ab initio molecular dynamics at finite temperature. We found that at low coverage, the adsorption geometry is determined by electronic interactions, depending sensitively on the adsorption sites and coverages, and the effect of temperature on geometries is significant. At coverage of 0.25 monolayer (ML), adsorbed NO at hollow sites prefer an upright configuration, while NO adsorbed at top sites prefer a tilting configuration. With increase in the coverage up to 0.50 ML, the enhanced steric repulsion lead to the tilting of hollow NO. We found that the tilting was enhanced by the thermal effects. At coverage of 0.75 ML with p(2 x 2)-3NO(fcc+hcp+top) structure, we found that there was no preferential orientation for tilted top NO. The interplay of the orbital hybridization, thermal effects, steric repulsion, and their effects on the adsorption geometries were highlighted at the end.

The role of electron localization in the atomic structure of transition-metal 13-atom clusters: the example of Co(13), Rh(13), and Hf(13)

The crystalline structure of transition-metals (TM) has been widely known for several decades, however, our knowledge on the atomic structure of TM clusters is still far from satisfactory, which compromises an atomistic understanding of the reactivity of TM clusters. For example, almost all density functional theory (DFT) calculations for TM clusters have been based on local (local density approximation-LDA) and semilocal (generalized gradient approximation-GGA) exchange-correlation functionals, however, it is well known that plain DFT fails to correct the self-interaction error, which affects the properties of several systems. To improve our basic understanding of the atomic and electronic properties of TM clusters, we report a DFT study within two nonlocal functionals, namely, the hybrid HSE (Heyd, Scuseria, and Ernzerhof) and GGA + U functionals, of the structural and electronic properties of the Co(13), Rh(13), and Hf(13) clusters. For Co(13) and Rh(13), we found that improved exchange-correlation functionals decrease the stability of open structures such as the hexagonal bilayer (HBL) and double simple-cubic (DSC) compared with the compact icosahedron (ICO) structure, however, DFT-GGA, DFT-GGA + U, and DFT-HSE yield very similar results for Hf(13). Thus, our results suggest that the DSC structure obtained by several plain DFT calculations for Rh(13) can be improved by the use of improved functionals. Using the sd hybridization analysis, we found that a strong hybridization favors compact structures, and hence, a correct description of the sd hybridization is crucial for the relative energy stability. For example, the sd hybridization decreases for HBL and DSC and increases for ICO in the case of Co(13) and Rh(13), while for Hf(13), the sd hybridization decreases for all configurations, and hence, it does not affect the relative stability among open and compact configurations.

Identification of COL6A1 as a differentially expressed gene in human astrocytomas

Diffuse infiltrating gliomas are the most common tumors of the central nervous system. Gliomas are classified by the WHO according to their histopathological and clinical characteristics into four classes: grade I (pilocytic astrocytoma), grade II (diffuse astrocytoma), grade III (anaplastic astrocytoma), and grade IV (glioblastoma multiforme). Several genes have already been correlated with astrocytomas, but many others are yet to be uncovered. By analyzing the public SAGE data from 21 patients, comprising low malignant grade astrocytomas and glioblastomas, we found COL6A1 to be differentially expressed, confirming this finding by real time RT-PCR in 66 surgical samples. To the best of our knowledge, COL6A1 has never been described in gliomas. The expression of this gene has significantly different means when normal glia is compared with low-grade astrocytomas (grades I and II) and high-grade astrocytomas (grades III and IV), with a tendency to be greater in higher grade samples, thus rendering it a powerful tumor marker.

Haplotype distribution of five nuclear genes based on network genealogies and Bayesian inference indicates that Trypanosoma cruzi hybrid strains are polyphyletic

Chagas disease is still a major public health problem in Latin America. Its causative agent, Trypanosoma cruzi, can be typed into three major groups, T. cruzi I, T. cruzi II and hybrids. These groups each have specific genetic characteristics and epidemiological distributions. Several highly virulent strains are found in the hybrid group; their origin is still a matter of debate. The null hypothesis is that the hybrids are of polyphyletic origin, evolving independently from various hybridization events. The alternative hypothesis is that all extant hybrid strains originated from a single hybridization event. We sequenced both alleles of genes encoding EF-1 alpha, actin and SSU rDNA of 26 T. cruzi strains and DHFR-TS and TR of 12 strains. This information was used for network genealogy analysis and Bayesian phylogenies. We found T. cruzi I and T. cruzi II to be monophyletic and that all hybrids had different combinations of T. cruzi I and T. cruzi II haplotypes plus hybrid-specific haplotypes. Bootstrap values (networks) and posterior probabilities (Bayesian phylogenies) of clades supporting the monophyly of hybrids were far below the 95% confidence interval, indicating that the hybrid group is polyphyletic. We hypothesize that T. cruzi I and T. cruzi II are two different species and that the hybrids are extant representatives of independent events of genome hybridization, which sporadically have sufficient fitness to impact on the epidemiology of Chagas disease.