Next station in microarray data analysis: GEPAS

Sex-Specific Associations between Particulate Matter Exposure and Gene Expression in Independent Discovery and Validation Cohorts of Middle-Aged Men and Women

BACKGROUND

Particulate matter (PM) exposure leads to premature death, mainly due to respiratory and cardiovascular diseases.

OBJECTIVES

Identification of transcriptomic biomarkers of air pollution exposure and effect in a healthy adult population.

METHODS

Microarray analyses were performed in 98 healthy volunteers (48 men, 50 women). The expression of eight sex-specific candidate biomarker genes (significantly associated with PM₁₀ in the discovery cohort and with a reported link to air pollution-related disease) was measured with qPCR in an independent validation cohort (75 men, 94 women). Pathway analysis was performed using Gene Set Enrichment Analysis. Average daily PM_2.5 and PM₁₀ exposures over 2-years were estimated for each participant's residential address using spatiotemporal interpolation in combination with a dispersion model.

RESULTS

Average long-term PM₁₀ was 25.9 (± 5.4) and 23.7 (± 2.3) μg/m³ in the discovery and validation cohorts, respectively. In discovery analysis, associations between PM₁₀ and the expression of individual genes differed by sex. In the validation cohort, long-term PM₁₀ was associated with the expression of DNAJB5 and EAPP in men and ARHGAP4 (p = 0.053) in women. AKAP6 and LIMK1 were significantly associated with PM₁₀ in women, although associations differed in direction between the discovery and validation cohorts. Expression of the eight candidate genes in the discovery cohort differentiated between validation cohort participants with high versus low PM₁₀ exposure (area under the receiver operating curve = 0.92; 95% CI: 0.85, 1.00; p = 0.0002 in men, 0.86; 95% CI: 0.76, 0.96; p = 0.004 in women).

CONCLUSIONS

Expression of the sex-specific candidate genes identified in the discovery population predicted PM₁₀ exposure in an independent cohort of adults from the same area. Confirmation in other populations may further support this as a new approach for exposure assessment, and may contribute to the discovery of molecular mechanisms for PM-induced health effects.

A Method for the Annotation of Functional Similarities of Coding DNA Sequences: the Case of a Populated Cluster of Transmembrane Proteins

The analysis of a large number of human and mouse genes codifying for a populated cluster of transmembrane proteins revealed that some of the genes significantly vary in their primary nucleotide sequence inter-species and also intra-species. In spite of that divergence and of the fact that all these genes share a common parental function we asked the question of whether at DNA level they have some kind of common compositional structure, not evident from the analysis of their primary nucleotide sequence. To reveal the existence of gene clusters not based on primary sequence relationships we have analyzed 13574 human and 14047 mouse genes by the composon-clustering methodology. The data presented show that most of the genes from each one of the samples are distributed in 18 clusters sharing the common compositional features between the particular human and mouse clusters. It was observed, in addition, that between particular human and mouse clusters having similar composon-profiles large variations in gene population were detected as an indication that a significant amount of orthologs between both species differs in compositional features. A gene cluster containing exclusively genes codifying for transmembrane proteins, an important fraction of which belongs to the Rhodopsin G-protein coupled receptor superfamily, was also detected. This indicates that even though some of them display low sequence similarity, all of them, in both species, participate with similar compositional features in terms of composons. We conclude that in this family of transmembrane proteins in general and in the Rhodopsin G-protein coupled receptor in particular, the composon-clustering reveals the existence of a type of common compositional structure underlying the primary nucleotide sequence closely correlated to function.

Lmx1b is required for the glutamatergic fates of a subset of spinal cord neurons

Background

Alterations in neurotransmitter phenotypes of specific neurons can cause imbalances in excitation and inhibition in the central nervous system (CNS), leading to diseases. Therefore, the correct specification and maintenance of neurotransmitter phenotypes is vital. As with other neuronal properties, neurotransmitter phenotypes are often specified and maintained by particular transcription factors. However, the specific molecular mechanisms and transcription factors that regulate neurotransmitter phenotypes remain largely unknown.

Methods

In this paper we use single mutant, double mutant and transgenic zebrafish embryos to elucidate the functions of Lmx1ba and Lmx1bb in the regulation of spinal cord interneuron neurotransmitter phenotypes.

Results

We demonstrate that lmx1ba and lmx1bb are both expressed in zebrafish spinal cord and that lmx1bb is expressed by both V0v cells and dI5 cells. Our functional analyses demonstrate that these transcription factors are not required for neurotransmitter fate specification at early stages of development, but that in embryos with at least two lmx1ba and/or lmx1bb mutant alleles there is a reduced number of excitatory (glutamatergic) spinal interneurons at later stages of development. In contrast, there is no change in the numbers of V0v or dI5 cells. These data suggest that lmx1b-expressing spinal neurons still form normally, but at least a subset of them lose, or do not form, their normal excitatory fates. As the reduction in glutamatergic cells is only seen at later stages of development, Lmx1b is probably required either for the maintenance of glutamatergic fates or to specify glutamatergic phenotypes of a subset of later forming neurons. Using double labeling experiments, we also show that at least some of the cells that lose their normal glutamatergic phenotype are V0v cells. Finally, we also establish that Evx1 and Evx2, two transcription factors that are required for V0v cells to acquire their excitatory neurotransmitter phenotype, are also required for lmx1ba and lmx1bb expression in these cells, suggesting that Lmx1ba and Lmx1bb act downstream of Evx1 and Evx2 in V0v cells.

Conclusions

Lmx1ba and Lmx1bb function at least partially redundantly in the spinal cord and three functional lmx1b alleles are required in zebrafish for correct numbers of excitatory spinal interneurons at later developmental stages. Taken together, our data significantly enhance our understanding of how spinal cord neurotransmitter fates are regulated.

Bayesian Network Inference Enables Unbiased Phenotypic Anchoring of Transcriptomic Responses to Cigarette Smoke in Humans

Microarray-based transcriptomic analysis has been demonstrated to hold the opportunity to study the effects of human exposure to, e.g., chemical carcinogens at the whole genome level, thus yielding broad-ranging molecular information on possible carcinogenic effects. Since genes do not operate individually but rather through concerted interactions, analyzing and visualizing networks of genes should provide important mechanistic information, especially upon connecting them to functional parameters, such as those derived from measurements of biomarkers for exposure and carcinogenic risk. Conventional methods such as hierarchical clustering and correlation analyses are frequently used to address these complex interactions but are limited as they do not provide directional causal dependence relationships. Therefore, our aim was to apply Bayesian network inference with the purpose of phenotypic anchoring of modified gene expressions. We investigated a use case on transcriptomic responses to cigarette smoking in humans, in association with plasma cotinine levels as biomarkers of exposure and aromatic DNA-adducts in blood cells as biomarkers of carcinogenic risk. Many of the genes that appear in the Bayesian networks surrounding plasma cotinine, and to a lesser extent around aromatic DNA-adducts, hold biologically relevant functions in inducing severe adverse effects of smoking. In conclusion, this study shows that Bayesian network inference enables unbiased phenotypic anchoring of transcriptomics responses. Furthermore, in all inferred Bayesian networks several dependencies are found which point to known but also to new relationships between the expression of specific genes, cigarette smoke exposure, DNA damaging-effects, and smoking-related diseases, in particular associated with apoptosis, DNA repair, and tumor suppression, as well as with autoimmunity.

