Recent technological advances and remarkable successes have led to genome-wide association studies (GWAS) becoming a tool of choice for investigating the genetic basis of common complex human diseases. These studies typically involve samples from thousands of individuals, scanning their DNA at up to a million loci along the genome to discover genetic variants that affect disease risk. Hundreds of such variants are now known for common diseases, nearly all discovered by GWAS over the last three years. As a result, many new studies are planned for the future or are already underway. In this thesis, I present analysis results from actual studies and some developments in theory and methodology. The Wellcome Trust Case Control Consortium (WTCCC) published one of the first large-scale GWAS in 2007. I describe my contribution to this study and present the results from some of my follow-up analyses. I also present results from a GWAS of a bipolar disorder sub-phenotype, and a recent and on-going fine mapping experiment. Building on methods developed as part of the WTCCC, I describe a Bayesian approach to GWAS analysis and compare it to widely used frequentist approaches. I do so both theoretically, by interpreting each approach from the perspective of the other, and empirically, by comparing their performance in the context of replicated GWAS findings. I discuss the implications of these comparisons on the interpretation and analysis of GWAS generally, highlighting the advantages of the Bayesian approach. Finally, I examine the effect of linkage disequilibrium on the detection and estimation of various types of genetic effects, particularly non-additive effects. I derive a theoretical result showing how the power to detect a departure from an additive model at a marker locus decays faster than the power to detect an association.

En 1975, Wilson et King ont proposé que l'évolution opère non seulement via des changements affectant la structure des protéines, mais aussi via des mutations qui modifient la régulation génétique. L'étude des éléments régulateurs de l'expression génétique a un rôle important dans la compréhension de l'expression de différentes maladies et de la réponse thérapeutique. Nous avons développé un algorithme bio- informatique qui nous permet rapidement de trouver des sites de régulation génétique à travers tout le génome et pour une grande quantité de gènes. Notre approche consiste à trouver des sites polymorphes (SNPs) qui sont en déséquilibre de liaison avec le débalancement allélique (AI) afin de cartographier la région régulatrice et le site responsable. Notre méthode est avantageuse par rapport à d'autres méthodes, car elle n'a pas besoin des données « phasées». De plus, les données de débalancement allélique ne sont pas affectées par des facteurs externes étant donné qu'ils sont mesurés dans la même cellule. Nous avons démontré que notre approche est fiable et qu'elle peut détecter des sites loin du gène. De plus, il peut être appliqué à des données de génotypage sans avoir besoin de les « phaser » . / Wilson and King (1975) proposed that evolution frequently operates through mutations affecting genetic regulation. Likewise, it is expected that genetic variation responsible for inter-individual differences will be due to variation in regulatory sites. Identifying such sites is thus important in the genetic and medical research. We have developed a new bioinformatics algorithm to find genome-wide regulatory sites for a big number of genes. Individuals carrying different alleles at a regulatory site will exhibit allelic imbalance(AI) due to differential expression of the two copies the same locus. Our approach consists of searching polymorphic sites (SNPs) in linkage disequilibrium with AI in order to map regulatory regions. We have detected many SNPs associated to the regulation of different genes pointed in previous studies. We have also found regulatory regions far from the transcription start site (TSS). The major advantage of this method is that phased data is not needed. In addition, AI data has the benefit of not being affected by external factors since it is measured in the same cell. The results show that our approach is reliable and it can detect sites far from the gene.

