Genetic association studies under the population stratification, family pedigree and application to genome-wide association studies

This dissertation has three separate parts: the first part deals with the general pedigree association testing incorporating continuous covariates; the second part deals with the association tests under population stratification using the conditional likelihood tests; the third part deals with the genome-wide association studies based on the real rheumatoid arthritis (RA) disease data sets from Genetic Analysis Workshop 16 (GAW16) problem 1. Many statistical tests are developed to test the linkage and association using either case-control status or phenotype covariates for family data structure, separately. Those univariate analyses might not use all the information coming from the family members in practical studies. On the other hand, the human complex disease do not have a clear inheritance pattern, there might exist the gene interactions or act independently. In part I, the new proposed approach MPDT is focused on how to use both the case control information as well as the phenotype covariates. This approach can be applied to detect multiple marker effects. Based on the two existing popular statistics in family studies for case-control and quantitative traits respectively, the new approach could be used in the simple family structure data set as well as general pedigree structure. The combined statistics are calculated using the two statistics; A permutation procedure is applied for assessing the p-value with adjustment from the Bonferroni for the multiple markers. We use simulation studies to evaluate the type I error rates and the powers of the proposed approach. Our results show that the combined test using both case-control information and phenotype covariates not only has the correct type I error rates but also is more powerful than the other existing methods. For multiple marker interactions, our proposed method is also very powerful. Selective genotyping is an economical strategy in detecting and mapping quantitative trait loci in the genetic dissection of complex disease. When the samples arise from different ethnic groups or an admixture population, all the existing selective genotyping methods may result in spurious association due to different ancestry distributions. The problem can be more serious when the sample size is large, a general requirement to obtain sufficient power to detect modest genetic effects for most complex traits. In part II, I describe a useful strategy in selective genotyping while population stratification is present. Our procedure used a principal component based approach to eliminate any effect of population stratification. The paper evaluates the performance of our procedure using both simulated data from an early study data sets and also the HapMap data sets in a variety of population admixture models generated from empirical data. There are one binary trait and two continuous traits in the rheumatoid arthritis dataset of Problem 1 in the Genetic Analysis Workshop 16 (GAW16): RA status, AntiCCP and IgM. To allow multiple traits, we suggest a set of SNP-level F statistics by the concept of multiple-correlation to measure the genetic association between multiple trait values and SNP-specific genotypic scores and obtain their null distributions. Hereby, we perform 6 genome-wide association analyses using the novel one- and two-stage approaches which are based on single, double and triple traits. Incorporating all these 6 analyses, we successfully validate the SNPs which have been identified to be responsible for rheumatoid arthritis in the literature and detect more disease susceptibility SNPs for follow-up studies in the future. Except for chromosome 13 and 18, each of the others is found to harbour susceptible genetic regions for rheumatoid arthritis or related diseases, i.e., lupus erythematosus. This topic is discussed in part III.

Cell death in immune thrombocytopenia: novel insights and perspectives

Immune thrombocytopenia (ITP) is a complex disease. The pathogenic and clinical heterogeneity of ITP is reflected by reports on variability in patient history and treatment response, in concert with recent evidence from mechanistic studies. Programmed cell death (PCD) pathways are thought to play a peculiar role in the megakaryocyte lineage in terms of hemostasis and the generation and function of megakaryocytes and platelets; unbalanced genetic or environmental disturbances of these tightly regulated pathways may cause thrombocytopenia. Dysregulated PCD has also been linked to peripheral platelet destruction, intramedullary apoptosis, and inefficient thrombopoiesis in ITP. In this article, we discuss novel and controversial findings on the role of PCD in the megakaryocyte lineage and their potential implications in terms of pathogenesis, diagnosis, and treatment of ITP.

Analysis of Variation in Clubfoot Candidate Genes

