Bayesian methodology for genetics of complex diseases

Genetic research of complex diseases is a challenging, but exciting, area of research. The early development of the research was limited, however, until the completion of the Human Genome and HapMap projects, along with the reduction in the cost of genotyping, which paves the way for understanding the genetic composition of complex diseases. In this thesis, we focus on the statistical methods for two aspects of genetic research: phenotype definition for diseases with complex etiology and methods for identifying potentially associated Single Nucleotide Polymorphisms (SNPs) and SNP-SNP interactions. With regard to phenotype definition for diseases with complex etiology, we firstly investigated the effects of different statistical phenotyping approaches on the subsequent analysis. In light of the findings, and the difficulties in validating the estimated phenotype, we proposed two different methods for reconciling phenotypes of different models using Bayesian model averaging as a coherent mechanism for accounting for model uncertainty. In the second part of the thesis, the focus is turned to the methods for identifying associated SNPs and SNP interactions. We review the use of Bayesian logistic regression with variable selection for SNP identification and extended the model for detecting the interaction effects for population based case-control studies. In this part of study, we also develop a machine learning algorithm to cope with the large scale data analysis, namely modified Logic Regression with Genetic Program (MLR-GEP), which is then compared with the Bayesian model, Random Forests and other variants of logic regression.

Genome-wide interaction analysis reveals replicated epistatic effects on brain structure

The discovery of several genes that affect the risk for Alzheimer's disease ignited a worldwide search for single-nucleotide polymorphisms (SNPs), common genetic variants that affect the brain. Genome-wide search of all possible SNP-SNP interactions is challenging and rarely attempted because of the complexity of conducting approximately 1011 pairwise statistical tests. However, recent advances in machine learning, for example, iterative sure independence screening, make it possible to analyze data sets with vastly more predictors than observations. Using an implementation of the sure independence screening algorithm (called EPISIS), we performed a genome-wide interaction analysis testing all possible SNP-SNP interactions affecting regional brain volumes measured on magnetic resonance imaging and mapped using tensor-based morphometry. We identified a significant SNP-SNP interaction between rs1345203 and rs1213205 that explains 1.9% of the variance in temporal lobe volume. We mapped the whole brain, voxelwise effects of the interaction in the Alzheimer's Disease Neuroimaging Initiative data set and separately in an independent replication data set of healthy twins (Queensland Twin Imaging). Each additional loading in the interaction effect was associated with approximately 5% greater brain regional brain volume (a protective effect) in both Alzheimer's Disease Neuroimaging Initiative and Queensland Twin Imaging samples.

Genetic determinants of bone density and fracture risk: State of the art and future directions

Context: Osteoporosis is a common, highly heritable condition that causes substantial morbidity and mortality, the etiopathogenesis of which is poorly understood. Genetic studies are making increasingly rapid progress in identifying the genes involved. Evidence Acquisition and Synthesis: In this review, we will summarize the current understanding of the genetics of osteoporosis based on publications from PubMed from the year 1987 onward. Conclusions: Most genes involved in osteoporosis identified to date encode components of known pathways involved in bone synthesis or resorption, but as the field progresses, new pathways are being identified. Only a small proportion of the total genetic variation involved in osteoporosis has been identified, and new approaches will be required to identify most of the remaining genes.

Functional evaluation of genetic variants associated with endometriosis near GREB1

