Construction of a global map of human gene expression : the process, tools and analysis

This thesis studies human gene expression space using high throughput gene expression data from DNA microarrays. In molecular biology, high throughput techniques allow numerical measurements of expression of tens of thousands of genes simultaneously. In a single study, this data is traditionally obtained from a limited number of sample types with a small number of replicates. For organism-wide analysis, this data has been largely unavailable and the global structure of human transcriptome has remained unknown. This thesis introduces a human transcriptome map of different biological entities and analysis of its general structure. The map is constructed from gene expression data from the two largest public microarray data repositories, GEO and ArrayExpress. The creation of this map contributed to the development of ArrayExpress by identifying and retrofitting the previously unusable and missing data and by improving the access to its data. It also contributed to creation of several new tools for microarray data manipulation and establishment of data exchange between GEO and ArrayExpress. The data integration for the global map required creation of a new large ontology of human cell types, disease states, organism parts and cell lines. The ontology was used in a new text mining and decision tree based method for automatic conversion of human readable free text microarray data annotations into categorised format. The data comparability and minimisation of the systematic measurement errors that are characteristic to each lab- oratory in this large cross-laboratories integrated dataset, was ensured by computation of a range of microarray data quality metrics and exclusion of incomparable data. The structure of a global map of human gene expression was then explored by principal component analysis and hierarchical clustering using heuristics and help from another purpose built sample ontology. A preface and motivation to the construction and analysis of a global map of human gene expression is given by analysis of two microarray datasets of human malignant melanoma. The analysis of these sets incorporate indirect comparison of statistical methods for finding differentially expressed genes and point to the need to study gene expression on a global level.

Cellular Membranes as a Playground for Semliki Forest Virus Replication Complex

All positive-strand RNA viruses utilize cellular membranes for the assembly of their replication complexes, which results in extensive membrane modification in infected host cells. These alterations act as structural and functional scaffolds for RNA replication, providing protection for the viral double-stranded RNA against host defences. It is known that different positive-strand RNA viruses alter different cellular membranes. However, the origin of the targeted membranes, the mechanisms that direct replication proteins to specific membranes and the steps in the formation of the membrane bound replication complex are not completely understood. Alphaviruses (including Semliki Forest virus, SFV), members of family Togaviridae, replicate their RNA in association with membranes derived from the endosomal and lysosomal compartment, inducing membrane invaginations called spherules. Spherule structures have been shown to be the specific sites for RNA synthesis. Four replication proteins, nsP1-nsP4, are translated as a polyprotein (P1234) which is processed autocatalytically and gives rise to a membrane-bound replication complex. Membrane binding is mediated via nsP1 which possesses an amphipathic α-helix (binding peptide) in the central region of the protein. The aim of this thesis was to characterize the association of the SFV replication complex with cellular membranes and the modification of the membranes during virus infection. Therefore, it was necessary to set up the system for determining which viral components are needed for inducing the spherules. In addition, the targeting of the replication complex, the formation site of the spherules and their intracellular trafficking were studied in detail. The results of current work demonstrate that mutations in the binding peptide region of nsP1 are lethal for virus replication and change the localization of the polyprotein precursor P123. The replication complex is first targeted to the plasma membrane where membrane invaginations, spherules, are induced. Using a specific regulated endocytosis event the spherules are internalized from the plasma membrane in neutral carrier vesicles and transported via an actin-and microtubule-dependent manner to the pericentriolar area. Homotypic fusions and fusions with pre-existing acidic organelles lead to the maturation of previously described cytopathic vacuoles with hundreds of spherules on their limiting membranes. This work provides new insights into the membrane binding mechanism of SFV replication complex and its role in the virus life cycle. Development of plasmid-driven system for studying the formation of the replication complex described in this thesis allows various applications to address different steps in SFV life cycle and virus-host interactions in the future. This trans-replication system could be applied for many different viruses. In addition, the current work brings up new aspects of membranes and cellular components involved in SFV replication leading to further understanding in the formation and dynamics of the membrane-associated replication complex.

Electron and proton transfer in NADH:ubiquinone oxidoreductase (Complex I) from Escherichia coli

