Genes involved in thoracic exoskeleton formation during the pupal-to-adult molt in a social insect model, Apis mellifera

Background: The insect exoskeleton provides shape, waterproofing, and locomotion via attached somatic muscles. The exoskeleton is renewed during molting, a process regulated by ecdysteroid hormones. The holometabolous pupa transforms into an adult during the imaginal molt, when the epidermis synthe3sizes the definitive exoskeleton that then differentiates progressively. An important issue in insect development concerns how the exoskeletal regions are constructed to provide their morphological, physiological and mechanical functions. We used whole-genome oligonucleotide microarrays to screen for genes involved in exoskeletal formation in the honeybee thoracic dorsum. Our analysis included three sampling times during the pupal-to-adult molt, i.e., before, during and after the ecdysteroid-induced apolysis that triggers synthesis of the adult exoskeleton. Results: Gene ontology annotation based on orthologous relationships with Drosophila melanogaster genes placed the honeybee differentially expressed genes (DEGs) into distinct categories of Biological Process and Molecular Function, depending on developmental time, revealing the functional elements required for adult exoskeleton formation. Of the 1,253 unique DEGs, 547 were upregulated in the thoracic dorsum after apolysis, suggesting induction by the ecdysteroid pulse. The upregulated gene set included 20 of the 47 cuticular protein (CP) genes that were previously identified in the honeybee genome, and three novel putative CP genes that do not belong to a known CP family. In situ hybridization showed that two of the novel genes were abundantly expressed in the epidermis during adult exoskeleton formation, strongly implicating them as genuine CP genes. Conserved sequence motifs identified the CP genes as members of the CPR, Tweedle, Apidermin, CPF, CPLCP1 and Analogous-to-Peritrophins families. Furthermore, 28 of the 36 muscle-related DEGs were upregulated during the de novo formation of striated fibers attached to the exoskeleton. A search for cis-regulatory motifs in the 5′-untranslated region of the DEGs revealed potential binding sites for known transcription factors. Construction of a regulatory network showed that various upregulated CP- and muscle-related genes (15 and 21 genes, respectively) share common elements, suggesting co-regulation during thoracic exoskeleton formation. Conclusions: These findings help reveal molecular aspects of rigid thoracic exoskeleton formation during the ecdysteroid-coordinated pupal-to-adult molt in the honeybee.

Cellular interference in craniofrontonasal syndrome: males mosaic for mutations in the X-linked EFNB1 gene are more severely affected than true hemizygotes

Craniofrontonasal syndrome (CFNS), an X-linked disorder caused by loss-of-function mutations of EFNB1, exhibits a paradoxical sex reversal in phenotypic severity: females characteristically have frontonasal dysplasia, craniosynostosis and additional minor malformations, but males are usually more mildly affected with hypertelorism as the only feature. X-inactivation is proposed to explain the more severe outcome in heterozygous females, as this leads to functional mosaicism for cells with differing expression of EPHRIN-B1, generating abnormal tissue boundaries-a process that cannot occur in hemizygous males. Apparently challenging this model, males occasionally present with a more severe female-like CFNS phenotype. We hypothesized that such individuals might be mosaic for EFNB1 mutations and investigated this possibility in multiple tissue samples from six sporadically presenting males. Using denaturing high performance liquid chromatography, massively parallel sequencing and multiplex-ligation-dependent probe amplification (MLPA) to increase sensitivity above standard dideoxy sequencing, we identified mosaic mutations of EFNB1 in all cases, comprising three missense changes, two gene deletions and a novel point mutation within the 5' untranslated region (UTR). Quantification by Pyrosequencing and MLPA demonstrated levels of mutant cells between 15 and 69%. The 5' UTR variant mutates the stop codon of a small upstream open reading frame that, using a dual-luciferase reporter construct, was demonstrated to exacerbate interference with translation of the wild-type protein. These results demonstrate a more severe outcome in mosaic than in constitutionally deficient males in an X-linked dominant disorder and provide further support for the cellular interference mechanism, normally related to X-inactivation in females.

Pair approximation for a model of vertical and horizontal transmission of parasites.

We apply Stochastic Dynamics method for a differential equations model, proposed by Marc Lipsitch and collaborators (Proc. R. Soc. Lond. B 260, 321, 1995), for which the transmission dynamics of parasites occurs from a parent to its offspring (vertical transmission), and by contact with infected host (horizontal transmission). Herpes, Hepatitis and AIDS are examples of diseases for which both horizontal and vertical transmission occur simultaneously during the virus spreading. Understanding the role of each type of transmission in the infection prevalence on a susceptible host population may provide some information about the factors that contribute for the eradication and/or control of those diseases. We present a pair mean-field approximation obtained from the master equation of the model. The pair approximation is formed by the differential equations of the susceptible and infected population densities and the differential equations of pairs that contribute to the former ones. In terms of the model parameters, we obtain the conditions that lead to the disease eradication, and set up the phase diagram based on the local stability analysis of fixed points. We also perform Monte Carlo simulations of the model on complete graphs and Erdös-Rényi graphs in order to investigate the influence of population size and neighborhood on the previous mean-field results; by this way, we also expect to evaluate the contribution of vertical and horizontal transmission on the elimination of parasite. Pair Approximation for a Model of Vertical and Horizontal Transmission of Parasites.

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.