STUDY QUESTION: Do DNA variants in the growth regulation by estrogen in breast cancer 1 (GREB1) region regulate endometrial GREB1 expression and increase the risk of developing endometriosis in women? SUMMARY ANSWER: We identified new single nucleotide polymorphisms (SNPs) with strong association with endometriosis at the GREB1 locus although we did not detect altered GREB1 expression in endometriosis patients with defined genotypes. WHAT IS ALREADY KNOWN: Genome-wide association studies have identified the GREB1 region on chromosome 2p25.1 for increasing endometriosis risk. The differential expression of GREB1 has also been reported by others in association with endometriosis disease phenotype. STUDY DESIGN, SIZE, DURATION: Fine mapping studies comprehensively evaluated SNPs within the GREB1 region in a large-scale data set (>2500 cases and >4000 controls). Publicly available bioinformatics tools were employed to functionally annotate SNPs showing the strongest association signal with endometriosis risk. Endometrial GREB1 mRNA and protein expression was studied with respect to phases of the menstrual cycle (n = 2-45 per cycle stage) and expression quantitative trait loci (eQTL) analysis for significant SNPs were undertaken for GREB1 [mRNA (n = 94) and protein (n = 44) in endometrium]. PARTICIPANTS/MATERIALS, SETTING, METHODS: Participants in this study are females who provided blood and/or endometrial tissue samples in a hospital setting. The key SNPs were genotyped using Sequenom MassARRAY. The functional roles and regulatory annotations for identified SNPs are predicted by various publicly available bioinformatics tools. Endometrial GREB1 expression work employed qRT-PCR, western blotting and immunohistochemistry studies. MAIN RESULTS AND THE ROLE OF CHANCE: Fine mapping results identified a number of SNPs showing stronger association (0.004 < P < 0.032) with endometriosis risk than the original GWAS SNP (rs13394619) (P = 0.034). Some of these SNPs were predicted to have functional roles, for example, interaction with transcription factor motifs. The haplotype (a combination of alleles) formed by the risk alleles from two common SNPs showed significant association (P = 0.026) with endometriosis and epistasis analysis showed no evidence for interaction between the two SNPs, suggesting an additive effect of SNPs on endometriosis risk. In normal human endometrium, GREB1 protein expression was altered depending on the cycle stage (significantly different in late proliferative versus late secretory, P < 0.05) and cell type (glandular epithelium, not stromal cells). However, GREB1 expression in endometriosis cases versus controls and eQTL analyses did not reveal any significant changes. LIMITATIONS, REASONS FOR CAUTION: In silico prediction tools are generally based on cell lines different to our tissue and disease of interest. Functional annotations drawn from these analyses should be considered with this limitation in mind. We identified cell-specific and hormone-specific changes in GREB1 protein expression. The lack of a significant difference observed following our GREB1 expression studies may be the result of moderate power on mixed cell populations in the endometrial tissue samples. WIDER IMPLICATIONS OF THE FINDINGS: This study further implicates the GREB1 region on chromosome 2p25.1 and the GREB1 gene with involvement in endometriosis risk. More detailed functional studies are required to determine the role of the novel GREB1 transcripts in endometriosis pathophysiology. STUDY FUNDING/COMPETING INTERESTS: Funding for this work was provided by NHMRC Project Grants APP1012245, APP1026033, APP1049472 and APP1046880. There are no competing interests.

Novel Molecular Mechanisms of Arabidopsis Disease Resistance