Les diagrammes de transitions d'états ont été réalisés avec le logiciel Latex. / Le processus ancestral est appliqué pour étudier la variabilité génétique et la mesure de déséquilibre de liaison de séquences d’ADN, et faire de l’inférence statistique sur les divers facteurs responsables de cette variabilité. En tenant compte, en premier lieu, des facteurs de dérive génétique, de mutation, et de recombinaison, les calculs exacts de la mesure de déséquilibre de liaison de deux loci sont retrouvés. De plus, une approximation du processus exact, SMC (sequentially Markov chain), est utilisée pour trouver la mesure d’association à deux loci, et une formule de covariance pour calculer cette mesure est corrigée. En intégrant le facteur de conversion dans le modèle de Moran, on trouve l’espérance des mesures de polymorphisme exprimées par les espérances des mesures de variation intra-locus et inter-locus. Celles-ci sont calculées à l’aide de temps espérés dans les états ancestraux. De plus, l’espérance du déséquilibre de liaison est trouvée et il est montré qu’elle diminue quand le taux de recombinaison augmente. En utilisant ces résultats théoriques, on présente une méthode pour estimer les paramètres de mutation, de recombinaison, et de conversion. / The ancestral process is applied to investigate the amount of DNA variation and the amount of linkage disequilibrium ; it is also applied to make statistical inference about the multiple factors responsible for this variation. Considering genetic drift, mutation, and recombination events, the exact solutions for linkage disequilibrium between two loci are obtained. Furthermore, the association measure between two loci is obtained by using an approximation of the exact process, SMC (sequentially Markov chain), and correcting a covariance formula. After introducing intrachromosomal gene conversion under the Moran model, the expected amounts of variation within and between two loci are obtained using expected times spent in the ancestral states. Furthermore, the expectation of linkage disequilibrium is obtained and it is shown to decrease as the recombination rate is increased. Using these theoretical results, a method for estimating the mutation, recombination and gene conversion parameters is presented.

The selection on genetic values predicted from markers could substantially increase the rate of genetic gain in animals by increasing accuracy of prediction and reducing generation interval, especially for difficult to measure traits, such as feed efficiency. Feed efficiency is the most important trait in animal production due to its impacts on cost of production and environmental factors. Many metrics measure the feed efficiency, such as ratio of gain to feed (FER), the ratio of feed to gain (FCR) and residual feed intake (RFI). Nevertheless, in ovine, no study with the aim of understand the genetic variants or the accuracy of genomic estimated breeding value (GEBV) for feed efficiency traits was published yet. Moreover, before to apply the genomic information, it is necessary to understand and characterized the population structure, for instance, by linkage disequilibrium (LD). Both genome-wide association studies (GWAS) and genomic selection (GS) leverage LD between marker and causal mutation. Based on the above considerations, the aim of this study was to map LD in ovine, characterized by Brazilian Santa Inês sheep; to search genetic variants for feed efficiency traits (FER, FCR and RFI) through GWAS; and to verify the accuracy of GEBV for RFI. In total, 396 samples (animals) of Longissimus dorsi muscle were collect. A high-density panel of SNP (Illumina High-Density Ovine SNP BeadChip®) comprising 54,241 SNPs was used to obtain the genotyping data. The phenotype data was comprised of 387 animals. The average LD between adjacent markers for two LD metrics, r² and |D\'|, were 0.166 and 0.617, respectively. The degree of LD estimated was lower than reported in other species and it was characterized by short haplotype blocks. Consequently, for genomic analyses, high-density panels of marker are recommended. Many markers were associated to feed efficiency traits in GWAS, mainly to RFI trait. Few candidate genes were reported in this study, highlighting NRF-1 (nuclear respiratory factor 1), which controls mitochondrial biosynthesis, the most important process responsible by a great fraction of the produced energy. Finally, we verified the accuracy of GEBV for RFI using few Bayesian regression models, and we found low accuracy, ranging from 0.033 (BayesB with π=0.9912) to 0.036 (BayesA), which might be explained by the low relationship among animals and small training population. / A seleção com base nos valores genéticos genômicos preditos pode aumentar substancialmente a taxa de ganho genético em animais por meio do aumento da acurácia de predição e redução do intervalo de gerações, especialmente para características de difícil e/ou onerosa mensuração, como eficiência alimentar. A eficiência alimentar é uma das características mais importantes na produção animal devido principalmente aos seus impactos econômicos e ambientais. Muitas métricas representam a eficiência alimentar, por exemplo: a relação do ganho de peso e consumo alimentar (EA), a proporção do consumo alimentar e ganho de peso (CA) e o consumo alimentar residual (CAR). Em ovinos, nenhum estudo com o objetivo de buscar variantes genéticas ou verificar a acurácia do valor genético genômico estimado para eficiência alimentar foi publicado. Adicionalmente, antes de aplicar a informação genômica, é necessário compreender e caracterizar a estrutura da população, como por meio do desequilíbrio de ligação (LD). O estudo de associação genômica (GWAS) e seleção genômica (GS) consideram o LD entre marcador e a mutação causal. Com base nas considerações acima, o objetivo deste estudo foi mapear o LD em ovinos, caracterizado pela raça ovina Santa Inês; localizar variantes genéticas para as características de eficiência alimentar (EA, CA e CAR) utilizando a abordagem GWAS; e verificar a acurácia da estimação dos valores genéticos genômico para o CAR. No total, foram coletadas 396 amostras (animais) do músculo Longissimus dorsi, para posterior genotipagem utilizando o painel de alta densidade (Illumina High-Density Ovine SNP BeadChip®), compreendendo 54.241 SNPs. O banco fenotípico é composto por 387 animais. O LD médio entre marcadores adjacentes para duas métricas de LD, r² e |D\'|, foram 0,166 e 0,617, respectivamente. O grau de LD estimado foi menor que o relatado em outras espécies e foi caracterizado por blocos de haplótipos curtos. Consequentemente, para as análises genômicas são recomendados painéis de marcadores de alta densidade. No GWAS, foram encontrados muitos marcadores associados aos fenótipos, em especial, à característica CAR. Alguns genes candidatos foram relatados neste estudo, destacando-se o NRF-1 (fator respiratório nuclear 1), que controla a biossíntese mitocondrial, o processo mais importante responsável por grande parte da produção de energia. Finalmente, verificamos a acurácia do valor genético genômico estimado para o CAR usando modelos de regressão Bayesiana, e encontramos baixos valores para acurácia (0,033 a 0,036) o que pode ser explicado pelo baixo grau de relacionamento entre os indivíduos e tamanho reduzido da população de treinamento.