Exome analysis of HIV patients submitted to dendritic cells therapeutic vaccine reveals an association of CNOT1 gene with response to the treatment

INTRODUCTION: With the aim of searching genetic factors associated with the response to an immune treatment based on autologous monocyte-derived dendritic cells pulsed with autologous inactivated HIV, we performed exome analysis by screening more than 240,000 putative functional exonic variants in 18 HIV-positive Brazilian patients that underwent the immune treatment. METHODS: Exome analysis has been performed using the ILLUMINA Infinium HumanExome BeadChip. zCall algorithm allowed us to recall rare variants. Quality control and SNP-centred analysis were done with GenABEL R package. An in-house implementation of the Wang method permitted gene-centred analysis. RESULTS: CCR4-NOT transcription complex, subunit 1 (CNOT1) gene (16q21), showed the strongest association with the modification of the response to the therapeutic vaccine (p=0.00075). CNOT1 SNP rs7188697 A/G was significantly associated with DC treatment response. The presence of a G allele indicated poor response to the therapeutic vaccine (p=0.0031; OR=33.00; CI=1.74-624.66), and the SNP behaved in a dominant model (A/A vs. A/G+G/G p=0.0009; OR=107.66; 95% CI=3.85-3013.31), being the A/G genotype present only in weak/transient responders, conferring susceptibility to poor response to the immune treatment. DISCUSSION: CNOT1 is known to be involved in the control of mRNA deadenylation and mRNA decay. Moreover, CNOT1 has been recently described as being involved in the regulation of inflammatory processes mediated by tristetraprolin (TTP). The TTP-CCR4-NOT complex (CNOT1 in the CCR4-NOT complex is the binding site for TTP) has been reported as interfering with HIV replication, through post-transcriptional control. Therefore, we can hypothesize that genetic variation occurring in the CNOT1 gene could impair the TTP-CCR4-NOT complex, thus interfering with HIV replication and/or host immune response. CONCLUSIONS: Being aware that our findings are exclusive to the 18 patients studied with a need for replication, and that the genetic variant of CNOT1 gene, localized at intron 3, has no known functional effect, we propose a novel potential candidate locus for the modulation of the response to the immune treatment, and open a discussion on the necessity to consider the host genome as another potential variant to be evaluated when designing an immune therapy study

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.

Angiogenesis and angioregression gene expression analyses in swine corpus luteum