Two distinct molecular subtypes of chronic lymphocytic leukemia give new insights on the pathogenesis of the disease and identify novel therapeutic targets

Biopsy samples of lymph nodes from 38 patients with CLL were analyzed. We found differential expression in 1092 genes in two different subgroups: 418 overexpressed in one subgroup and 674 in another. Molecular pathways identified in one subgroup appear to be characterized by greater dependence of signaling by cytokines and activation of the NFkB pathway, while in the other seem to depend on cell cycle. Despite having found a differential expression between both subgroups, none of these genes reached FDR < 0.25. We have not found significant association with survival or any prognostic factors. Analysis of the differences between normal lymph node and CLL in 253 genes with difference in the intensity of expression revealed upregulated genes different to BCR: CD40, TCL1, IL-7, and PAX5. Using large-scale molecular analysis, we may obtain information about molecular mechanisms of CLL pathogenesis and may contribute to the identification of new therapeutic targets.

Babelomics 5.0: functional interpretation for new generations of genomic data

Babelomics has been running for more than one decade offering a user-friendly interface for the functional analysis of gene expression and genomic data. Here we present its fifth release, which includes support for Next Generation Sequencing data including gene expression (RNA-seq), exome or genome resequencing. Babelomics has simplified its interface, being now more intuitive. Improved visualization options, such as a genome viewer as well as an interactive network viewer, have been implemented. New technical enhancements at both, client and server sides, makes the user experience faster and more dynamic. Babelomics offers user-friendly access to a full range of methods that cover: (i) primary data analysis, (ii) a variety of tests for different experimental designs and (iii) different enrichment and network analysis algorithms for the interpretation of the results of such tests in the proper functional context. In addition to the public server, local copies of Babelomics can be downloaded and installed. Babelomics is freely available at: http://www.babelomics.org.

Identification of stage-specific markers during differentiation of hair cells from mouse inner ear stem cells or progenitor cells in vitro

The induction of inner ear hair cells from stem cells or progenitor cells in the inner ear proceeds through a committed inner ear sensory progenitor cell stage prior to hair cell differentiation. To increase the efficacy of inducing inner ear hair cell differentiation from the stem cells or progenitor cells, it is essential to identify comprehensive markers for the stem cells/progenitor cells from the inner ear, the committed inner ear sensory progenitor cells and the differentiating hair cells to optimize induction conditions. Here, we report that we efficiently isolated and expanded the stem cells or progenitor cells from postnatal mouse cochleae, and induced the generation of inner ear progenitor cells and subsequent differentiation of hair cells. We profiled the gene expression of the stem cells or progenitor cells, the inner ear progenitor cells, and hair cells using aRNA microarray analysis. The pathway and gene ontology (GO) analysis of differentially expressed genes was performed. Analysis of genes exclusively detected in one particular cellular population revealed 30, 38, and 31 genes specific for inner ear stem cells, inner ear progenitor cells, and hair cells, respectively. We further examined the expression of these genes in vivo and determined that Gdf10+Ccdc121, Tmprss9+Orm1, and Chrna9+Espnl are marker genes specific for inner ear stem cells, inner ear progenitor cells, and differentiating hair cells, respectively. The identification of these marker genes will likely help the effort to increase the efficacy of hair cell induction from the stem cells or progenitor cells.

Up-regulation of anti-apoptotic genes confers resistance to the novel anti-leukaemic compound PEP005 in primary AML cells

We showed previously that PEP005 induced apoptosis in leukaemic cell lines and blasts from patients with acute myeloid leukaemia (AML). Here we assess the anti-leukeamic effects of PEP005 in vivo and determine the mechanism of resistance of PEP005 non-responsive cells. We used 2 human xenograft mouse models of AML to assess the anti-leukaemic effects of PEP005 in vivo. Expression microarray analysis of primary AML blasts following treatment with PEP005 was used to determine patterns of gene expression that conferred resistance. PEP005 significantly reduced tumour burden in two human leukaemia mouse xenograft models. We also assessed responsiveness of 33 AML samples to PEP005, with 78% of the samples entering apoptosis at 100nM. Resistance to PEP005 was not restricted to a particular AML subtype. Expression microarray analysis of resistant samples following treatment with PEP005 revealed a significant up regulation of the anti-apoptotic genes Bcl-2A1, Mcl-1, and PHLDA1 which was verified using RT-PCR. We conclude that PEP005 shows broad efficacy against AML subtypes and that up regulation of anti-apoptotic genes underlies resistance to this agent and could be used to screen for patients unlikely to benefit from a therapeutic regime involving PEP005.

Uncaria rhynchophylla and rhynchophylline improved kainic acid-induced epileptic seizures via IL-1β and brain-derived neurotrophic factor

Uncaria rhynchophylla (UR) has been used for the treatment of convulsions and epilepsy in traditional Chinese medicine. This study reported the major anti-convulsive signaling pathways and effective targets of UR and rhynchophylline (RP) using genomic and immunohistochemical studies. Epileptic seizure model was established by intraperitoneal injection of kainic acid (KA) in rats. Electroencephalogram and electromyogram recordings indicated that UR and RP improved KA-induced epileptic seizures. Toll-like receptor (TLR) and neurotrophin signaling pathways were regulated by UR in both cortex and hippocampus of KA-treated rats. KA upregulated the expression levels of interleukin-1β (IL-1β) and brain-derived neurotrophin factor (BDNF), which were involved in TLR and neurotrophin signaling pathways, respectively. However, UR and RP downregulated the KA-induced IL-1β and BDNF gene expressions. Our findings suggested that UR and RP exhibited anti-convulsive effects in KA-induced rats via the regulation of TLR and neurotrophin signaling pathways, and the subsequent inhibition of IL-1β and BDNF gene expressions.

Genomic analysis of clonal eosinophils by CGH arrays reveals new genetic regions involved in chronic eosinophilia

To assess the presence of genetic imbalances in patients with myeloproliferative neoplasms (MPNs), 38 patients with chronic eosinophilia were studied by array comparative genomic hybridization (aCGH): seven had chronic myelogenous leukaemia (CML), BCR-ABL1 positive, nine patients had myeloproliferative neoplasia Ph- (MPN-Ph-), three had a myeloid neoplasm associated with a PDGFRA rearrangement, and the remaining two cases were Lymphoproliferative T neoplasms associated with eosinophilia. In addition, 17 patients had a secondary eosinophilia and were used as controls. Eosinophilic enrichment was carried out in all cases. Genomic imbalances were found in 76% of all MPN patients. Losses on 20q were the most frequent genetic abnormality in MPNs (32%), affected the three types of MPN studied. This study also found losses at 11q13.3 in 26% of patients with MPN-Ph- and in 19p13.11 in two of the three patients with an MPN associated with a PDGFRA rearrangement. In addition, 29% of patients with CML had losses on 8q24. In summary, aCGH revealed clonality in eosinophils in most MPNs, suggesting that it could be a useful technique for defining clonality in these diseases. The presence of genetic losses in new regions could provide new insights into the knowledge of these MPN associated with eosinophilia.

MiR-210 is induced by Oct-2, regulates B cells, and inhibits autoantibody production