Isolated clubfoot, a common birth defect occurring in more than 135,000 livebirths worldwide each year, is associated with significant health care and financial burdens. Clubfoot is defined by forefoot adduction, hindfoot varus, midfoot cavus and hindfoot equinus. Isolated clubfoot, which is the focus of these studies, is distinct from syndromic clubfoot because there are no other associated malformations. Population, family, twin and segregation analysis studies provide evidence that genetic and environmental factors play an etiologic role in isolated clubfoot. The studies described in this thesis were performed to define the role of genetic variation in isolated clubfoot. Interrogation of a deletion region associated with syndromic clubfoot, suggested that CASP8 and CASP10, two apoptotic genes, play a role in isolated clubfoot. To explore the role of apoptotic genes in clubfoot, SNPs spanning genes involved in the apoptotic pathway in the six chromosomal deletion regions, and limb patterning genes, HOXD and HOXA, were interrogated. SNPs in mitochondrial mediated apoptotic genes and several SNPs in HOXA and HOXD genes were modestly associated with clubfoot with the most significant SNP, rs3801776, located in the basal promoter of HOXA9. Several significant associations were found with SNPs in NFAT2 and TNIP2. Significant gene interactions were detected between SNPs in HOX and apoptotic genes. These findings suggest a model for clubfoot in which variation in one gene is not sufficient to cause the malformation but requires variation several genes to perturb protein expression sufficiently to alter muscle and foot development. These results significantly impact our knowledge base by delineating underlying mechanisms causing clubfoot.

Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

MODELING SPORADIC TUMOR FORMATION DRIVEN BY TELOMERE DYSFUNCTION IN THE GASTROINTESTINAL TRACT

Colorectal cancer is a complex disease that is thought to arise when cells accumulate mutations that allow for uncontrolled growth. There are several recognized mechanisms for generating such mutations in sporadic colon cancer; one of which is chromosomal instability (CIN). One hypothesized driver of CIN in cancer is the improper repair of dysfunctional telomeres. Telomeres comprise the linear ends of chromosomes and play a dual role in cancer. Its length is maintained by the ribonucleoprotein, telomerase, which is not a normally expressed in somatic cells and as cells divide, telomeres continuously shorten. Critically shortened telomeres are considered dysfunctional as they are recognized as sites of DNA damage and cells respond by entering into replicative senescence or apoptosis, a process that is p53-dependent and the mechanism for telomere-induced tumor suppression. Loss of this checkpoint and improper repair of dysfunctional telomeres can initiate a cycle of fusion, bridge and breakage that can lead to chromosomal changes and genomic instability, a process that can lead to transformation of normal cells to cancer cells. Mouse models of telomere dysfunction are currently based on knocking out the telomerase protein or RNA component; however, the naturally long telomeres of mice require multiple generational crosses of telomerase null mice to achieve critically short telomeres. Shelterin is a complex of six core proteins that bind to telomeres specifically. Pot1a is a highly conserved member of this complex that specifically binds to the telomeric single-stranded 3’ G-rich overhang. Previous work in our lab has shown that Pot1a is essential for chromosomal end protection as deletion of Pot1a in murine embryonic fibroblasts (MEFs) leads to open telomere ends that initiate a DNA damage response mediated by ATR, resulting in p53-dependent cellular senescence. Loss of Pot1a in the background of p53 deficiency results in increased aberrant homologous recombination at telomeres and elevated genomic instability, which allows Pot1a-/-, p53-/- MEFs to form tumors when injected into SCID mice. These phenotypes are similar to those seen in cells with critically shortened telomeres. In this work, we created a mouse model of telomere ysfunction in the gastrointestinal tract through the conditional deletion of Pot1a that recapitulates the microscopic features seen in severe telomere attrition. Combined intestinal loss of Pot1a and p53 lead to formation of invasive adenocarcinomas in the small and large intestines. The tumors formed with long latency, low multiplicity and had complex genomes due to chromosomal instability, features similar to those seen in sporadic human colorectal cancers. Taken together, we have developed a novel mouse model of intestinal tumorigenesis based on genomic instability driven by telomere dysfunction.

Non-coronary atherosclerosis

