Natural gene-expression variation in Down syndrome modulates the outcome of gene-dosage imbalance.

Down syndrome (DS) is characterized by extensive phenotypic variability, with most traits occurring in only a fraction of affected individuals. Substantial gene-expression variation is present among unaffected individuals, and this variation has a strong genetic component. Since DS is caused by genomic-dosage imbalance, we hypothesize that gene-expression variation of human chromosome 21 (HSA21) genes in individuals with DS has an impact on the phenotypic variability among affected individuals. We studied gene-expression variation in 14 lymphoblastoid and 17 fibroblast cell lines from individuals with DS and an equal number of controls. Gene expression was assayed using quantitative real-time polymerase chain reaction on 100 and 106 HSA21 genes and 23 and 26 non-HSA21 genes in lymphoblastoid and fibroblast cell lines, respectively. Surprisingly, only 39% and 62% of HSA21 genes in lymphoblastoid and fibroblast cells, respectively, showed a statistically significant difference between DS and normal samples, although the average up-regulation of HSA21 genes was close to the expected 1.5-fold in both cell types. Gene-expression variation in DS and normal samples was evaluated using the Kolmogorov-Smirnov test. According to the degree of overlap in expression levels, we classified all genes into 3 groups: (A) nonoverlapping, (B) partially overlapping, and (C) extensively overlapping expression distributions between normal and DS samples. We hypothesize that, in each cell type, group A genes are the most dosage sensitive and are most likely involved in the constant DS traits, group B genes might be involved in variable DS traits, and group C genes are not dosage sensitive and are least likely to participate in DS pathological phenotypes. This study provides the first extensive data set on HSA21 gene-expression variation in DS and underscores its role in modulating the outcome of gene-dosage imbalance.

Genome-wide prediction of matrix attachment regions that increase gene expression in mammalian cells.

Gene transfer in eukaryotic cells and organisms suffers from epigenetic effects that result in low or unstable transgene expression and high clonal variability. Use of epigenetic regulators such as matrix attachment regions (MARs) is a promising approach to alleviate such unwanted effects. Dissection of a known MAR allowed the identification of sequence motifs that mediate elevated transgene expression. Bioinformatics analysis implied that these motifs adopt a curved DNA structure that positions nucleosomes and binds specific transcription factors. From these observations, we computed putative MARs from the human genome. Cloning of several predicted MARs indicated that they are much more potent than the previously known element, boosting the expression of recombinant proteins from cultured cells as well as mediating high and sustained expression in mice. Thus we computationally identified potent epigenetic regulators, opening new strategies toward high and stable transgene expression for research, therapeutic production or gene-based therapies.

Minimization of biosynthetic costs in adaptive gene expression responses of yeast to environmental changes

Yeast successfully adapts to an environmental stress by altering physiology and fine-tuning metabolism. This fine-tuning is achieved through regulation of both gene expression and protein activity, and it is shaped by various physiological requirements. Such requirements impose a sustained evolutionary pressure that ultimately selects a specific gene expression profile, generating a suitable adaptive response to each environmental change. Although some of the requirements are stress specific, it is likely that others are common to various situations. We hypothesize that an evolutionary pressure for minimizing biosynthetic costs might have left signatures in the physicochemical properties of proteins whose gene expression is fine-tuned during adaptive responses. To test this hypothesis we analyze existing yeast transcriptomic data for such responses and investigate how several properties of proteins correlate to changes in gene expression. Our results reveal signatures that are consistent with a selective pressure for economy in protein synthesis during adaptive response of yeast to various types of stress. These signatures differentiate two groups of adaptive responses with respect to how cells manage expenditure in protein biosynthesis. In one group, significant trends towards downregulation of large proteins and upregulation of small ones are observed. In the other group we find no such trends. These results are consistent with resource limitation being important in the evolution of the first group of stress responses.

Gene expression following induction of regeneration in Drosophila wing imaginal discs