MicroRNAs (MiRs) are small, noncoding RNAs that regulate gene expression posttranscriptionally. In this study, we show that MiR-210 is induced by Oct-2, a key transcriptional mediator of B cell activation. Germline deletion of MiR-210 results in the development of autoantibodies from 5 mo of age. Overexpression of MiR-210 in vivo resulted in cell autonomous expansion of the B1 lineage and impaired fitness of B2 cells. Mice overexpressing MiR-210 exhibited impaired class-switched Ab responses, a finding confirmed in wild-type B cells transfected with a MiR-210 mimic. In vitro studies demonstrated defects in cellular proliferation and cell cycle entry, which were consistent with the transcriptomic analysis demonstrating downregulation of genes involved in cellular proliferation and B cell activation. These findings indicate that Oct-2 induction of MiR-210 provides a novel inhibitory mechanism for the control of B cells and autoantibody production.

Gender-specific transcriptomic response to environmental exposure in Flemish adults

Flanders, Belgium, is one of the most densely populated areas in Europe. The Flemish Environment and Health Survey (2002-2006) aimed at determining exposure to pollutants of neonates, adolescents, and older adults and to assess associated biological and health effects. This study investigated genome wide gene expression changes associated with a range of environmental pollutants, including cadmium, lead, PCBs, dioxin, hexachlorobenzene, p,p'-DDE, benzene, and PAHs. Gene expression levels were measured in peripheral blood cells of 20 adults with relatively high and 20 adults with relatively low combined internal exposure levels, all non-smokers aged 50-65. Pearson correlation was used to analyze associations between pollutants and gene expression levels, separately for both genders. Pollutant- and gender-specific correlation analysis results were obtained. For organochlorine pollutants, analysis within genders revealed that genes were predominantly regulated in opposite directions in males and females. Significantly modulated pathways were found to be associated with each of the exposure biomarkers measured. Pathways and/or genes related to estrogen and STAT5 signaling were correlated to organochlorine exposures in both genders. Our work demonstrates that gene expression in peripheral blood is influenced by environmental pollutants. In particular, gender-specific changes are associated with organochlorine pollutants, including gender-specific modulation of endocrine related pathways and genes. These pathways and genes have previously been linked to endocrine disruption related disorders, which in turn have been associated with organochlorine exposure. Based on our results, we recommend that males and females be considered separately when analyzing gene expression changes associated with exposures that may include chemicals with endocrine disrupting properties.

Clinical applications of BAC array-CGH to the study of diffuse large B-cell lymphomas

BAC array-CGH is a powerful method to identify DNA copy number changes (gains, amplifications and deletions) on a genome-wide scale, and to map these changes to genomic sequence. It is based on the analysis of genomic DNA isolated from test and reference cell populations, the differential labelling with fluorescent dyes and the co-hybridization with a genomic array. BAC array-CGH has proven to be a specific, sensitive, and reliable technique, with considerable advantages compared to other methods used for the analysis of DNA copy number changes. The application of genome scanning technologies and the recent advances in bioinformatics tools that enable us to perform a robust and highly sensitive analysis of array-CGH data, useful not only for genome scanning of tumor cells but also in the identification of novel cancer related genes, oncogenes and suppressor genes. Cytogenetic analysis provides essential information for diagnosis and prognosis in patients with hematologic malignancies such as lymphomas. However, the chromosomal interpretation in non-Hodgkin lymphoma (NHL) is sometimes inconclusive. Copy number aberrations identified by BAC array-CGH analyses could be a complementary methodology to chromosomal analysis. In NHL the genomic imbalances might have a prognostic rather than a diagnostic value. In fact, the diagnosis of NHL is based on pathological and molecular cytogenetics data. Furthermore genetic variations and their association with specific types of lymphoma development, and elucidation of the variable genetic pathways leading to lymphoma development, are important directions for future cancer research. Array-CGH, along with FISH and PCR, will be used for routine diagnostic purposes in near future.

Differential gene expression of medullary thyroid carcinoma reveals specific markers associated with genetic conditions

Medullary thyroid carcinoma accounts for 2% to 5% of thyroid malignancies, of which 75% are sporadic and the remaining 25% are hereditary and related to multiple endocrine neoplasia type 2 syndrome. Despite a genotype-phenotype correlation with specific germline RET mutations, knowledge of pathways specifically associated with each mutation and with non-RET-mutated sporadic MTC remains lacking. Gene expression patterns have provided a tool for identifying molecular events related to specific tumor types and to different clinical features that could help identify novel therapeutic targets. Using transcriptional profiling of 49 frozen MTC specimens classified as RET mutation, we identified PROM1, LOXL2, GFRA1, and DKK4 as related to RET(M918T) and GAL as related to RET(634) mutation. An independent series of 19 frozen and 23 formalin-fixed, paraffin-embedded (FFPE) MTCs was used for validation by RT-qPCR. Two tissue microarrays containing 69 MTCs were available for IHC assays. According to pathway enrichment analysis and gene ontology biological processes, genes associated with the MTC(M918T) group were involved mainly in proliferative, cell adhesion, and general malignant metastatic effects and with Wnt, Notch, NFκB, JAK/Stat, and MAPK signaling pathways. Assays based on silencing of PROM1 by siRNAs performed in the MZ-CRC-1 cell line, harboring RET(M918T), caused an increase in apoptotic nuclei, suggesting that PROM1 is necessary for survival of these cells. This is the first report of PROM1 overexpression among primary tumors.

Unifying the integration, analysis and interpretation of multi-omic datasets: exploration of the disease networks of Obstructive Nephropathy in children

The wealth of data amassed by the utilization of various high-throughput techniques, in various layers of molecular dissection, stresses the critical role of the unification of the computational methodologies applied in biological data handling, storage, analysis and visualization. In this article, a generic workflow is showcased in a multi-omic dataset that is used to study Obstructive Nephropathy (ON) in children, integrating microarray data from several biological layers (transcriptomic, post-transcriptomic, proteomic). The workflow exploits raw measurements and through several analytical stages (preprocessing, statistical and functional), which entail various parsing steps, reaches the visualization stage of the heterogeneous, broader, molecular interacting network derived. This network, where the interconnected entities are exploiting the knowledge stored in public repositories, represents a systems level interpretation of the pathological state probed.

Genes are often sheltered from the global histone hyperacetylation induced by HDAC inhibitors

Histone deacetylase inhibitors (HDACi) are increasingly used as therapeutic agents, but the mechanisms by which they alter cell behaviour remain unclear. Here we use microarray expression analysis to show that only a small proportion of genes (∼9%) have altered transcript levels after treating HL60 cells with different HDACi (valproic acid, Trichostatin A, suberoylanilide hydroxamic acid). Different gene populations respond to each inhibitor, with as many genes down- as up-regulated. Surprisingly, HDACi rarely induced increased histone acetylation at gene promoters, with most genes examined showing minimal change, irrespective of whether genes were up- or down-regulated. Many genes seem to be sheltered from the global histone hyperacetyation induced by HDACi.

Identification of novel mechanisms of silymarin on the carbon tetrachloride-induced liver fibrosis in mice by nuclear factor-κB bioluminescent imaging-guided transcriptomic analysis