The corpus luteum (CL) lifespan is characterized by a rapid growth, differentiation and controlled regression of the luteal tissue, accompanied by an intense angiogenesis and angioregression. Indeed, the CL is one of the most highly vascularised tissue in the body with a proliferation rate of the endothelial cells 4- to 20-fold more intense than in some of the most malignant human tumours. This angiogenic process should be rigorously controlled to allow the repeated opportunities of fertilization. After a first period of rapid growth, the tissue becomes stably organized and prepares itself to switch to the phenotype required for its next apoptotic regression. In pregnant swine, the lifespan of the CLs must be extended to support embryonic and foetal development and vascularisation is necessary for the maintenance of luteal function. Among the molecules involved in the angiogenesis, Vascular Endothelial Growth Factor (VEGF) is the main regulator, promoting endothelial cells proliferation, differentiation and survival as well as vascular permeability and vessel lumen formation. During vascular invasion and apoptosis process, the remodelling of the extracellular matrix is essential for the correct evolution of the CL, particularly by the action of specific class of proteolytic enzymes known as matrix metalloproteinases (MMPs). Another important factor that plays a role in the processes of angiogenesis and angioregression during the CL formation and luteolysis is the isopeptide Endothelin-1 (ET-1), which is well-known to be a potent vasoconstrictor and mitogen for endothelial cells. The goal of the present thesis was to study the role and regulation of vascularisation in an adult vascular bed. For this purpose, using a precisely controlled in vivo model of swine CL development and regression, we determined the levels of expression of the members of VEGF system (VEGF total and specific isoforms; VEGF receptor-1, VEGFR-1; VEGF receptor-2, VEGFR-2) and ET- 1 system (ET-1; endothelin converting enzyme-1, ECE-1; endothelin receptor type A, ET-A) as well as the activity of the Ca++/Mg++-dependent endonucleases and gelatinases (MMP-2 and MMP-9). Three experiments were conducted to reach such objectives in CLs isolated from ovaries of cyclic, pregnant or fasted gilts. In the Experiment I, we evaluated the influence of acute fasting on VEGF production and VEGF, VEGFR-2, ET-1, ECE-1 and ET-A mRNA expressions in CLs collected on day 6 after ovulation (midluteal phase). The results indicated a down-regulation of VEGF, VEGFR-2, ET-1 and ECE-1 mRNA expression, although no change was observed for VEGF protein. Furthermore, we observed that fasting stimulated steroidogenesis by luteal cells. On the basis of the main effects of VEGF (stimulation of vessel growth and endothelial permeability) and ET-1 (stimulation of endothelial cell proliferation and vasoconstriction, as well as VEGF stimulation), we concluded that feed restriction possibly inhibited luteal vessel development. This could be, at least in part, compensated by a decrease of vasal tone due to a diminution of ET-1, thus ensuring an adequate blood flow and the production of steroids by the luteal cells. In the Experiment II, we investigated the relationship between VEGF, gelatinases and Ca++/Mg++-dependent endonucleases activities with the functional CL stage throughout the oestrous cycle and at pregnancy. The results demonstrated differential patterns of expression of those molecules in correspondence to the different phases of the oestrous cycle. Immediately after ovulation, VEGF mRNA/protein levels and MMP-9 activity are maximal. On days 5–14 after ovulation, VEGF expression and MMP-2 and -9 activities are at basal levels, while Ca++/Mg++-dependent endonuclease levels increased significantly in relation to day 1. Only at luteolysis (day 17), Ca++/Mg++-dependent endonuclease and MMP-2 spontaneous activity increased significantly. At pregnancy, high levels of MMP-9 and VEGF were observed. These results suggested that during the very early luteal phase, high MMPs activities coupled with high VEGF levels drive the tissue to an angiogenic phenotype, allowing CL growth under LH (Luteinising Hormone) stimulus, while during the late luteal phase, low VEGF and elevate MMPs levels may play a role in the apoptotic tissue and extracellular matrix remodelling during structural luteolysis. In the Experiment III, we described the expression patterns of all distinct VEGF isoforms throughout the oestrous cycle. Furthermore, the mRNA expression and protein levels of both VEGF receptors were also evaluated. Four novel VEGF isoforms (VEGF144, VEGF147, VEGF182, and VEGF164b) were found for the first time in swine and the seven identified isoforms presented four different patterns of expression. All isoforms showed their highest mRNA levels in newly formed CLs (day 1), followed by a decrease during mid-late luteal phase (days 10–17), except for VEGF182, VEGF188 and VEGF144 that showed a differential regulation during late luteal phase (day 14) or at luteolysis (day 17). VEGF protein levels paralleled the most expressed and secreted VEGF120 and VEGF164 isoforms. The VEGF receptors mRNAs showed a different pattern of expression in relation to their ligands, increasing between day 1 and 3 and gradually decreasing during the mid-late luteal phase. The differential regulation of some VEGF isoforms principally during the late luteal phase and luteolysis suggested a specific role of VEGF during tissue remodelling process that occurs either for CL maintenance in case of pregnancy or for noncapillary vessel development essential for tissue removal during structural luteolysis. In summary, our findings allow us to determine relationships among factors involved in the angiogenesis and angioregression mechanisms that take place during the formation and regression of the CL. Thus, CL provides a very interesting model for studying such factors in different fields of the basic research.

Inibizione selettiva del gene MYCN mediante PNA (acidi peptido nucleici) anti-gene nel rabdomiosarcoma umano