During the last decades, the clinical and research interest in atherosclerosis has been mostly focused on coronary arteries. After the publications of the European Society Guidelines and AHA/ACC Guidelines on Peripheral artery diseases, and of the Registry REduction in Atherothrombosis for Continued Health Registry, there has been an increased interest in atherosclerosis of the lower extremity arteries and its presence in multifocal disease. However, awareness in the general population and the medical community of non-coronary artery diseases, and of its major prognostic implications remain relatively low. The aim of this general review stemming out of an ESC Working Group on Peripheral Circulation meeting in 2011 is to enhance awareness of this complex disease highlighting the importance of the involvement of atherosclerosis at different levels with respect to clinical presentation, diagnosis, and co-existence of the disease in the distinct arterial territories. We also emphasize the need of an interdisciplinary approach to face the broad and complex spectrum of multifocal disease, and try to propose a series of tentative recommendations and measures to be implemented in non-coronary atherosclerosis.

Evaluation of candidate genes for susceptibility to spina bifida

Neural tube defects (NTDs) are the most common severely disabling birth defects in the United States, with a frequency of approximately 1–2 of every 1,000 births. This text includes the identification and evaluation of candidate susceptibility genes that confer risk for the development of neural tube defects (NTDs). The project focused on isolated meningomyelocele, also termed spina bifida (SB). ^ Spina bifida is a complex disease with multifactorial inheritance, therefore the subject population (consisting of North American Caucasians and Hispanics of Mexicali-American descent) was composed of 459 simplex SB families who were tested for genetic associations utilizing the transmission disequilibrium test (TDT), a nonparametric linkage technique. Three categories of candidate genes were studied, including (1) human equivalents of genes determined in mouse models to cause NTDs, (2) HOX and PAX genes, and (3) the MTHFR gene involved in the metabolic pathway of folate. ^ The C677T variant of the 5,10-methylenetetrahydrofolate reductase (MTHFR) gene was the first mutation in this gene to be implicated as a risk factor for NTDs. Our evaluation of the MTHFR gene provides evidence that maternal C677T homozygosity is a risk factor for upper level spina bifida defects in Hispanics [OR = 2.3, P = 0.02]. This observed risk factor is of great importance due to the high prevalence of this homozygous genotype in the Hispanic population. Additionally, maternal C677T/A1298C compound heterozygosity is a risk factor for upper level spina bifida defects in non-Hispanic whites [OR = 3.6, P = 0.03]. ^ For TDT analysis, our total population of 1128 subjects were genotyped for 54 markers from within and/or flanking the 20 candidate genes/gene regions of interest. Significant TDT findings were obtained for 3 of the 54 analyzed markers: d20s101 flanking the PAX1 gene (P = 0.019), d1s228 within the PAX7 gene (P = 0.011), and d2s110 within the PAX8 gene (P = 0.013). These results were followed-up by testing the genes directly for mutations utilizing single-strand conformational analysis (SSCA) and direct sequencing. Multiple variations were detected in each of these PAX genes; however, these variations were not passed from parent to child in phase with the positively transmitted alleles. Therefore, these variations do not contribute to the susceptibility of spina bifida, but rather are previously unreported single nucleotide polymorphisms. ^

Delineating the molecular basis of human genetic diseases: Epigenetic and functional studies