In this study, we applied bioluminescent imaging-guided transcriptomic analysis to evaluate and identify the therapeutic potentials and novel mechanisms of silymarin on carbon tetrachloride (CCl(4))-induced liver fibrosis. Transgenic mice, carrying the luciferase genes driven by nuclear factor-κB (NF-κB), were given with CCl(4) and/or silymarin. In vivo NF-κB activity was evaluated by bioluminescent imaging, liver fibrosis was judged by Sirius red staining and immunohistochemistry, and gene expression profiles of silymarin-treated livers were analyzed by DNA microarray. CCl(4) enhanced the NF-κB-dependent hepatic luminescence and induced hepatic fibrosis, while silymarin reduced the CCl(4)-induced hepatic luminescence and improved CCl(4)-induced liver fibrosis. Microarray analysis showed that silymarin altered the transforming growth factor-β-mediated pathways, which play pivotal roles in the progression of liver fibrosis. Moreover, we newly identified that silymarin downregulated the expression levels of cytoskeleton organization genes and mitochondrion electron-transfer chain genes, such as cytochrome c oxidase Cox6a2, Cox7a1, and Cox8b genes. In conclusion, the correlation of NF-κB-dependent luminescence and liver fibrosis suggested the feasibility of NF-κB bioluminescent imaging for the evaluation of liver fibrosis progression and therapeutic potentials. Moreover, our findings suggested that silymarin might exhibit anti-fibrotic effects in vivo via altering the expression of genes involved in cytoskeleton organization and mitochondrion electron-transfer chain.

Transcriptome profiling of sexual maturation and mating in the Mediterranean fruit fly, Ceratitis capitata

Sexual maturation and mating in insects are generally accompanied by major physiological and behavioural changes. Many of these changes are related to the need to locate a mate and subsequently, in the case of females, to switch from mate searching to oviposition behaviour. The prodigious reproductive capacity of the Mediterranean fruit fly, Ceratitis capitata, is one of the factors that has led to its success as an invasive pest species. To identify the molecular changes related to maturation and mating status in male and female medfly, a microarray-based gene expression approach was used to compare the head transcriptomes of sexually immature, mature virgin, and mated individuals. Attention was focused on the changes in abundance of transcripts related to reproduction, behaviour, sensory perception of chemical stimulus, and immune system processes. Broad transcriptional changes were recorded during female maturation, while post-mating transcriptional changes in females were, by contrast, modest. In male medfly, transcriptional changes were consistent both during maturation and as a consequence of mating. Of particular note was the lack of the mating-induced immune responses that have been recorded for Drosophila melanogaster, that may be due to the different reproductive strategies of these species. This study, in addition to increasing our understanding of the molecular machinery behind maturation and mating in the medfly, has identified important gene targets that might be useful in the future management of this pest.

Histone modifications and expression of DAM6 gene in peach are modulated during bud dormancy release in a cultivar-dependent manner

• Bud dormancy release in many woody perennial plants responds to the seasonal accumulation of chilling stimulus. MADS-box transcription factors encoded by DORMANCY ASSOCIATED MADS-box (DAM) genes in peach (Prunus persica) are implicated in this pathway, but other regulatory factors remain to be identified. In addition, the regulation of DAM gene expression is not well known at the molecular level. • A microarray hybridization approach was performed to identify genes whose expression correlates with the bud dormancy-related behaviour in 10 different peach cultivars. Histone modifications in DAM6 gene were investigated by chromatin immunoprecipitation in two different cultivars. • The expression of DAM4-DAM6 and several genes related to abscisic acid and drought stress response correlated with the dormancy behaviour of peach cultivars. The trimethylation of histone H3 at K27 in the DAM6 promoter, coding region and the second large intron was preceded by a decrease in acetylated H3 and trimethylated H3K4 in the region of translation start, coinciding with repression of DAM6 during dormancy release. • Analysis of chromatin modifications reinforced the role of epigenetic mechanisms in DAM6 regulation and bud dormancy release, and highlighted common features with the vernalization process in Arabidopsis thaliana and cereals.

Thyroid paraganglioma. Report of 3 cases and description of an immunohistochemical profile useful in the differential diagnosis with medullary thyroid carcinoma, based on complementary DNA array results

Thyroid paraganglioma is a rare disorder that sometimes poses problems in differential diagnosis with medullary thyroid carcinoma. So far, differential diagnosis is solved with the help of some markers that are frequently expressed in medullary thyroid carcinoma (thyroid transcription factor 1, calcitonin, and carcinoembryonic antigen). However, some of these markers are not absolutely specific of medullary thyroid carcinoma and may be expressed in other tumors. Here we report 3 new cases of thyroid paraganglioma and describe our strategy to design a diagnostic immunohistochemical battery. First, we performed a comparative analysis of the expression profile of head and neck paragangliomas and medullary thyroid carcinoma, obtained after complementary DNA array analysis of 2 series of fresh-frozen samples of paragangliomas and medullary thyroid carcinoma, respectively. Seven biomarkers showing differential expression were selected (nicotinamide adenine dinucleotide dehydrogenase 1 alpha subcomplex, 4-like 2, NDUFA4L2; cytochrome c oxidase subunit IV isoform 2; vesicular monoamine transporter 2; calcitonin gene-related protein/calcitonin; carcinoembryonic antigen; and thyroid transcription factor 1) for immunohistochemical analysis. Two tissue microarrays were constructed from 2 different series of paraffin-embedded samples of paragangliomas and medullary thyroid carcinoma. We provide a classifying rule for differential diagnosis that combines negativity or low staining for calcitonin gene-related protein (histologic score, <10) or calcitonin (histologic score, <50) together with positivity of any of NADH dehydrogenase 1 alpha subcomplex, 4-like 2; cytochrome c oxidase subunit IV isoform 2; or vesicular monoamine transporter 2 to predict paragangliomas, showing a prediction error of 0%. Finally, the immunohistochemical battery was checked in paraffin-embedded blocks from 4 examples of thyroid paraganglioma (1 previously reported case and 3 new cases), showing also a prediction error of 0%. Our results suggest that the comparative expression profile, obtained by complementary DNA arrays, seems to be a good tool to design immunohistochemical batteries used in differential diagnosis.

Transcriptomic profiling of citrus fruit peel tissues reveals fundamental effects of phenylpropanoids and ethylene on induced resistance

Penicillium spp. are the major postharvest pathogens of citrus fruit in Mediterranean climatic regions. The induction of natural resistance constitutes one of the most promising alternatives to avoid the environmental contamination and health problems caused by chemical fungicides. To understand the bases of the induction of resistance in citrus fruit against Penicillium digitatum, we have used a 12k citrus cDNA microarray to study transcriptional changes in the outer and inner parts of the peel (flavedo and albedo, respectively) of elicited fruits. The elicitor treatment led to an over-representation of biological processes associated with secondary metabolism, mainly phenylpropanoids and cellular amino acid biosynthesis and methionine metabolism, and the down-regulation of genes related to biotic and abiotic stresses. Among phenylpropanoids, we detected the over-expression of a large subset of genes important for the synthesis of flavonoids, coumarins and lignin, especially in the internal tissue. Furthermore, these genes and those of ethylene biosynthesis showed the highest induction. The involvement of both phenylpropanoid and ethylene pathways was confirmed by examining changes in gene expression and ethylene production in elicited citrus fruit. Therefore, global results indicate that secondary metabolism, mainly phenylpropanoids, and ethylene play important roles in the induction of resistance in citrus fruit.

Defence responses of Arabidopsis thaliana to infection by Pseudomonas syringae are regulated by the circadian clock