MYCN oncogene amplification/expression is a feature of many childhood tumors, and some adult tumors, and it is associated with poor prognosis. While MYC expression is ubiquitary, MYCN has a restricted expression after birth and it is an ideal target for an effective therapy. PNAs belong to the latest class of nucleic acid-based therapeutics, and they can bind chromosomal DNA and block gene transcription (anti-gene activity). We have developed an anti-gene PNA that targets specifically the MYCN gene to block its transcription. We report for the first time MYCN targeted inhibition in Rhabdomyosarcoma (RMS) by the anti-MYCN-PNA in RMS cell lines (four ARMS and four ERMS) and in a xenograft RMS mouse model. Rhabdomyosarcoma is the most common pediatric soft-tissue sarcoma, comprising two main subgroups [Alveolar (ARMS) and Embryonal (ERMS)]. ARMS is associated with a poorer prognosis. MYCN amplification is a feature of both the ERMS and ARMS, but the MYCN amplification and expression levels shows a significant correlation and are greater in ARMS, in which they are associated with adverse outcome. We found that MYCN mRNA and protein levels were higher in the four ARMS (RH30, RH4, RH28 and RMZ-RC2) than in the four ERMS (RH36, SMS-CTR, CCA and RD) cell lines. The potent inhibition of MYCN transcription was highly specific, it did not affect the MYC expression, it was followed by cell-growth inhibition in the RMS cell lines which correlated with the MYCN expression rate, and it led to complete cell-growth inhibition in ARMS cells. We used a mutated- PNA as control. MYCN silencing induced apoptosis. Global gene expression analysis (Affymetrix microarrays) in ARMS cells treated with the anti-MYCN-PNA revealed genes specifically induced or repressed, with both genes previously described as targets of N-myc or Myc, and new genes undescribed as targets of N-myc or Myc (mainly involved in cell cycle, apoptosis, cell motility, metastasis, angiogenesis and muscle development). The changes in the expression of the most relevant genes were confirmed by Real-Time PCR and western blot, and their expression after the MYCN silencing was evaluated in the other RMS cell lines. The in vivo study, using an ARMS xenograft murine model evaluated by micro-PET, showed a complete elimination of the metabolic tumor signal in most of the cases (70%) after anti-MYCN-PNA treatment (without toxicity), whereas treatment with the mutated-PNA had no effect. Our results strongly support the development of MYCN anti-gene therapy for the treatment of RMS, particularly for poor prognosis ARMS, and of other MYCN-expressing tumors.

Transcriptional regulation of human mu-opioid receptor gene: functional characterization of activating and inhibitory transcription factors