Plants are capable of recognizing phytopathogens through the perception of pathogen-derived molecules or plant cell-wall degradation products due to the activities of pathogen-secreted enzymes. Such elicitor recognition events trigger an array of inducible defense responses involving signal transduction networks and massive transcriptional re-programming. The outcome of a pathogen infection relies on the balance between different signaling pathways, which are integrated by regulatory proteins. This thesis characterized two key regulatory components: a damage control enzyme, chlorophyllase 1 (AtCHL1), and a transcription factor, WRKY70. Their roles in defense signaling were then investigated. The Erwinia-derived elicitors rapidly activated the expression of AtCLH1 and WRKY70 through different signaling pathways. The expression of the AtCHL1 gene was up-regulated by jasmonic acid (JA) but down-regulated by salicylic acid (SA), whereas WRKY70 was activated by SA and repressed by JA. In order to elucidate the functions of AtCLH1 and WRKY70 in plant defense, stable transgenic lines were produced where these genes were overexpressed or silenced. Additionally, independent knockout lines were also characterized. Bacterial and fungal pathogens were then used to assess the contribution of these genes to the Arabidopsis disease resistance. The transcriptional modulation of AtCLH1 by either the constitutive over-expression or RNAi silencing caused alterations in the chlorophyll-to-chlorophyllide ratio, supporting the claim that chlorophyllase 1 has a role in the chlorophyll degradation pathway. Silencing of this gene led to light-dependent over-accumulation of the reactive oxygen species (ROS) in response to infection by Erwinia carotovora subsp. carotovora SCC1. This was followed by an enhanced induction of SA-dependent defense genes and an increased resistance to this pathogen. Interestingly, little effect on the pathogen-induced SA accumulation at the early infection was observed, suggesting that action of ROS might potentiate SA signaling. In contrast, the pathogen-induced JA production was significantly reduced in the RNAi silenced plants. Moreover, JA signaling and resistance to Alternaria brassicicola were impaired. These observations provide support for the argument that the ROS generated in chloroplasts might have a negative impact on JA signaling. The over-expression of WRKY70 resulted in an enhanced resistance to E. carotovora subsp. carotovora SCC1, Pseudomonas syringae pv. tomato DC3000 and Erysiphe cichoracearum UCSC1, whilst an antisense suppression or an insertional inactivation of WRKY70 led to a compromised resistance to E. carotovora subsp. carotovora SCC1 and to E. cichoracearum UCSC1 but not to P. syringae pv. tomato DC3000. Gene expression analysis revealed that WRKY70 activated many known defense-related genes associated with the SAR response but suppressed a subset of the JA-responsive genes. In particular, I was able to show that both the basal and the induced expression of AtCLH1 was enhanced by the antisense silencing or the insertional inactivation of WRKY70, whereas a reduction in AtCLH1 expression was observed in the WRKY70 over-expressors following an MeJA application or an A. brassicicola infection. Moreover, the SA-induced suppression of AtCLH1 was relieved in wrky70 mutants. These results indicate that WRKY70 down-regulates AtCLH1. An epistasis analysis suggested that WRKY70 functions downstream of the NPR1 in an SA-dependent signaling pathway. When challenged with A. brassicicola, WRKY70 over-expressing plants exhibited a compromised disease resistance while wrky70 mutants had the opposite effect. These results confirmed the WRKY70-mediated inhibitory effects on JA signaling. Furthermore, the WRKY70-controlled suppression of A. brassicicola resistance was mainly through an NPR1-dependent mechanism. Taking all the data together, I suggest that the pathogen-responsive transcription factor WRKY70 is a common component in both SA- and JA-dependent pathways and plays a crucial role in the SA-mediated suppression of JA signaling.

Variational principles in evolution

For a one-locus selection model, Svirezhev introduced an integral variational principle by defining a Lagrangian which remained stationary on the trajectory followed by the population undergoing selection. It is shown here (i) that this principle can be extended to multiple loci in some simple cases and (ii) that the Lagrangian is defined by a straightforward generalization of the one-locus case, but (iii) that in two-locus or more general models there is no straightforward extension of this principle if linkage and epistasis are present. The population trajectories can be constructed as trajectories of steepest ascent in a Riemannian metric space. A general method is formulated to find the metric tensor and the surface-in the metric space on which the trajectories, which characterize the variations in the gene structure of the population, lie. The local optimality principle holds good in such a space. In the special case when all possible linkage disequilibria are zero, the phase point of the n-locus genetic system moves on the surface of the product space of n higher dimensional unit spheres in a certain Riemannian metric space of gene frequencies so that the rate of change of mean fitness is maximum along the trajectory. In the two-locus case the corresponding surface is a hyper-torus.

The LIN-3 gene of the nematode C. elegans is the vulval-inducing signal

The development of the vulva of the nematode Caenorhabditis elegans is induced by a signal from the anchor cell of the somatic gonad. Activity of the gene lin-3 is required for the Vulval Precursor Cells (VPCs) to assume vulval fates. It is shown here that lin-3 encodes the vulval-inducing signal.

