Modeling gene expression regulatory networks with the sparse vector autoregressive model

Abstract Background To understand the molecular mechanisms underlying important biological processes, a detailed description of the gene products networks involved is required. In order to define and understand such molecular networks, some statistical methods are proposed in the literature to estimate gene regulatory networks from time-series microarray data. However, several problems still need to be overcome. Firstly, information flow need to be inferred, in addition to the correlation between genes. Secondly, we usually try to identify large networks from a large number of genes (parameters) originating from a smaller number of microarray experiments (samples). Due to this situation, which is rather frequent in Bioinformatics, it is difficult to perform statistical tests using methods that model large gene-gene networks. In addition, most of the models are based on dimension reduction using clustering techniques, therefore, the resulting network is not a gene-gene network but a module-module network. Here, we present the Sparse Vector Autoregressive model as a solution to these problems. Results We have applied the Sparse Vector Autoregressive model to estimate gene regulatory networks based on gene expression profiles obtained from time-series microarray experiments. Through extensive simulations, by applying the SVAR method to artificial regulatory networks, we show that SVAR can infer true positive edges even under conditions in which the number of samples is smaller than the number of genes. Moreover, it is possible to control for false positives, a significant advantage when compared to other methods described in the literature, which are based on ranks or score functions. By applying SVAR to actual HeLa cell cycle gene expression data, we were able to identify well known transcription factor targets. Conclusion The proposed SVAR method is able to model gene regulatory networks in frequent situations in which the number of samples is lower than the number of genes, making it possible to naturally infer partial Granger causalities without any a priori information. In addition, we present a statistical test to control the false discovery rate, which was not previously possible using other gene regulatory network models.

Multimodality and flexibility of stochastic gene expression

We consider a general class of mathematical models for stochastic gene expression where the transcription rate is allowed to depend on a promoter state variable that can take an arbitrary (finite) number of values. We provide the solution of the master equations in the stationary limit, based on a factorization of the stochastic transition matrix that separates timescales and relative interaction strengths, and we express its entries in terms of parameters that have a natural physical and/or biological interpretation. The solution illustrates the capacity of multiple states promoters to generate multimodal distributions of gene products, without the need for feedback. Furthermore, using the example of a three states promoter operating at low, high, and intermediate expression levels, we show that using multiple states operons will typically lead to a significant reduction of noise in the system. The underlying mechanism is that a three-states promoter can change its level of expression from low to high by passing through an intermediate state with a much smaller increase of fluctuations than by means of a direct transition.

Structural synaptic elements are differentially regulated in superior temporal cortex of schizophrenia patients

Inaccurate wiring and synaptic pathology appear to be major hallmarks of schizophrenia. A variety of gene products involved in synaptic neurotransmission and receptor signaling are differentially expressed in brains of schizophrenia patients. However, synaptic pathology may also develop by improper expression of intra- and extra-cellular structural elements weakening synaptic stability. Therefore, we have investigated transcription of these elements in the left superior temporal gyrus of 10 schizophrenia patients and 10 healthy controls by genome-wide microarrays (Illumina). Fourteen up-regulated and 22 downregulated genes encoding structural elements were chosen from the lists of differentially regulated genes for further qRT-PCR analysis. Almost all genes confirmed by this method were downregulated. Their gene products belonged to vesicle-associated proteins, that is, synaptotagmin 6 and syntaxin 12, to cytoskeletal proteins, like myosin 6, pleckstrin, or to proteins of the extracellular matrix, such as collagens, or laminin C3. Our results underline the pivotal roles of structural genes that control formation and stabilization of pre- and post-synaptic elements or influence axon guidance in schizophrenia. The glial origin of collagen or laminin highlights the close interrelationship between neurons and glial cells in establishment and maintenance of synaptic strength and plasticity. It is hypothesized that abnormal expression of these and related genes has a major impact on the pathophysiology of schizophrenia.

Efficient RNAi-induced Protein Knockdown in Somatic Cells Using Diced or Chemically Produced Small Interfering RNAs (siRNA)

