Functional impact of genetic variation on gene expression

Numerous links between genetic variants and phenotypes are known and genome-wide association studies dramatically increased the number of genetic variants associated with traits during the last decade. However, how changes in the DNA perturb the molecular mechanisms and impact on the phenotype of an organism remains elusive. Studies suggest that many traitassociated variants are in the non-coding region of the genome and probably act through regulation of gene expression. During my thesis I investigated how genetic variants affect gene expression through gene regulatory mechanisms. The first chapter was a collaborative project with a pharmaceutical company, where we investigated genome-wide copy number variation (CNVs) among Cynomolgus monkeys (Macaca fascicularis) used in pharmaceutical studies, and associated them to changes in gene expression. We found substantial copy number variation and identified CNVs linked to tissue-specific expression changes of proximal genes. The second and third chapters focus on genetic variation in humans and its effects on gene regulatory mechanisms and gene expression. The second chapter studies two human trios, where the allelic effects of genetic variation on genome-wide gene expression, protein-DNA binding and chromatin modifications were investigated. We found abundant allele specific activity across all measured molecular phenotypes and show extended coordinated behavior among them. In the third chapter, we investigated the impact of genetic variation on these phenotypes in 47 unrelated individuals. We found that chromatin phenotypes are organized into local variable modules, often linked to genetic variation and gene expression. Our results suggest that chromatin variation emerges as a result of perturbations of cis-regulatory elements by genetic variants, leading to gene expression changes. The work of this thesis provides novel insights into how genetic variation impacts gene expression by perturbing regulatory mechanisms. -- De nombreux liens entre variations génétiques et phénotypes sont connus. Les études d'association pangénomique ont considérablement permis d'augmenter le nombre de variations génétiques associées à des phénotypes au cours de la dernière décennie. Cependant, comprendre comment ces changements perturbent les mécanismes moléculaires et affectent le phénotype d'un organisme nous échappe encore. Des études suggèrent que de nombreuses variations, associées à des phénotypes, sont situées dans les régions non codantes du génome et sont susceptibles d'agir en modifiant la régulation d'expression des gènes. Au cours de ma thèse, j'ai étudié comment les variations génétiques affectent les niveaux d'expression des gènes en perturbant les mécanismes de régulation de leur expression. Le travail présenté dans le premier chapitre est un projet en collaboration avec une société pharmaceutique. Nous avons étudié les variations en nombre de copies (CNV) présentes chez le macaque crabier (Macaca fascicularis) qui est utilisé dans les études pharmaceutiques, et nous les avons associées avec des changements d'expression des gènes. Nous avons découvert qu'il existe une variabilité substantielle du nombre de copies et nous avons identifié des CNVs liées aux changements d'expression des gènes situés dans leur voisinage. Ces associations sont présentes ou absentes de manière spécifique dans certains tissus. Les deuxième et troisième chapitres se concentrent sur les variations génétiques dans les populations humaines et leurs effets sur les mécanismes de régulation des gènes et leur expression. Le premier se penche sur deux trios humains, père, mère, enfant, au sein duquel nous avons étudié les effets alléliques des variations génétiques sur l'expression des gènes, les liaisons protéine-ADN et les modifications de la chromatine. Nous avons découvert que l'activité spécifique des allèles est abondante abonde dans tous ces phénotypes moléculaires et nous avons démontré que ces derniers ont un comportement coordonné entre eux. Dans le second, nous avons examiné l'impact des variations génétiques de ces phénotypes moléculaires chez 47 individus, sans lien de parenté. Nous avons observé que les phénotypes de la chromatine sont organisés en modules locaux, qui sont liés aux variations génétiques et à l'expression des gènes. Nos résultats suggèrent que la variabilité de la chromatine est due à des variations génétiques qui perturbent des éléments cis-régulateurs, et peut conduire à des changements dans l'expression des gènes. Le travail présenté dans cette thèse fournit de nouvelles pistes pour comprendre l'impact des différentes variations génétiques sur l'expression des gènes à travers les mécanismes de régulation.

Developmental and tissue-specific involvement of peroxisome proliferator-activated receptor-alpha in the control of mouse uncoupling protein-3 gene expression.