lin-3 was molecularly cloned by transposon-tagging and shown to encode a nematode member ofthe Epidermal Growth Factor (EGF) family. Genetic epistasis experiments indicate that lin-3 acts upstream of let-23, which encodes a homologue of the EGF-Receptor.

lin-3 transgenes that contain multiple copies of wild-type lin-3 genomic DNA clones confer a dominant multivulva phenotype in which up to all six of the VPCs assume vulval fates. The properties of these trans genes suggest that lin-3 can act in the anchor cell to induce vulval fates. Ablation of the gonadal precursors, which prevents the development of the AC, strongly reduces the ability of lin-3 transgenes to stimulate vulval development. A lin-3 recorder transgene that retains the ability to stimulate vulval development is expressed specifically in the anchor cell at the time of vulval induction.

Expression of an obligate secreted form of the EGF domain of Lin-S from a heterologous promoter is sufficient to induce vulval fates in the absence of the normal source of the inductive signal. This result suggests that Lin-S may act as a secreted factor, and that Lin-S may be the sole vulval-inducing signal made by the anchor cell.

lin-3 transgenes can cause adjacent VPCs to assume the 1° vulval fate and thus can override the action of the lateral signal mediated by lin-12 that normally prevents adjacent 1° fates. This indicates that the production of Lin-3 by the anchor cell must be limited to allow the VPCs to assume the proper pattern of fates of so 3° 3° 2° 1° 2° 3°.

Estudio del receptor 5HT2A y de la vía Akt/GSK3 como mecanismos moleculares de la psicosis inducida por el consumo crónico de cannabis

La esquizofrenia es una enfermedad mental de carácter crónico que afecta aproximadamente al 1% de la población , principalmente con edades comprendidas entre los 15 y 4 5 años . Se trata de una enfermedad de inicio en la adolescencia o comienzo de la madurez y cuyo potencial discapacitante persiste y se agrava a lo largo de la vida. La etiología de la esquizofrenia es multifactorial con evidencias de factores genéticos y ambientales. El carácter hereditario de la esquizofr e nia se ha estimado en un 73 - 90% y los factores genéticos que predisponen a la enfermedad son múltiples y heterogéneos. Los estudios epidemiológicos parecen sugerir que el mayor riesgo de desarrollar esquizofrenia es prod ucto de interacciones genético/ ambi entales (predisposición genética a agresiones ambientales) y del fenómeno de epistasis (fenotipo dependiente de la interacción entre diversos genes). Las manifestaciones clínicas de la esquizofrenia se dividen en síntomas positivos, síntomas negativos y dé ficits cognitivos. A causa de los déficits asociados y de su carácter crónico, la esquizofrenia se encuentra entre las diez causas principales de discapacidad por enfermedad del mundo y según la World Health Organization (2001), se estima que es la quinta enfermedad más costosa para la s ociedad en términos de atención requ erida y pérdida de productividad, con un coste anual en la Unión Europea que supera los 35 mil millones de €, según Andlin - Sobocki y Rössler (2005) . A nivel nacional, en el año 2009 los an tipsicóticos atípicos más empleados (risperidona y olanzapina) se encontraron entre los 8 primeros fármacos que más gasto generaron, con un coste económico superior a 350 millones de euros (Sistema Nacional de Salud (2010)) . Además, según el estudio de Pal mer et al. (2005), el riesgo de muerte prematura en la población esquizofrénica es aproximadamente dos veces mayor en la población general, siendo el suicidio la principal causa de este exceso de mortalidad, con una prevalencia estimada entre 5 - 10%. En est e sentido, se ha estimado que entre el 25% y 50% de los pacientes con esquizofrenia cometen por lo menos un intento de suicidio en su vida , según el estudio llevado a cabo por Meltzer (2002) . Dada su prevalencia, la tendencia a la cronicidad y el riesgo su icida de pacientes con esquizofrenia así como los elevados costes socio - sanitarios, la investigación de esta enfermedad es un objetivo prioritario de los sistemas de salud mundiales