The organization of the nervous and immune systems is characterized by obvious differences and striking parallels. Both systems need to relay information across very short and very long distances. The nervous system communicates over both long and short ranges primarily by means of more or less hardwired intercellular connections, consisting of axons, dendrites, and synapses. Longrange communication in the immune system occurs mainly via the ordered and guided migration of immune cells and systemically acting soluble factors such as antibodies, cytokines, and chemokines. Its short-range communication either is mediated by locally acting soluble factors or transpires during direct cell–cell contact across specialized areas called “immunological synapses” (Kirschensteiner et al., 2003). These parallels in intercellular communication are complemented by a complex array of factors that induce cell growth and differentiation: these factors in the immune system are called cytokines; in the nervous system, they are called neurotrophic factors. Neither the cytokines nor the neurotrophic factors appear to be completely exclusive to either system (Neumann et al., 2002). In particular, mounting evidence indicates that some of the most potent members of the neurotrophin family, for example, nerve growth factor (NGF) and brainderived neurotrophic factor (BDNF), act on or are produced by immune cells (Kerschensteiner et al., 1999) There are, however, other neurotrophic factors, for example the insulin-like growth factor-1 (IGF-1), that can behave similarly (Kermer et al., 2000). These factors may allow the two systems to “cross-talk” and eventually may provide a molecular explanation for the reports that inflammation after central nervous system (CNS) injury has beneficial effects (Moalem et al., 1999). In order to shed some more light on such a cross-talk, therefore, transcription factors modulating mu-opioid receptor (MOPr) expression in neurons and immune cells are here investigated. More precisely, I focused my attention on IGF-I modulation of MOPr in neurons and T-cell receptor induction of MOPr expression in T-lymphocytes. Three different opioid receptors [mu (MOPr), delta (DOPr), and kappa (KOPr)] belonging to the G-protein coupled receptor super-family have been cloned. They are activated by structurallyrelated exogenous opioids or endogenous opioid peptides, and contribute to the regulation of several functions including pain transmission, respiration, cardiac and gastrointestinal functions, and immune response (Zollner and Stein 2007). MOPr is expressed mainly in the central nervous system where it regulates morphine-induced analgesia, tolerance and dependence (Mayer and Hollt 2006). Recently, induction of MOPr expression in different immune cells induced by cytokines has been reported (Kraus et al., 2001; Kraus et al., 2003). The human mu-opioid receptor gene (OPRM1) promoter is of the TATA-less type and has clusters of potential binding sites for different transcription factors (Law et al. 2004). Several studies, primarily focused on the upstream region of the OPRM1 promoter, have investigated transcriptional regulation of MOPr expression. Presently, however, it is still not completely clear how positive and negative transcription regulators cooperatively coordinate cellor tissue-specific transcription of the OPRM1 gene, and how specific growth factors influence its expression. IGF-I and its receptors are widely distributed throughout the nervous system during development, and their involvement in neurogenesis has been extensively investigated (Arsenijevic et al. 1998; van Golen and Feldman 2000). As previously mentioned, such neurotrophic factors can be also produced and/or act on immune cells (Kerschenseteiner et al., 2003). Most of the physiologic effects of IGF-I are mediated by the type I IGF surface receptor which, after ligand binding-induced autophosphorylation, associates with specific adaptor proteins and activates different second messengers (Bondy and Cheng 2004). These include: phosphatidylinositol 3-kinase, mitogen-activated protein kinase (Vincent and Feldman 2002; Di Toro et al. 2005) and members of the Janus kinase (JAK)/STAT3 signalling pathway (Zong et al. 2000; Yadav et al. 2005). REST plays a complex role in neuronal cells by differentially repressing target gene expression (Lunyak et al. 2004; Coulson 2005; Ballas and Mandel 2005). REST expression decreases during neurogenesis, but has been detected in the adult rat brain (Palm et al. 1998) and is up-regulated in response to global ischemia (Calderone et al. 2003) and induction of epilepsy (Spencer et al. 2006). Thus, the REST concentration seems to influence its function and the expression of neuronal genes, and may have different effects in embryonic and differentiated neurons (Su et al. 2004; Sun et al. 2005). In a previous study, REST was elevated during the early stages of neural induction by IGF-I in neuroblastoma cells. REST may contribute to the down-regulation of genes not yet required by the differentiation program, but its expression decreases after five days of treatment to allow for the acquisition of neural phenotypes. Di Toro et al. proposed a model in which the extent of neurite outgrowth in differentiating neuroblastoma cells was affected by the disappearance of REST (Di Toro et al. 2005). The human mu-opioid receptor gene (OPRM1) promoter contains a DNA sequence binding the repressor element 1 silencing transcription factor (REST) that is implicated in transcriptional repression. Therefore, in the fist part of this thesis, I investigated whether insulin-like growth factor I (IGF-I), which affects various aspects of neuronal induction and maturation, regulates OPRM1 transcription in neuronal cells in the context of the potential influence of REST. A series of OPRM1-luciferase promoter/reporter constructs were transfected into two neuronal cell models, neuroblastoma-derived SH-SY5Y cells and PC12 cells. In the former, endogenous levels of human mu-opioid receptor (hMOPr) mRNA were evaluated by real-time PCR. IGF-I upregulated OPRM1 transcription in: PC12 cells lacking REST, in SH-SY5Y cells transfected with constructs deficient in the REST DNA binding element, or when REST was down-regulated in retinoic acid-differentiated cells. IGF-I activates the signal transducer and activator of transcription-3 (STAT3) signaling pathway and this transcription factor, binding to the STAT1/3 DNA element located in the promoter, increases OPRM1 transcription. T-cell receptor (TCR) recognizes peptide antigens displayed in the context of the major histocompatibility complex (MHC) and gives rise to a potent as well as branched intracellular signalling that convert naïve T-cells in mature effectors, thus significantly contributing to the genesis of a specific immune response. In the second part of my work I exposed wild type Jurkat CD4+ T-cells to a mixture of CD3 and CD28 antigens in order to fully activate TCR and study whether its signalling influence OPRM1 expression. Results were that TCR engagement determined a significant induction of OPRM1 expression through the activation of transcription factors AP-1, NF-kB and NFAT. Eventually, I investigated MOPr turnover once it has been expressed on T-cells outer membrane. It turned out that DAMGO induced MOPr internalisation and recycling, whereas morphine did not. Overall, from the data collected in this thesis we can conclude that that a reduction in REST is a critical switch enabling IGF-I to up-regulate human MOPr, helping these findings clarify how human MOPr expression is regulated in neuronal cells, and that TCR engagement up-regulates OPRM1 transcription in T-cells. My results that neurotrophic factors a and TCR engagement, as well as it is reported for cytokines, seem to up-regulate OPRM1 in both neurons and immune cells suggest an important role for MOPr as a molecular bridge between neurons and immune cells; therefore, MOPr could play a key role in the cross-talk between immune system and nervous system and in particular in the balance between pro-inflammatory and pro-nociceptive stimuli and analgesic and neuroprotective effects.

New approaches to open problems in gene expression microarray data