Identifying and characterizing the genes responsible for inherited human diseases will ultimately lead to a more holistic understanding of disease pathogenesis, catalyze new diagnostic and treatment modalities, and provide insights into basic biological processes. This dissertation presents research aimed at delineating the genetic and molecular basis of human diseases through epigenetic and functional studies and can be divided into two independent areas of research. The first area of research describes the development of two high-throughput melting curve based methods to assay DNA methylation, referred to as McMSP and McCOBRA. The goal of this project was to develop DNA methylation methods that can be used to rapidly determine the DNA methylation status at a specific locus in a large number of samples. McMSP and McCOBRA provide several advantages over existing methods, as they are simple, accurate, robust, and high-throughput making them applicable to large-scale DNA methylation studies. McMSP and McCOBRA were then used in an epigenetic study of the complex disease Ankylosing spondylitis (AS). Specifically, I tested the hypothesis that aberrant patterns of DNA methylation in five AS candidate genes contribute to disease susceptibility. While no statistically significant methylation differences were observed between cases and controls, this is the first study to investigate the hypothesis that epigenetic variation contributes to AS susceptibility and therefore provides the conceptual framework for future studies. ^ In the second area of research, I performed experiments to better delimit the function of aryl hydrocarbon receptor-interacting protein-like 1 (AIPL1), which when mutated causes various forms of inherited blindness such as Leber congenital amaurosis. A yeast two-hybrid screen was performed to identify putative AIPL1-interacting proteins. After screening 2 × 106 bovine retinal cDNA library clones, 6 unique putative AIPL1-interacting proteins were identified. While these 6 AIPL1 protein-protein interactions must be confirmed, their identification is an important step in understanding the functional role of AIPL1 within the retina and will provide insight into the molecular mechanisms underlying inherited blindness. ^

Candidate gene identification following linkage: Search for hypertension susceptibility genes

Following up genetic linkage studies to identify the underlying susceptibility gene(s) for complex disease traits is an arduous yet biologically and clinically important task. Complex traits, such as hypertension, are considered polygenic with many genes influencing risk, each with small effects. Chromosome 2 has been consistently identified as a genomic region with genetic linkage evidence suggesting that one or more loci contribute to blood pressure levels and hypertension status. Using combined positional candidate gene methods, the Family Blood Pressure Program has concentrated efforts in investigating this region of chromosome 2 in an effort to identify underlying candidate hypertension susceptibility gene(s). Initial informatics efforts identified the boundaries of the region and the known genes within it. A total of 82 polymorphic sites in eight positional candidate genes were genotyped in a large hypothesis-generating sample consisting of 1640 African Americans, 1339 whites, and 1616 Mexican Americans. To adjust for multiple comparisons, resampling-based false discovery adjustment was applied, extending traditional resampling methods to sibship samples. Following this adjustment for multiple comparisons, SLC4A5, a sodium bicarbonate transporter, was identified as a primary candidate gene for hypertension. Polymorphisms in SLC4A5 were subsequently genotyped and analyzed for validation in two populations of African Americans (N = 461; N = 778) and two of whites (N = 550; N = 967). Again, SNPs within SLC4A5 were significantly associated with blood pressure levels and hypertension status. While not identifying a single causal DNA sequence variation that is significantly associated with blood pressure levels and hypertension status across all samples, the results further implicate SLC4A5 as a candidate hypertension susceptibility gene, validating previous evidence for one or more genes on chromosome 2 that influence hypertension related phenotypes in the population-at-large. The methodology and results reported provide a case study of one approach for following up the results of genetic linkage analyses to identify genes influencing complex traits. ^

Expression study of the atherosclerosis mouse aorta reveals significant disturbance of calcium signaling pathway

Atherosclerosis is a complex disease resulting from interactions of genetic and environmental risk factors leading to heart failure and stroke. Using an atherosclerotic mouse model (ldlr-/-, apobec1-/- designated as LDb), we performed microarray analysis to identify candidate genes and pathways, which are most perturbed in changes in the following risk factors: genetics (control C57BL/6 vs. LDb mice), shearstress (lesion-prone vs. lesion-resistant regions in LDb mice), diet (chow vs. high fat fed LDb mice) and age (2-month-old vs. 8-month old LDb mice). ^ Atherosclerotic lesion quantification and lipid profile studies were performed to assess the disease phenotype. A microarray study was performed on lesion-prone and lesion-resistant regions of each aorta. Briefly, 32 male C57BL/6 and LDb mice (n =16/each) were fed on either chow or high fat diet, sacrificed at 2- and 8-months old, and RNA isolated from the aortic lesion-prone and aortic lesion-resistant segments. Using 64 Affymetrix Murine 430 2.0 chips, we profiled differentially expressed genes with the cut off value of FDR ≤ 0.15 for t-test, and q <0.0001 for the ANOVA. The data were normalized using two normalization methods---invariant probe sets (Loess) and Quantile normalization, the statistical analysis was performed using t-tests and ANOVA, and pathway characterization was done using Pathway Express (Wayne State). The result identified the calcium signaling pathway as the most significant overrepresented pathway, followed by focal adhesion. In the calcium signaling pathway, 56 genes were found to be significantly differentially expressed out of 180 genes listed in the KEGG calcium signaling pathway. Nineteen of these genes were consistently identified by both statistical tests, 11 of which were unique to the test, and 26 were unique to the ANOVA test, using the cutoffs noted above. ^ In conclusion, this finding suggested that hypercholesterolemia drives the disease progression by altering the expression of calcium channels and regulators which subsequently results in cell differentiation, growth, adhesion, cytoskeletal change and death. Clinically, this pathway may serve as an important target for future therapeutic intervention, and thus the calcium signaling pathway may serve as an important target for future diagnostic and therapeutic intervention. ^

