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.

Myb-binding protein 1a augments AhR-dependent gene expression

We have studied the mechanism by which an acidic domain (amino acids 515-583) of the aromatic hydrocarbon receptor (AhR) transactivates a target gene. Studies with glutathione S-transferase fusion proteins demonstrate that the wild-type acidic domain associates in vitro with Myb-binding protein la, whereas a mutant domain (F542A, 1569A) does not. AhR-defective cells reconstituted with an AhR containing the wild-type acidic domain exhibit normal AhR function; however, cells reconstituted with an AhR containing the mutant acidic domain do not function normally. Transient transfection of Myb-binding protein la into mouse hepatoma cells is associated with augmentation of AhR-dependent gene expression. Such augmentation does not occur when Myb-binding protein la is transfected into AhR-defective cells that have been reconstituted with an AhR that lacks the acidic domain. We infer that 1) Myb-binding protein la associates with AhR, thereby enhancing transactivation, and 2) the presence of AhR's acidic domain is both necessary and sufficient for Myb-binding protein la to increase AhR-dependent gene expression.

ASPM gene expression in medulloblastoma

Medulloblastomas are the most common malignant tumors of the central nervous system in childhood. The incidence is about 19-20% between children younger than 16 years old with peak incidence between 4 and 7 years. Despite its sensibility to no specific therapeutic means like chemotherapy and radiotherapy, the treatment is very aggressive and frequently results in regression, growth deficit, and endocrine dysfunction. From this point of view, new treatment approaches are needed such as molecular targeted therapies. Studies in glioblastoma demonstrated that ASPM gene was overexpressed when compared to normal brain and ASPM inhibition by siRNA-mediated inhibits tumor cell proliferation and neural stem cell proliferation, supporting ASPM gene as a potential molecular target in glioblastoma. The aim of this work was to evaluate ASPM expression in medulloblastoma fragment samples, and to compare the results with the patient clinical features. Analysis of gene expression was performed by quantitative PCR real time using SYBR Green system in tumor samples from 37 children. The t test was used to analyze the gene expression, and Mann-Whitney test was performed to analyze the relationship between gene expressions and clinical characteristics. Kaplan-Meier test evaluated curve survival. All samples overexpressed ASPM gene more than 40-fold. However, we did not find any association between the overexpressed samples and the clinical parameters. ASPM overexpression may modify the ability of stem cells to differentiate during the development of the central nervous system, contributing to the development of medulloblastoma, a tumor of embryonic origin from cerebellar progenitor cells.

Differential gene expression of peripheral blood mononuclear cells from rheumatoid arthritis patients may discriminate immunogenetic, pathogenic and treatment features

This study aimed to evaluate the association between the differential gene expression profiling of peripheral blood mononuclear cells of rheumatoid arthritis patients with their immunogenetic (human leucocyte antigen shared-epitope, HLA-SE), autoimmune response [anti-cyclic citrullinated peptide (CCP) antibodies], disease activity score (DAS-28) and treatment (disease-modifying antirheumatic drugs and tumour necrosis factor blocker) features. Total RNA samples were copied into Cy3-labelled complementary DNA probes, hybridized onto a glass slide microarray containing 4500 human IMAGE complementary DNA target sequences. The Cy3-monocolour microarray images from patients were quantified and normalized. Analysis of the data using the significance analysis of microarrays algorithm together with a Venn diagram allowed the identification of shared and of exclusively modulated genes, according to patient features. Thirteen genes were exclusively associated with the presence of HLA-SE alleles, whose major biological function was related to signal transduction, phosphorylation and apoptosis. Ninety-one genes were associated with disease activity, being involved in signal transduction, apoptosis, response to stress and DNA damage. One hundred and one genes were associated with the presence of anti-CCP antibodies, being involved in signal transduction, cell proliferation and apoptosis. Twenty-eight genes were associated with tumour necrosis factor blocker treatment, being involved in intracellular signalling cascade, phosphorylation and protein transport. Some of these genes had been previously associated with rheumatoid arthritis pathogenesis, whereas others were unveiled for future research.

Shared and Unique Gene Expression in Systemic Lupus Erythematosus Depending on Disease Activity

Patients presenting with active Systemic lupus erythematosus (SLE) manifestations may exhibit distinct pathogenetic features in relation to inactive SLE. Also, cDNA microarrays may potentially discriminate the gene expression profile of a disease or disease variant. Therefore, we evaluated the expression profile of 4500 genes in peripheral blood lymphocytes (PBL) of SLE patients. We studied 11 patients with SLE (seven with active SLE and four with inactive SLE) and eight healthy controls. Total RNA was isolated from PBL, reverse transcribed into cDNA, and postlabeled with Cy3 fluorochrome. These probes were then hybridized to a glass slide cDNA microarray containing 4500 human IMAGE cDNA target sequences. An equimolar amount of total RNA from human cell lines served as reference. The microarray images were quantified, normalized, and analyzed using the R environment (ANOVA, significant analysis of microarrays, and cluster-tree view algorithms). Disease activity was assessed by the SLE disease activity index. Compared to the healthy controls, 104 genes in active SLE patients (80 repressed and 24 induced) and 52 genes in nonactive SLE patients (31 induced and 21 repressed) were differentially expressed. The modulation of 12 genes, either induced or repressed, was found in both disease variants; however, each disease variant had differential expression of different genes. Taken together, these results indicate that the two lupus variants studied have common and unique differentially expressed genes. Although the biological significance of the differentially expressed genes discussed above has not been completely understood, they may serve as a platform to further explore the molecular basis of immune deregulation in SLE.