In the past decade, the advent of efficient genome sequencing tools and high-throughput experimental biotechnology has lead to enormous progress in the life science. Among the most important innovations is the microarray tecnology. It allows to quantify the expression for thousands of genes simultaneously by measurin the hybridization from a tissue of interest to probes on a small glass or plastic slide. The characteristics of these data include a fair amount of random noise, a predictor dimension in the thousand, and a sample noise in the dozens. One of the most exciting areas to which microarray technology has been applied is the challenge of deciphering complex disease such as cancer. In these studies, samples are taken from two or more groups of individuals with heterogeneous phenotypes, pathologies, or clinical outcomes. these samples are hybridized to microarrays in an effort to find a small number of genes which are strongly correlated with the group of individuals. Eventhough today methods to analyse the data are welle developed and close to reach a standard organization (through the effort of preposed International project like Microarray Gene Expression Data -MGED- Society [1]) it is not unfrequant to stumble in a clinician's question that do not have a compelling statistical method that could permit to answer it.The contribution of this dissertation in deciphering disease regards the development of new approaches aiming at handle open problems posed by clinicians in handle specific experimental designs. In Chapter 1 starting from a biological necessary introduction, we revise the microarray tecnologies and all the important steps that involve an experiment from the production of the array, to the quality controls ending with preprocessing steps that will be used into the data analysis in the rest of the dissertation. While in Chapter 2 a critical review of standard analysis methods are provided stressing most of problems that In Chapter 3 is introduced a method to adress the issue of unbalanced design of miacroarray experiments. In microarray experiments, experimental design is a crucial starting-point for obtaining reasonable results. In a two-class problem, an equal or similar number of samples it should be collected between the two classes. However in some cases, e.g. rare pathologies, the approach to be taken is less evident. We propose to address this issue by applying a modified version of SAM [2]. MultiSAM consists in a reiterated application of a SAM analysis, comparing the less populated class (LPC) with 1,000 random samplings of the same size from the more populated class (MPC) A list of the differentially expressed genes is generated for each SAM application. After 1,000 reiterations, each single probe given a "score" ranging from 0 to 1,000 based on its recurrence in the 1,000 lists as differentially expressed. The performance of MultiSAM was compared to the performance of SAM and LIMMA [3] over two simulated data sets via beta and exponential distribution. The results of all three algorithms over low- noise data sets seems acceptable However, on a real unbalanced two-channel data set reagardin Chronic Lymphocitic Leukemia, LIMMA finds no significant probe, SAM finds 23 significantly changed probes but cannot separate the two classes, while MultiSAM finds 122 probes with score >300 and separates the data into two clusters by hierarchical clustering. We also report extra-assay validation in terms of differentially expressed genes Although standard algorithms perform well over low-noise simulated data sets, multi-SAM seems to be the only one able to reveal subtle differences in gene expression profiles on real unbalanced data. In Chapter 4 a method to adress similarities evaluation in a three-class prblem by means of Relevance Vector Machine [4] is described. In fact, looking at microarray data in a prognostic and diagnostic clinical framework, not only differences could have a crucial role. In some cases similarities can give useful and, sometimes even more, important information. The goal, given three classes, could be to establish, with a certain level of confidence, if the third one is similar to the first or the second one. In this work we show that Relevance Vector Machine (RVM) [2] could be a possible solutions to the limitation of standard supervised classification. In fact, RVM offers many advantages compared, for example, with his well-known precursor (Support Vector Machine - SVM [3]). Among these advantages, the estimate of posterior probability of class membership represents a key feature to address the similarity issue. This is a highly important, but often overlooked, option of any practical pattern recognition system. We focused on Tumor-Grade-three-class problem, so we have 67 samples of grade I (G1), 54 samples of grade 3 (G3) and 100 samples of grade 2 (G2). The goal is to find a model able to separate G1 from G3, then evaluate the third class G2 as test-set to obtain the probability for samples of G2 to be member of class G1 or class G3. The analysis showed that breast cancer samples of grade II have a molecular profile more similar to breast cancer samples of grade I. Looking at the literature this result have been guessed, but no measure of significance was gived before.

Morphogenesis of follicular epithelium in Drosophila melanogaster: function of von Hippel-Lindau tumour suppressor gene

During my PhD I have been involved in several projects regarding the morphogenesis of the follicular epithelium, such as the analysis of the pathways that correlate follicular epithelium patterning and eggshell genes expression. Moreover, I used the follicular epithelium as a model system to analyze the function of the Drosophila homolog of the human von Hippel-Lindau (d-VHL) during oogenesis, in order to gain insight into the role of h-VHL for the pathogenesis of VHL disease. h-VHL is implicated in a variety of processes and there is now a greater appreciation of HIF-independent h-VHL functions that are relevant to tumour development, including maintenance and organization of the primary cilium, maintenance of the differentiated phenotype in renal cells and regulation of epithelial-mesenchymal transition. However, the function of h-VHL gene during development has not been fully understood. It was previously shown that d-VHL down-regulates the motility of tubular epithelial cells (tracheal cells) during embryogenesis. Epithelial morphogenesis is important for organogenesis and pivotal for carcinogenesis, but mechanisms that control it are poorly understood. The Drosophila follicular epithelium is a genetically tractable model to understand these mechanisms in vivo. Therefore, to examine whether d-VHL has a role in epithelial morphogenesis and maintenance, I performed genetic and molecular analyses by using in vivo and in vitro approaches. From my analysis, I determined that d-VHL binds to and stabilizes microtubules. Loss of d-VHL depolymerizes the microtubule network during oogenesis, leading to a possible deregulation in the subcellular trafficking transport of polarity markers from Golgi apparatus to the different domains in which follicle cells are divided. The analysis carried out has allowed to establish a significant role of d-VHL in the maintenance of the follicular epithelium integrity.