Cells of every living organism on our planet − bacterium, plant or animal − are organized in such a way that despite differences in structure and function they utilize the same metabolic energy represented by electrochemical proton gradient across a membrane. This gradient of protons is generated by the series of membrane bound multisubunit proteins, Complex I, II, III and IV, organized in so-called respiratory or electron transport chain. In the eukaryotic cell it locates in the inner mitochondrial membrane while in the bacterial cell it locates in the cytoplasmic membrane. The function of the respiratory chain is to accept electrons from NADH and ubiquinol and transfer them to oxygen resulting in the formation of water. The free energy released upon these redox reactions is converted by respiratory enzymes into an electrochemical proton gradient, which is used for synthesis of ATP as well as for many other energy dependent processes. This thesis is focused on studies of the first member of the respiratory chain − NADH:ubiquinone oxidoreductase or Complex I. This enzyme has a boot-shape structure with hydrophilic and hydrophobic domains, the former of which has all redox groups of the protein, the flavin and eight to nine iron-sulfur clusters. Complex I serves as a proton pump coupling transfer of two electrons from NADH to ubiquinone to the translocation of four protons across the membrane. So far the mechanism of energy transduction by Complex I is unknown. In the present study we applied a set of different methods to study the electron and proton transfer reactions in Complex I from Escherichia coli. The main achievement was the experiment that showed that the electron transfer through the hydrophilic domain of Complex I is unlikely to be coupled to proton transfer directly or to conformational changes in the protein. In this work for the first time properties of all redox centers of Complex I were characterized in the intact purified bacterial enzyme. We also probed the role of several conserved amino acid residues in the electron transfer of Complex I. Finally, we found that highly conserved amino acid residues in several membrane subunits form a common pattern with a very prominent feature – the presence of a few lysines within the membrane. Based on the experimental data, we suggested a tentative principle which may govern the redox-coupled proton pumping in Complex I.

The range expansion of the European map butterfly in Finland

Climate change will influence the living conditions of all life on Earth. For some species the change in the environmental conditions that has occurred so far has already increased the risk of extinction, and the extinction risk is predicted to increase for large numbers of species in the future. Some species may have time to adapt to the changing environmental conditions, but the rate and magnitude of the change are too great to allow many species to survive via evolutionary changes. Species responses to climate change have been documented for some decades. Some groups of species, like many insects, respond readily to changes in temperature conditions and have shifted their distributions northwards to new climatically suitable regions. Such range shifts have been well documented especially in temperate zones. In this context, butterflies have been studied more than any other group of species, partly for the reason that their past geographical ranges are well documented, which facilitates species-climate modelling and other analyses. The aim of the modelling studies is to examine to what extent shifts in species distributions can be explained by climatic and other factors. Models can also be used to predict the future distributions of species. In this thesis, I have studied the response to climate change of one species of butterfly within one geographically restricted area. The study species, the European map butterfly (Araschnia levana), has expanded rapidly northwards in Finland during the last two decades. I used statistical and dynamic modelling approaches in combination with field studies to analyse the effects of climate warming and landscape structure on the expansion. I studied possible role of molecular variation in phosphoglucose isomerase (PGI), a glycolytic enzyme affecting flight metabolism and thereby flight performance, in the observed expansion of the map butterfly at two separate expansion fronts in Finland. The expansion rate of the map butterfly was shown to be correlated with the frequency of warmer than average summers during the study period. The result is in line with the greater probability of occurrence of the second generation during warm summers and previous results on this species showing greater mobility of the second than first generation individuals. The results of a field study in this thesis indicated low mobility of the first generation butterflies. Climatic variables alone were not sufficient to explain the observed expansion in Finland. There are also problems in transferring the climate model to new regions from the ones from which data were available to construct the model. The climate model predicted a wider distribution in the south-western part of Finland than what has been observed. Dynamic modelling of the expansion in response to landscape structure suggested that habitat and landscape structure influence the rate of expansion. In southern Finland the landscape structure may have slowed down the expansion rate. The results on PGI suggested that allelic variation in this enzyme may influence flight performance and thereby the rate of expansion. Genetic differences of the populations at the two expansion fronts may explain at least partly the observed differences in the rate of expansion. Individuals with the genotype associated with high flight metabolic rate were most frequent in eastern Finland, where the rate of range expansion has been highest.

Major histocompatibility complex and coronary artery disease: Special emphasis on Chlamydia pneumoniae, and periodontitis