The circadian clock allows plants to anticipate predictable daily changes in abiotic stimuli, such as light; however, whether the clock similarly allows plants to anticipate interactions with other organisms is unknown. Here we show that Arabidopsis thaliana (Arabidopsis) has circadian clock-mediated variation in resistance to the virulent bacterial pathogen Pseudomonas syringae pv. tomato DC3000 (Pst DC3000), with plants being least susceptible to infection in the subjective morning. We suggest that the increased resistance to Pst DC3000 observed in the morning in Col-0 plants results from clock-mediated modulation of pathogen associated molecular pattern (PAMP)-triggered immunity. Analysis of publicly available microarray data revealed that a large number of Arabidopsis defence-related genes showed both diurnal- and circadian-regulation, including genes involved in the perception of the PAMP flagellin which exhibit a peak in expression in the morning. Accordingly, we observed that PAMP-triggered callose deposition was significantly higher in wild-type plants inoculated with Pst DC3000 hrpA in the subjective morning than in the evening, while no such temporal difference was evident in arrhythmic plants. Our results suggest that PAMP-triggered immune responses are modulated by the circadian clock and that temporal regulation allows plants to anticipate and respond more effectively to pathogen challenges in the daytime.

Molecular characterization of the region 7q22.1 in splenic marginal zone lymphomas

Splenic marginal zone lymphomas (SMZL) are an uncommon type of B-cell non-Hodgkin's lymphoma (NHL-B) in which no specific chromosomal translocations have been described. In contrast, the most frequent cytogenetic abnormality is the loss of the long arm of chromosome 7 (7q). Previous reports have located this loss in the 7q32 region. In order to better characterize the genomic imbalances in SMZL, molecular studies were carried out in 73 patients with SMZL. To gain insight into the mapping at 7q a tiling array was also used. The results confirmed the loss of 7q as the most frequent change. In addition, several abnormalities, including 4q22.1, 1q21.3–q22, 6q25.3, 20q13.33, 3q28, 2q23.3–q24.1 and 17p13, were also present. A loss of 7q22.1 at 99925039–101348479 bp was observed in half of the cases. The region of 7q22.1 has not previously been characterised in SMZL. Our results confirmed the presence of a new region of loss on chromosome 7 in these NHL.

Analysis of normal-tumour tissue interaction in tumours: prediction of prostate cancer features from the molecular profile of adjacent normal cells

Statistical modelling, in combination with genome-wide expression profiling techniques, has demonstrated that the molecular state of the tumour is sufficient to infer its pathological state. These studies have been extremely important in diagnostics and have contributed to improving our understanding of tumour biology. However, their importance in in-depth understanding of cancer patho-physiology may be limited since they do not explicitly take into consideration the fundamental role of the tissue microenvironment in specifying tumour physiology. Because of the importance of normal cells in shaping the tissue microenvironment we formulate the hypothesis that molecular components of the profile of normal epithelial cells adjacent the tumour are predictive of tumour physiology. We addressed this hypothesis by developing statistical models that link gene expression profiles representing the molecular state of adjacent normal epithelial cells to tumour features in prostate cancer. Furthermore, network analysis showed that predictive genes are linked to the activity of important secreted factors, which have the potential to influence tumor biology, such as IL1, IGF1, PDGF BB, AGT, and TGFβ.

Gene expression profiling of subcutaneous adipose tissue in morbid obesity using a focused microarray: distinct expression of cell-cycle- and differentiation-related genes

Background

Obesity results from an imbalance between food intake and energy expenditure, which leads to an excess of adipose tissue. The excess of adipose tissue and adipocyte dysfunction associated with obesity are linked to the abnormal regulation of adipogenesis. The objective of this study was to analyze the expression profile of cell-cycle- and lipid-metabolism-related genes of adipose tissue in morbid obesity.

Methods

We used a custom-made focused cDNA microarray to determine the adipose tissue mRNA expression profile. Gene expression of subcutaneous abdominal fat samples from 15 morbidly obese women was compared with subcutaneous fat samples from 10 nonobese control patients. The findings were validated in an independent population of 31 obese women and 9 obese men and in an animal model of obesity (Lep^ob/ob mice) by real-time RT-PCR.

Results

Microarray analysis revealed that transcription factors that regulate the first stages of adipocyte differentiation, such as CCAAT/enhancer binding protein beta (C/EBPβ) and JUN, were upregulated in the adipose tissues of morbidly obese patients. The expression of peroxisome proliferator-activated receptor gamma (PPARγ), a transcription factor which controls lipid metabolism and the final steps of preadipocyte conversion into mature adipocytes, was downregulated. The expression of three cyclin-dependent kinase inhibitors that regulate clonal expansion and postmitotic growth arrest during adipocyte differentiation was also altered in obese subjects: p18 and p27 were downregulated, and p21 was upregulated. Angiopoietin-like 4 (ANGPTL4), which regulates angiogenesis, lipid and glucose metabolism and it is know to increase dramatically in the early stages of adipocyte differentiation, was upregulated. The expression of C/EBPβ, p18, p21, JUN, and ANGPTL4 presented similar alterations in subcutaneous adipose tissue of Lep^ob/ob mice.

Conclusions

Our microarray gene profiling study revealed that the expression of genes involved in adipogenesis is profoundly altered in the subcutaneous adipose tissue of morbidly obese subjects. This expression pattern is consistent with an immature adipocyte phenotype that could reflect the expansion of the adipose tissue during obesity.

Introns form compositional clusters in parallel with the compositional clusters of the coding sequences to which they pertain

This report deals with the study of compositional properties of human gene sequences evaluating similarities and differences among functionally distinct sectors of the gene independently of the reading frame. To retrieve the compositional information of DNA, we present a neighbor base dependent coding system in which the alphabet of 64 letters (DNA triplets) is compressed to an alphabet of 14 letters here termed triplet composons. The triplets containing the same set of distinct bases in whatever order and number form a triplet composon. The reading of the DNA sequence is performed starting at any letter of the initial triplet and then moving, triplet-to-triplet, until the end of the sequence. The readings were made in an overlapping way along the length of the sequences. The analysis of the compositional content in terms of the composon usage frequencies of the gene sequences shows that: (i) the compositional content of the sequences is far from that of random sequences, even in the case of non-protein coding sequences; (ii) coding sequences can be classified as components of compositional clusters; and (iii) intron sequences in a cluster have the same composon usage frequencies, even as their base composition differs notably from that of their home coding sequences. A comparison of the composon usage frequencies between human and mouse homologous genes indicated that two clusters found in humans do not have their counterpart in mouse whereas the others clusters are stable in both species with respect to their composon usage frequencies in both coding and noncoding sequences.

Research resource: Transcriptional profiling reveals different pseudohypoxic signatures in SDHB and VHL-related pheochromocytomas

The six major genes involved in hereditary susceptibility for pheochromocytoma (PCC)/paraganglioma (PGL) (RET, VHL, NF1, SDHB, SDHC, and SDHD) have been recently integrated into the same neuronal apoptotic pathway where mutations in any of these genes lead to cell death. In this model, prolyl hydroxylase 3 (EglN3) abrogation plays a pivotal role, but the molecular mechanisms underlying its inactivation are currently unknown. The aim of the study was to decipher specific alterations associated with the different genetic classes of PCCs/PGLs. With this purpose, 84 genetically characterized tumors were analyzed by means of transcriptional profiling. The analysis revealed a hypoxia-inducible factor (HIF)-related signature common to succinate dehydrogenase (SDH) and von Hippel-Lindau (VHL) tumors, that differentiated them from RET and neurofibromatosis type 1 cases. Both canonical HIF-1α and HIF-2α target genes were overexpressed in the SDH/VHL cluster, suggesting that a global HIF deregulation accounts for this common profile. Nevertheless, when we compared VHL tumors with SDHB cases, which often exhibit a malignant behavior, we found that HIF-1α target genes showed a predominant activation in the VHL PCCs. Expression data from 67 HIF target genes was sufficient to cluster SDHB and VHL tumors into two different groups, demonstrating different pseudo-hypoxic signatures. In addition, VHL-mutated tumors showed an unexpected overexpression of EglN3 mRNA that did not lead to significantly different EglN3 protein levels. These findings pave the way for more specific therapeutic approaches for malignant PCCs/PGLs management based on the patient's genetic alteration.