Background: Regeneration is the ability of an organism to rebuild a body part that has been damaged or amputated, and can be studied at the molecular level using model organisms. Drosophila imaginal discs, which are the larval primordia of adult cuticular structures, are capable of undergoing regenerative growth after transplantation and in vivo culture into the adult abdomen. Results: Using expression profile analyses, we studied the regenerative behaviour of wing discs at 0, 24 and 72 hours after fragmentation and implantation into adult females. Based on expression level, we generated a catalogue of genes with putative role in wing disc regeneration, identifying four classes: 1) genes with differential expression within the first 24 hours; 2) genes with differential expression between 24 and 72 hours; 3) genes that changed significantly in expression levels between the two time periods; 4) genes with a sustained increase or decrease in their expression levels throughout regeneration. Among these genes, we identified members of the JNK and Notch signalling pathways and chromatin regulators. Through computational analysis, we recognized putative binding sites for transcription factors downstream of these pathways that are conserved in multiple Drosophilids, indicating a potential relationship between members of the different gene classes. Experimental data from genetic mutants provide evidence of a requirement of selected genes in wing disc regeneration. Conclusions: We have been able to distinguish various classes of genes involved in early and late steps of the regeneration process. Our data suggests the integration of signalling pathways in the promoters of regulated genes.

Gene expression following induction of regeneration in Drosophila wing imaginal discs

Background: Regeneration is the ability of an organism to rebuild a body part that has been damaged or amputated, and can be studied at the molecular level using model organisms. Drosophila imaginal discs, which are the larval primordia of adult cuticular structures, are capable of undergoing regenerative growth after transplantation and in vivo culture into the adult abdomen. Results: Using expression profile analyses, we studied the regenerative behaviour of wing discs at 0, 24 and 72 hours after fragmentation and implantation into adult females. Based on expression level, we generated a catalogue of genes with putative role in wing disc regeneration, identifying four classes: 1) genes with differential expression within the first 24 hours; 2) genes with differential expression between 24 and 72 hours; 3) genes that changed significantly in expression levels between the two time periods; 4) genes with a sustained increase or decrease in their expression levels throughout regeneration. Among these genes, we identified members of the JNK and Notch signalling pathways and chromatin regulators. Through computational analysis, we recognized putative binding sites for transcription factors downstream of these pathways that are conserved in multiple Drosophilids, indicating a potential relationship between members of the different gene classes. Experimental data from genetic mutants provide evidence of a requirement of selected genes in wing disc regeneration. Conclusions: We have been able to distinguish various classes of genes involved in early and late steps of the regeneration process. Our data suggests the integration of signalling pathways in the promoters of regulated genes.

Metabolic gene expression changes in astrocytes in Multiple Sclerosis cerebral cortex are indicative of immune-mediated signaling.

Emerging as an important correlate of neurological dysfunction in Multiple Sclerosis (MS), extended focal and diffuse gray matter abnormalities have been found and linked to clinical manifestations such as seizures, fatigue and cognitive dysfunction. To investigate possible underlying mechanisms we analyzed the molecular alterations in histopathological normal appearing cortical gray matter (NAGM) in MS. By performing a differential gene expression analysis of NAGM of control and MS cases we identified reduced transcription of astrocyte specific genes involved in the astrocyte-neuron lactate shuttle (ANLS) and the glutamate-glutamine cycle (GGC). Additional quantitative immunohistochemical analysis demonstrating a CX43 loss in MS NAGM confirmed a crucial involvement of astrocytes and emphasizes their importance in MS pathogenesis. Concurrently, a Toll-like/IL-1β signaling expression signature was detected in MS NAGM, indicating that immune-related signaling might be responsible for the downregulation of ANLS and GGC gene expression in MS NAGM. Indeed, challenging astrocytes with immune stimuli such as IL-1β and LPS reduced their ANLS and GGC gene expression in vitro. The detected upregulation of IL1B in MS NAGM suggests inflammasome priming. For this reason, astrocyte cultures were treated with ATP and ATP/LPS as for inflammasome activation. This treatment led to a reduction of ANLS and GGC gene expression in a comparable manner. To investigate potential sources for ANLS and GGC downregulation in MS NAGM, we first performed an adjuvant-driven stimulation of the peripheral immune system in C57Bl/6 mice in vivo. This led to similar gene expression changes in spinal cord demonstrating that peripheral immune signals might be one source for astrocytic gene expression changes in the brain. IL1B upregulation in MS NAGM itself points to a possible endogenous signaling process leading to ANLS and GGC downregulation. This is supported by our findings that, among others, MS NAGM astrocytes express inflammasome components and that astrocytes are capable to release Il-1β in-vitro. Altogether, our data suggests that immune signaling of immune- and/or central nervous system origin drives alterations in astrocytic ANLS and GGC gene regulation in the MS NAGM. Such a mechanism might underlie cortical brain dysfunctions frequently encountered in MS patients.