Cytokinesis proteins Tum and Pav have a nuclear role in Wnt regulation.

Wg/Wnt signals specify cell fates in both invertebrate and vertebrate embryos and maintain stem-cell populations in many adult tissues. Deregulation of the Wnt pathway can transform cells to a proliferative fate, leading to cancer. We have discovered that two Drosophila proteins that are crucial for cytokinesis have a second, largely independent, role in restricting activity of the Wnt pathway. The fly homolog of RacGAP1, Tumbleweed (Tum)/RacGAP50C, and its binding partner, the kinesin-like protein Pavarotti (Pav), negatively regulate Wnt activity in fly embryos and in cultured mammalian cells. Unlike many known regulators of the Wnt pathway, these molecules do not affect stabilization of Arm/beta-catenin (betacat), the principal effector molecule in Wnt signal transduction. Rather, they appear to act downstream of betacat stabilization to control target-gene transcription. Both Tum and Pav accumulate in the nuclei of interphase cells, a location that is spatially distinct from their cleavage-furrow localization during cytokinesis. We show that this nuclear localization is essential for their role in Wnt regulation. Thus, we have identified two modulators of the Wnt pathway that have shared functions in cell division, which hints at a possible link between cytokinesis and Wnt activity during tumorigenesis.

Environmental and genetic determinants of colony morphology in yeast.

Nutrient stresses trigger a variety of developmental switches in the budding yeast Saccharomyces cerevisiae. One of the least understood of such responses is the development of complex colony morphology, characterized by intricate, organized, and strain-specific patterns of colony growth and architecture. The genetic bases of this phenotype and the key environmental signals involved in its induction have heretofore remained poorly understood. By surveying multiple strain backgrounds and a large number of growth conditions, we show that limitation for fermentable carbon sources coupled with a rich nitrogen source is the primary trigger for the colony morphology response in budding yeast. Using knockout mutants and transposon-mediated mutagenesis, we demonstrate that two key signaling networks regulating this response are the filamentous growth MAP kinase cascade and the Ras-cAMP-PKA pathway. We further show synergistic epistasis between Rim15, a kinase involved in integration of nutrient signals, and other genes in these pathways. Ploidy, mating-type, and genotype-by-environment interactions also appear to play a role in the controlling colony morphology. Our study highlights the high degree of network reuse in this model eukaryote; yeast use the same core signaling pathways in multiple contexts to integrate information about environmental and physiological states and generate diverse developmental outputs.

Assessment of LD matrix measures for the analysis of biological pathway association.

Complex diseases will have multiple functional sites, and it will be invaluable to understand the cross-locus interaction in terms of linkage disequilibrium (LD) between those sites (epistasis) in addition to the haplotype-LD effects. We investigated the statistical properties of a class of matrix-based statistics to assess this epistasis. These statistical methods include two LD contrast tests (Zaykin et al., 2006) and partial least squares regression (Wang et al., 2008). To estimate Type 1 error rates and power, we simulated multiple two-variant disease models using the SIMLA software package. SIMLA allows for the joint action of up to two disease genes in the simulated data with all possible multiplicative interaction effects between them. Our goal was to detect an interaction between multiple disease-causing variants by means of their linkage disequilibrium (LD) patterns with other markers. We measured the effects of marginal disease effect size, haplotype LD, disease prevalence and minor allele frequency have on cross-locus interaction (epistasis). In the setting of strong allele effects and strong interaction, the correlation between the two disease genes was weak (r=0.2). In a complex system with multiple correlations (both marginal and interaction), it was difficult to determine the source of a significant result. Despite these complications, the partial least squares and modified LD contrast methods maintained adequate power to detect the epistatic effects; however, for many of the analyses we often could not separate interaction from a strong marginal effect. While we did not exhaust the entire parameter space of possible models, we do provide guidance on the effects that population parameters have on cross-locus interaction.