Genetic and physical interactions between Microphthalmia transcription factor and PU.1 are necessary for osteoclast gene expression and differentiation

The microphthalmia transcription factor (MITF), a basic-helix-loop-helix zipper factor, regulates distinct target genes in several cell types. We hypothesized that interaction with the Ets family factor PU.1, whose expression is limited to hematopoietic cells, might be necessary for activation of target genes like tartrate-resistant acid phosphatase (TRAP) in osteoclasts. Several lines of evidence were consistent with this model. The combination of MITF and PU.1 synergistically activated the TRAP promoter in transient assays. This activation was dependent on intact binding sites for both factors in the TRAP promoter. MITF and PU.1 physically interacted when coexpressed in COS cells or in vitro when purified recombinant proteins were studied. The minimal regions of MITF and PU.1 required for the interaction were the basic-helix-loop-helix zipper domain and the Ets DNA binding domain, respectively. Significantly, mice heterozygous for both the mutant mi allele and a PU.1 null allele developed osteopetrosis early in life which resolved with age. The size and number of osteoclasts were not altered in the double heterozygous mutant mice, indicating that the defect lies in mature osteoclast function. Taken in total, the results afford an example of how lineage-specific gene regulation can be achieved by the combinatorial action of two broadly expressed transcription factors.

Fsp27/CIDEC is a CREB target gene induced during fasting and regulated by fatty acid oxidation rate

Dissertação para obtenção do Grau de Mestre em Biotecnologia

Controlling gene expression in enterobacteriaceae: studies on sRNAs and strategies for synthetic biology

Dissertation presented to obtain the Master Degree in Molecular, Genetics and Biomedicine

Reconstruction of the regulatory network for Bacillus subtilis and reconciliation with gene expression data

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb. 2016.00275

Endothelial Cell-Derived Nitric Oxide Enhances Aerobic Glycolysis in Astrocytes via HIF-1α-Mediated Target Gene Activation.

Astrocytes exhibit a prominent glycolytic activity, but whether such a metabolic profile is influenced by intercellular communication is unknown. Treatment of primary cultures of mouse cortical astrocytes with the nitric oxide (NO) donor DetaNONOate induced a time-dependent enhancement in the expression of genes encoding various glycolytic enzymes as well as transporters for glucose and lactate. Such an effect was shown to be dependent on the hypoxia-inducible factor HIF-1α, which is stabilized and translocated to the nucleus to exert its transcriptional regulation. NO action was dependent on both the PI3K/Akt/mTOR and MEK signaling pathways and required the activation of COX, but was independent of the soluble guanylate cyclase pathway. Furthermore, as a consequence of NO treatment, an enhanced lactate production and release by astrocytes was evidenced, which was prevented by downregulating HIF-1α. Several brain cell types represent possible sources of NO. It was found that endothelial cells, which express the endothelial NO synthase (eNOS) isoform, constitutively produced the largest amount of NO in culture. When astrocytes were cocultured with primary cultures of brain vascular endothelial cells, stabilization of HIF-1α and an enhancement in glucose transporter-1, hexokinase-2, and monocarboxylate transporter-4 expression as well as increased lactate production was found in astrocytes. This effect was inhibited by the NOS inhibitor l-NAME and was not seen when astrocytes were cocultured with primary cultures of cortical neurons. Our findings suggest that endothelial cell-derived NO participates to the maintenance of a high glycolytic activity in astrocytes mediated by astrocytic HIF-1α activation.

The Interleukin-1 receptor antagonist is a direct target gene of PPARalpha in liver.

BACKGROUND/AIMS: The Peroxisome Proliferator-Activated Receptor (PPAR) alpha belongs to the superfamily of Nuclear Receptors and plays an important role in numerous cellular processes, including lipid metabolism. It is known that PPARalpha also has an anti-inflammatory effect, which is mainly achieved by down-regulating pro-inflammatory genes. The objective of this study was to further characterize the role of PPARalpha in inflammatory gene regulation in liver. RESULTS: According to Affymetrix micro-array analysis, the expression of various inflammatory genes in liver was decreased by treatment of mice with the synthetic PPARalpha agonist Wy14643 in a PPARalpha-dependent manner. In contrast, expression of Interleukin-1 receptor antagonist (IL-1ra), which was acutely stimulated by LPS treatment, was induced by PPARalpha. Up-regulation of IL-1ra by LPS was lower in PPARalpha -/- mice compared to Wt mice. Transactivation and chromatin immunoprecipitation studies identified IL-1ra as a direct positive target gene of PPARalpha with a functional PPRE present in the promoter. Up-regulation of IL-1ra by PPARalpha was conserved in human HepG2 hepatoma cells and the human monocyte/macrophage THP-1 cell line. CONCLUSIONS: In addition to down-regulating expression of pro-inflammatory genes, PPARalpha suppresses the inflammatory response by direct up-regulation of genes with anti-inflammatory properties.