Impact of a cis-associated gene expression SNP on chromosome 20q11.22 on bipolar disorder susceptibility, hippocampal structure and cognitive performance.

BackgroundBipolar disorder is a highly heritable polygenic disorder. Recent enrichment analyses suggest that there may be true risk variants for bipolar disorder in the expression quantitative trait loci (eQTL) in the brain.AimsWe sought to assess the impact of eQTL variants on bipolar disorder risk by combining data from both bipolar disorder genome-wide association studies (GWAS) and brain eQTL.MethodTo detect single nucleotide polymorphisms (SNPs) that influence expression levels of genes associated with bipolar disorder, we jointly analysed data from a bipolar disorder GWAS (7481 cases and 9250 controls) and a genome-wide brain (cortical) eQTL (193 healthy controls) using a Bayesian statistical method, with independent follow-up replications. The identified risk SNP was then further tested for association with hippocampal volume (n = 5775) and cognitive performance (n = 342) among healthy individuals.ResultsIntegrative analysis revealed a significant association between a brain eQTL rs6088662 on chromosome 20q11.22 and bipolar disorder (log Bayes factor = 5.48; bipolar disorder P = 5.85×10(-5)). Follow-up studies across multiple independent samples confirmed the association of the risk SNP (rs6088662) with gene expression and bipolar disorder susceptibility (P = 3.54×10(-8)). Further exploratory analysis revealed that rs6088662 is also associated with hippocampal volume and cognitive performance in healthy individuals.ConclusionsOur findings suggest that 20q11.22 is likely a risk region for bipolar disorder; they also highlight the informative value of integrating functional annotation of genetic variants for gene expression in advancing our understanding of the biological basis underlying complex disorders, such as bipolar disorder.

Comparison of differential gene expression to water stress among bacteria with relevant pollutant-degradation properties.

Resistance to semi-dry environments has been considered a crucial trait for superior growth and survival of strains used for bioaugmentation in contaminated soils. In order to compare water stress programmes, we analyse differential gene expression among three phylogenetically different strains capable of aromatic compound degradation: Arthrobacter chlorophenolicus A6, Sphingomonas wittichii RW1 and Pseudomonas veronii 1YdBTEX2. Standardized laboratory-induced water stress was imposed by shock exposure of liquid cultures to water potential decrease, induced either by addition of solutes (NaCl, solute stress) or by addition of polyethylene glycol (matric stress), both at absolute similar stress magnitudes and at those causing approximately similar decrease of growth rates. Genome-wide differential gene expression was recorded by micro-array hybridizations. Growth of P. veronii 1YdBTEX2 was the most sensitive to water potential decrease, followed by S. wittichii RW1 and A. chlorophenolicus A6. The number of genes differentially expressed under decreasing water potential was lowest for A. chlorophenolicus A6, increasing with increasing magnitude of the stress, followed by S. wittichii RW1 and P. veronii 1YdBTEX2. Gene inspection and gene ontology analysis under stress conditions causing similar growth rate reduction indicated that common reactions among the three strains included diminished expression of flagellar motility and increased expression of compatible solutes (which were strain-specific). Furthermore, a set of common genes with ill-defined function was found between all strains, including ABC transporters and aldehyde dehydrogenases, which may constitute a core conserved response to water stress. The data further suggest that stronger reduction of growth rate of P. veronii 1YdBTEX2 under water stress may be an indirect result of the response demanding heavy NADPH investment, rather than the presence or absence of a suitable stress defence mechanism per se.