Uncoupling protein-3 (UCP3) is a member of the mitochondrial carrier family expressed preferentially in skeletal muscle and heart. It appears to be involved in metabolic handling of fatty acids in a way that minimizes excessive production of reactive oxygen species. Fatty acids are powerful regulators of UCP3 gene transcription. We have found that the role of peroxisome proliferator-activated receptor-α (PPARα) on the control of UCP3 gene expression depends on the tissue and developmental stage. In adults, UCP3 mRNA expression is unaltered in skeletal muscle from PPARα-null mice both in basal conditions and under the stimulus of starvation. In contrast, UCP3 mRNA is down-regulated in adult heart both in fed and fasted PPARα-null mice. This occurs despite the increased levels of free fatty acids caused by fasting in PPARα-null mice. In neonates, PPARα-null mice show impaired UCP3 mRNA expression in skeletal muscle in response to milk intake, and this is not a result of reduced free fatty acid levels. The murine UCP3 promoter is activated by fatty acids through either PPARα or PPARδ but not by PPARγ or retinoid X receptor alone. PPARδ-dependent activation could be a potential compensatory mechanism to ensure appropriate expression of UCP3 gene in adult skeletal muscle in the absence of PPARα. However, among transcripts from other PPARα and PPARδ target genes, only those acutely induced by milk intake in wild-type neonates were altered in muscle or heart from PPARα-null neonates. Thus, PPARα-dependent regulation is required for appropriate gene regulation of UCP3 as part of the subset of fatty-acid-responsive genes in neonatal muscle and heart.

Individual differences in gene expression : insights from transcriptional control of the CSL gene and from human population studies

CSL is a key transcription factor, mostly acting as a repressor, which has been shown to have a highly context-dependent function. While known as the main effector of Notch signaling, it can also exert Notch-independent functions. The downstream effects of the Notch/CSL signaling pathway and its involvement in several biological processes have been intensively studied. We recently showed that CSL is important to maintain skin homeostasis, as its specific deletion in mouse dermal fibroblasts -or downmodulation in human stromal fibroblasts- creates an inducing environment for squamous cell carcinoma (SCC) development, possibly due to the conversion of stromal fibroblasts into cancer associated fibroblasts (CAFs). Despite the wide interest in CSL as a transcriptional regulator, the mechanism of its own regulation has so far been neglected. We show here that CSL expression levels differ between individuals, and correlate among others with genes involved in DNA damage response. Starting from this finding we show that in dermal fibroblasts CSL is under transcriptional control of stress inducers such as UVA irradiation and Reactive Oxygen Species (ROS) induction, and that a main player in CSL transcriptional regulation is the transcription factor p53. In a separate line of work, we focused on individual variability, studying the differences in gene expression between human populations in various cancer types, particularly focusing on the Caucasian and African populations. It is indeed widely known that these populations have different incidences and mortalities for various cancers, and response to cancer treatment may also vary between them. We show here several genes that are differentially expressed and could be of interest in the study of population differences in cancer. -- CSL est un facteur de transcription agissant essentiellement comme répresseur, et qui a une fonction hautement dépendant du contexte. C'est l'effecteur principal de la voie de signalisation de Notch, mais il peut également exercer ses fonctions dans une façon Notch- indépendante. Nous avons récemment montré que CSL est important pour maintenir l'homéostasie de la peau. Sa suppression spécifique dans les fibroblastes dermiques de la souris ou dans les fibroblastes stromales humaines crée un environnement favorable pour le développement du carcinome épidermoïde (SCC), probablement en raison de la conversion des fibroblastes en fibroblastes associé au cancer (CAF). Malgré le grand intérêt de CSL comme régulateur transcriptionnel, le mécanisme de sa propre régulation a été jusqu'ici négligée. Nous montrons ici que dans les fibroblastes dermiques CSL est sous le contrôle transcriptionnel de facteurs de stress tels que l'irradiation UVA et l'induction des ROS dont p53 est l'acteur principal de cette régulation. Nous montrons aussi que les niveaux d'expression de CSL varient selon les individus, en corrélation avec d'autres gènes impliqués dans la réponse aux dommages de l'ADN. Dans une autre axe de recherche, concernant la variabilité individuelle, nous avons étudié les différences dans l'expression des gènes dans différents types de cancer entre les populations humaines, en se concentrant particulièrement sur les populations africaines et caucasiennes. Il est en effet bien connu que ces populations montrent des variations dans l'incidence des cancers, la mortalité, ainsi que pour les réponses au traitement. Nous montrons ici plusieurs gènes qui sont exprimés différemment et pourraient être digne d'intérêt dans l'étude du cancer au sein de différentes populations.