Transcriptional characteristics of familial non-BRCA1/BRCA2 breast tumors

To better understand the alterations present in the group of the so-called BRCAX tumors, we have used a cDNA microarray containing genes related to tumorigenesis and analyzed a series of 49 tumors consisting of 13 BRCA1, 14 BRCAX and 22 sporadic. We have confirmed that the BRCAX tumors are heterogeneous and can be divided in at least two main subgroups, so-called A and B, transcriptionally distinguishable and with different altered pathways within each of the groups. We have found that BRCAX-A and B subgroups, can be classified as Luminal A and Luminal B, respectively, taking into account the intrinsic phenotypes defined for the sporadic breast tumors. We have found that, at the somatic level, the BRCAX-B tumors are identical to their sporadic Luminal B counterparts, whereas BRCAX-A, despite having a Luminal A phenotype, shows additional genomic alterations. We have found 21 deregulated genes in the BRCAX-A group that we have called "the BRCAX susceptibility pathway" and suggested it as a candidate to search for new genes involved in the inherited susceptibility underlying the disease in this group.

There is a steady-state transcriptome in exponentially growing yeast cells

The growth of yeast cells in batches in glucose-based media is a standard condition in most yeast laboratories. Most gene expression experiments are done by taking this condition as a reference. Presumably, cells are in a stable physiological condition that can be easily reproduced in other laboratories. With this assumption, however, it is necessary to consider that the average amount of the mRNAs per cell for most genes does not change during exponential growth. That is to say, there is a steady-state condition for the transcriptome. However, this has not been rigorously demonstrated to date. In this work we take several cell samples during the exponential phase growth to perform a kinetic study using the genomic run-on (GRO) technique, which allows simultaneous measurement of the amount of mRNA and transcription rate variation at the genomic level. We show here that the steady-state condition is fulfilled for almost all the genes during most exponential growth in yeast extract-peptone-dextrose medium (YPD) and, therefore, that simultaneous measures of the transcription rates and the amounts of mRNA can be used for indirect mRNA stability calculations. With this kinetic approach, we were also able to determine the relative influence of the transcription rate and the mRNA stability changes for the mRNA variation for those genes that deviate from the steady state.

Seasonal reorganization of the xylem transcriptome at different tree ages reveals novel insights into wood formation in Pinus radiata

*Seasonal wood development produces earlywood (EW) and latewood (LW) with distinct properties. The molecular mechanisms controlling EW and LW formation at different tree ages are poorly understood. *Seasonal reorganization of the xylem transcriptome was investigated in Pinus radiata at four tree ages using cDNA microarrays. Transcriptome profiles were compared with seasonal wood variation measured by SilviScan (CSIRO, Clayton, Australia). *The xylem transcriptome was considerably reorganized during seasonal change, and this reorganization showed a maturation-related pattern. The greater reorganization occurred at the transition (30%) and juvenile (21%) stages, but it declined with tree maturity (11-13%). However, this pattern does not correlate well with maturation-related patterns of seasonal wood variation. In total, 319 seasonal-responsive xylem candidate genes were identified. Many transcripts involved in primary and secondary wall biosynthesis were preferentially accumulated in EW and LW, respectively. A large proportion (45-81%) of the candidate genes are preferentially regulated at a single age and their transcript abundance may influence maturation-related patterns of seasonal wood variation. *Seasonal reorganization of the xylem transcriptome is significantly affected by tree age. Physiological changes at the transition stage may contribute to its greater seasonal transcriptome reorganization. Identified stage-preferential xylem transcripts could influence seasonal wood variation at different tree ages.

Deregulated expression of the polycomb-group protein SUZ12 target genes characterizes mantle cell lymphoma

Polycomb proteins are known to be of great importance in human cancer pathogenesis. SUZ12 is a component of the Polycomb PRC2 complex that, along with EZH2, is involved in embryonic stem cell differentiation. EZH2 plays an essential role in many cancer types, but an equivalent involvement of SUZ12 has not been as thoroughly demonstrated. Here we show that SUZ12 is anomalously expressed in human primary tumors, especially in mantle cell lymphoma (MCL), pulmonary carcinomas and melanoma, and is associated with gene locus amplification in some cases. Using MCL as a model, functional and genomic studies demonstrate that SUZ12 loss compromises cell viability, increases apoptosis, and targets genes involved in central oncogenic pathways associated with MCL pathogenesis. Our results support the hypothesis that the abnormal expression of SUZ12 accounts for some of the unexplained features of MCL, such as abnormal DNA repair and increased resistance to apoptosis.

Babelomics: an integrative platform for the analysis of transcriptomics, proteomics and genomic data with advanced functional profiling

Babelomics is a response to the growing necessity of integrating and analyzing different types of genomic data in an environment that allows an easy functional interpretation of the results. Babelomics includes a complete suite of methods for the analysis of gene expression data that include normalization (covering most commercial platforms), pre-processing, differential gene expression (case-controls, multiclass, survival or continuous values), predictors, clustering; large-scale genotyping assays (case controls and TDTs, and allows population stratification analysis and correction). All these genomic data analysis facilities are integrated and connected to multiple options for the functional interpretation of the experiments. Different methods of functional enrichment or gene set enrichment can be used to understand the functional basis of the experiment analyzed. Many sources of biological information, which include functional (GO, KEGG, Biocarta, Reactome, etc.), regulatory (Transfac, Jaspar, ORegAnno, miRNAs, etc.), text-mining or protein–protein interaction modules can be used for this purpose. Finally a tool for the de novo functional annotation of sequences has been included in the system. This provides support for the functional analysis of non-model species. Mirrors of Babelomics or command line execution of their individual components are now possible. Babelomics is available at http://www.babelomics.org.

Coordinated modular functionality and prognostic potential of a heart failure biomarker-driven interaction network

Background

The identification of potentially relevant biomarkers and a deeper understanding of molecular mechanisms related to heart failure (HF) development can be enhanced by the implementation of biological network-based analyses. To support these efforts, here we report a global network of protein-protein interactions (PPIs) relevant to HF, which was characterized through integrative bioinformatic analyses of multiple sources of "omic" information.

Results

We found that the structural and functional architecture of this PPI network is highly modular. These network modules can be assigned to specialized processes, specific cellular regions and their functional roles tend to partially overlap. Our results suggest that HF biomarkers may be defined as key coordinators of intra- and inter-module communication. Putative biomarkers can, in general, be distinguished as "information traffic" mediators within this network. The top high traffic proteins are encoded by genes that are not highly differentially expressed across HF and non-HF patients. Nevertheless, we present evidence that the integration of expression patterns from high traffic genes may support accurate prediction of HF. We quantitatively demonstrate that intra- and inter-module functional activity may be controlled by a family of transcription factors known to be associated with the prevention of hypertrophy.

Conclusion