The effects of maternal hyperglycemia on human umbilical vascular gene expression and neonatal rat lung development

Background: Maternal diabetes affects many fetal organ systems, including the vasculature and the lungs. The offspring of diabetic mothers have respiratory adaptation problems after birth. The mechanisms are multifactorial and the effects are prolonged during the postnatal period. An increasing incidence of diabetic pregnancies accentuates the importance of identifying the pathological mechanisms, which cause the metabolic and genetic changes that occur in offspring, born to diabetic mothers. Aims and methods: The aim of this thesis was to determine changes both in human umbilical cord exposed to maternal type 1 diabetes and in neonatal rat lungs after streptozotocin-induced maternal hyperglycemia, during pregnancy. Rat lungs were used as a model for the potential disease mechanisms. Gene expression alterations were determined in human umbilical cords at birth and in rat pup lungs at two week of age. During the first two postnatal weeks, rat lung development was studied morphologically and histologically. Further, the effect of postnatal hyperoxia on hyperglycemia-primed rat lungs was investigated at one week of age to mimic the clinical situation of supplemental oxygen treatment. Results: In the umbilical cord, maternal diabetes had a major negative effect on the expression of genes involved in blood vessel development. The genes regulating vascular tone were also affected. In neonatal rat lungs, intrauterine hyperglycemia had a prolonged effect on gene expression during late alveolarization. The most affected pathway was the upregulation of extracellular matrix proteins. Newborn rat lungs exposed to intrauterine hyperglycemia had thinner saccular walls without changes in airspace size, a smaller relative lung weight and lung total tissue area, and increased cellular apoptosis and proliferation compared to control lungs, possibly reflecting an aberrant maturational adaptation. At one and two weeks of age, cell proliferation and secondary crest formation were accelerated in hyperglycemia-exposed lungs. Postnatal hyperoxic exposure, alone caused arrested alveolarization with thin-walled and enlarged alveoli. In contrast, the dual exposure of intrauterine hyperglycemia and postnatal hyperoxia resulted in the phenotype of thick septa together with arrested alveolarization and decreased number of small pulmonary arteries. Conclusions: Maternal diabetic environment seems to alter the umbilical cord gene expression profile of the regulation of vascular development and function. Fetal hyperglycemia may additionally affect the genetic regulation of the postnatal lung development and may actually induce prolonged structural alterations in neonatal lungs together with a modifying effect on the deleterious pulmonary exposure of postnatal hyperoxia. This, combined with the novel human umbilical cord gene data could serve as stepping stones for future therapies to curb developmental aberrations.

Algorithmic Techniques in Gene Expression Processing. From Imputation to Visualization

The amount of biological data has grown exponentially in recent decades. Modern biotechnologies, such as microarrays and next-generation sequencing, are capable to produce massive amounts of biomedical data in a single experiment. As the amount of the data is rapidly growing there is an urgent need for reliable computational methods for analyzing and visualizing it. This thesis addresses this need by studying how to efficiently and reliably analyze and visualize high-dimensional data, especially that obtained from gene expression microarray experiments. First, we will study the ways to improve the quality of microarray data by replacing (imputing) the missing data entries with the estimated values for these entries. Missing value imputation is a method which is commonly used to make the original incomplete data complete, thus making it easier to be analyzed with statistical and computational methods. Our novel approach was to use curated external biological information as a guide for the missing value imputation. Secondly, we studied the effect of missing value imputation on the downstream data analysis methods like clustering. We compared multiple recent imputation algorithms against 8 publicly available microarray data sets. It was observed that the missing value imputation indeed is a rational way to improve the quality of biological data. The research revealed differences between the clustering results obtained with different imputation methods. On most data sets, the simple and fast k-NN imputation was good enough, but there were also needs for more advanced imputation methods, such as Bayesian Principal Component Algorithm (BPCA). Finally, we studied the visualization of biological network data. Biological interaction networks are examples of the outcome of multiple biological experiments such as using the gene microarray techniques. Such networks are typically very large and highly connected, thus there is a need for fast algorithms for producing visually pleasant layouts. A computationally efficient way to produce layouts of large biological interaction networks was developed. The algorithm uses multilevel optimization within the regular force directed graph layout algorithm.