Identification of molecular subtypes and gene expression patterns of breast cancer analysing RNA-seq data

Breast cancer is the most common diagnosed cancer and the leading cause of cancer death among females worldwide. It is considered a highly heterogeneous disease and it must be classified into more homogeneous groups. Hence, the purpose of this study was to classify breast tumors based on variations in gene expression patterns derived from RNA sequencing by using different class discovery methods. 42 breast tumors paired-samples were sequenced by Illumine Genome Analyzer and the data was analyzed and prepared by TopHat2 and htseq-count. As reported previously, breast cancer could be grouped into five main groups known as basal epithelial-like group, HER2 group, normal breast-like group and two Luminal groups with a distinctive expression profile. Classifying breast tumor samples by using PAM50 method, the most common subtype was Luminal B and was significantly associated with ESR1 and ERBB2 high expression. Luminal A subtype had ESR1 and SLC39A6 significant high expression, whereas HER2 subtype had a high expression of ERBB2 and CNNE1 genes and low luminal epithelial gene expression. Basal-like and normal-like subtypes were associated with low expression of ESR1, PgR and HER2, and had significant high expression of cytokeratins 5 and 17. Our results were similar compared with TGCA breast cancer data results and with known studies related with breast cancer classification. Classifying breast tumors could add significant prognostic and predictive information to standard parameters, and moreover, identify marker genes for each subtype to find a better therapy for patients with breast cancer.

Regulation of B Cell Gene Expression and Function by Ikaros, Helios and Bcl6

B lymphocytes constitute a key branch of adaptive immunity by providing specificity to recognize a vast variety of antigens by B cell antigen receptors (BCR) and secreted antibodies. Antigen recognition activates the cells and can produce antibody secreting plasma cells via germinal center reaction that leads to the maturation of antigen recognition affinity and switching of antibody effector class. The specificity of antigen recognition is achieved through a multistep developmental pathway that is organized by interplay of transcription factors and signals through BCR. Lymphoid malignancies arise from different stages of development in abnormal function of transcriptional regulation. To understand the B cell development and the function of B cells, a thorough understanding of the regulation of gene expression is important. The transcription factors of the Ikaros family and Bcl6 are frequently associated with lymphoma generation. The aim of this study was to reveal the targets of Ikaros, Helios and Bcl6 mediated gene regulation and to find out the function of Ikaros and Helios in B cells. This study uses gene targeted DT40 B cell lines and establishes a role for Ikaros family factors Ikaros and Helios in the regulation of BCR signaling that is important at developmental checkpoints, for cell survival and in activation. Ikaros and Helios had opposing roles in the regulation of BCR signals. Ikaros was found to directly repress the SHIP gene that encodes a signaling lipid-metabolizing enzyme, whereas Helios had activating effect on SHIP expression. The findings demonstrate a balancing function for these two Ikaros family transcription factors in the regulation of BCR signaling as well as in the regulation of gene expression. Bcl6 was found to repress plasma cell gene expression program while maintaining gene expression profile of B cells. Analysis of direct Bcl6 target genes suggested novel mechanisms for Bcl6-mediated suppression of plasma cell differentiation and promoting germinal center phenotype.

Cytokine gene expression and molecular detection of Mycobacterium avium subspecies paratuberculosisin organs of experimentally infected mice