Molecular bases, pathogenic mechanisms and possible therapeutic approach in Leber's Hereditary Optic Neuropathy

Leber’s hereditary optic neuropathy (LHON) is a mitochondrial disease characterized by a rapid loss of central vision and optic atrophy, due to the selective degeneration of retinal ganglion cells. The age of onset is around 20, and the degenerative process is fast and usually the second eye becomes affected in weeks or months. Even if this pathology is well known and has been well characterized, there are still open questions on its pathophysiology, such as the male prevalence, the incomplete penetrance and the tissue selectivity. This maternally inherited disease is caused by mutations in mitochondrial encoded genes of NADH ubiquinone oxidoreductase (complex I) of the respiratory chain. The 90% of LHON cases are caused by one of the three common mitochondrial DNA mutations (11778/ND4, 14484/ND6 and 3460/ND1) and the remaining 10% is caused by rare pathogenic mutations, reported in literature in one or few families. Moreover, there is also a small subset of patients reported with new putative pathogenic nucleotide changes, which awaits to be confirmed. We here clarify some molecular aspects of LHON, mainly the incomplete penetrance and the role of rare mtDNA mutations or variants on LHON expression, and attempt a possible therapeutic approach using the cybrids cell model. We generated novel structural models for mitochondrial encoded complex I subunits and a conservation analysis and pathogenicity prediction have been carried out for LHON reported mutations. This in-silico approach allowed us to locate LHON pathogenic mutations in defined and conserved protein domains and can be a useful tool in the analysis of novel mtDNA variants with unclear pathogenic/functional role. Four rare LHON pathogenic mutations have been identified, confirming that the ND1 and ND6 genes are mutational hot spots for LHON. All mutations were previously described at least once and we validated their pathogenic role, suggesting the need for their screening in LHON diagnostic protocols. Two novel mtDNA variants with a possible pathogenic role have been also identified in two independent branches of a large pedigree. Functional studies are necessary to define their contribution to LHON in this family. It also been demonstrated that the combination of mtDNA rare polymorphic variants is relevant in determining the maternal recurrence of myoclonus in unrelated LHON pedigrees. Thus, we suggest that particular mtDNA backgrounds and /or the presence of specific rare mutations may increase the pathogenic potential of the primary LHON mutations, thereby giving rise to the extraocular clinical features characteristic of the LHON “plus” phenotype. We identified the first molecular parameter that clearly discriminates LHON affected individuals from asymptomatic carriers, the mtDNA copy number. This provides a valuable mechanism for future investigations on variable penetrance in LHON. However, the increased mtDNA content in LHON individuals was not correlated to the functional polymorphism G1444A of PGC-1 alpha, the master regulator of mitochondrial biogenesis, but may be due to gene expression of genes involved in this signaling pathway, such as PGC-1 alpha/beta and Tfam. Future studies will be necessary to identify the biochemical effects of rare pathogenic mutations and to validate the novel candidate mutations here described, in terms of cellular bioenergetic characterization of these variants. Moreover, we were not able to induce mitochondrial biogenesis in cybrids cell lines using bezafibrate. However, other cell line models are available, such as fibroblasts harboring LHON mutations, or other approaches can be used to trigger the mitochondrial biogenesis.

Myc-mediated control of gene transcription in cancer cells