Transcriptional profiling reveals divergent roles of PPARalpha and PPARbeta/delta in regulation of gene expression in mouse liver.

Little is known about the role of the transcription factor peroxisome proliferator-activated receptor (PPAR) beta/delta in liver. Here we set out to better elucidate the function of PPARbeta/delta in liver by comparing the effect of PPARalpha and PPARbeta/delta deletion using whole genome transcriptional profiling and analysis of plasma and liver metabolites. In fed state, the number of genes altered by PPARalpha and PPARbeta/delta deletion was similar, whereas in fasted state the effect of PPARalpha deletion was much more pronounced, consistent with the pattern of gene expression of PPARalpha and PPARbeta/delta. Minor overlap was found between PPARalpha- and PPARbeta/delta-dependent gene regulation in liver. Pathways upregulated by PPARbeta/delta deletion were connected to innate immunity and inflammation. Pathways downregulated by PPARbeta/delta deletion included lipoprotein metabolism and various pathways related to glucose utilization, which correlated with elevated plasma glucose and triglycerides and reduced plasma cholesterol in PPARbeta/delta-/- mice. Downregulated genes that may underlie these metabolic alterations included Pklr, Fbp1, Apoa4, Vldlr, Lipg, and Pcsk9, which may represent novel PPARbeta/delta target genes. In contrast to PPARalpha-/- mice, no changes in plasma free fatty acid, plasma beta-hydroxybutyrate, liver triglycerides, and liver glycogen were observed in PPARbeta/delta-/- mice. Our data indicate that PPARbeta/delta governs glucose utilization and lipoprotein metabolism and has an important anti-inflammatory role in liver. Overall, our analysis reveals divergent roles of PPARalpha and PPARbeta/delta in regulation of gene expression in mouse liver.

Hitchhiking Trypanosoma cruzi minicircle DNA affects gene expression in human host cells via LINE-1 retrotransposon

The horizontal transfer of Trypanosoma cruzi mitochondrial minicircle DNA to the genomes of naturally infected humans may play an important role in the pathogenesis of Chagas disease. Minicircle integrations within LINE-1 elements create the potential for foreign DNA mobility within the host genome via the machinery associated with this retrotransposon. Here we document integration of minicircle DNA fragments in clonal human macrophage cell lines and their mobilization over time. The movement of an integration event in a clonal transfected cell line was tracked at three months and three years post-infection. The minicircle sequence integrated into a LINE-1 retrotransposon; one such foreign fragment subsequently relocated to another genomic location in association with associated LINE-1 elements. The p15 locus was altered at three years as a direct effect of minicircle/LINE-1 acquisition, resulting in elimination of p15 mRNA. Here we show for the first time a molecular pathology stemming from mobilization of a kDNA/LINE-1 mutation. These genomic changes and detected transcript variations are consistent with our hypothesis that minicircle integration is a causal component of parasite-independent, autoimmune-driven lesions seen in the heart and other target tissues associated with Chagas disease.

Steroid hormone pulsing drives cyclic gene expression.

Transcriptional cycling of activated glucocorticoid receptor (GR) and ultradian glucocorticoid secretion are well established processes. Ultradian hormone release is now shown to result in pulsatile gene transcription through dynamic exchange of GR with the target-gene promoter and GR cycling through the chaperone machinery.

Zinc-alpha 2-glycoprotein gene expression in adipose tissue is related with insulin resistance and lipolytic genes in morbidly obese patients.

OBJECTIVE Zinc-α(2) glycoprotein (ZAG) stimulates lipid loss by adipocytes and may be involved in the regulation of adipose tissue metabolism. However, to date no studies have been made in the most extreme of obesity. The aims of this study are to analyze ZAG expression levels in adipose tissue from morbidly obese patients, and their relationship with lipogenic and lipolytic genes and with insulin resistance (IR). METHODS mRNA expression levels of PPARγ, IRS-1, IRS-2, lipogenic and lipolytic genes and ZAG were quantified in visceral (VAT) and subcutaneous adipose tissue (SAT) of 25 nondiabetic morbidly obese patients, 11 with low IR and 14 with high IR. Plasma ZAG was also analyzed. RESULTS The morbidly obese patients with low IR had a higher VAT ZAG expression as compared with the patients with high IR (p = 0.023). In the patients with low IR, the VAT ZAG expression was greater than that in SAT (p = 0.009). ZAG expression correlated between SAT and VAT (r = 0.709, p<0.001). VAT ZAG expression was mainly predicted by insulin, HOMA-IR, plasma adiponectin and expression of adiponectin and ACSS2. SAT ZAG expression was only predicted by expression of ATGL. CONCLUSIONS ZAG could be involved in modulating lipid metabolism in adipose tissue and is associated with insulin resistance. These findings suggest that ZAG may be a useful target in obesity and related disorders, such as diabetes.