Most of the genes in the MHC region are involveed in adaptive and innate immunity, with essential function in inflammatory reactions and in protection against infections. These genes might serve as a candidate region for infection and inflammation associated diseases. CAD is an inflammatory disease. The present set of studies was performed to assess whether the MHC region harbors genetic markers for CAD, and whether these genetic markers explain the CAD risk factors: e.g., C. pneumoniae, periodontitis, and periodontal pathogens. Study I was performed using two separate patient materials and age- and sex-matched healthy controls, categorizing them into two independent studies: the HTx and ACS studies. Both studies consistently showed the HLA-A3– B35– DR1 (35 ancestral haplotype) haplotype as a susceptible MHC genetic marker for CAD. HLA-DR1 alone was associated not only with CAD, but also with CAD risk factor diseases, e.g., diabetes mellitus, and hyperlipidemia. The ACS study further showed the HLA-B*07 and -DRB1*15 -related haplotype as a protective MHC haplotype for CAD. Study II showed that patients with CAD showed signs of chronic C. pneumoniae infection when compared to age- and sex-matched healthy controls. HLA-B*35 or -related haplotypes associated with the C. pneumoniae infection markers. Among these haplotype carriers, males and smokers associated with elevated C. pneumoniae infection markers. Study III showed that CAD patients with periodontitis had elevated serum markers of P. gingivalis and occurrence of the pathogen in saliva. LTA+496C strongly associated with periodontitis, while HLA-DRB1*01 with periodontitis and with the elevated serum antibodies of P. gingivalis. Study IV showed that the increased level of C3/C4 ratio was a new risk factor and was associated with recurrent cardiovascular end-points. The increased C3 and decreased C4 concentrations in serum explained the increased level of the C3/C4 ratio. Both the higher than cut-off value (4.53) and the highest quartile of the C3/C4 ratio were also associated with worst survival, increased end-points, and C4 null alleles. The presence of C4 null alleles associated with decreased serum C4 concentration, and increased C3/C4 ratio. In conclusion, the present studies show that the CAD susceptibility haplotype (HLA-A3− B35− DR1 -related haplotypes, Study I) partially explains the development of CAD in patients possessing several recognized and novel risk factors: diabetes mellitus, increased LDL, smoking, C4B*Q0, C. pneumnoiae, periodontitis, P. gingivalis, and complement C3/C4 ratio (Study II, III, and IV).

Complex trait of acute pancreatitis : Studies of candidate genes and adipokines

Acute pancreatitis (AP) is a common disease. Mild disease resolves spontaneously in a few days. Severe forms of the disease can lead to local complications, necrosis, and abscesses in and around the pancreas. Systemic inflammation in severe AP is associated with distant organ failures. The aim of this study is to identify genetically determined prognostic factors involved in the clinical features of AP. The study employs a candidate-gene approach, and the genes are involved in trysinogen activation in the initiation phase of the disease, as well as in the systemic inflammation as the disease proceeds. The last study examines adipokines, fat-derived hormones characterized with the capacity to modify inflammation. SPINK 1 is a gene coding trypsin activation inhibitor. Mutations N34S and P55N were determined by minisequencing methods in 371 AP patients and in 459 controls. The mutation N34S was more common in AP patients (7.8%) than in controls (2.6%). This suggests that SPINK 1 gene mutation N34S is a risk factor for AP. In the fourth study, in 12 matched pairs of patients with severe and mild AP, levels of adipokines, adiponectin, and leptin were evaluated. Plasma adipokine levels did not differ between patients with mild and severe AP. The results suggest that in AP, adipokine plasma levels are not factors predisposing to organ failures. This study identified the SPINK 1 mutation N34S to be a risk factor for AP in the general population. As AP is a multifactorial disease, and extensive genetic heterogeneity is likely, further identification of genetic factors in the disease requires larger future studies with more advanced genetic study models. Further identification of the patient characteristics associated with organ failures offers another direction of the study to achieve more detailed understanding of the severe form of AP.

Two views on interstellar dust: near-infrared scattering and polarized thermal dust emission

Interstellar clouds are not featureless, but show quite complex internal structures of filaments and clumps when observed with high enough resolution. These structures have been generated by 1) turbulent motions driven mainly by supernovae, 2) magnetic fields working on the ions and, through neutral-ion collisions, on neutral gas as well, and 3) self-gravity pulling a dense clump together to form a new star. The study of the cloud structure gives us information on the relative importance of each of these mechanisms, and helps us to gain a better understanding of the details of the star formation process. Interstellar dust is often used as a tracer for the interstellar gas which forms the bulk of the interstellar matter. Some of the methods that are used to derive the column density are summarized in this thesis. A new method, which uses the scattered light to map the column density in large fields with high spatial resolution, is introduced. This thesis also takes a look at the grain alignment with respect to the magnetic fields. The aligned grains give rise to the polarization of starlight and dust emission, thus revealing the magnetic field. The alignment mechanisms have been debated for the last half century. The strongest candidate at present is the radiative torques mechanism. In the first four papers included in this thesis, the scattered light method of column density estimation is formulated, tested in simulations, and finally used to obtain a column density map from observations. They demonstrate that the scattered light method is a very useful and reliable tool in column density estimation, and is able to provide higher resolution than the near-infrared color excess method. These two methods are complementary. The derived column density maps are also used to gain information on the dust emissivity within the observed cloud. The two final papers present simulations of polarized thermal dust emission assuming that the alignment happens by the radiative torques mechanism. We show that the radiative torques can explain the observed decline of the polarization degree towards dense cores. Furthermore, the results indicate that the dense cores themselves might not contribute significantly to the polarized signal, and hence one needs to be careful when interpreting the observations and deriving the magnetic field.