Die Photoparoxysmale Reaktion (PPR) oder Photosensitivität ist elektrophysiologisch durch eine synchronisierte Erregungssteigerung des visuellen Kortex infolge intermittierender Lichtreize gekennzeichnet. Familienstudien belegen, dass die PPR eine genetisch komplexe Disposition aufweist. Visuell induzierte Anfälle treten bei bis zu 10% der Kinder mit Epilepsien auf. Besonders häufig (bis zu 30%) ist die PPR mit den idiopathisch generalisierten Epilepsien (IGE) assoziiert. Die PPR repräsentiert somit einen geeigneten Endophänotyp, um die pathogenetische Komplexität der IGE zu reduzieren und auf einfachere Phänotyp-Genotyp-Beziehungen zurückzuführen. Die Aufgabenstellung der vorliegenden Arbeit war die molekulargenetische Identifizierung von PPR-Determinanten und die Untersuchung des genetischen Zusammenhanges dieser Determinanten mit den häufigen IGE. Dazu wurden die Familien in zwei distinkte, phänotypisch homogenere Subgruppen unterteilt. Bei 25 Familien mit photosensitiver IGE fand sich ein signifikanter IGE/PPR-Locus in der Region 13q31, während sich bei 19 Familien mit PPR und visuell induzierten Anfällen ohne IGE ein signifikanter PPR-Locus im Bereich 6p21 eingrenzen ließ. Ein positionell und funktionell interessantes Kandidatengen für den 6p21-Locus ist das Gen BRD2 (Bromodomain-containing Protein 2). Für die juvenile myoklonische Epilepsie (JME) wurde in der Region 6p21.3 ein Genort kartiert und ein Kopplungsungleichgewicht zwischen BRD2 Polymorphismen und JME nachgewiesen. Ca. 30% der JME Patienten sind photosensitiv. Unsere Studie ergab eine signifikante Assoziation von sechs BRD2-Polymorphismen mit PPR, die auf eine Beteiligung des BRD2 Gens an der Pathogenese der PPR hinweisen. / Photosensitivity or photoparoxysmal response (PPR) is an electroencephalographic (EEG) trait characterized by an abnormal visual sensitivity of the brain in reaction to intermittent photic stimulation. Family and twin studies provide unequivocal evidence that PPR is genetically determined. Seizures triggered by visual stimuli occur in up to 10% of patients with epilepsies in childhood. PPR is frequently (30%) associated with idiopathic generalized epilepsies (IGE). By focusing on epilepsy-related EEG patterns such as PPR, which can be evoked in the EEG laboratory independent of obvious clinical expression of seizures, PPR represents a suitable endophenotype to search for susceptibility loci involved in cortical hyperexcitability and epileptogenesis. We carried out a genome wide linkage scan to map susceptibility loci for PPR and to explore their genetic relationship with idiopathic generalized epilepsy (IGE). To dissect PPR-related traits, we defined two distinct subgroups of families differing in their familial background on IGE. Our analyses provided significant evidence for linkage to the region 13q31 in PPR/IGE families, in which PPR was strongly associated with IGE, and significant evidence for linkage to the region 6p21 in PPR-families with predominantly pure PPR. The locus on 13q31 seems to confer susceptibility to photosensitive IGE, whereas the locus on 6p21 may predispose to PPR itself. The linkage region on 6p21 harbors a high-ranking candidate gene: BRD2 (bromodomain containing protein2). Linkage disequilibrium mapping in the candidate region 6p21.3 revealed strong associations between polymorphisms of the gene BRD2 and juvenile myoclonic epilepsy (JME), which exhibits PPR in 30% of JME patients. We found associations between PPR and six BRD2 polymorphisms. Considering the strong clinical association of JME and PPR, the present results support evidence that PPR and JME share epileptogenic pathways, for which BRD2 might be an underlying susceptibility gene.

