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

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

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.

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

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

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.

Establishing Zebrafish as a Novel Exercise Model: Swimming Economy, Swimming-enhanced Growth and Muscle Growth Marker Gene Expression

Zebrafish has been largely accepted as a vertebrate multidisciplinary model but its usefulness as a model for exercise physiology has been hampered by the scarce knowledge on its swimming economy, optimal swimming speeds and cost of transport. Therefore, we have performed individual and group-wise swimming experiments to quantify swimming economy and to demonstrate the exercise effects on growth in adult zebrafish.

Cytoplasmic accumulation of NCoR in malignant melanoma: Consequences of altered gene repression and prognostic significance

Invasive malignant melanoma (MM) is an aggressive tumor with no curative therapy available in advanced stages. Nuclear corepressor (NCoR) is an essential regulator of gene transcription, and its function has been found deregulated in different types of cancer. In colorectal cancer cells, loss of nuclear NCoR is induced by Inhibitor of kappa B kinase (IKK) through the phosphorylation of specific serine residues. We here investigate whether NCoR function impacts in MM, which might have important diagnostic and prognostic significance. By IHC, we here determined the subcellular distribution of NCoR in a cohort of 63 primary invasive MM samples, and analyzed its possible correlation with specific clinical parameters. We therefore used a microarray-based strategy to determine global gene expression differences in samples with similar tumor stage, which differ in the presence of cytoplasmic or nuclear NCoR. We found that loss of nuclear NCoR results in upregulation of a specific cancer-related genetic signature, and is significantly associated with MM progression. Inhibition of IKK activity in melanoma cells reverts NCoR nuclear distribution and specific NCoR-regulated gene transcription. Analysis of public database demonstrated that inactivating NCoR mutations are highly prevalent in MM, showing features of driver oncogene.

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.

Gene expression profile of ABC transporters and cytotoxic effect of ibuprofen and acetaminophen in an epithelial ovarian cancer cell line in vitro

PURPOSES: To determine the basic expression of ABC transporters in an epithelial ovarian cancer cell line, and to investigate whether low concentrations of acetaminophen and ibuprofen inhibited the growth of this cell line in vitro. METHODS: TOV-21 G cells were exposed to different concentrations of acetaminophen (1.5 to 15 μg/mL) and ibuprofen (2.0 to 20 μg/mL) for 24 to 48 hours. The cellular growth was assessed using a cell viability assay. Cellular morphology was determined by fluorescence microscopy. The gene expression profile of ABC transporters was determined by assessing a panel including 42 genes of the ABC transporter superfamily. RESULTS: We observed a significant decrease in TOV-21 G cell growth after exposure to 15 μg/mL of acetaminophen for 24 (p=0.02) and 48 hours (p=0.01), or to 20 μg/mL of ibuprofen for 48 hours (p=0.04). Assessing the morphology of TOV-21 G cells did not reveal evidence of extensive apoptosis. TOV-21 G cells had a reduced expression of the genes ABCA1, ABCC3, ABCC4, ABCD3, ABCD4 and ABCE1 within the ABC transporter superfamily. CONCLUSIONS: This study provides in vitro evidence of inhibitory effects of growth in therapeutic concentrations of acetaminophen and ibuprofen on TOV-21 G cells. Additionally, TOV-21 G cells presented a reduced expression of the ABCA1, ABCC3, ABCC4, ABCD3, ABCD4 and ABCE1 transporters.

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.

Gene Regulation in The Human Immune System. Gene Expression Signatures of Th2 Cell Differentiation, Type 1 Diabetes, and Intrauterine Immune Adaptation

The human immune system is constantly interacting with the surrounding stimuli and microorganisms. However, when directed against self or harmless antigens, these vital defense mechanisms can cause great damage. In addition, the understanding the underlying mechanism of several human diseases caused by aberrant immune cell functions, for instance type 1 diabetes and allergies, remains far from being complete. In this Ph.D. study these questions were addressed using genome-wide transcriptomic analyses. Asthma and allergies are characterized by a hyperactive response of the T helper 2 (Th2) immune cells. In this study, the target genes of the STAT6 transcription factor in naïve human T cells were identified with RNAi for the first time. STAT6 was shown to act as a central activator of the genes expression upon IL-4 signaling, with both direct and indirect effects on Th2 cell transcriptome. The core transcription factor network induced by IL-4 was identified from a kinetic analysis of the transcriptome. Type 1 diabetes is an autoimmune disease influenced by both the genetic susceptibility of an individual and the disease-triggering environmental factors. To improve understanding of the autoimmune processes driving pathogenesis in the prediabetic phase in humans, a unique series of prospective whole-blood RNA samples collected from HLA-susceptible children in the Finnish Type 1 Diabetes Prediction and Prevention (DIPP) study was studied. Changes in different timewindows of the pathogenesis process were identified, and especially the type 1 interferon response was activated early and throughout the preclinical T1D. The hygiene hypothesis states that allergic diseases, and lately also autoimmune diseases, could be prevented by infections and other microbial contacts acquired in early childhood, or even prenatally. To study the effects of the standard of hygiene on the development of neonatal immune system, cord blood samples from children born in Finland (high standard of living), Estonia (rapid economic growth) and Russian Karelia (low standard of living) were compared. Children born in Russian Karelia deviated from Finnish and Estonian children in many aspects of the neonatal immune system, which was developmentally more mature in Karelia, resembling that of older infants. The results of this thesis offer significant new information on the regulatory networks associated with immune-mediated diseases in human. The results will facilitate understanding and further research on the role of the identified target genes and mechanisms driving the allergic inflammation and type 1 diabetes, hopefully leading to a new era of drug development.