Proton magnetic resonance spectroscopy of a boron neutron capture therapy 10B-carrier, L-p-boronophenylalanine-fructose complex

Boron neutron capture therapy (BNCT) is a radiotherapy that has mainly been used to treat malignant brain tumours, melanomas, and head and neck cancer. In BNCT, the patient receives an intravenous infusion of a 10B-carrier, which accumulates in the tumour area. The tumour is irradiated with epithermal or thermal neutrons, which result in a boron neutron capture reaction that generates heavy particles to damage tumour cells. In Finland, boronophenylalanine fructose (BPA-F) is used as the 10B-carrier. Currently, the drifting of boron from blood to tumour as well as the spatial and temporal accumulation of boron in the brain, are not precisely known. Proton magnetic resonance spectroscopy (1H MRS) could be used for selective BPA-F detection and quantification as aromatic protons of BPA resonate in the spectrum region, which is clear of brain metabolite signals. This study, which included both phantom and in vivo studies, examined the validity of 1H MRS as a tool for BPA detection. In the phantom study, BPA quantification was studied at 1.5 and 3.0 T with single voxel 1H MRS, and at 1.5 T with magnetic resonance imaging (MRSI). The detection limit of BPA was determined in phantom conditions at 1.5 T and 3.0 T using single voxel 1H MRS, and at 1.5 T using MRSI. In phantom conditions, BPA quantification accuracy of ± 5% and ± 15% were achieved with single voxel MRS using external or internal (internal water signal) concentration references, respectively. For MRSI, a quantification accuracy of <5% was obtained using an internal concentration reference (creatine). The detection limits of BPA in phantom conditions for the PRESS sequence were 0.7 (3.0 T) and 1.4 mM (1.5 T) mM with 20 × 20 × 20 mm3 single voxel MRS, and 1.0 mM with acquisition-weighted MRSI (nominal voxel volume 10(RL) × 10(AP) × 7.5(SI) mm3), respectively. In the in vivo study, an MRSI or single voxel MRS or both was performed for ten patients (patients 1-10) on the day of BNCT. Three patients had glioblastoma multiforme (GBM), and five patients had a recurrent or progressing GBM or anaplastic astrocytoma gradus III, and two patients had head and neck cancer. For nine patients (patients 1-9), MRS/MRSI was performed 70-140 min after the second irradiation field, and for one patient (patient 10), the MRSI study began 11 min before the end of the BPA-F infusion and ended 6 min after the end of the infusion. In comparison, single voxel MRS was performed before BNCT, for two patients (patients 3 and 9), and for one patient (patient 9), MRSI was performed one month after treatment. For one patient (patient 10), MRSI was performed four days before infusion. Signals from the tumour spectrum aromatic region were detected on the day of BNCT in three patients, indicating that in favourable cases, it is possible to detect BPA in vivo in the patient’s brain after BNCT treatment or at the end of BPA-F infusion. However, because the shape and position of the detected signals did not exactly match the BPA spectrum detected in the in vitro conditions, assignment of BPA is difficult. The opportunity to perform MRS immediately after the end of BPA-F infusion for more patients is necessary to evaluate the suitability of 1H MRS for BPA detection or quantification for treatment planning purposes. However, it could be possible to use MRSI as criteria in selecting patients for BNCT.

Molecular line and continuum studies of the early stages of star formation