AbstractMycobacterium avium subspecies paratuberculosis (MAP) can infect ruminants and remain subclinical for long periods within herds. The identification of organs that are more susceptible to infection and the evaluation of cytokine expression at the site of infection are important to understand the pathogenesis of MAP. In this study, the probability of detection of MAP-DNA and the expression of cytokines in organs of C57BL/6 mice infected intraperitoneally for 120 days were evaluated. Among the evaluated organs, the spleen (85%), colon (75%) and liver (60%) had the highest frequency of positivity. When compared these frequencies between organs, it has been found that the spleen had 1.54 times as likely to be positive in relation to the ileum, and 2.0 times more likely in relation to the Peyer's patches. In addition, at 60 days post-infection, the spleen and the liver were responsible for upregulation of IFN-γ , and the ileum by TNF-α and IL-4. The results indicate that the spleen is the best organ for evaluating an experimental infection by MAP, especially in the initial stages of the infection. Moreover, it showed that the spleen, liver and ileum have a direct role in the inflammatory response in experimental models.

Effects of Cinnamon extract on biochemical enzymes, TNF-α and NF-κB gene expression levels in liver of broiler chickens inoculated with Escherichia coli

Abstract: Infection with Escherichia coli (E. coli) is a common disease in poultry industry. The use of antibiotics to treat diseases is facing serious criticism and concerns. The medicinal plants may be effective alternatives because of their multiplex activities. The aim of this study was to investigate the effects of cinnamon extract on the levels of liver enzymes, tumor necrosis factor-alpha (TNF-α) and nuclear factor-kappa B (NF-κB) gene expressions in liver of broiler chickens infected with E. coli. Ninety Ross-308 broilers were divided into healthy or E. coli-infected groups, receiving normal or cinnamon extract (in concentrations of 100 or 200mg/kg of food) supplemented diets. E. coli suspension (108cfu) was injected subcutaneously after 12 days cinnamon administration. Seventy-two hours after E. coli injection, the blood samples were taken for biochemical analysis of liver enzymes in serum (spectrophotometrically), and liver tissue samples were obtained for detection of gene expression of inflammatory markers TNF-α and NF-κB, using real-time PCR. Infection with E. coli significantly increased the levels of TNF-α and NF-κB gene expressions as well as some liver enzymes including creatine-kinase (CK), lactate-dehydrogenase (LDH), alanine-transferase (ALT) and aspartate-transferase (AST) as compared with control group (P<0.05). Pre-administration of cinnamon extract in broilers diet (in both concentrations) significantly reduced the tissue levels of TNF-α and NF-κB gene expressions and enzymes CK and ALT in serum of broiler chickens inoculated with E. coli in comparison with E. coli group (P<0.05 and P<0.01). The levels of LDH and AST were significantly decreased only by 200mg/kg cinnamon extract in infected broilers. The level of alkaline-phosphatase (ALP) was not affected in any groups. Pre-administration of cinnamon extract in diets of broiler chickens inoculated with E. coli could significantly reduce the gene expression levels of pro-inflammatory mediators and liver enzymes activities, thereby protecting the liver against this pathologic condition.

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.

Control of gene expression in Trypanosomatidae

The study of mechanisms which control gene expression in trypanosomatids has developed at an increasing rate since 1989 when the first successful DNA transfection experiments were reported. Using primarily Trypanosoma brucei as a model, several groups have begun to elucidate the basic control mechanisms and to define the cellular factors involved in mRNA transcription, processing and translation in these parasites. This review focuses on the most recent studies regarding a subset of genes that are expressed differentially during the life cycle of three groups of parasites. In addition to T. brucei, I will address studies on gene regulation in a few species of Leishmania and the results obtained by a much more limited group of laboratories studying gene expression in Trypanosoma cruzi. It is becoming evident that the regulatory strategies chosen by different species of trypanosomatids are not similar, and that for these very successful parasites it is probably advantageous to employ multiple mechanisms simultaneously. In addition, with the increasing numbers of parasite genes that have now been submitted to molecular dissection, it is also becoming evident that, among the various strategies for gene expression control, there is a predominance of regulatory pathways acting at the post-transcriptional level.

Pituitary glycoprotein hormone a-subunit secretion by cirrhotic patients