Background: RNA interference (RNAi) is a post-transcriptional gene silencing process in which double-stranded RNA (dsRNA) directs the degradation of a specific corresponding target mRNA. The mediators of this process are small dsRNAs of approximately 21 to 23 bp in length, called small interfering RNAs (siRNAs), which can be prepared in vitro and used to direct the degradation of specific mRNAs inside cells. Hence, siRNAs represent a powerful tool to study and control gene and cell function. Rapid progress has been made in the use of siRNA as a means to attenuate the expression of any protein for which the cDNA sequence is known. Individual siRNAs can be chemically synthesized, in vitro-transcribed, or expressed in cells from siRNA expression vectors. However, screening for the most efficient siRNAs for post-transcriptional gene silencing in cells in culture is a laborious and expensive process. In this study, the effectiveness of two siRNA production strategies for the attenuation of abundant proteins for DNA repair were compared in human cells: (a) the in vitro production of siRNA mixtures by the Dicer enzyme (Diced siRNAs); and (b) the chemical synthesis of very specific and unique siRNA sequences (Stealth RNai (TM)). Materials, Methods & Results: For in vitro-produced siRNAs, two segments of the human Ku70 (167 bp in exon 5; and 249 bp in exon 13; NM001469) and Xrcc4 (172 bp in exon 2; and 108 bp in exon 6; NM003401) genes were chosen to generate dsRNA for subsequent "Dicing" to create mixtures of siRNAs. The Diced fragments of siRNA for each gene sequence were pooled and stored at -80 degrees C. Alternatively, chemically synthesized Stealth siRNAs were designed and generated to match two very specific gene sequence regions for each target gene of interest (Ku70 and Xrcc4). HCT116 cells were plated at 30% confluence in 24- or 6-well culture plates. The next day, cells were transfected by lipofection with either Diced or Stealth siRNAs for Ku70 or Xrcc4, in duplicate, at various doses, with blank and sham transfections used as controls. Cells were harvested at 0, 24, 48, 72 and 96 h post-transfection for protein determination. The knockdown of specific targeted gene products was quantified by Western blot using GAPDH as control. Transfection of gene-specific siRNA to either Ku70 or Xrcc4 with both Diced and Stealth siRNAs resulted in a down regulation of the targeted proteins to approximately 10 to 20% of control levels 48 h after transfection, with recovery to pre-treatment levels by 96 h. Discussion: By transfecting cells with Diced or chemically synthesized Stealth siRNAs, Ku70 and Xrcc4, two highly expressed proteins in cells, were effectively attenuated, demonstrating the great potential for the use of both siRNA production strategies as tools to perform loss of function experiments in mammalian cells. In fact, down-regulation of Ku70 and Xrcc4 has been shown to reduce the activity of the non-homologous end joining DNA pathway, a very desirable approach for the use of homologous recombination technology for gene targeting or knockout studies. Stealth RNAi (TM) was developed to achieve high specificity and greater stability when compared with mixtures of enzymatically-produced (Diced) siRNA fragments. In this study, both siRNA approaches inhibited the expression of Ku70 and Xrcc4 gene products, with no detectable toxic effects to the cells in culture. However, similar knockdown effects using Diced siRNAs were only attained at concentrations 10-fold higher than with Stealth siRNAs. The application of RNAi technology will expand and continue to provide new insights into gene regulation and as potential applications for new therapies, transgenic animal production and basic research.

DM-1, sodium 4-[5-(4-hydroxy-3-methoxyphenyl)-3-oxo-penta-1,4-dienyl]-2-methoxy-phenolate: a curcumin analog with a synergic effect in combination with paclitaxel in breast cancer treatment

This paper describes a new method for the preparation of sodium 4-[5-(4-hydroxy-3-methoxyphenyl)-3-oxo-penta-1,4-dienyl]-2-methoxy-phenolate, DM-1, and 3-oxo-penta-1,4-dienyl-bis (2-methoxy-phenolate), DM-2. The aim of this work was to evaluate the antitumor effects of DM-1 in adjuvant chemotherapy for breast cancer treatment. Mice bearing mammary adenocarcinomas (Ehrlich ascites tumors) were treated with paclitaxel alone, DM-1 alone, and paclitaxel + DM-1. Tumor samples were used to perform cytological analysis by the Papanicolaou method and apoptosis analysis by annexin V and phosphorylated caspase 3. The paclitaxel + DM-1 group had decreased tumor areas and tumor volumes, and the frequency of metastasis was significantly reduced. This caused a decrease in cachexia, which is usually caused by the tumor. Furthermore, treatment with paclitaxel + DM-1 and DM-1 alone increased the occurrence of apoptosis up to 40% in tumor cells, which is 35% more than in the group treated with paclitaxel alone. This cell death was mainly caused through phosphorylated caspase 3 (11% increase in paclitaxel + DM-1 compared to the paclitaxel group), as confirmed by reduced malignancy criteria in the ascitic fluid. DM-1 emerges as a potential treatment for breast cancer and may act as an adjuvant in chemotherapy, enhancing antitumor drug activity with reduced side effects.

Cloning and characterization of a novel alternatively spliced transcript of the human CHD7 putative helicase

Abstract Background The CHD7 (Chromodomain Helicase DNA binding protein 7) gene encodes a member of the chromodomain family of ATP-dependent chromatin remodeling enzymes. Mutations in the CHD7 gene are found in individuals with CHARGE, a syndrome characterized by multiple birth malformations in several tissues. CHD7 was identified as a binding partner of PBAF complex (Polybromo and BRG Associated Factor containing complex) playing a central role in the transcriptional reprogramming process associated to the formation of multipotent migratory neural crest, a transient cell population associated with the genesis of various tissues. CHD7 is a large gene containing 38 annotated exons and spanning 200 kb of genomic sequence. Although genes containing such number of exons are expected to have several alternative transcripts, there are very few evidences of alternative transcripts associated to CHD7 to date indicating that alternative splicing associated to this gene is poorly characterized. Findings Here, we report the cloning and characterization by experimental and computational studies of a novel alternative transcript of the human CHD7 (named CHD7 CRA_e), which lacks most of its coding exons. We confirmed by overexpression of CHD7 CRA_e alternative transcript that it is translated into a protein isoform lacking most of the domains displayed by the canonical isoform. Expression of the CHD7 CRA_e transcript was detected in normal liver, in addition to the DU145 human prostate carcinoma cell line from which it was originally isolated. Conclusions Our findings indicate that the splicing event associated to the CHD7 CRA_e alternative transcript is functional. The characterization of the CHD7 CRA_e novel isoform presented here not only sets the basis for more detailed functional studies of this isoform, but, also, contributes to the alternative splicing annotation of the CHD7 gene and the design of future functional studies aimed at the elucidation of the molecular functions of its gene products.

Horizontal gene transfer confers fermentative metabolism in the respiratory-deficient plant trypanosomatid Phytomonas serpens

Among trypanosomatids, the genus Phytomonas is the only one specifically adapted to infect plants. These hosts provide a particular habitat with a plentiful supply of carbohydrates. Phytomonas sp. lacks a cytochrome-mediated respiratory chain and Krebs cycle, and ATP production relies predominantly on glycolysis. We have characterised the complete gene encoding a putative pyruvate/indolepyruvate decarboxylase (PDC/IPDC) (548 amino acids) of P. serpens, that displays high amino acid sequence similarity with phytobacteria and Leishmania enzymes. No orthologous PDC/IPDC genes were found in Trypanosoma cruzi or T. brucei. Conservation of the PDC/IPDC gene sequence was verified in 14 Phytomonas isolates. A phylogenetic analysis shows that Phytomonas protein is robustly monophyletic with Leishmania spp. and C. fasciculata enzymes. In the trees this clade appears as a sister group of indolepyruvate decarboxylases of gamma-proteobacteria. This supports the proposition that a horizontal gene transfer event from a donor phytobacteria to a recipient ancestral trypanosome has occurred prior to the separation between Phytomonas. Leishmania and Crithidia. We have measured the PDC activity in P. serpens cell extracts. The enzyme has a Km value for pyruvate of 1.4 mM. The acquisition of a PDC, a key enzyme in alcoholic fermentation, explains earlier observations that ethanol is one of the major end-products of glucose catabolism under aerobic and anaerobic conditions. This represents an alternative and necessary route to reoxidise part of the NADH produced in the highly demanding glycolytic pathway and highlights the importance of this type of event in metabolic adaptation. (C) 2012 Elsevier B.V. All rights reserved.

Formaldehyde inhalation reduces respiratory mechanics in a rat model with allergic lung inflammation by altering the nitric oxide/cyclooxygenase-derived products relationship

Bronchial hyperresponsiveness is a hallmark of asthma and many factors modulate bronchoconstriction episodes. A potential correlation of formaldehyde (FA) inhalation and asthma has been observed; however, the exact role of FA remains controversial. We investigated the effects of FA inhalation on Ovalbumin (OVA) sensitisation using a parameter of respiratory mechanics. The involvement of nitric oxide (NO) and cyclooxygenase-derived products were also evaluated. The rats were submitted, or not, to FA inhalation (1%, 90 min/day, 3 days) and were OVA-sensitised and challenged 14 days later. Our data showed that previous FA exposure in allergic rats reduced bronchial responsiveness, respiratory resistance (Rrs) and elastance (Ers) to methacholine. FA exposure in allergic rats also increased the iNOS gene expression and reduced COX-1. L-NAME treatment exacerbated the bronchial hyporesponsiveness and did not modify the Ers and Rrs, while Indomethacin partially reversed all of the parameters studied. The L-NAME and Indomethacin treatments reduced leukotriene B4 levels while they increased thromboxane B2 and prostaglandin E2. In conclusion, FA exposure prior to OVA sensitisation reduces the respiratory mechanics and the interaction of NO and PGE2 may be representing a compensatory mechanism in order to protect the lung from bronchoconstriction effects.