New stars form in dense interstellar clouds of gas and dust called molecular clouds. The actual sites where the process of star formation takes place are the dense clumps and cores deeply embedded in molecular clouds. The details of the star formation process are complex and not completely understood. Thus, determining the physical and chemical properties of molecular cloud cores is necessary for a better understanding of how stars are formed. Some of the main features of the origin of low-mass stars, like the Sun, are already relatively well-known, though many details of the process are still under debate. The mechanism through which high-mass stars form, on the other hand, is poorly understood. Although it is likely that the formation of high-mass stars shares many properties similar to those of low-mass stars, the very first steps of the evolutionary sequence are unclear. Observational studies of star formation are carried out particularly at infrared, submillimetre, millimetre, and radio wavelengths. Much of our knowledge about the early stages of star formation in our Milky Way galaxy is obtained through molecular spectral line and dust continuum observations. The continuum emission of cold dust is one of the best tracers of the column density of molecular hydrogen, the main constituent of molecular clouds. Consequently, dust continuum observations provide a powerful tool to map large portions across molecular clouds, and to identify the dense star-forming sites within them. Molecular line observations, on the other hand, provide information on the gas kinematics and temperature. Together, these two observational tools provide an efficient way to study the dense interstellar gas and the associated dust that form new stars. The properties of highly obscured young stars can be further examined through radio continuum observations at centimetre wavelengths. For example, radio continuum emission carries useful information on conditions in the protostar+disk interaction region where protostellar jets are launched. In this PhD thesis, we study the physical and chemical properties of dense clumps and cores in both low- and high-mass star-forming regions. The sources are mainly studied in a statistical sense, but also in more detail. In this way, we are able to examine the general characteristics of the early stages of star formation, cloud properties on large scales (such as fragmentation), and some of the initial conditions of the collapse process that leads to the formation of a star. The studies presented in this thesis are mainly based on molecular line and dust continuum observations. These are combined with archival observations at infrared wavelengths in order to study the protostellar content of the cloud cores. In addition, centimetre radio continuum emission from young stellar objects (YSOs; i.e., protostars and pre-main sequence stars) is studied in this thesis to determine their evolutionary stages. The main results of this thesis are as follows: i) filamentary and sheet-like molecular cloud structures, such as infrared dark clouds (IRDCs), are likely to be caused by supersonic turbulence but their fragmentation at the scale of cores could be due to gravo-thermal instability; ii) the core evolution in the Orion B9 star-forming region appears to be dynamic and the role played by slow ambipolar diffusion in the formation and collapse of the cores may not be significant; iii) the study of the R CrA star-forming region suggests that the centimetre radio emission properties of a YSO are likely to change with its evolutionary stage; iv) the IRDC G304.74+01.32 contains candidate high-mass starless cores which may represent the very first steps of high-mass star and star cluster formation; v) SiO outflow signatures are seen in several high-mass star-forming regions which suggest that high-mass stars form in a similar way as their low-mass counterparts, i.e., via disk accretion. The results presented in this thesis provide constraints on the initial conditions and early stages of both low- and high-mass star formation. In particular, this thesis presents several observational results on the early stages of clustered star formation, which is the dominant mode of star formation in our Galaxy.

Functional role of the Mso1p-Sec1p complex in membrane fusion regulation

Sec1/Munc18 (SM) protein family members are evolutionary conserved proteins. They perform an essential, albeit poorly understood function in SNARE complex formation in membrane fusion. In addition to the SNARE complex components, only a few SM protein binding proteins are known. Typically, their binding modes to SM proteins and their contribution to the membrane fusion regulation is poorly characterised. We identified Mso1p as a novel Sec1p interacting partner. It was shown that Mso1p and Sec1p interact at sites of polarised secretion and that this localisation is dependent on the Rab GTPase Sec4p and its GEF Sec2p. Using targeted mutagenesis and N- and C-terminal deletants, it was discovered that the interaction between an N-terminal peptide of Mso1p and the putative Syntaxin N-peptide binding area in Sec1p domain 1 is important for membrane fusion regulation. The yeast Syntaxin homologues Sso1p and Sso2p lack the N-terminal peptide. Our results show that in addition to binding to the putative N-peptide binding area in Sec1p, Mso1p can interact with Sso1p and Sso2p. This result suggests that Mso1p can mimic the N-peptide binding to facilitate membrane fusion. In addition to Mso1p, a novel role in membrane fusion regulation was revealed for the Sec1p C-terminal tail, which is missing in its mammalian homologues. Deletion of the Sec1p-tail results in temperature sensitive growth and reduced sporulation. Using in vivo and in vitro experiments, it was shown that the Sec1p-tail mediates SNARE complex binding and assembly. These results propose a regulatory role for the Sec1p-tail in SNARE complex formation. Furthermore, two novel interaction partners for Mso1p, the Rab GTPase Sec4p and plasma membrane phospholipids, were identified. The Sec4p link was identified using Bimolecular Fluorescence Complementation assays with Mso1p and the non-SNARE binding Sec1p(1-657). The assay revealed that Mso1p can target Sec1p(1-657) to sites of secretion. This effect is mediated via the Mso1p C-terminus, which previously has been genetically linked to Sec4p. These results and in vitro binding experiments suggest that Mso1p acts in cooperation with the GTP-bound form of Sec4p on vesicle-like structures prior to membrane fusion. Mso1p shares homology with the PIP2 binding domain of the mammalian Munc18 binding Mint proteins. It was shown both in vivo and in vitro that Mso1p is a phospholipid inserting protein and that this insertion is mediated by the conserved Mso1p amino terminus. In vivo, the Mso1p phospholipid binding is needed for sporulation and Mso1p-Sec1p localisation at the sites of secretion at the plasma membrane. The results reveal a novel layer of membrane fusion regulation in exocytosis and propose a coordinating role for Mso1p in connection with membrane lipids, Sec1p, Sec4p and SNARE complexes in this process.