The systems-driven analysis reported here provides the basis for the identification of potentially novel biomarkers and understanding HF-related mechanisms in a more comprehensive and integrated way.

ETS transcription factors control transcription of EZH2 and epigenetic silencing of the tumor suppressor gene Nkx3.1 in prostate cancer

Background

ETS transcription factors regulate important signaling pathways involved in cell differentiation and development in many tissues and have emerged as important players in prostate cancer. However, the biological impact of ETS factors in prostate tumorigenesis is still debated.

Methodology/Principal Findings

We performed an analysis of the ETS gene family using microarray data and real-time PCR in normal and tumor tissues along with functional studies in normal and cancer cell lines to understand the impact in prostate tumorigenesis and identify key targets of these transcription factors. We found frequent dysregulation of ETS genes with oncogenic (i.e., ERG and ESE1) and tumor suppressor (i.e., ESE3) properties in prostate tumors compared to normal prostate. Tumor subgroups (i.e., ERG^high, ESE1^high, ESE3^low and NoETS tumors) were identified on the basis of their ETS expression status and showed distinct transcriptional and biological features. ERG^high and ESE3^low tumors had the most robust gene signatures with both distinct and overlapping features. Integrating genomic data with functional studies in multiple cell lines, we demonstrated that ERG and ESE3 controlled in opposite direction transcription of the Polycomb Group protein EZH2, a key gene in development, differentiation, stem cell biology and tumorigenesis. We further demonstrated that the prostate-specific tumor suppressor gene Nkx3.1 was controlled by ERG and ESE3 both directly and through induction of EZH2.

Conclusions/Significance

These findings provide new insights into the role of the ETS transcriptional network in prostate tumorigenesis and uncover previously unrecognized links between aberrant expression of ETS factors, deregulation of epigenetic effectors and silencing of tumor suppressor genes. The link between aberrant ETS activity and epigenetic gene silencing may be relevant for the clinical management of prostate cancer and design of new therapeutic strategies.

Integrated protein network and microarray analysis to identify potential biomarkers after myocardial infarction

A significant proportion of patients develop left ventricular (LV) dysfunction and heart failure (HF) after acute myocardial infarction (MI). Existing biomarkers of HF provide limited information after MI. To identify new prognostic biomarkers in MI patients, we designed an approach combining protein interaction networks and microarray analysis of blood cells. Blood samples for RNA and protein analysis were taken from 127 acute MI patients. Echocardiography was performed at one month. Assuming that angiogenesis is related to cardiac repair after MI, a protein-protein interaction network of angiogenesis was constructed and analyzed. Among the 556 proteins and 686 interactions of this network, a cluster of 53 proteins highly specialized in regulation of cell growth was identified. Of these 53 proteins, 38 were found differentially expressed by microarrays between low (< or = 40%) and high (>40%) LV ejection fraction (EF) patients (n = 32). Among these 38 genes, prediction analysis identified a set of three genes able to predict significant LV dysfunction (EF < or = 40%) with an area under the receiver operating characteristic curve (AUC) of 0.82. These three genes-vascular endothelial growth factor B, thrombospondin-1 and placental growth factor-had a stronger predictive value than brain natriuretic peptide and troponin T (AUC of 0.63). Independent validations on protein expression and quantitative PCR datasets confirmed the results. In conclusion, a new strategy is described that allows identifying new potential biomarkers. The three specific biomarkers described here remain to be validated in a larger patient population.

Bone progenitor dysfunction induces myelodysplasia and secondary leukaemia

Mesenchymal cells contribute to the 'stroma' of most normal and malignant tissues, with specific mesenchymal cells participating in the regulatory niches of stem cells. By examining how mesenchymal osteolineage cells modulate haematopoiesis, here we show that deletion of Dicer1 specifically in mouse osteoprogenitors, but not in mature osteoblasts, disrupts the integrity of haematopoiesis. Myelodysplasia resulted and acute myelogenous leukaemia emerged that had acquired several genetic abnormalities while having intact Dicer1. Examining gene expression altered in osteoprogenitors as a result of Dicer1 deletion showed reduced expression of Sbds, the gene mutated in Schwachman-Bodian-Diamond syndrome-a human bone marrow failure and leukaemia pre-disposition condition. Deletion of Sbds in mouse osteoprogenitors induced bone marrow dysfunction with myelodysplasia. Therefore, perturbation of specific mesenchymal subsets of stromal cells can disorder differentiation, proliferation and apoptosis of heterologous cells, and disrupt tissue homeostasis. Furthermore, primary stromal dysfunction can result in secondary neoplastic disease, supporting the concept of niche-induced oncogenesis.

Different proliferative potential and migratory characteristics of human CD4+ regulatory T cells that express either CD45RA or CD45RO

Although human naturally occurring regulatory T cells (Tregs) may express either CD45RA or CD45RO, we find in agreement with previous reports that the ( approximately 80%) majority of natural Tregs in adults are CD45RO(+). The proportion of CD45RA(+) Tregs decreases, whereas CD45RO(+) Tregs increase significantly with age. Nevertheless, a small proportion of CD45RA(+) Tregs are found even in old (>80 y) adults and a proportion of these express CD31, a marker for recent thymic emigrants. We found that CD45RO(+) Tregs were highly proliferative compared with their CD45RA(+) counterparts. This was due in part to the conversion of CD45RA Tregs to CD45RO expression after activation. Another difference between these two Treg populations was their preferential migration to different tissues in vivo. Whereas CD45RA(+) Tregs were preferentially located in the bone marrow, associated with increased CXCR4 expression, CD45RO(+) Tregs were preferentially located in the skin, and this was associated with their increased expression of CLA and CCR4. Our studies therefore show that proliferation features strongly in maintenance of the adult Treg pool in humans and that the thymus may make a minor contribution to the maintenance of the peripheral pool of these cells, even in older adults. Furthermore, the different tissue compartmentalization of these cells suggests that different Treg niches exist in vivo, which may have important roles for their maturation and function.

Integrative pathway-centric modeling of ventricular dysfunction after myocardial infarction

Background

A significant proportion of myocardial infarction (MI) patients undergo complex, coordinated perturbations at the molecular level that may eventually drive the occurrence of ventricular dysfunction and heart failure. Despite advances in the elucidation of key processes implicated in this condition, traditional methods relying on gene expression data and the identification of individual biomarkers in isolation pose major limitations not only for improving prediction power, but also for model interpretability. Mechanisms underlying clinical responses after MI remain elusive and there is no biomarker with the capacity to accurately predict ventricular dysfunction after MI. This calls for the exploration of system-level modeling of ventricular dysfunction in post-MI patients. Within this discovery framework key perturbations and predictive patterns are characterized by the integrated biological activity levels observed in pathways, rather than in individual genes.

Methodology/Principal Findings

Here we report an integrative approach to identifying pathways related with ventricular dysfunction post MI with potential prognostic and therapeutic value. We found that a diversity of pathway-level perturbations can be profiled in samples of patients with ventricular dysfunction post MI, most of which represent major reductions of gene expression. Highly perturbed pathways included those implicated in antigen-dependent B-cell activation and the synthesis of leucine. By analyzing patient-specific samples encoded with information derived from highly-perturbed pathways, it is possible to visualize differential prognostic patterns and to perform computational classification of patients with areas under the receiver operating characteristic curve above 0.75. We also demonstrate how the integration of the outcomes generated by different pathway-based analysis models may improve ventricular dysfunction prediction performance.