Gene expression in IFN-gamma-activated murine macrophages

Macrophages are critical for natural immunity and play a central role in specific acquired immunity. The IFN-gamma activation of macrophages derived from A/J or BALB/c mice yielded two different patterns of antiviral state in murine hepatitis virus 3 infection, which were related to a down-regulation of the main virus receptor. Using cDNA hybridization to evaluate mRNA accumulation in the cells, we were able to identify several genes that are differently up- or down-regulated by IFN-gamma in A/J (267 and 266 genes, respectively, up- and down-regulated) or BALB/c (297 and 58 genes, respectively, up- and down-regulated) mouse macrophages. Macrophages from mice with different genetic backgrounds behave differently at the molecular level and comparison of the patterns of non-activated and IFN-gamma-activated A/J or BALB/c mouse macrophages revealed, for instance, an up-regulation and a down-regulation of genes coding for biological functions such as enzymatic reactions, nucleic acid synthesis and transport, protein synthesis, transport and metabolism, cytoskeleton arrangement and extracellular matrix, phagocytosis, resistance and susceptibility to infection and tumors, inflammation, and cell differentiation or activation. The present data are reported in order to facilitate future correlation of proteomic/transcriptomic findings as well as of results obtained from a classical approach for the understanding of biological phenomena. The possible implication of the role of some of the gene products relevant to macrophage biology can now be further scrutinized. In this respect, a down-regulation of the main murine hepatitis virus 3 receptor gene was detected only in IFN-gamma-activated macrophages of resistant mice.

Pim-1 kinase inhibits the activation of reporter gene expression in Elk-1 and c-Fos reporting systems but not the endogenous gene expression: an artifact of the reporter gene assay by transient co-transfection

We have studied the molecular mechanism and signal transduction of pim-1, an oncogene encoding a serine-threonine kinase. This is a true oncogene which prolongs survival and inhibits apoptosis of hematopoietic cells. In order to determine whether the effects of Pim-1 occur by regulation of the mitogen-activated protein kinase pathway, we used a transcriptional reporter assay by transient co-transfection as a screening method. In this study, we found that Pim-1 inhibited the Elk-1 and NFkappaB transcriptional activities induced by activation of the mitogen-activated protein kinase cascade in reporter gene assays. However, Western blots showed that the induction of Elk-1-regulated expression of endogenous c-Fos was not affected by Pim-1. The phosphorylation and activation of neither Erk1/2 nor Elk-1 was influenced by Pim-1. Also, in the gel shift assay, the pattern of endogenous NFkappaB binding to its probe was not changed in any manner by Pim-1. These data indicate that Pim-1 does not regulate the activation of Erk1/2, Elk-1 or NFkappaB. These contrasting results suggest a pitfall of the transient co-transfection reporter assay in analyzing the regulation of transcription factors outside of the chromosome context. It ensures that results from reporter gene expression assay should be verified by study of endogenous gene expression.

Cardiac gene expression and systemic cytokine profile are complementary in a murine model of post-ischemic heart failure

After myocardial infarction (MI), activation of the immune system and inflammatory mechanisms, among others, can lead to ventricular remodeling and heart failure (HF). The interaction between these systemic alterations and corresponding changes in the heart has not been extensively examined in the setting of chronic ischemia. The main purpose of this study was to investigate alterations in cardiac gene and systemic cytokine profile in mice with post-ischemic HF. Plasma was tested for IgM and IgG anti-heart reactive repertoire and inflammatory cytokines. Heart samples were assayed for gene expression by analyzing hybridization to AECOM 32k mouse microarrays. Ischemic HF significantly increased the levels of total serum IgM (by 5.2-fold) and total IgG (by 3.6-fold) associated with a relatively high content of anti-heart specificity. A comparable increase was observed in the levels of circulating pro-inflammatory cytokines such as IL-1β (3.8X) and TNF-α (6.0X). IFN-γ was also increased by 3.1-fold in the MI group. However, IL-4 and IL-10 were not significantly different between the MI and sham-operated groups. Chemokines such as MCP-1 and IL-8 were 1.4- and 13-fold increased, respectively, in the plasma of infarcted mice. We identified 2079 well annotated unigenes that were significantly regulated by post-ischemic HF. Complement activation and immune response were among the most up-regulated processes. Interestingly, 21 of the 101 quantified unigenes involved in the inflammatory response were significantly up-regulated and none were down-regulated. These data indicate that post-ischemic heart remodeling is accompanied by immune-mediated mechanisms that act both systemically and locally.