The Myc oncoproteins belong to a family of transcription factors composed by Myc, N-Myc and L-Myc. The most studied components of this family are Myc and N-Myc because their expressions are frequently deregulated in a wide range of cancers. These oncoproteins can act both as activators or repressors of gene transcription. As activators, they heterodimerize with Max (Myc associated X-factor) and the heterodimer recognizes and binds a specific sequence elements (E-Box) onto gene promoters recruiting histone acetylase and inducing transcriptional activation. Myc-mediated transcriptional repression is a quite debated issue. One of the first mechanisms defined for the Myc-mediated transcriptional repression consisted in the interaction of Myc-Max complex Sp1 and/or Miz1 transcription factors already bound to gene promoters. This interaction may interfere with their activation functions by recruiting co-repressors such as Dnmt3 or HDACs. Moreover, in the absence of , Myc may interfere with the Sp1 activation function by direct interaction and subsequent recruitment of HDACs. More recently the Myc/Max complex was also shown to mediate transcriptional repression by direct binding to peculiar E-box. In this study we analyzed the role of Myc overexpression in Osteosarcoma and Neuroblastoma oncogenesis and the mechanisms underling to Myc function. Myc overexpression is known to correlate with chemoresistance in Osteosarcoma cells. We extended this study by demonstrating that c-Myc induces transcription of a panel of ABC drug transporter genes. ABCs are a large family trans-membrane transporter deeply involved in multi drug resistance. Furthermore expression levels of Myc, ABCC1, ABCC4 and ABCF1 were proved to be important prognostic tool to predict conventional therapy failure. N-Myc amplification/overexpression is the most important prognostic factor for Neuroblastoma. Cyclin G2 and Clusterin are two genes often down regulated in neuroblastoma cells. Cyclin G2 is an atypical member of Cyclin family and its expression is associated with terminal differentiation and apoptosis. Moreover it blocks cell cycle progression and induces cell growth arrest. Instead, CLU is a multifunctional protein involved in many physiological and pathological processes. Several lines of evidences support the view that CLU may act as a tumour suppressor in Neuroblastoma. In this thesis I showed that N-Myc represses CCNG2 and CLU transcription by different mechanisms. • N-Myc represses CCNG2 transcription by directly interacting with Sp1 bound in CCNG2 promoter and recruiting HDAC2. Importantly, reactivation of CCNG2 expression through epigenetic drugs partially reduces N-Myc and HDAC2 mediated cell proliferation. • N-Myc/Max complex represses CLU expression by direct binding to a peculiar E-box element on CLU promoter and by recruitment of HDACs and Polycomb Complexes, to the CLU promoter. Overall our findings strongly support the model in which Myc overexpression/amplification may contribute to some aspects of oncogenesis by a dual action: i) transcription activation of genes that confer a multidrug resistant phenotype to cancer cells; ii), transcription repression of genes involved in cell cycle inhibition and cellular differentiation.

Genetic Manipulation, Gene Silencing and Biomarker Development in Multiple Experimental Models.

The thesis is set in three different parts, according to the relative experimental models. First, the domestic pig (Sus scrofa) is part of the study on reproductive biotechnologies: the transgenesis technique of Sperm Mediated Gene Transfer is widely studied starting from the quality of the semen, through the study of multiple uptakes of exogenous DNA and lastly used in the production of multi-transgenic blastocysts. Finally we managed to couple the transgenesis pipeline with sperm sorting and therefore produced transgenic embryos of predetermined sex. In the second part of the thesis the attention is on the fruit fly (Drosophila melanogaster) and on its derived cell line: the S2 cells. The in vitro and in vivo models are used to develop and validate an efficient way to knock down the myc gene. First an efficient in vitro protocol is described, than we demonstrate how the decrease in myc transcript remarkably affects the ribosome biogenesis through the study of Polysome gradients, rRNA content and qPCR. In vivo we identified two optimal drivers for the conditional silencing of myc, once the flies are fed with RU486: the first one is throughout the whole body (Tubulin), while the second is a head fat body driver (S32). With these results we present a very efficient model to study the role of myc in multiple aspects of translation. In the third and last part, the focus is on human derived lung fibroblasts (hLF-1), mouse tail fibroblasts and mouse tissues. We developed an efficient assay to quantify the total protein content of the nucleus on a single cell level via fluorescence. We coupled the protocol with classical immunofluorescence so to have at the same time general and particular information, demonstrating that during senescence nuclear proteins increase by 1.8 fold either in human cells, mouse cells and mouse tissues.

Stochastic models and dynamic measures for the characterization of bistable circuits

During the last few years, a great deal of interest has risen concerning the applications of stochastic methods to several biochemical and biological phenomena. Phenomena like gene expression, cellular memory, bet-hedging strategy in bacterial growth and many others, cannot be described by continuous stochastic models due to their intrinsic discreteness and randomness. In this thesis I have used the Chemical Master Equation (CME) technique to modelize some feedback cycles and analyzing their properties, including experimental data. In the first part of this work, the effect of stochastic stability is discussed on a toy model of the genetic switch that triggers the cellular division, which malfunctioning is known to be one of the hallmarks of cancer. The second system I have worked on is the so-called futile cycle, a closed cycle of two enzymatic reactions that adds and removes a chemical compound, called phosphate group, to a specific substrate. I have thus investigated how adding noise to the enzyme (that is usually in the order of few hundred molecules) modifies the probability of observing a specific number of phosphorylated substrate molecules, and confirmed theoretical predictions with numerical simulations. In the third part the results of the study of a chain of multiple phosphorylation-dephosphorylation cycles will be presented. We will discuss an approximation method for the exact solution in the bidimensional case and the relationship that this method has with the thermodynamic properties of the system, which is an open system far from equilibrium.In the last section the agreement between the theoretical prediction of the total protein quantity in a mouse cells population and the observed quantity will be shown, measured via fluorescence microscopy.