Identifying Genetic Risk Factors in Canine Autoimmune Disorders

Autoimmune diseases are more common in dogs than in humans and are already threatening the future of some highly predisposed dog breeds. Susceptibility to autoimmune diseases is controlled by environmental and genetic factors, especially the major histocompatibility complex (MHC) gene region. Dogs show a similar physiology, disease presentation and clinical response as humans, making them an excellent disease model for autoimmune diseases common to both species. The genetic background of canine autoimmune disorders is largely unknown, but recent annotation of the dog genome and subsequent development of new genomic tools offer a unique opportunity to map novel autoimmune genes in various breeds. Many autoimmune disorders show breed-specific enrichment, supporting a strong genetic background. Furthermore, the presence of hundreds of breeds as genetic isolates facilitates gene mapping in complex autoimmune disorders. Identification of novel predisposing genes establishes breeds as models and may reveal novel candidate genes for the corresponding human disorders. Genetic studies will eventually shed light on common biological functions and interactions between genes and the environment. This study aimed to identify genetic risk factors in various autoimmune disorders, including systemic lupus erythematosus (SLE)-related diseases, comprising immune-mediated rheumatic disease (IMRD) and steroid-responsive meningitis arteritis (SMRA) as well as Addison s disease (AD) in Nova Scotia Duck Tolling Retrievers (NSDTRs) and chronic superficial keratitis (CSK) in German Shepherd dogs (GSDs). We used two different approaches to identify genetic risk factors. Firstly, a candidate gene approach was applied to test the potential association of MHC class II, also known as a dog leukocyte antigen (DLA) in canine species. Secondly, a genome-wide association study (GWAS) was performed to identify novel risk loci for SLE-related disease and AD in NSDTRs. We identified DLA risk haplotypes for an IMRD subphenotype of SLE-related disease, AD and CSK, but not in SMRA, and show that the MHC class II gene region is a major genetic risk factor in canine autoimmune diseases. An elevated risk was found for IMRD in dogs that carried the DLA-DRB1*00601/DQA1*005011/DQB1*02001 haplotype (OR = 2.0, 99% CI = 1.03-3.95, p = 0.01) and for ANA-positive IMRD dogs (OR = 2.3, 99% CI = 1.07-5.04, p-value 0.007). We also found that DLA-DRB1*01502/DQA*00601/DQB1*02301 haplotype was significantly associated with AD in NSDTRs (OR = 2.1, CI = 1.0-4.4, P = 0.044) and the DLA-DRB1*01501/DQA1*00601/DQB1*00301 haplotype with the CSK in GSDs (OR=2.67, CI=1.17-6.44, p= 0.02). In addition, we found that homozygosity for the risk haplotype increases the risk for each disease phenotype and that an overall homozygosity for the DLA region predisposes to CSK and AD. Our results have enabled the development of genetic tests to improve breeding practices by avoiding the production of puppies homozygous for risk haplotypes. We also performed the first successful GWAS for a complex disease in dogs. With less than 100 cases and 100 controls, we identified five risk loci for SLE-related disease and AD and found strong candidate genes involved in a novel T-cell activation pathway. We show that an inbred dog population has fewer risk factors, but each of them has a stronger genetic risk. Ongoing studies aim to identify the causative mutations and bring new knowledge to help diagnostics, treatment and understanding of the aetiology of SLE-related diseases.