Single nucleotide polymorphisms (SNPs) associated with TGF-beta pathway and their significance in systemic sclerosis - A multilevel analysis

Systemic sclerosis (SSc) or Scleroderma is a complex disease and its etiopathogenesis remains unelucidated. Fibrosis in multiple organs is a key feature of SSc and studies have shown that transforming growth factor-β (TGF-β) pathway has a crucial role in fibrotic responses. For a complex disease such as SSc, expression quantitative trait loci (eQTL) analysis is a powerful tool for identifying genetic variations that affect expression of genes involved in this disease. In this study, a multilevel model is described to perform a multivariate eQTL for identifying genetic variation (SNPs) specifically associated with the expression of three members of TGF-β pathway, CTGF, SPARC and COL3A1. The uniqueness of this model is that all three genes were included in one model, rather than one gene being examined at a time. A protein might contribute to multiple pathways and this approach allows the identification of important genetic variations linked to multiple genes belonging to the same pathway. In this study, 29 SNPs were identified and 16 of them located in known genes. Exploring the roles of these genes in TGF-β regulation will help elucidate the etiology of SSc, which will in turn help to better manage this complex disease. ^

Biomedical informatics applications for epilepsy management: A systematic review

Epilepsy is a very complex disease which can have a variety of etiologies, co-morbidities, and a long list of psychosocial factors4. Clinical management of epilepsy patients typically includes serological tests, EEG's, and imaging studies to determine the single best antiepileptic drug (AED). Self-management is a vital component of achieving optimal health when living with a chronic disease. For patients with epilepsy self-management includes any necessary actions to control seizures and cope with any subsequent effects of the condition9; including aspects of treatment, seizure, and lifestyle. The use of computer-based applications can allow for more effective use of clinic visits and ultimately enhance the patient-provider relationship through focused discussion of determinants affecting self-management. ^ The purpose of this study is to conduct a systematic literature review on informatics application in epilepsy self-management in an effort to describe current evidence for informatics applications and decision support as an adjunct to successful clinical management of epilepsy. Each publication was analyzed for the type of study design utilized. ^ A total of 68 publications were included and categorized by the study design used, development stage, and clinical domain. Descriptive study designs comprised of three-fourths of the publications and indicate an underwhelming use of prospective studies. The vast majority of prospective studies also focused on clinician use to increase knowledge in treating patients with epilepsy. ^ Due to the chronic nature of epilepsy and the difficulty that both clinicians and patients can experience in managing epilepsy, more prospective studies are needed to evaluate applications that can effectively increase management activities. Within the last two decades of epilepsy research, management studies have employed the use of biomedical informatics applications. While the use of computer applications to manage epilepsy has increased, more progress is needed.^

A new class of obesity genes encodes leukocyte adhesion receptors

