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.

Constitutive gene expression profile segregates toxicity in locally advanced breast cancer patients treated with high-dose hyperfractionated radical radiotherapy

[EN] Breast cancer patients show a wide variation in normal tissue reactions after radiotherapy. The individual sensitivity to x-rays limits the efficiency of the therapy. Prediction of individual sensitivity to radiotherapy could help to select the radiation protocol and to improve treatment results. The aim of this study was to assess the relationship between gene expression profiles of ex vivo un-irradiated and irradiated lymphocytes and the development of toxicity due to high-dose hyperfractionated radiotherapy in patients with locally advanced breast cancer. Raw data from microarray experiments were uploaded to the Gene Expression Omnibus Database http://www.ncbi.nlm.nih.gov/geo/ (GEO accession GSE15341). We obtained a small group of 81 genes significantly regulated by radiotherapy, lumped in 50 relevant pathways. Using ANOVA and t-test statistical tools we found 20 and 26 constitutive genes (0 Gy) that segregate patients with and without acute and late toxicity, respectively. Non-supervised hierarchical clustering was used for the visualization of results. Six and 9 pathways were significantly regulated respectively. Concerning to irradiated lymphocytes (2 Gy), we founded 29 genes that separate patients with acute toxicity and without it. Those genes were gathered in 4 significant pathways. We could not identify a set of genes that segregates patients with and without late toxicity. In conclusion, we have found an association between the constitutive gene expression profile of peripheral blood lymphocytes and the development of acute and late toxicity in consecutive, unselected patients. These observations suggest the possibility of predicting normal tissue response to irradiation in high-dose non-conventional radiation therapy regimens. Prospective studies with higher number of patients are needed to validate these preliminary results.

Variation in peach fruit susceptibility to Monilinia laxa and gene expression