Manipulation of the periovulatory sex steroidal milieu affects endometrial but not luteal gene expression in early diestrus Nelore cows

In beef cattle, the ability to conceive has been associated positively with size of the preovulatory follicle (POF). Proestrus estradiol and subsequent progesterone concentrations can regulate the endometrium to affect receptivity and fertility. The aim of the present study was to verify the effect of the size of the POF on luteal and endometrial gene expression during subsequent early diestrus in beef cattle. Eighty-three multiparous, nonlactating, presynchronized Nelore cows received a progesterone-releasing device and estradiol benzoate on Day–10 (D 10). Animals received cloprostenol (large follicle-large CL group; LF-LCL; N ¼ 42) or not (small follicle-small CL group; SF-SCL; N ¼ 41) on D 10. Progesterone devices were withdrawn and cloprostenol administered 42 to 60 hours (LF-LCL) or 30 to 36 hours (SF-SCL) before GnRH treatment (D0). Tissues were collected at slaughter on D7. The LF-LCL group had larger (P < 0.0001) POF (13.24 0.33 mm vs. 10.76 0.29 mm), greater (P < 0.0007) estradiol concentrations on D0 (2.94 0.28 pg/mL vs. 1.27 0.20 pg/mL), and greater (P < 0.01) progesterone concentrations on D7 (3.71 0.25 ng/mL vs. 2.62 0.26 ng/mL) compared with the SF-SCL group. Luteal gene expression of vascular endothelial growth factor A, kinase insert domain receptor, fms-related tyrosine kinase 1, steroidogenic acute regulatory protein, cytochrome P450, family 11, subfamily A, polypeptide 1, and hydroxy-delta-5-steroid dehydrogenase, 3 beta- and steroid deltaisomerase 7 was similar between groups. Endometrial gene expression of oxytocin receptor and peptidase inhibitor 3, skin-derived was reduced, and estrogen receptor alpha 2, aldo-keto reductase family 1, member C4, and lipoprotein lipase expression was increased in LF-LCL versus SF-SCL. Results support the hypothesis that the size of the POF alters the periovulatory endocrine milieu (i.e., proestrus estradiol and diestrus progesterone concentrations) and acts on the uterus to alter endometrial gene expression. It is proposed that the uterine environment and receptivity might also be modulated. Additionally, it is suggested that increased progesterone secretion of cows ovulating larger follicles is likely due to increased CL size rather than increased luteal expression of steroidogenic genes.

Cytoprotective role of heme oxygenase-1 and heme degradation derived end products in liver injury

The activation of heme oxygenase-1 (HO-1) appears to be an endogenous defensive mechanism used by cells to reduce inflammation and tissue damage in a number of injury models. HO-1, a stress-responsive enzyme that catabolizes heme into carbon monoxide (CO), biliverdin and iron, has previously been shown to protect grafts from ischemia/reperfusion and rejection. In addition, the products of the HO-catalyzed reaction, particularly CO and biliverdin/bilirubin, have been shown to exert protective effects in the liver against a number of stimuli, as in chronic hepatitis C and in transplanted liver grafts. Furthermore, the induction of HO-1 expression can protect the liver against damage caused by a number of chemical compounds. More specifically, the CO derived from HO-1-mediated heme catabolism has been shown to be involved in the regulation of inflammation; furthermore, administration of low concentrations of exogenous CO has a protective effect against inflammation. Both murine and human HO-1 deficiencies have systemic manifestations associated with iron metabolism, such as hepatic overload (with signs of a chronic hepatitis) and iron deficiency anemia (with paradoxical increased levels of ferritin). Hypoxia induces HO-1 expression in multiple rodent, bovine and monkey cell lines, but interestingly, hypoxia represses expression of the human HO-1 gene in a variety of human cell types (endothelial cells, epithelial cells, T cells). These data suggest that HO-1 and CO are promising novel therapeutic molecules for patients with inflammatory diseases. In this review, we present what is currently known regarding the role of HO-1 in liver injuries and in particular, we focus on the implications of targeted induction of HO-1 as a potential therapeutic strategy to protect the liver against chemically induced injury.