Interaction between interleukin 3 and dystrobrevin-binding protein 1 in schizophrenia

Schizophrenia is a common psychotic mental disorder that is believed to result from the effects of multiple genetic and environmental factors. In this study, we explored gene-gene interactions and main effects in both case-control (657 cases and 411 controls) and family-based (273 families, 1350 subjects) datasets of English or Irish ancestry. Fifty three markers in 8 genes were genotyped in the family sample and 44 markers in 7 genes were genotyped in the case-control sample. The Multifactor Dimensionality Reduction Pedigree Disequilibrium Test (MDR-PDT) was used to examine epistasis in the family dataset and a 3-locus model was identified (permuted p=0.003). The 3-locus model involved the IL3 (rs2069803), RGS4 (rs2661319), and DTNBP1 (rs21319539) genes. We used MDR to analyze the case-control dataset containing the same markers typed in the RGS4, IL3 and DTNBP1 genes and found evidence of a joint effect between IL3 (rs31400) and DTNBP1 (rs760761) (cross-validation consistency 4/5, balanced prediction accuracy=56.84%, p=0.019). While this is not a direct replication, the results obtained from both the family and case-control samples collectively suggest that IL3 and DTNBP1 are likely to interact and jointly contribute to increase risk for schizophrenia. We also observed a significant main effect in DTNBP1, which survived correction for multiple comparisons, and numerous nominally significant effects in several genes. (C) 2008 Elsevier B.V. All rights reserved.

Genetic mapping of the rice ionome in leaves and grain:Identification of QTLs for 17 elements including arsenic, cadmium, iron and selenium

Research into the composition of cereal grains is motivated by increased interest in food quality. Here multi-element analysis is conducted on leaves and grain of the Bala x Azucena rice mapping population grown in the field. Quantitative trait loci (QTLs) for the concentration of 17 elements were detected, revealing 36 QTLs for leaves and 41 for grains. Epistasis was detected for most elements. There was very little correlation between leaf and grain element concentrations. For selenium, lead, phosphorus and magnesium QTLs were detected in the same location for both tissues. In general, there were no major QTL clusters, suggesting separate regulation of each element. QTLs for grain iron, zinc, molybdenum and selenium are potential targets for marker assisted selection to improve seed nutritional quality. An epistatic interaction for grain arsenic also looks promising to decrease the concentration of this carcinogenic element. © Springer Science+Business Media B.V. 2009.

The protein phosphatase 7 regulates phytochrome signaling in Arabidopsis.

The psi2 mutant of Arabidopsis displays amplification of the responses controlled by the red/far red light photoreceptors phytochrome A (phyA) and phytochrome B (phyB) but no apparent defect in blue light perception. We found that loss-of-function alleles of the protein phosphatase 7 (AtPP7) are responsible for the light hypersensitivity in psi2 demonstrating that AtPP7 controls the levels of phytochrome signaling. Plants expressing reduced levels of AtPP7 mRNA display reduced blue-light induced cryptochrome signaling but no noticeable deficiency in phytochrome signaling. Our genetic analysis suggests that phytochrome signaling is enhanced in the AtPP7 loss of function alleles, including in blue light, which masks the reduced cryptochrome signaling. AtPP7 has been found to interact both in yeast and in planta assays with nucleotide-diphosphate kinase 2 (NDPK2), a positive regulator of phytochrome signals. Analysis of ndpk2-psi2 double mutants suggests that NDPK2 plays a critical role in the AtPP7 regulation of the phytochrome pathway and identifies NDPK2 as an upstream element involved in the modulation of the salicylic acid (SA)-dependent defense pathway by light. Thus, cryptochrome- and phytochrome-specific light signals synchronously control their relative contribution to the regulation of plant development. Interestingly, PP7 and NDPK are also components of animal light signaling systems.

Utilisation de triades cas-parents dans la régression logique : exploration d'interaction génétique

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal