Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

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.

Characterization of global transcription profile of normal and HPV-immortalized keratinocytes and their response to TNF treatment

Background: Persistent infection by high risk HPV types (e.g. HPV-16, -18, -31, and -45) is the main risk factor for development of cervical intraepithelial neoplasia and cervical cancer. Tumor necrosis factor (TNF) is a key mediator of epithelial cell inflammatory response and exerts a potent cytostatic effect on normal or HPV16, but not on HPV18 immortalized keratinocytes. Moreover, several cervical carcinoma-derived cell lines are resistant to TNF anti-proliferative effect suggesting that the acquisition of TNF-resistance may constitute an important step in HPV-mediated carcinogenesis. In the present study, we compared the gene expression profiles of normal and HPV16 or 18 immortalized human keratinocytes before and after treatment with TNF for 3 or 60 hours. Methods: In this study, we determined the transcriptional changes 3 and 60 hours after TNF treatment of normal, HPV16 and HPV18 immortalized keratinocytes by microarray analysis. The expression pattern of two genes observed by microarray was confirmed by Northern Blot. NF-kappa B activation was also determined by electrophoretic mobility shift assay (EMSA) using specific oligonucleotides and nuclear protein extracts. Results: We observed the differential expression of a common set of genes in two TNF-sensitive cell lines that differs from those modulated in TNF-resistant ones. This information was used to define genes whose differential expression could be associated with the differential response to TNF, such as: KLK7 (kallikrein 7), SOD2 (superoxide dismutase 2), 100P (S100 calcium binding protein P), PI3 (protease inhibitor 3, skin-derived), CSTA (cystatin A), RARRES1 (retinoic acid receptor responder 1), and LXN (latexin). The differential expression of the KLK7 and SOD2 transcripts was confirmed by Northern blot. Moreover, we observed that SOD2 expression correlates with the differential NF-kappa B activation exhibited by TNF-sensitive and TNF-resistant cells. Conclusion: This is the first in depth analysis of the differential effect of TNF on normal and HPV16 or HPV18 immortalized keratinocytes. Our findings may be useful for the identification of genes involved in TNF resistance acquisition and candidate genes which deregulated expression may be associated with cervical disease establishment and/or progression.

Construct and Compare Gene Coexpression Networks with DAPfinder and DAPview

Background: DAPfinder and DAPview are novel BRB-ArrayTools plug-ins to construct gene coexpression networks and identify significant differences in pairwise gene-gene coexpression between two phenotypes. Results: Each significant difference in gene-gene association represents a Differentially Associated Pair (DAP). Our tools include several choices of filtering methods, gene-gene association metrics, statistical testing methods and multiple comparison adjustments. Network results are easily displayed in Cytoscape. Analyses of glioma experiments and microarray simulations demonstrate the utility of these tools. Conclusions: DAPfinder is a new friendly-user tool for reconstruction and comparison of biological networks.

The impact of measurement errors in the identification of regulatory networks

Background: There are several studies in the literature depicting measurement error in gene expression data and also, several others about regulatory network models. However, only a little fraction describes a combination of measurement error in mathematical regulatory networks and shows how to identify these networks under different rates of noise. Results: This article investigates the effects of measurement error on the estimation of the parameters in regulatory networks. Simulation studies indicate that, in both time series (dependent) and non-time series (independent) data, the measurement error strongly affects the estimated parameters of the regulatory network models, biasing them as predicted by the theory. Moreover, when testing the parameters of the regulatory network models, p-values computed by ignoring the measurement error are not reliable, since the rate of false positives are not controlled under the null hypothesis. In order to overcome these problems, we present an improved version of the Ordinary Least Square estimator in independent (regression models) and dependent (autoregressive models) data when the variables are subject to noises. Moreover, measurement error estimation procedures for microarrays are also described. Simulation results also show that both corrected methods perform better than the standard ones (i.e., ignoring measurement error). The proposed methodologies are illustrated using microarray data from lung cancer patients and mouse liver time series data. Conclusions: Measurement error dangerously affects the identification of regulatory network models, thus, they must be reduced or taken into account in order to avoid erroneous conclusions. This could be one of the reasons for high biological false positive rates identified in actual regulatory network models.

KeaA, a Dictyostelium kelch-domain protein that regulates the response to stress and development

Background: The protein kinase YakA is responsible for the growth arrest and induction of developmental processes that occur upon starvation of Dictyostelium cells. yakA-cells are aggregation deficient, have a faster cell cycle and are hypersensitive to oxidative and nitrosoative stress. With the aim of isolating members of the YakA pathway, suppressors of the death induced by nitrosoative stress in the yakA-cells were identified. One of the suppressor mutations occurred in keaA, a gene identical to DG1106 and similar to Keap1 from mice and the Kelch protein from Drosophila, among others that contain Kelch domains. Results: A mutation in keaA suppresses the hypersensitivity to oxidative and nitrosoative stresses but not the faster growth phenotype of yakA-cells. The growth profile of keaA deficient cells indicates that this gene is necessary for growth. keaA deficient cells are more resistant to nitrosoative and oxidative stress and keaA is necessary for the production and detection of cAMP. A morphological analysis of keaA deficient cells during multicellular development indicated that, although the mutant is not absolutely deficient in aggregation, cells do not efficiently participate in the process. Gene expression analysis using cDNA microarrays of wild-type and keaA deficient cells indicated a role for KeaA in the regulation of the cell cycle and pre-starvation responses. Conclusions: KeaA is required for cAMP signaling following stress. Our studies indicate a role for kelch proteins in the signaling that regulates the cell cycle and development in response to changes in the environmental conditions.

Identification of protein-coding and non-coding RNA expression profiles in CD34(+) and in stromal cells in refractory anemia with ringed sideroblasts

Background: Myelodysplastic syndromes (MDS) are a group of clonal hematological disorders characterized by ineffective hematopoiesis with morphological evidence of marrow cell dysplasia resulting in peripheral blood cytopenia. Microarray technology has permitted a refined high-throughput mapping of the transcriptional activity in the human genome. Non-coding RNAs (ncRNAs) transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression, and in the regulation of exon-skipping and intron retention. Characterization of ncRNAs in progenitor cells and stromal cells of MDS patients could be strategic for understanding gene expression regulation in this disease. Methods: In this study, gene expression profiles of CD34(+) cells of 4 patients with MDS of refractory anemia with ringed sideroblasts (RARS) subgroup and stromal cells of 3 patients with MDS-RARS were compared with healthy individuals using 44 k combined intron-exon oligoarrays, which included probes for exons of protein-coding genes, and for non-coding RNAs transcribed from intronic regions in either the sense or antisense strands. Real-time RT-PCR was performed to confirm the expression levels of selected transcripts. Results: In CD34(+) cells of MDS-RARS patients, 216 genes were significantly differentially expressed (q-value <= 0.01) in comparison to healthy individuals, of which 65 (30%) were non-coding transcripts. In stromal cells of MDS-RARS, 12 genes were significantly differentially expressed (q-value <= 0.05) in comparison to healthy individuals, of which 3 (25%) were non-coding transcripts. Conclusions: These results demonstrated, for the first time, the differential ncRNA expression profile between MDS-RARS and healthy individuals, in CD34(+) cells and stromal cells, suggesting that ncRNAs may play an important role during the development of myelodysplastic syndromes.

Environmental stresses inhibit splicing in the aquatic fungus Blastocladiella emersonii

Background: Exposure of cells to environmental stress conditions can lead to the interruption of several intracellular processes, in particular those performed by macromolecular complexes such as the spliceosome. Results: During nucleotide sequencing of cDNA libraries constructed using RNA isolated from B. emersonii cells submitted to heat shock and cadmium stress, a large number of ESTs with retained introns was observed. Among the 6,350 ESTs obtained through sequencing of stress cDNA libraries, 181 ESTs presented putative introns (2.9%), while sequencing of cDNA libraries from unstressed B. emersonii cells revealed only 0.2% of ESTs containing introns. These data indicate an enrichment of ESTs with introns in B. emersonii stress cDNA libraries. Among the 85 genes corresponding to the ESTs that retained introns, 19 showed more than one intron and three showed three introns, with intron length ranging from 55 to 333 nucleotides. Canonical splicing junctions were observed in most of these introns, junction sequences being very similar to those found in introns from genes previously characterized in B. emersonii, suggesting that inhibition of splicing during stress is apparently a random process. Confirming our observations, analyses of gpx3 and hsp70 mRNAs by Northern blot and S1 protection assays revealed a strong inhibition of intron splicing in cells submitted to cadmium stress. Conclusion: In conclusion, data indicate that environmental stresses, particularly cadmium treatment, inhibit intron processing in B. emersonii, revealing a new adaptive response to cellular exposure to this heavy metal.

Interference with Hemozoin Formation Represents an Important Mechanism of Schistosomicidal Action of Antimalarial Quinoline Methanols

Background: The parasitic trematode Schistosoma mansoni is one of the major causative agents of human schistosomiasis, which afflicts 200 million people worldwide. Praziquantel remains the main drug used for schistosomiasis treatment, and reliance on the single therapy has been prompting the search for new therapeutic compounds against this disease. Our group has demonstrated that heme crystallization into hemozoin (Hz) within the S. mansoni gut is a major heme detoxification route with lipid droplets involved in this process and acting as a potential chemotherapeutical target. In the present work, we investigated the effects of three antimalarial compounds, quinine (QN), quinidine (QND) and quinacrine (QCR) in a murine schistosomiasis model by using a combination of biochemical, cell biology and molecular biology approaches. Methodology/Principal Findings: Treatment of S. mansoni-infected female Swiss mice with daily intraperitoneal injections of QN, and QND (75 mg/kg/day) from the 11(th) to 17(th) day after infection caused significant decreases in worm burden (39%-61%) and egg production (42%-98%). Hz formation was significantly inhibited (40%-65%) in female worms recovered from QN- and QND-treated mice and correlated with reduction in the female worm burden. We also observed that QN treatment promoted remarkable ultrastructural changes in male and female worms, particularly in the gut epithelium and reduced the granulomatous reaction to parasite eggs trapped in the liver. Microarray gene expression analysis indicated that QN treatment increased the expression of transcripts related to musculature, protein synthesis and repair mechanisms. Conclusions: The overall significant reduction in several disease burden parameters by the antimalarial quinoline methanols indicates that interference with Hz formation in S. mansoni represents an important mechanism of schistosomicidal action of these compounds and points out the heme crystallization process as a valid chemotherapeutic target to treat schistosomiasis.

LATENT RANK CHANGE DETECTION FOR ANALYSIS OF SPLICE-JUNCTION MICROARRAYS WITH NONLINEAR EFFECTS

Alternative splicing of gene transcripts greatly expands the functional capacity of the genome, and certain splice isoforms may indicate specific disease states such as cancer. Splice junction microarrays interrogate thousands of splice junctions, but data analysis is difficult and error prone because of the increased complexity compared to differential gene expression analysis. We present Rank Change Detection (RCD) as a method to identify differential splicing events based upon a straightforward probabilistic model comparing the over-or underrepresentation of two or more competing isoforms. RCD has advantages over commonly used methods because it is robust to false positive errors due to nonlinear trends in microarray measurements. Further, RCD does not depend on prior knowledge of splice isoforms, yet it takes advantage of the inherent structure of mutually exclusive junctions, and it is conceptually generalizable to other types of splicing arrays or RNA-Seq. RCD specifically identifies the biologically important cases when a splice junction becomes more or less prevalent compared to other mutually exclusive junctions. The example data is from different cell lines of glioblastoma tumors assayed with Agilent microarrays.

New genes of Xanthomonas citri subsp citri involved in pathogenesis and adaptation revealed by a transposon-based mutant library

Background: Citrus canker is a disease caused by the phytopathogens Xanthomonas citri subsp. citri, Xanthomonas fuscans subsp. aurantifolli and Xanthomonas alfalfae subsp. citrumelonis. The first of the three species, which causes citrus bacterial canker type A, is the most widely spread and severe, attacking all citrus species. In Brazil, this species is the most important, being found in practically all areas where citrus canker has been detected. Like most phytobacterioses, there is no efficient way to control citrus canker. Considering the importance of the disease worldwide, investigation is needed to accurately detect which genes are related to the pathogen-host adaptation process and which are associated with pathogenesis. Results: Through transposon insertion mutagenesis, 10,000 mutants of Xanthomonas citri subsp. citri strain 306 (Xcc) were obtained, and 3,300 were inoculated in Rangpur lime (Citrus limonia) leaves. Their ability to cause citrus canker was analyzed every 3 days until 21 days after inoculation; a set of 44 mutants showed altered virulence, with 8 presenting a complete loss of causing citrus canker symptoms. Sequencing of the insertion site in all 44 mutants revealed that 35 different ORFs were hit, since some ORFs were hit in more than one mutant, with mutants for the same ORF presenting the same phenotype. An analysis of these ORFs showed that some encoded genes were previously known as related to pathogenicity in phytobacteria and, more interestingly, revealed new genes never implicated with Xanthomonas pathogenicity before, including hypothetical ORFs. Among the 8 mutants with no canker symptoms are the hrpB4 and hrpX genes, two genes that belong to type III secretion system (TTSS), two hypothetical ORFS and, surprisingly, the htrA gene, a gene reported as involved with the virulence process in animal-pathogenic bacteria but not described as involved in phytobacteria virulence. Nucleic acid hybridization using labeled cDNA probes showed that some of the mutated genes are differentially expressed when the bacterium is grown in citrus leaves. Finally, comparative genomic analysis revealed that 5 mutated ORFs are in new putative pathogenicity islands. Conclusion: The identification of these new genes related with Xcc infection and virulence is a great step towards the understanding of plant-pathogen interactions and could allow the development of strategies to control citrus canker.

Searching for Moniliophthora perniciosa pathogenicity genes

The basidiomycete Moniliophthora perniciosa is the causal agent of witches` broom disease of Theobroma cacao (cacao). Pathogenesis mechanisms of this hemibiotrophic fungus are largely unknown. An approach to identify putative pathogenicity genes is searching for sequences induced in mycelia grown under in vitro conditions. Using this approach, genes from M. perniciosa induced under limiting nitrogen and light were identified from a cDNA library enriched by suppression subtractive hybridization as potential putative pathogenicity genes. From the 159 identified unique sequences, 59 were annotated and classified by gene ontology. Two sequences were categorized as ""Defence genes, Virulence, and Cell response"" presumably coding for allergenic proteins, whose homologues from other fungi are inducers of animal or plant defences. Differential gene expression was evaluated by quantitative amplification of reversed transcripts (RT-qPCR) of the putative identified genes coding for the two allergenic proteins (Aspf13 and 88KD), and for the enzymes Arylsulfatase (AS); Aryl-Alcohol Oxidase; Aldo-Keto Reductase (AK); Cytochrome P450 (P450); Phenylalanine Ammonia-Lyase; and Peroxidase from mycelia grown under contrasting N concentrations. All genes were validated for differential expression, except for the putative Peroxidase. The same eight genes were analysed for expression in susceptible plants inoculated with M. perniciosa, and six were induced during the early asymptomatic stage of the disease. In infected host tissues, transcripts of 88KD and AS were found more abundant at the biotrophic phase, while those from Aspf13, AK, PAL, and P450 accumulated at the necrotrophic phase, enabling to suggest that mycelia transition from biotrophic to necrotrophic might occur earlier than currently considered. These sequences appeared to be virulence life-style genes, which encode factors or enzymes that enable invasion, colonization or intracellular survival, or manipulate host factors to benefit the pathogen`s own survival in the hostile environment. (C) 2010 The British Mycological Society. Published by Elsevier Ltd. All rights reserved.

MicroRNAs 29 are involved in the improvement of ventricular compliance promoted by aerobic exercise training in rats

Soci UPR, Fernandes T, Hashimoto NY, Mota GF, Amadeu MA, Rosa KT, Irigoyen MC, Phillips MI, Oliveira EM. MicroRNAs 29 are involved in the improvement of ventricular compliance promoted by aerobic exercise training in rats. Physiol Genomics 43: 665-673, 2011. First published March 29, 2011; doi:10.1152/physiolgenomics.00145.2010.-MiRNAs regulate cardiac development, hypertrophy, and angiogenesis, but their role in cardiac hypertrophy (CH) induced by aerobic training has not previously been studied. Aerobic training promotes physiological CH preserving cardiac function. This study assessed involvement of miRNAs-29 in CH of trained rats. Female Wistar rats (n = 7/group) were randomized into three groups: sedentary (S), training 1 (T1), training 2 (T2). T1: swimming sessions of 60 min/5 days/wk/10 wk. T2: similar to T1 until 8th wk. On the 9th wk rats swam 2x/day, and on the 10th wk 3x/day. MiRNAs analysis was performed by miRNA microarray and confirmed by real-time PCR. We assessed: markers of training, CH by ratio of left ventricle (LV) weight/body wt and cardiomyocytes diameter, pathological markers of CH (ANF, skeletal alpha-actin, alpha/beta-MHC), collagen I and III (COLIAI and COLIIIAI) by real-time PCR, protein collagen by hydroxyproline (OH-proline) concentration, CF and CH by echocardiography. Training improved aerobic capacity and induced CH. MiRNAs-1, 133a, and 133b were downregulated as observed in pathological CH, however, without pathological markers. MiRNA-29c expression increased in T1 (52%) and T2 (123%), correlated with a decrease in COLIAI and COLIIIAI expression in T1 (27%, 38%) and T2 (33%, 48%), respectively. MiRNA-29c was inversely correlated to OH-proline concentration (r = 0.61, P = 0.05). The E/A ratio increased in T2, indicating improved LV compliance. Thus, these results show that aerobic training increase miR-29 expression and decreased collagen gene expression and concentration in the heart, which is relevant to the improved LV compliance and beneficial cardiac effects, associated with aerobic high performance training.