Brown rot caused by Monilinia laxa and Monilinia fructigena is considered one of the most important diseases affecting Prunus species. Although some losses can result from the rotten fruits in the orchard, most of the damage is caused to fruits during the post-harvest phase. Several studies reported that brown rot incidence during fruit development highly varies; it was found that at a period corresponding to the the pit hardening stage, fruit susceptibility drastically decreases, to be quickly restored afterwards. However the molecular basis of this phenomenon is still not well understood. Furthermore, no difference in the rot incidence was found between wound and un-wound fruits, suggesting that resistance associated more to a specifc biochemical response of the fruit, rather than to a higher mechanical resistance. So far, the interaction Monilinia-peach was analyzed through chemical approaches. In this study, a bio-molecular approach was undertaken in order to reveal alteration in gene expression associated to the variation of susceptibility. In this thesis three different methods for gene expression analysis were used to analyze the alterations in gene expression occurring in peach fruits during the pit hardening stage, in a period encompassing the temporary change in Monilinia susceptibility: real time PCR, microarray and cDNA AFLP techniques. In 2005, peach fruits (cv.K2) were weekly harvested during a 19-week long-period, starting from the fourth week after full bloom, until full maturity. At each sampling time, three replicates of 5 fruits each were dipped in the M.laxa conidial suspension or in distilled water, as negative control. The fruits were maintained at room temperature for 3 hours; afterwards, they were peeled with a scalpel; the peel was immediately frozen in liquid nitrogen and transferred to -80 °C until use. The degree of susceptibility of peach fruit to the pathogen was determined on 3 replicates of 20 fruits each, as percentage of infected fruits, after one week at 20 °C. Real time PCR analysis was performed to study the variation in expression of those genes encoding for the enzymes of the phenylpropanoid pathway (phenylalanine ammonia lyase (PAL), chalcone synthase (CHS), cinnamate 4-hydroxylase (C4H), leucoanthocyanidine reductase (LAR), hydroxycinnamoyl CoA quinate hydroxycinnamoyl transferase (HQT) and of the jasmonate pathway, such as lipoxygenase (LOX), both involved in the production of important defense compounds. Alteration in gene expression was monitored on fruit samples of a period encompassing the pit hardening stage and the corresponding temporary resistance to M.laxa infections, weekly, from the 6thto the 12th week after full bloom (AFB) inoculated with M. laxa or mock-inoculated. The data suggest a critical change in the expression level of the phenylpropanoid pathway from the 7th to the 8th week AFB; such change could be directly physiologically associated to the peach growth and it could indirectly determine the decrease of susceptibility of peach fruit to Monilinia rot during the subsequent weeks. To investigate on the transcriptome variation underneath the temporary loss of susceptibility of peach fruits to Monilinia rot, the microarray and the cDNA AFLP techniques were used. The samples harvested on the 8th week AFB (named S, for susceptible ones) and on the 12th week AFB (named R, for resistant ones) were compared, both inoculated or mock-inoculated. The microarray experiments were carried out at the University of Padua (Dept. of Environmental Agronomy and Crop Science), using the μPEACH1.0 microarray together with the suited protocols. The analysis showed that 30 genes (corresponding to the 0.6% of the total sequences (4806) contained in the μPeach1.0 microarray) were found up-regulated and 31 ( 0.6%) down regulated in RH vs. SH fruits. On the other hand, 20 genes (0.4%) were shown to be up-regulated and 13 (0.3%) down-regulated in the RI vs. SI fruit. No genes were found differentially expressed in the mock-inoculated resistant fruits (RH) vs. the inoculated resistant ones (RI). Among the up-regulated genes an ATP sulfurylase, an heat shock protein 70, the major allergen Pru P1, an harpin inducing protein and S-adenosylmethionine decarboxylase were found, conversely among the down-regulated ones, cinnamyl alcohol dehydrogenase, an histidine- containing phosphotransfer protein and the ferritin were found. The microarray experimental results and the data indirectly derived, were tested by Real Time PCR analysis. cDNA AFLP analysis was also performed on the same samples. 339 transcript derived fragments considered significant for Monilinia resistance, were selected, sequenced and classified. Genes potentially involved in cell rescue and defence were well represented (8%); several genes (12.1%) involved in the protein folding, post-transductional modification and genes (9.2%) involved in cellular transport were also found. A further 10.3% of genes were classified as involved in the metabolism of aminoacid, carbohydrate and fatty acid. On the other hand, genes involved in the protein synthesis (5.7%) and in signal transduction and communication (5.7%) were found. Among the most interesting genes found differentially expressed between susceptible and resistant fruits, genes encoding for pathogenesis related (PR) proteins were found. To investigate on the association of Monilinia resistance and PR biological function, the major allergen Pru P1 (GenBank accession AM493970) and its isoform (here named Pru P2), were expressed in heterologous system and in vitro assayed for their anti-microbial activity. The ribonuclease activity of the recombinant Pru P1 and Pru P2 proteins was assayed against peach total RNA. As the other PR10 proteins, they showed a ribonucleolytic activity, that could be important to contrast pathogen penetration. Moreover Pru P1 and Pru P2 recombinant proteins were checked for direct antimicrobial activity. No inhibitory effect of Pru P1 or Pru P2 was detected against the selected fungi.

Streptococcus agalactiae adapts to glucose stress conditions by modulating gene expression profile

Diabetes mellitus is considered a risk factor for Group B Streptococcus (GBS) infections. Typically, this pathology is associated to high glucose levels in the bloodstream. Although clinical evidences support this notion, the physiological mechanisms underlying GBS adaptation to such conditions are not yet defined. In the attempt to address this issue, we performed comparative global gene expression analysis of GBS grown under glucose-stress conditions and observed that a number of metabolic and virulence genes was differentially regulated. Of importance, we also demonstrated that by knocking-out the csrRS locus the transcription profile of GBS grown in high-glucose conditions was profoundly affected, with more than a third of glucose-dependent genes, including the virulence factor bibA, found to be controlled by this two-component system. Furthermore, in vitro molecular analysis showed that CsrR specifically binds to the bibA promoter and the phosphorilation increases the affinity of the regulator to this promoter region. Moreover, we demonstrated that CsrR acts as a repressor of bibA expression by binding to its promoter in vivo. In conclusion, this work by elucidating both the response of GBS to pathological glucose conditions and the underlined molecular mechanisms will set the basis for a better understanding of GBS pathogenesis.

Feeding entrainment of locomotor activity and clock gene expression in Senegalese sole, Solea senegalensis

The present study was conducted to investigate the influence of restricted food access on Solea senegalensis behaviour and daily expression of clock genes in central (diencephalon and optic tectum) and pheripheral (liver) tissues. The Senegalese sole is a marine teleost fish belonging to the Class of Actinopterygii, Order Pleuronectiformes and Family Soleidae. Its geographical distribution in the Mediterranean sea is fairly broad, covering the south and east of the Iberian Peninsula, the North of Africa and Middle East until the coast of Turkey. From a commercial perspective Solea senegalensis has acquired in recent years, a key role in aquacolture industry of the Iberian Peninsula. The Senegalese sole is also acquiring an important relevance in chronobiological studies as the number of published works focused on the sole circadian system has increased in the last few years. The molecular mechanisms underlying sole circadian rhythms has also been explored recently, both in adults and developing sole. Moreover, the consideration of the Pleuronectiformes Order as one of the most evolved teleost groups make the Senegalese sole a species of high interest under a comparative and phylogenetic point of view. All these facts have reinforced the election of Senegalese sole as model species for the present study. The animals were kept under 12L:12D photoperiod conditions and divided into three experimental groups depending on the feeding time: fed at midlight (ML), middark (MD) or random (RND) times. Throughout the experiment, the existence of a daily activity rhythm and it synchronization to the light-dark and feeding cycles was checked. To this end locomotor activity was registred by means of two infrared photocells placed in pvc tube 10 cm below the water surface (upper photocell) and the other one was located 10 cm above the bottom of the tank (bottom photocell). The photocell were connected to a computer so that every time a fish interrupted the infrared light beam, it produced an output signal that was recorded. The number of light beam interruptions was stored every 10 minutes by specialized software for data acquisition.

Genus Bifidobacterium: taxonomy studies and gene expression analysis on folate pathway

Folates (vitamin B9) are essential water soluble vitamins, whose deficiency in humans may contribute to the onset of several diseases, such as anaemia, cancer, cardiovascular diseases, neurological problems as well as defects in embryonic development. Human and other mammals are unable to synthesize ex novo folate obtaining it from exogenous sources, via intestinal absorption. Recently the gut microbiota has been identified as an important source of folates and the selection and use of folate producing microorganisms represents an innovative strategy to increase human folate levels. The aim of this thesis was to gain a fundamental understanding of folate metabolism in Bifidobacterium adolescentis. The work was subdivided in three main phases, also aimed to solve different problems encountered working with Bifidobacterium strains. First, a new identification method (based on PCR-RFLP of hsp60 gene) was specifically developed to identify Bifidobacterium strains. Secondly, Bifidobacterium adolescentis biodiversity was explored in order to recognize representing strains of this species to be screened for their folate production ability. Results showed that this species is characterized by a wide variability and support the idea that a possible new taxonomic re-organization would be required. Finally B. adolescentis folate metabolism was studied using a double approach. A quantitative analysis of folate content was complemented by the examination of expression levels of genes involved in folate related pathways. For the normalization process, required to increase the robustness of the qRT-PCR analysis, an appropriate set of reference genes was tested using two different algorithms. Results demonstrate that B.adolescentis strains may represent an endogenous source of natural folate and they could be used to fortify fermented dairy products. This bio-fortification strategy presents many advantages for the consumer, providing native folate forms more bio-available, and not implicated in the discussed controversy concerning the safety of high intake of synthetic folic acid.

Mechanistic study of the effects of 8-oxo-7,8-dihydroguanine (8-oxoG) on reporter gene expression in mammalian cells

DNA damage causes replication errors, leading to genetic instability or cell death. Besides that, many types of DNA base modifications have been shown to interfere with transcriptional elongation if they are located in the transcribed DNA strand of active genes, acting as roadblocks for RNA polymerases. It is widely assumed that transcription blockage by endogenous DNA damage is responsible for the early cell senescence in organs and accelerated ageing observed in individuals with compromised nucleotide excision repair.rnThe aims of this work were to design new experimental systems for testing transcription blocking potentials of DNA base modifications in an individual gene and to apply these test systems to the investigation of the effects of a frequent endogenously generated base modification, namely 8-oxo-7,8-hydroxyguanine (8-oxoG), on the gene transcription in cells. Several experimental strategies were employed for this purpose. First, I constructed an episomal vector encoding for a short-lived EGFP-ODC fusion protein and measured expression of the reporter gene in permanently transfected clonal cell lines exposed to DNA damaging agents. Second, the expression of plasmid-borne EGFP gene damaged with photosensitisers to obtain one or several oxidative purine modifications per plasmid molecule was determined in transiently transfected human and mouse host cells in an approach known as “host cell reactivation”. As a prerequisite for these experiments, a robust method of precise quantitative measurement of the EGFP gene expression in transiently transfected cells by flow cytometry was developed and validated. Third, I elaborated a very efficient procedure for insertion of synthetic oligonucleotides carrying 8-oxoG into plasmid DNA, avoiding any unwanted base damage and strand breaks. The consequences of 8-oxoG placed in defined positions in opposing DNA strands of the EGFP gene for transcription were measured by host cell reactivation in cells with functional 8-oxoguanine DNA glycosylase (OGG1) gene and in OGG1 null cells.rnThe results obtained in Ogg1-/- cells demonstrated that unrepaired 8-oxoG, even if situated in the transcribed DNA strand, does not have any negative effect on the reporter gene transcription. On the other hand, as few as one 8-oxoG was sufficient to cause a significant decrease of the gene expression in OGG1-proficient cell lines, i.e. in the presence of base excision repair. For two analysed positions of 8-oxoG in the plasmid DNA, the inhibition of gene transcription by the base modification correlated with the efficiency of its excision by purified OGG1 protein under cell-free conditions. Based on these findings, it has to be concluded that the observed decrease of transcription is mediated by excision of the base modification by OGG1 and probably caused by the repair-induced single-strand breaks. The mechanism of transcription inhibition by 8-oxoG is therefore clearly distinct from stalling of elongating RNA polymerase II complexes at the modified base.

Gene expression responses to acute and chronic heat stress in the common reef-building coral Pocillopora verrucosa

Global climate change is impacting coral reefs worldwide, with approximately 19% of reefs being permanently degraded, 15% showing symptoms of imminent collapse, and 20% at risk of becoming critically affected in the next few decades. This alarming level of reef degradation is mainly due to an increase in frequency and intensity of natural and anthropogenic disturbances. Recent evidence has called into question whether corals have the capacity to acclimatize or adapt to climate changes and some groups of corals showed inherent physiological tolerance to environmental stressors. The aim of the present study was to evaluate mRNA expression patterns underlying differences in thermal tolerance in specimen of the common reef-building coral Pocillopora verrucosa collected at different locations in Bangka Island waters (North Sulawesi, Indonesia). Part of the experimental work was carried out at the CoralEye Reef Research Outpost (Bangka Island). This includes sampling of corals at selected sites and at different depths (3 and 12 m) as well as their experimental exposure to an increased water temperature under controlled conditions for 3 and 7 days. Levels of mRNAs encoding ATP synthase (ATPs) NADH dehydrogenase (NDH) and a 70kDa Heat Shock Protein (HSP70) were evaluated by quantitative real time PCR. Transcriptional profiles evaluated under field conditions suggested an adaptation to peculiar local environmental conditions in corals collected at different sites and at the low depth. Nevertheless, high–depth collected corals showed a less pronounced site-to-site separation suggesting more homogenous environmental conditions. Exposure to an elevated temperature under controlled conditions pointed out that corals adapted to the high depth are more sensitive to the effects of thermal stress, so that reacted to thermal challenge by significantly over-expressing the selected gene products. Being continuously exposed to fluctuating environmental conditions, low-depth adapted corals are more resilient to the stress stimulus, and indeed showed unaffected or down-regulated mRNA expression profiles. Overall these results highlight that transcriptional profiles of selected genes involved in cellular stress response are modulated by natural seasonal temperature changes in P. verrucosa. Moreover, specimens living in more variable habitats (low-depth) exhibit higher basal HSP70 mRNA levels, possibly enhancing physiological tolerance to environmental stressors.

Hepatic gene expression profile in mice perorally infected with Echinococcus multilocularis eggs

Alveolar echinococcosis (AE) is a severe chronic hepatic parasitic disease currently emerging in central and eastern Europe. Untreated AE presents a high mortality (>90%) due to a severe hepatic destruction as a result of parasitic metacestode proliferation which behaves like a malignant tumor. Despite this severe course and outcome of disease, the genetic program that regulates the host response leading to organ damage as a consequence of hepatic alveolar echinococcosis is largely unknown.

Pheochromocytoma in rats with multiple endocrine neoplasia (MENX) shares gene expression patterns with human pheochromocytoma

Pheochromocytomas are rare neoplasias of neural crest origin arising from chromaffin cells of the adrenal medulla and sympathetic ganglia (extra-adrenal pheochromocytoma). Pheochromocytoma that develop in rats homozygous for a loss-of-function mutation in p27Kip1 (MENX syndrome) show a clear progression from hyperplasia to tumor, offering the possibility to gain insight into tumor pathobiology. We compared the gene-expression signatures of both adrenomedullary hyperplasia and pheochromocytoma with normal rat adrenal medulla. Hyperplasia and tumor show very similar transcriptome profiles, indicating early determination of the tumorigenic signature. Overrepresentation of developmentally regulated neural genes was a feature of the rat lesions. Quantitative RT-PCR validated the up-regulation of 11 genes, including some involved in neural development: Cdkn2a, Cdkn2c, Neurod1, Gal, Bmp7, and Phox2a. Overexpression of these genes precedes histological changes in affected adrenal glands. Their presence at early stages of tumorigenesis indicates they are not acquired during progression and may be a result of the lack of functional p27Kip1. Adrenal and extra-adrenal pheochromocytoma development clearly follows diverged molecular pathways in MENX rats. To correlate these findings to human pheochromocytoma, we studied nine genes overexpressed in the rat lesions in 46 sporadic and familial human pheochromocytomas. The expression of GAL, DGKH, BMP7, PHOX2A, L1CAM, TCTE1, EBF3, SOX4, and HASH1 was up-regulated, although with different frequencies. Immunohistochemical staining detected high L1CAM expression selectively in 27 human pheochromocytomas but not in 140 nonchromaffin neuroendocrine tumors. These studies reveal clues to the molecular pathways involved in rat and human pheochromocytoma and identify previously unexplored biomarkers for clinical use.