Significance

This research offers an alternative, comprehensive view of key relationships and perturbations that may trigger the emergence or prevention of ventricular dysfunction post-MI.

Shrinkage-based diagonal discriminant analysis and its applications in high-dimensional data

High-dimensional data such as microarrays have brought us new statistical challenges. For example, using a large number of genes to classify samples based on a small number of microarrays remains a difficult problem. Diagonal discriminant analysis, support vector machines, and k-nearest neighbor have been suggested as among the best methods for small sample size situations, but none was found to be superior to others. In this article, we propose an improved diagonal discriminant approach through shrinkage and regularization of the variances. The performance of our new approach along with the existing methods is studied through simulations and applications to real data. These studies show that the proposed shrinkage-based and regularization diagonal discriminant methods have lower misclassification rates than existing methods in many cases.

Microarray analysis reveals the inhibition of nuclear factor-kappa B signaling by aristolochic acid in normal human kidney (HK-2) cells

AIM

To study the molecular mechanism underlying the effect of aristolochic acid (AA), a major active component of plants from the Aristolochiaceae family using microarray analysis.

METHODS

Human kidney (HK-2) cells were treated with AA (0, 10, 30, and 90 micromol/L) for 24 h, and the cell viability was measured by a 3-(4,5-dimethylthiazol-2yl)-2,5-diphenyltetrazolium bromide assay. Complementary DNA microarrays were used to investigate the gene expression pattern of HK-2 cells exposed to AA in triplicate. A quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) assay was used to verify the microarray data for selected nuclear factor kappa B (NF-kappaB)-regulated genes. Furthermore, the subcellular localization of NF-kappaB p65 was visualized by immunofluorescence confocal microscopy in HK-2 cells. The NF-kappaB activity was examined by a luciferase reporter assay in HK-2/NF-kappaB transgenic cells.

RESULTS

AA exhibited a dose-dependent cytotoxic effect in HK-2 cells and induced alterations in the gene expression profiles related to the DNA damage response, DNA repair, macromolecule metabolic process, carbohydrate metabolic process, DNA metabolic process, apoptosis, cell cycle, and transcription. In addition, 9 biological pathways associated with immunomodulatory functions were down-regulated in AA-treated HK-2 cells. A network analysis revealed that NF-kappaB played a central role in the network topology. Among NF-kappaB-regulated genes, 8 differentially expressed genes were verified by qRT-PCR. The inhibition of NF-kappaB activity by AA was further confirmed by immunofluorescence confocal microscopy and by NF-kappaB luciferase reporter assay.

CONCLUSION

Our data revealed that AA could suppress NF-kappaB activity in normal human cells, perhaps partially accounting for the reported anti-inflammatory effects of some plants from the genus Aristolochia.

Metaanalysis of ChIP-chip data

Genome-wide analysis of histone modifications via ChIP-chip (chromatin immunoprecipitation followed by whole genome tiling array hybridization) may generate lists of up to several thousand potential target genes. In the case of the model organism Arabidopsis thaliana, several databases are available to alleviate further characterization and classification of genomic data sets. The term metaanalysis has been coined for this type of multidatabase comparison. In this chapter, we describe open source software and web tools that perform transcriptional and functional analysis of target genes. Sources of transcription data and clustering tools to subdivide genes according to their expression pattern are described. The user is guided through all necessary steps, including data download and formatting. In addition, the Gene Ontology (GO) vocabulary and methods to uncover over- or underrepresented functions among target genes are introduced. Genomic targets of the histone H3K27me3 modification are presented as a case study to demonstrate that metaanalysis can uncover novel functions that were hidden in genomic data sets.

Changes in the pattern of DNA methylation associate with twin discordance in systemic lupus erythematosus

Monozygotic (MZ) twins are partially concordant for most complex diseases, including autoimmune disorders. Whereas phenotypic concordance can be used to study heritability, discordance suggests the role of non-genetic factors. In autoimmune diseases, environmentally driven epigenetic changes are thought to contribute to their etiology. Here we report the first high-throughput and candidate sequence analyses of DNA methylation to investigate discordance for autoimmune disease in twins. We used a cohort of MZ twins discordant for three diseases whose clinical signs often overlap: systemic lupus erythematosus (SLE), rheumatoid arthritis, and dermatomyositis. Only MZ twins discordant for SLE featured widespread changes in the DNA methylation status of a significant number of genes. Gene ontology analysis revealed enrichment in categories associated with immune function. Individual analysis confirmed the existence of DNA methylation and expression changes in genes relevant to SLE pathogenesis. These changes occurred in parallel with a global decrease in the 5-methylcytosine content that was concomitantly accompanied with changes in DNA methylation and expression levels of ribosomal RNA genes, although no changes in repetitive sequences were found. Our findings not only identify potentially relevant DNA methylation markers for the clinical characterization of SLE patients but also support the notion that epigenetic changes may be critical in the clinical manifestations of autoimmune disease.

PHOENIX, a web interface for (re)analysis of microarray data

Microarrays are tools to study the expression profile of an entire genome. Technology, statistical tools and biological knowledge in general have evolved over the past ten years and it is now possible to improve analysis of previous datasets. We have developed a web interface called PHOENIX that automates the analysis of microarray data from preprocessing to the evaluation of significance through manual or automated parameterization. At each analytical step, several methods are possible for (re)analysis of data. PHOENIX evaluates a consensus score from several methods and thus determines the performance level of the best methods (even if the best performing method is not known). With an estimate of the true gene list, PHOENIX can evaluate the performance of methods or compare the results with other experiments. Each method used for differential expression analysis and performance evaluation has been implemented in the PEGASE back-end package, along with additional tools to further improve PHOENIX. Future developments will involve the addition of steps (CDF selection, geneset analysis, meta-analysis), methods (PLIER, ANOVA, Limma), benchmarks (spike-in and simulated datasets), and illustration of the results (automatically generated report).

Gene expression profiling integrated into network modelling reveals heterogeneity in the mechanisms of BRCA1 tumorigenesis

BACKGROUND

Gene expression profiling has distinguished sporadic breast tumour classes with genetic and clinical differences. Less is known about the molecular classification of familial breast tumours, which are generally considered to be less heterogeneous. Here, we describe molecular signatures that define BRCA1 subclasses depending on the expression of the gene encoding for oestrogen receptor, ESR1.

METHODS

For this purpose, we have used the Oncochip v2, a cancer-related cDNA microarray to analyze 14 BRCA1-associated breast tumours.

RESULTS

Signatures were found to be molecularly associated with different biological processes and transcriptional regulatory programs. The signature of ESR1-positive tumours was mainly linked to cell proliferation and regulated by ER, whereas the signature of ESR1-negative tumours was mainly linked to the immune response and possibly regulated by transcription factors of the REL/NFkappaB family. These signatures were then verified in an independent series of familial and sporadic breast tumours, which revealed a possible prognostic value for each subclass. Over-expression of immune response genes seems to be a common feature of ER-negative sporadic and familial breast cancer and may be associated with good prognosis. Interestingly, the ESR1-negative tumours were substratified into two groups presenting slight differences in the magnitude of the expression of immune response transcripts and REL/NFkappaB transcription factors, which could be dependent on the type of BRCA1 germline mutation.

CONCLUSION

This study reveals the molecular complexity of BRCA1 breast tumours, which are found to display similarities to sporadic tumours, and suggests possible prognostic implications.