Recombinant rabies virus glycoprotein synthesis in bioreactor by transfected Drosophila melanogaster S2 cells carrying a constitutive or an inducible promoter

S2 cell populations (S2AcRVGP2K and S2MtRVGP-Hy) were selected after transfection of gene expression vectors carrying the cDNA encoding the rabies virus glycoprotein (RVGP) gene under the control of the constitutive (actin) or inductive (metallothionein) promoters. These cell populations were cultivated in a 1 L bioreactor mimicking a large scale bioprocess. Cell cultures were carried out at 90 rpm and monitored/controlled for temperature (28 degrees C) and dissolved oxygen (10 or 50% air saturation). Cell growth attained similar to 1.5-3 x 10(7) cells/mL after 3-4 clays of cultivation. The constitutive synthesis of RVGP in S2AcRVGP2K cells led to values of 0.76 mu g/10(7) cells at day 4 of culture. The RVGP synthesis in S2MtRVGP-Hy cell fraction increased upon CuSO(4) induction attaining specific productivities of 1.5-2 mu g/10(7) cells at clays 4-5. RVGP values in supernatant as a result of cell lysis were always very low (<0.2 mu g/mL) indicating good integrity of cells in culture. Overall the RVGP productivity was of 1.5-3 mg/L. Our data showed an important influence of dissolved oxygen on RVGP synthesis allowing a higher and sustained productivity by S2MtRVGP-Hy cells when cultivated with a DO of 10% air saturation. The RVGP productivity in bioreactors shown here mirrors those previously observed for T-flasks and shaker bottles and allow the preparation of the large RVGP quantities required for studies of structure and function. (C) 2010 Elsevier B.V. All rights reserved.

Time course proteomic profiling of human myocardial infarction plasma samples: An approach to new biomarker discovery

Background: The aim of this study was to identify novel candidate biomarker proteins differentially expressed in the plasma of patients with early stage acute myocardial infarction (AMI) using SELDI-TOF-MS as a high throughput screening technology. Methods: Ten individuals with recent acute ischemic-type chest pain (< 12 h duration) and ST-segment elevation AMI (1STEMI) and after a second AMI (2STEMI) were selected. Blood samples were drawn at six times after STEMI diagnosis. The first stage (T(0)) was in Emergency Unit before receiving any medication, the second was just after primary angioplasty (T(2)), and the next four stages occurred at 12 h intervals after T(0). Individuals (n = 7) with similar risk factors for cardiovascular disease and normal ergometric test were selected as a control group (CG). Plasma proteomic profiling analysis was performed using the top-down (i.e. intact proteins) SELDI-TOF-MS, after processing in a Multiple Affinity Removal Spin Cartridge System (Agilent). Results: Compared with the CG, the 1STEMI group exhibited 510 differentially expressed protein peaks in the first 48 h after the AMI (p < 0.05). The 2STEMI group, had similar to 85% fewer differently expressed protein peaks than those without previous history of AMI (76, p < 0.05). Among the 16 differentially-regulated protein peaks common to both STEMI cohorts (compared with the CG at T(0)), 6 peaks were persistently down-regulated at more than one time-stage, and also were inversed correlated with serum protein markers (cTnI, CK and CKMB) during 48 h-period after IAM. Conclusions: Proteomic analysis by SELDI-TOF-MS technology combined with bioinformatics tools demonstrated differential expression during a 48 h time course suggests a potential role of some of these proteins as biomarkers for the very early stages of AMI, as well as for monitoring early cardiac ischemic recovery. (C) 2011 Elsevier B.V. All rights reserved.