Submitted by Luciana Ferreira (lucgeral@gmail.com) on 2017-01-04T13:00:39Z No. of bitstreams: 2 Dissertação - Kaisson Ernane dos Santos - 2013.pdf: 2738475 bytes, checksum: 6c79ab9177dd126fa5cb677127debc2c (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Approved for entry into archive by Luciana Ferreira (lucgeral@gmail.com) on 2017-01-04T13:01:13Z (GMT) No. of bitstreams: 2 Dissertação - Kaisson Ernane dos Santos - 2013.pdf: 2738475 bytes, checksum: 6c79ab9177dd126fa5cb677127debc2c (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Made available in DSpace on 2017-01-04T13:01:13Z (GMT). No. of bitstreams: 2 Dissertação - Kaisson Ernane dos Santos - 2013.pdf: 2738475 bytes, checksum: 6c79ab9177dd126fa5cb677127debc2c (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) Previous issue date: 2013-11-25 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / The Major Histocompatibility Complex is mainly composed by genes of the adaptive immune response. In humans, part of this complex is known as the Human Leukocyte Antigens (HLA), whose genes are responsible for specific antigen presentation to effector immune cells. The classical class I HLA genes (HLA-A, -B and -C) are responsible for antigen presentation to T CD8+ cells and they constitute the most polymorphic genes in the human genome. This variability is maintained by selection mediated by microorganisms. In contrast to their classical counterparts, the non classical class I genes (HLA-G, -E and -F) present low variability and are associated with immune tolerance due to the interaction with NK and T cells inhibitor receptors. HLA-G is the most studied non classical gene, which is associated with immune response modulation, mainly during pregnancy. Considering that natural selection is acting on the HLA-G regulatory regions maintaining high heterozigosity in this region, we evaluated a nearby Alu insertion (AluyHG) correlating this Alu element with coding and 3’UTR HLA-G polymorphisms. The AluyHG insertion was particularly associated with the HLA-G haplotype known as G*01:01:01:01/UTR-1, considered a high-expressing HLA-G haplotype. The G*01:01:01:01/UTR-1/AluyHG haplotype would be the most recent HLA-G haplotypes, in spite of its high frequency in worldwide populations. / O Complexo Principal de Histocompatibilidade (MHC) é formado principalmente por genes que participam da resposta imunológica adaptativa. Entre esses genes encontramos o grupo denominado de Antígenos Leucocitários Humanos (HLA), que são responsáveis pela apresentação de antígenos específicos às células efetoras do sistema imunológico. Os genes HLA de classe I clássicos (HLA-A, -B e -C), responsáveis pela apresentação antigênica aos linfócitos T citotóxicos, são considerado como os mais polimórficos do genoma humano e de outros vertebrados. A variabilidade desses genes e elevada heterozigose é mantida por seleção mediada por microrganismos. Diferentemente dos genes clássicos, os genes HLA de classe I não clássicos (HLA-G, -E e -F) apresentam variabilidade reduzida e como função principal a tolerância imunológica, por meio de sua interação com receptores inibitórios presentes nas células NK e T. O HLA-G é o mais estudado entre esses genes e, devido sua importância como molécula imunomoduladora e sua importância em situações como gestação, e considerando evidências anteriores de seleção natural mantendo uma elevada heterozigose nas regiões regulatórias do HLA-G, avaliamos a presença de uma inserção Alu (AluyHG) próxima a este gene correlacionando os achados com a variabilidade contida nas suas regiões codificadora e 3’ não traduzida. A inserção AluyHG mostrou-se em desequilíbrio de ligação (LD) com os polimorfismos do gene HLA-G. Especificamente, o elemento inserido apresentou-se em LD com um haplótipo denominado G*01:01:01:01/UTR-1, considerado como um haplótipo de alta produção da molécula de HLA-G. Esse haplótipo aparentemente é o mais jovem entre humanos, apesar de sua elevada frequência nas populações estudadas até o momento.

<p>The breakdown of self-tolerance is the main driving force behind susceptibility to SLE. When this occurs, T and B cells are activated in an uncontrolled manner and produce autoantibodies against self fragmented DNA, RNA and sometimes other parts of the cell such as cardiolipin, phosphatidylserine, etc.</p><p>The mechanism behind the breakdown of self-tolerance may be genetic factors that are triggered by environmental factors. SLE is not caused by a single gene, but by many genes, and is thus a polygenic disease. So far only a few genes have been found to be associated with SLE including PDCD1, FcγRs, and PTPN22. The main aim of my thesis is to find susceptibility genes responsible for SLE.</p><p>Recently, a gene called IRF5 was found to be associated with SLE. In paper one, we performed a thorough study and confirmed its association to SLE. In addition, we found a few other SNPs in the gene that were associated to the disease. Among them, SNP rs2004640 is very strongly associated and was found to affect the splicing of the gene. Another SNP, rs2280714, correlated with overexpression of the gene, although SNP rs10954213 was much more highly correlated with expression adding to this, in paper two we found a few other SNPs that were associated to SLE and played crucial roles in gene function. An indel in exon 6, though not associated by itself, regulated which isoforms were expressed. Individuals with 2 repeats expressed isoforms V1 and V4, while individuals with 4 repeats expressed isoforms V5 and V6. SNP rs2070197 was also very strongly associated, but did not have a functional role. In paper three, the same polymorphisms were studied in a Mexican population, which showed an even stronger association when compared to a European population.</p><p>It is known that autoimmune diseases share susceptibility genes, therefore we wanted to see if the IRF5 gene is associated with any other autoimmune diseases. In papers four and five, we tested its association to RA (using three sets of patients and controls from Sweden, Argentina and Spain) and psoriasis (using a set of patients and controls from Sweden). Association was not found in either of the diseases. Therefore, we believe that this association may be SLE-specific.</p>

The breakdown of self-tolerance is the main driving force behind susceptibility to SLE. When this occurs, T and B cells are activated in an uncontrolled manner and produce autoantibodies against self fragmented DNA, RNA and sometimes other parts of the cell such as cardiolipin, phosphatidylserine, etc. The mechanism behind the breakdown of self-tolerance may be genetic factors that are triggered by environmental factors. SLE is not caused by a single gene, but by many genes, and is thus a polygenic disease. So far only a few genes have been found to be associated with SLE including PDCD1, FcγRs, and PTPN22. The main aim of my thesis is to find susceptibility genes responsible for SLE. Recently, a gene called IRF5 was found to be associated with SLE. In paper one, we performed a thorough study and confirmed its association to SLE. In addition, we found a few other SNPs in the gene that were associated to the disease. Among them, SNP rs2004640 is very strongly associated and was found to affect the splicing of the gene. Another SNP, rs2280714, correlated with overexpression of the gene, although SNP rs10954213 was much more highly correlated with expression adding to this, in paper two we found a few other SNPs that were associated to SLE and played crucial roles in gene function. An indel in exon 6, though not associated by itself, regulated which isoforms were expressed. Individuals with 2 repeats expressed isoforms V1 and V4, while individuals with 4 repeats expressed isoforms V5 and V6. SNP rs2070197 was also very strongly associated, but did not have a functional role. In paper three, the same polymorphisms were studied in a Mexican population, which showed an even stronger association when compared to a European population. It is known that autoimmune diseases share susceptibility genes, therefore we wanted to see if the IRF5 gene is associated with any other autoimmune diseases. In papers four and five, we tested its association to RA (using three sets of patients and controls from Sweden, Argentina and Spain) and psoriasis (using a set of patients and controls from Sweden). Association was not found in either of the diseases. Therefore, we believe that this association may be SLE-specific.

The dog was the first animal to be domesticated and the process started at least 15 000 years ago. Today it is the most morphologically diverse mammal, with a huge variation in size and shape. Dogs have always been useful to humans in several ways, from being a food source, hunting companion, guard, social companion and lately also a model for scientific research. This thesis describes some of the changes that have occurred in the dog’s genome, both during the domestication process and later through breed creation. To give a more comprehensive view, three genetic systems were studied: maternally inherited mitochondrial DNA, paternally inherited Y chromosome and biparental autosomal chromosomes. I also sequenced complete mitochondrial genomes to view the effect new living conditions might have had on dogs’ genes after domestication. Finally, knowledge of the genetic structure in purebred dogs was used to test analytic methods usable in other species or in natural populations where little information is available. The domestication process appears to have caused a relaxation of the selective constraint in the mitochondrial genome, leading to a faster rate of accumulation of nonsynonymous changes in the mitochondrial genes. Later, the process of breed creation resulted in genetically separated breed groups. Breeds are a result from an unequal contribution of males and females with only a few popular sires contributing and a larger amount of dams. However, modern breeder preferences might lead to disruptive selective forces within breeds, which can result in additional fragmentation of breeds. The increase in linkage disequilibrium that this represents increases the value of purebred dogs as model organisms for the identification and mapping of diseases and traits. Purebred dogs’ potential for these kinds of studies will probably increase the more we know about the dog’s genome.

Polymorphism in coding regions and regions controlling gene expression is the major determinant of adaptive differences in natural populations. Genes of the major histocompatibility complex (MHC) possess a high level of genetic variation, which is maintained by selection over long coalescence times. MHC genes encode antigen-presenting molecules in the adaptive immune system, which protects the host from infectious diseases. However, MHC molecules may also present self-peptides and for most autoimmune diseases there is a genetic factor associated with the MHC. MHC genes have been used to learn about the interplay of selection and historical population events. In domestic dogs and their progenitor, the wolf, I explored factors associated with domestication and breed formation and their influence not only on MHC coding regions but also on the haplotypic structure of the class II region. Polymorphism and strong selection was demonstrated in the proximal promoters of MHC genes in dogs and wolves. Hence, genetic variation associated with MHC gene expression may have at least equal importance for a well functioning immune system. Associations between promoter sequences and particular coding alleles suggested allele-specific expression patterns. SNP haplotypes of the MHC class II region revealed ancestral as well as convergent haplotypes, in which combinations of alleles are kept by selection. Interestingly, weaker allelic associations were found between different genes and between coding regions and promoters in dogs compared to wolves. Potentially, this could cause insufficient defense against infections and predispose dogs to autoimmune diseases. For example, I identified a site in the promoter region that showed a consistent difference between haplotypes conferring susceptibility and protection to diabetes in dogs, which should be investigated further. Furthermore, I investigated how selection and demographic changes associated with glacial and inter-glacial periods have affected MHC variation in European hedgehogs and extended the prevailing knowledge concerning their population history.