Secretion of the a-subunit of pituitary glycoprotein hormones usually follows the secretion of intact gonadotropins and is increased in gonadal failure and decreased in isolated gonadotropin deficiency. The aim of the present study was to determine the levels of the a-subunit in the serum of patients with cirrhosis of the liver and to compare the results obtained for eugonadal cirrhotic patients with those obtained for cirrhotic patients with hypogonadotropic hypogonadism. Forty-seven of 63 patients with cirrhosis (74.6%) presented hypogonadism (which was central in 45 cases and primary in 2), 7 were eugonadal, and 9 women were in normal menopause. The serum a-subunit was measured by the fluorimetric method using monoclonal antibodies. Cross-reactivity with LH, TSH, FSH and hCG was 6.5, 1.2, 4.3 and 1.1%, respectively, with an intra-assay coefficient of variation (CV) of less than 5% and an interassay CV of 5%, and sensitivity limit of 4 ng/l. The serum a-subunit concentration ranged from 36 to 6253 ng/l, with a median of 273 ng/l. The median was 251 ng/l for patients with central hypogonadism and 198 ng/l for eugonadal patients. The correlation between the a-subunit and basal LH levels was significant both in the total sample (r = 0.48, P<0.01) and in the cirrhotic patients with central hypogonadism (r = 0.33, P = 0.02). Among men with central hypogonadism there was a negative correlation between a-subunit levels and total testosterone levels (r = 0.54, P<0.01) as well as free testosterone levels (r = -0.53, P<0.01). In conclusion, although the a-subunit levels are correlated with LH levels, at present they cannot be used as markers for hypogonadism in patients with cirrhosis of the liver.

Control of CD4 gene expression: connecting signals to outcomes in T cell development

The control of CD4 gene expression is essential for proper T lymphocyte development. Signals transmitted from the T-cell antigen receptor (TCR) during the thymic selection processes are believed to be linked to the regulation of CD4 gene expression during specific stages of T cell development. Thus, a study of the factors that control CD4 gene expression may lead to further insight into the molecular mechanisms that drive thymic selection. In this review, we discuss the work conducted to date to identify and characterize the cis-acting transcriptional control elements in the CD4 locus and the DNA-binding factors that mediate their function. From these studies, it is becoming clear that the molecular mechanisms controlling CD4 gene expression are very complex and differ at each stage of development. Thus, the control of CD4 expression is subject to many different influences as the thymocyte develops.

The neuroimmune-endocrine axis: pathophysiological implications for the central nervous system cytokines and hypothalamus-pituitary-adrenal hormone dynamics

Cytokines are molecules that were initially discovered in the immune system as mediators of communication between various types of immune cells. However, it soon became evident that cytokines exert profound effects on key functions of the central nervous system, such as food intake, fever, neuroendocrine regulation, long-term potentiation, and behavior. In the 80's and 90's our group and others discovered that the genes encoding various cytokines and their receptors are expressed in vascular, glial, and neuronal structures of the adult brain. Most cytokines act through cell surface receptors that have one transmembrane domain and which transduce a signal through the JAK/STAT pathway. Of particular physiological and pathophysiological relevance is the fact that cytokines are potent regulators of hypothalamic neuropeptidergic systems that maintain neuroendocrine homeostasis and which regulate the body's response to stress. The mechanisms by which cytokine signaling affects the function of stress-related neuroendocrine systems are reviewed in this article.

Sm14 gene expression in different stages of the Schistosoma mansoni life cycle and immunolocalization of the Sm14 protein within the adult worm

Sm14 is a 14-kDa vaccine candidate antigen from Schistosoma mansoni that seems to be involved in cytoplasmic trafficking of fatty acids. Although schistosomes have a high requirement for lipids, they are not able to synthesize fatty acids and sterols de novo. Thus, they must acquire host lipids. In order to determine whether Sm14 is present in different stages of the life cycle of the parasite, we performed RT-PCR. Sm14 mRNA was identified in all stages of the life cycle studied, mainly schistosomulum, adult worm and egg. Additionally, we used a rabbit anti-Sm14 polyclonal antibody in an indirect immunofluorescence assay to localize Sm14 in adult worm sections. The basal lamella of the tegument and the gut epithelium were strongly labeled. These tissues have a high flow of and demand for lipids, a finding that supports the putative role of Sm14 as an intracellular transporter of fatty acids from host cells.