Hypothyroidism decreases proinsulin gene expression and the attachment of its mRNA and eEF1A protein to the actin cytoskeleton of INS-1E cells

The actions of thyroid hormone (TH) on pancreatic beta cells have not been thoroughly explored, with current knowledge being limited to the modulation of insulin secretion in response to glucose, and beta cell viability by regulation of pro-mitotic and pro-apoptotic factors. Therefore, the effects of TH on proinsulin gene expression are not known. This led us to measure: a) proinsulin mRNA expression, b) proinsulin transcripts and eEF1A protein binding to the actin cytoskeleton, c) actin cytoskeleton arrangement, and d) proinsulin mRNA poly(A) tail length modulation in INS-1E cells cultured in different media containing: i) normal fetal bovine serum - FBS (control); ii) normal FBS plus 1 µM or 10 nM T3, for 12 h, and iii) FBS depleted of TH for 24 h (Tx). A decrease in proinsulin mRNA content and attachment to the cytoskeleton were observed in hypothyroid (Tx) beta cells. The amount of eEF1A protein anchored to the cytoskeleton was also reduced in hypothyroidism, and it is worth mentioning that eEF1A is essential to attach transcripts to the cytoskeleton, which might modulate their stability and rate of translation. Proinsulin poly(A) tail length and cytoskeleton arrangement remained unchanged in hypothyroidism. T3 treatment of control cells for 12 h did not induce any changes in the parameters studied. The data indicate that TH is important for proinsulin mRNA expression and translation, since its total amount and attachment to the cytoskeleton are decreased in hypothyroid beta cells, providing evidence that effects of TH on carbohydrate metabolism also include the control of proinsulin gene expression.

Arginine induces GH gene expression by activating NOS/NO signaling in rat isolated hemi-pituitaries

The amino acid arginine (Arg) is a recognized secretagogue of growth hormone (GH), and has been shown to induce GH gene expression. Arg is the natural precursor of nitric oxide (NO), which is known to mediate many of the effects of Arg, such as GH secretion. Arg was also shown to increase calcium influx in pituitary cells, which might contribute to its effects on GH secretion. Although the mechanisms involved in the effects of Arg on GH secretion are well established, little is known about them regarding the control of GH gene expression. We investigated whether the NO pathway and/or calcium are involved in the effects of Arg on GH gene expression in rat isolated pituitaries. To this end, pituitaries from approximately 170 male Wistar rats (~250 g) were removed, divided into two halves, pooled (three hemi-pituitaries) and incubated or not with Arg, as well as with different pharmacological agents. Arg (71 mM), the NO donor sodium nitroprusside (SNP, 1 and 0.1 mM) and a cyclic guanosine monophosphate (cGMP) analogue (8-Br-cGMP, 1 mM) increased GH mRNA expression 60 min later. The NO acceptor hemoglobin (0.3 µM) blunted the effect of SNP, and the combined treatment with Arg and L-NAME (a NO synthase (NOS) inhibitor, 55 mM) abolished the stimulatory effect of Arg on GH gene expression. The calcium channel inhibitor nifedipine (3 µM) also abolished Arg-induced GH gene expression. The present study shows that Arg directly induces GH gene expression in hemi-pituitaries isolated from rats, excluding interference from somatostatinergic neurons, which are supposed to be inhibited by Arg. Moreover, the data demonstrate that the NOS/NO signaling pathway and calcium mediate the Arg effects on GH gene expression.