Obesity is a complex disease, and multiple genes contribute to the trait. The description of five genes (ob, db, tub, Ay, and fat) responsible for distinct syndromes of spontaneous monogenic obesity in mice has advanced our knowledge of the genetics of obesity. However, many other genes involved in the expression of this disease remain to be determined. We report here the identification of an additional class of genes involved in the regulation of adipose tissue mass. These genes encode receptors mediating leukocyte adhesion. Mice deficient in intercellular adhesion molecule-1 became spontaneously obese in old age on normal mouse chow or at a young age when provided with a diet rich in fat. Mice deficient in the counterreceptor for intercellular adhesion molecule-1, the leukocyte integrin αMβ2 (Mac-1), showed a similar obesity phenotype. Since all mice consumed approximately the same amount of food as controls, the leukocyte function appears to be in regulating lipid metabolism and/or energy expenditure. Our results indicate that (i) leukocytes play a role in preventing excess body fat deposition and (ii) defects in leukocyte adhesion receptors can result in obesity.

Localization of atherosclerosis susceptibility loci to chromosomes 4 and 6 using the Ldlr knockout mouse model

Atherosclerosis is a complex disease resulting from the interaction of multiple genes. We have used the Ldlr knockout mouse model in an interspecific genetic cross to map atherosclerosis susceptibility loci. A total of 174 (MOLF/Ei × B6.129S7-Ldlrtm1Her) × C57BL/6J-Ldlrtm1Her backcross mice, homozygous for the Ldlr null allele, were fed a Western-type diet for 3 months and then killed for quantification of aortic lesions. A genome scan was carried out by using DNA pools and microsatellite markers spaced at ≈18-centimorgan intervals. Quantitative trait locus analysis of individual backcross mice confirmed linkages to chromosomes 4 (Athsq1, logarithm of odds = 6.2) and 6 (Athsq2, logarithm of odds = 6.7). Athsq1 affected lesions in females only whereas Athsq2 affected both sexes. Among females, the loci accounted for ≈50% of the total variance of lesion area. The susceptible allele at Athsq1 was derived from the MOLF/Ei genome whereas the susceptible allele at Athsq2 was derived from C57BL/6J. Inheritance of susceptible alleles at both loci conferred a 2-fold difference in lesion area, suggesting an additive effect of Athsq1 and Athsq2. No associations were observed between the quantitative trait loci and levels of plasma total cholesterol, high density lipoprotein cholesterol, non-high density lipoprotein cholesterol, insulin, or body weight. We provide strong evidence for complex inheritance of atherosclerosis in mice with elevated plasma low density lipoprotein cholesterol and show a major influence of nonlipoprotein-related factors on disease susceptibility. Athsq1 and Athsq2 represent candidate susceptibility loci for human atherosclerosis, most likely residing on chromosomes 1p36–32 and 12p13–12, respectively.

Stress management in mothers of children with Autism Spectrum Disorder : psychological characteristics, coping strategies and biological correlates

Un figlio con un Disturbo dello Spettro Autistico, caratterizzato da gravi difficoltà nelle relazioni, nei comportamenti e nella comunicazione, costringe tutto il sistema familiare a gestire un notevole stress dovuto alla gestione quotidiana di una patologia così complessa. Per questi motivi, i genitori necessitano di un sostegno il più possibile personalizzato rispetto alle caratteristiche del loro contesto familiare. Per fare questo sarebbe importante individuare quali siano i parametri correlati ai livelli di stress nei familiari di pazienti con autismo e che potrebbero avere un’influenza sul benessere familiare. Lo scopo di questo studio è quello di valutare quali caratteristiche di personalità, stili di coping e capacità di gestire le emozioni possano essere in relazione con la reattività individuale alle situazioni di stress, valutata attraverso alcuni correlati biologici, quali il livello di cortisolo (l’ormone dello stress) e la variabilità della frequenza cardiaca. L’ottica di ricerca applicata fa sì che gli obiettivi ultimi di questo lavoro siano anche quelli di diminuire l’accesso ai servizi per questi soggetti, considerando il fatto che progetti individualizzati di sostegno genitoriale costituiscono un fattore protettivo rispetto a conseguenze fisiche e psicologiche di disagio se implementati tenendo conto della variabilità individuale rispetto alle caratteristiche sopra citate.