THE NATURAL AND ORTHOGONAL INTERACTION (NOIA) MODELS FOR QUANTITATIVE TRAITS (QTs) AND COMPLEX DISEASES

My dissertation focuses on developing methods for gene-gene/environment interactions and imprinting effect detections for human complex diseases and quantitative traits. It includes three sections: (1) generalizing the Natural and Orthogonal interaction (NOIA) model for the coding technique originally developed for gene-gene (GxG) interaction and also to reduced models; (2) developing a novel statistical approach that allows for modeling gene-environment (GxE) interactions influencing disease risk, and (3) developing a statistical approach for modeling genetic variants displaying parent-of-origin effects (POEs), such as imprinting. In the past decade, genetic researchers have identified a large number of causal variants for human genetic diseases and traits by single-locus analysis, and interaction has now become a hot topic in the effort to search for the complex network between multiple genes or environmental exposures contributing to the outcome. Epistasis, also known as gene-gene interaction is the departure from additive genetic effects from several genes to a trait, which means that the same alleles of one gene could display different genetic effects under different genetic backgrounds. In this study, we propose to implement the NOIA model for association studies along with interaction for human complex traits and diseases. We compare the performance of the new statistical models we developed and the usual functional model by both simulation study and real data analysis. Both simulation and real data analysis revealed higher power of the NOIA GxG interaction model for detecting both main genetic effects and interaction effects. Through application on a melanoma dataset, we confirmed the previously identified significant regions for melanoma risk at 15q13.1, 16q24.3 and 9p21.3. We also identified potential interactions with these significant regions that contribute to melanoma risk. Based on the NOIA model, we developed a novel statistical approach that allows us to model effects from a genetic factor and binary environmental exposure that are jointly influencing disease risk. Both simulation and real data analyses revealed higher power of the NOIA model for detecting both main genetic effects and interaction effects for both quantitative and binary traits. We also found that estimates of the parameters from logistic regression for binary traits are no longer statistically uncorrelated under the alternative model when there is an association. Applying our novel approach to a lung cancer dataset, we confirmed four SNPs in 5p15 and 15q25 region to be significantly associated with lung cancer risk in Caucasians population: rs2736100, rs402710, rs16969968 and rs8034191. We also validated that rs16969968 and rs8034191 in 15q25 region are significantly interacting with smoking in Caucasian population. Our approach identified the potential interactions of SNP rs2256543 in 6p21 with smoking on contributing to lung cancer risk. Genetic imprinting is the most well-known cause for parent-of-origin effect (POE) whereby a gene is differentially expressed depending on the parental origin of the same alleles. Genetic imprinting affects several human disorders, including diabetes, breast cancer, alcoholism, and obesity. This phenomenon has been shown to be important for normal embryonic development in mammals. Traditional association approaches ignore this important genetic phenomenon. In this study, we propose a NOIA framework for a single locus association study that estimates both main allelic effects and POEs. We develop statistical (Stat-POE) and functional (Func-POE) models, and demonstrate conditions for orthogonality of the Stat-POE model. We conducted simulations for both quantitative and qualitative traits to evaluate the performance of the statistical and functional models with different levels of POEs. Our results showed that the newly proposed Stat-POE model, which ensures orthogonality of variance components if Hardy-Weinberg Equilibrium (HWE) or equal minor and major allele frequencies is satisfied, had greater power for detecting the main allelic additive effect than a Func-POE model, which codes according to allelic substitutions, for both quantitative and qualitative traits. The power for detecting the POE was the same for the Stat-POE and Func-POE models under HWE for quantitative traits.

BAYESIAN STATISTICAL METHODS IN GENE-ENVIRONMENT AND GENE-GENE INTERACTION STUDIES

Complex diseases such as cancer result from multiple genetic changes and environmental exposures. Due to the rapid development of genotyping and sequencing technologies, we are now able to more accurately assess causal effects of many genetic and environmental factors. Genome-wide association studies have been able to localize many causal genetic variants predisposing to certain diseases. However, these studies only explain a small portion of variations in the heritability of diseases. More advanced statistical models are urgently needed to identify and characterize some additional genetic and environmental factors and their interactions, which will enable us to better understand the causes of complex diseases. In the past decade, thanks to the increasing computational capabilities and novel statistical developments, Bayesian methods have been widely applied in the genetics/genomics researches and demonstrating superiority over some regular approaches in certain research areas. Gene-environment and gene-gene interaction studies are among the areas where Bayesian methods may fully exert its functionalities and advantages. This dissertation focuses on developing new Bayesian statistical methods for data analysis with complex gene-environment and gene-gene interactions, as well as extending some existing methods for gene-environment interactions to other related areas. It includes three sections: (1) Deriving the Bayesian variable selection framework for the hierarchical gene-environment and gene-gene interactions; (2) Developing the Bayesian Natural and Orthogonal Interaction (NOIA) models for gene-environment interactions; and (3) extending the applications of two Bayesian statistical methods which were developed for gene-environment interaction studies, to other related types of studies such as adaptive borrowing historical data. We propose a Bayesian hierarchical mixture model framework that allows us to investigate the genetic and environmental effects, gene by gene interactions (epistasis) and gene by environment interactions in the same model. It is well known that, in many practical situations, there exists a natural hierarchical structure between the main effects and interactions in the linear model. Here we propose a model that incorporates this hierarchical structure into the Bayesian mixture model, such that the irrelevant interaction effects can be removed more efficiently, resulting in more robust, parsimonious and powerful models. We evaluate both of the 'strong hierarchical' and 'weak hierarchical' models, which specify that both or one of the main effects between interacting factors must be present for the interactions to be included in the model. The extensive simulation results show that the proposed strong and weak hierarchical mixture models control the proportion of false positive discoveries and yield a powerful approach to identify the predisposing main effects and interactions in the studies with complex gene-environment and gene-gene interactions. We also compare these two models with the 'independent' model that does not impose this hierarchical constraint and observe their superior performances in most of the considered situations. The proposed models are implemented in the real data analysis of gene and environment interactions in the cases of lung cancer and cutaneous melanoma case-control studies. The Bayesian statistical models enjoy the properties of being allowed to incorporate useful prior information in the modeling process. Moreover, the Bayesian mixture model outperforms the multivariate logistic model in terms of the performances on the parameter estimation and variable selection in most cases. Our proposed models hold the hierarchical constraints, that further improve the Bayesian mixture model by reducing the proportion of false positive findings among the identified interactions and successfully identifying the reported associations. This is practically appealing for the study of investigating the causal factors from a moderate number of candidate genetic and environmental factors along with a relatively large number of interactions. The natural and orthogonal interaction (NOIA) models of genetic effects have previously been developed to provide an analysis framework, by which the estimates of effects for a quantitative trait are statistically orthogonal regardless of the existence of Hardy-Weinberg Equilibrium (HWE) within loci. Ma et al. (2012) recently developed a NOIA model for the gene-environment interaction studies and have shown the advantages of using the model for detecting the true main effects and interactions, compared with the usual functional model. In this project, we propose a novel Bayesian statistical model that combines the Bayesian hierarchical mixture model with the NOIA statistical model and the usual functional model. The proposed Bayesian NOIA model demonstrates more power at detecting the non-null effects with higher marginal posterior probabilities. Also, we review two Bayesian statistical models (Bayesian empirical shrinkage-type estimator and Bayesian model averaging), which were developed for the gene-environment interaction studies. Inspired by these Bayesian models, we develop two novel statistical methods that are able to handle the related problems such as borrowing data from historical studies. The proposed methods are analogous to the methods for the gene-environment interactions on behalf of the success on balancing the statistical efficiency and bias in a unified model. By extensive simulation studies, we compare the operating characteristics of the proposed models with the existing models including the hierarchical meta-analysis model. The results show that the proposed approaches adaptively borrow the historical data in a data-driven way. These novel models may have a broad range of statistical applications in both of genetic/genomic and clinical studies.

Functional analysis of GNA -1, a Galphai superfamily member found in the filamentous fungus Neurospora crassa

Heterotrimeric GTP-binding proteins, G proteins, are integral components of eukaryotic signaling systems linking extracellular signals to intracellular responses. Through coupling to seven-transmembrane helix receptors, G proteins convey primary signaling events into multi-leveled cascades of intracellular activity by regulating downstream enzymes, collectively called effectors. The effector enzymes regulated by G proteins include adenylyl cyclase, cAMP phosphodiesterase, phospolipase C-β, mitogen-activated protein kinases, and ion channels. ^ Neurospora crassa is a multicellular, filamentous fungus that is capable of both asexual and sexual reproduction by elaboration of specialized, developmentally controlled structures that give rise to either asexual or sexual spores, respectively. N. crassa possesses at least three heterotrimeric Gα proteins (GNA-1–3) and one Gβ subunit (GNB-1). GNA-1 was the first microbial protein that could be classified in the Gαi superfamily based on its amino acid identity and demonstration that it is a substrate for ADP-ribosylation by pertussis toxin. ^ Experiments were designed to identify the signal transduction pathways and the effector enzymes regulated by GNA-1. Targeted gene-replacement of gna-1 revealed that GNA-1 controls multiple developmental pathways including both asexual and sexual reproduction, maintenance of growth, and resistance to osmotic stress. The Gαi and Gαz members of the Gαi superfamily negatively regulate adenylyl cyclase activity in mammalian cells; therefore, adenylyl cyclase and cAMP levels were measured in Δgna-1 strains and also in strains that were deleted for both gna-1 and gna-2, a second Gα in N. crassa shown to have overlapping functions with GNA-1. Direct measurements of adenylyl cyclase activity revealed that GNA-1, but not GNA-2, was responsible for GTP-stimulated adenylyl cyclase activity in N. crassa. Furthermore, anti-GNA-1 IgG could specifically inhibit GTP-stimulated adenylyl cyclase activity in wild-type strain extracts. These studies also provided evidence that N. crassa possesses feedback mechanisms that control steady-state cAMP levels through indirect regulation of cAMP-phosphodiesterase activity; mutations in gna-1 and gna-2 were additive in their effect on lowering cAMP-phosphodiesterase activity under growth conditions where steady-state cAMP levels were normal but GTP-stimulated adenylyl cyclase activity was reduced 90% in comparison to control strains. ^ Genetic and biochemical epistasis experiments utilizing a Δ gna-1 cr-1 mutant suggest that GNA-1 is essential for female fertility in a cAMP-independent pathway. Furthermore, deletion of gna-1 in a cr-1 background exacerbated many of the defects already observed in the cr-1 strain including more severe growth restriction and developmental defects. However, deletion of gna-1 had no effect on the increased thermotolerance of cr-1, which has been attributed to loss of cAMP. cr-1 possesses GNA-1 protein, and crude membrane fractions from this strain reconstituted GTP-stimulated adenylyl cyclase activity in Δgna-1 membrane fractions. These studies provide direct evidence for the involvement of Gα proteins in the regulation of adenylyl cyclase activity in eukaryotic microbes. ^

The SasS -SasR two-component signal transduction system required for early Myxococcus xanthus multicellular *development

Initiation of Myxococcus xanthus multicellular development requires both nutrient limitation and high cell density. The extracellular signal, A signal, which consists of a set of amino acids at specific concentrations, serves as a cell density signal in M. xanthus early development. A reporter gene, designated 4521, that requires both starvation and A signal for developmental expression was used to identify mutations in the signal transduction pathways. A group of point mutations located in the chromosomal sasB locus that bypasses both requirements was previously isolated. One of these point mutations, sasB7, was mapped to the sasS gene, which is predicted to encode a transmembrane histidine protein kinase required for normal development. SasS is a positive regulator of 4521 and a candidate A signal sensor. This dissertation continues the characterization of the sasB locus, focusing on the sasR gene and the functional relationship of SasS and SasR. ^ The sasR gene is located 2.2-kb downstream of sasS. It is predicted to encode an NtrC-like response regulator, which belongs to the family of sigma54 transcriptional activators. SasR is a positive regulator of 4521 gene and is required for normal development. The sasR mutant displays phenotypes similar to that of sasS mutant. Both SasS and SasR are required for the A-signal-dependent 4521 expression. Genetic epistasis analysis indicates that SasR functions downstream of SasS. Biochemical studies show that SasS has autokinase activity, and phosphorylated SasS is able to transfer its phosphate to SasR. We propose that SasS and SasR form a two-component signal transduction system in the A signal transduction pathway. ^ To search for the genes regulated by SasS and SasR, expression patterns of a group of developmental genes were compared in wild-type and sasS null mutant backgrounds. SasS and SasR were found to positively regulate sasN and 4521. The sasN gene was previously identified as a negative regulator of 4521, located at about 170-bp downstream of sasR. It is required for normal fruiting body development. Based on the above data, a regulatory network consisting of sasS, sasR, sasN, and 4521 is hypothesized, and the interactions of the components in this network can now be further studied. ^

New insights into the regulation of NADPH oxidase dependent reactive oxygen species signaling during the plant immune response

Las NADPH oxidasas de plantas, denominadas “respiratory burst oxidase homologues” (RBOHs), producen especies reactivas del oxígeno (ROS) que median un amplio rango de funciones. En la célula vegetal, el ajuste preciso de la producción de ROS aporta la especificidad de señal para generar una respuesta apropiada ante las amenazas ambientales. RbohD y RbohF, dos de los diez genes Rboh de Arabidopsis, son pleiotrópicos y median diversos procesos fisiológicos en respuesta a patógenos. El control espacio-temporal de la expresión de los genes RbohD y RbohF podría ser un aspecto crítico para determinar la multiplicidad de funciones de estas oxidasas. Por ello, generamos líneas transgénicas de Arabidopsis con fusiones de los promoters de RbohD y RbohF a los genes delatores de la B-glucuronidasa y la luciferasa. Estas líneas fueron empleadas para revelar el patrón de expresión diferencial de RbohD y RbohF durante la respuesta inmune de Arabidopsis a la bacteria patógena Pseudomonas syringae pv. tomato DC3000, el hongo necrótrofo Plectosphaerella cucumerina y en respuesta a señales relacionadas con la respuesta inmune. Nuestros experimentos revelan un patrón de expresión diferencial de los promotores de RbohD y RbohF durante el desarrollo de la planta y en la respuesta inmune de Arabidopsis. Además hemos puesto de manifiesto que existe una correlación entre el nivel de actividad de los promotores de RbohD y RbohF con la acumulación de ROS y el nivel de muerte celular en respuesta a patógenos. La expression de RbohD y RbohF también es modulada de manera diferencial en respuesta a patrones moleculares asociados a patógenos (PAMPs) y por ácido abscísico (ABA). Cabe destacar que, mediante una estrategia de intercambio de promotores, hemos revelado que la región promotora de RbohD, es necesaria para dirigir la producción de ROS en respuesta a P. cucumerina. Adicionalmente, la activación del promotor de RbohD en respuesta al aislado de P. cucumerina no adaptado a Arabidopsis 2127, nos llevó a realizar ensayos de susceptibilidad con el doble mutante rbohD rbohF que han revelado un papel desconocido de estas oxidasas en resistencia no-huesped. La interacción entre la señalización dependiente de las RBOHs y otros componentes de la respuesta inmune de plantas podría explicar también las distintas funciones que median estas oxidasas en relación con la respuesta inmune. Entre la gran cantidad de señales coordinadas con la actividad de las RBOHs, existen evidencias genéticas y farmacológicas que indican que las proteínas G heterotriméricas están implicadas en algunas de las rutas de señalización mediadas por ROS derivadas de los RBOHs en respuesta a señales ambientales. Por ello hemos estudiado la relación entre estas RBOH-NADPH oxidasas y AGB1, la subunidad β de las proteínas G heterotriméricas en la respuesta inmune de Arabidopsis. Análisis de epistasis indican que las proteínas G heterotriméricas están implicadas en distintas rutas de señalización en defensa mediadas por las RBOHs. Nuestros resultados ilustran la relación compleja entre la señalización mediada por las RBOHs y las proteínas G heterotriméricas, que varía en función de la interacción planta-patógeno analizada. Además, hemos explorado la posible asociación entre AGB1 con RBOHD y RBOHF en eventos tempranos de la respuesta immune. Cabe señalar que experimentos de coímmunoprecipitación apuntan a una posible asociación entre AGB1 y la kinasa citoplasmática reguladora de RBOHD, BIK1. Esto indica un posible mecanismo de control de la función de esta NADPH oxidase por AGB1. En conjunto, estos datos aportan nuevas perspectivas sobre cómo, a través del control transcripcional o mediante la interacción con las proteínas G heterotriméricas, las NADPH oxidases de plantas median la producción de ROS y la señalización por ROS en la respuesta inmune. Nuestro trabajo ejemplifica cómo la regulación diferencial de dos miembros de una familia multigénica, les permite realizar distintas funciones fisiológicas especializadas usando un mismo mecanismo enzimático. ABSTRACT The plant NADPH oxidases, termed respiratory burst oxidase homologues (RBOHs), produce reactive oxygen species (ROS) which mediate a wide range of functions. Fine tuning this ROS production provides the signaling specificity to the plant cell to produce the appropriate response to environmental threats. RbohD and RbohF, two of the ten Rboh genes present in Arabidopsis, are pleiotropic and mediate diverse physiological processes in response to pathogens. One aspect that may prove critical to determine the multiplicity of functions of RbohD and RbohF is the spatio-temporal control of their gene expression. Thus, we generated Arabidopsis transgenic lines with RbohD- and RbohF-promoter fusions to the β-glucuronidase and the luciferase reporter genes. These transgenics were employed to reveal RbohD and RbohF promoter activity during Arabidopsis immune response to the pathogenic bacterium Pseudomonas syringae pv tomato DC3000, the necrotrophic fungus Plectosphaerella cucumerina and in response to immunity-related cues. Our experiments revealed a differential expression pattern of RbohD and RbohF throughout plant development and during Arabidopsis immune response. Moreover, we observed a correlation between the level of RbohD and RbohF promoter activity, the accumulation of ROS and the amount of cell death in response to pathogens. RbohD and RbohF gene expression was also differentially modulated by pathogen associated molecular patterns and abscisic acid. Interestingly, a promoter-swap strategy revealed the requirement for the promoter region of RbohD to drive the production of ROS in response to P. cucumerina. Additionally, since the RbohD promoter was activated during Arabidopsis interaction with a non-adapted P. cucumerina isolate 2127, we performed susceptibility tests to this fungal isolate that uncovered a new role of these oxidases on non-host resistance. The interplay between RBOH-dependent signaling with other components of the plant immune response might also explain the different immunity-related functions mediated by these oxidases. Among the plethora of signals coordinated with RBOH activity, pharmacological and genetic evidence indicates that heterotrimeric G proteins are involved in some of the signaling pathways mediated by RBOH–derived ROS in response to environmental cues. Therefore, we analysed the interplay between these RBOH-NADPH oxidases and AGB1, the Arabidopsis β-subunit of heterotrimeric G proteins during Arabidopsis immune response. We carried out epistasis studies that allowed us to test the implication of AGB1 in different RBOH-mediated defense signaling pathways. Our results illustrate the complex relationship between RBOH and heterotrimeric G proteins signaling, that varies depending on the type of plant-pathogen interaction. Furthermore, we tested the potential association between AGB1 with RBOHD and RBOHF during early immunity. Interestingly, our co-immunoprecipitation experiments point towards an association of AGB1 and the RBOHD regulatory kinase BIK1, thus providing a putative mechanism in the control of the NADPH oxidase function by AGB1. Taken all together, these studies provide further insights into the role that transcriptional control or the interaction with heterotrimeric G-proteins have on RBOH-NADPH oxidase-dependent ROS production and signaling in immunity. Our work exemplifies how, through a differential regulation, two members of a multigenic family achieve specialized physiological functions using a common enzymatic mechanism.

Identification and characterization of Saccharomyces cerevisiae EXO1, a gene encoding an exonuclease that interacts with MSH2

A two-hybrid screen was used to identify Saccharomyces cerevisiae genes encoding proteins that interact with MSH2. One gene was found to encode a homologue of Schizosaccharomyces pombe EXO1, a double-stranded DNA-specific 5′–3′ exonuclease. S. cerevisiae EXO1 interacted with both S. cerevisiae and human MSH2 in two-hybrid and coimmunoprecipitation experiments. exo1 mutants showed a mutator phenotype, and epistasis analysis was consistent with EXO1 functioning in the MSH2-dependent mismatch repair pathway. exo1 mutations were lethal in combination with rad27 mutations, and overexpression of EXO1 suppressed both the temperature sensitive and mutator phenotypes of rad27 mutants.

Mitochondrial Inheritance Is Delayed in Saccharomyces cerevisiae Cells Lacking the Serine/Threonine Phosphatase PTC1

In wild-type yeast mitochondrial inheritance occurs early in the cell cycle concomitant with bud emergence. Cells lacking the PTC1 gene initially produce buds without a mitochondrial compartment; however, these buds later receive part of the mitochondrial network from the mother cell. Thus, the loss of PTC1 causes a delay, but not a complete block, in mitochondrial transport. PTC1 encodes a serine/threonine phosphatase in the high-osmolarity glycerol response (HOG) pathway. The mitochondrial inheritance delay in the ptc1 mutant is not attributable to changes in intracellular glycerol concentrations or defects in the organization of the actin cytoskeleton. Moreover, epistasis experiments with ptc1Δ and mutations in HOG pathway kinases reveal that PTC1 is not acting through the HOG pathway to control the timing of mitochondrial inheritance. Instead, PTC1 may be acting either directly or through a different signaling pathway to affect the mitochondrial transport machinery in the cell. These studies indicate that the timing of mitochondrial transport in wild-type cells is genetically controlled and provide new evidence that mitochondrial inheritance does not depend on a physical link between the mitochondrial network and the incipient bud site.

Muscle LIM Proteins Are Associated with Muscle Sarcomeres and Require dMEF2 for Their Expression during Drosophila Myogenesis

A genetic hierarchy of interactions, involving myogenic regulatory factors of the MyoD and myocyte enhancer-binding 2 (MEF2) families, serves to elaborate and maintain the differentiated muscle phenotype through transcriptional regulation of muscle-specific target genes. Much work suggests that members of the cysteine-rich protein (CRP) family of LIM domain proteins also play a role in muscle differentiation; however, the specific functions of CRPs in this process remain undefined. Previously, we characterized two members of the Drosophila CRP family, the muscle LIM proteins Mlp60A and Mlp84B, which show restricted expression in differentiating muscle lineages. To extend our analysis of Drosophila Mlps, we characterized the expression of Mlps in mutant backgrounds that disrupt specific aspects of muscle development. We show a genetic requirement for the transcription factor dMEF2 in regulating Mlp expression and an ability of dMEF2 to bind, in vitro, to consensus MEF2 sites derived from those present in Mlp genomic sequences. These data suggest that the Mlp genes may be direct targets of dMEF2 within the genetic hierarchy controlling muscle differentiation. Mutations that disrupt myoblast fusion fail to affect Mlp expression. In later stages of myogenic differentiation, which are dedicated primarily to assembly of the contractile apparatus, we analyzed the subcellular distribution of Mlp84B in detail. Immunofluorescent studies revealed the localization of Mlp84B to muscle attachment sites and the periphery of Z-bands of striated muscle. Analysis of mutations that affect expression of integrins and α-actinin, key components of these structures, also failed to perturb Mlp84B distribution. In conclusion, we have used molecular epistasis analysis to position Mlp function downstream of events involving mesoderm specification and patterning and concomitant with terminal muscle differentiation. Furthermore, our results are consistent with a structural role for Mlps as components of muscle cytoarchitecture.

Functions of the DNA damage response pathway target Ho endonuclease of yeast for degradation via the ubiquitin-26S proteasome system

Ho endonuclease of Saccharomyces cerevisiae is a homing endonuclease that makes a site-specific double-strand break in the MAT gene in late G1. Here we show that Ho is rapidly degraded via the ubiquitin-26S proteasome system through two ubiquitin-conjugating enzymes UBC2Rad6 and UBC3Cdc34. UBC2Rad6 is complexed with the ring finger DNA-binding protein Rad18, and we find that Ho is stabilized in rad18 mutants. We show that the Ho degradation pathway involving UBC3Cdc34 goes through the Skp1/Cdc53/F-box (SCF) ubiquitin ligase complex and identify a F-box protein, Yml088w, that is required for Ho degradation. Components of a defined pathway of the DNA damage response, MEC1, RAD9, and CHK1, are also necessary for Ho degradation, whereas functions of the RAD24 epistasis group and the downstream effector RAD53 have no role in degradation of Ho. Our results indicate a link between the endonuclease function of Ho and its destruction.

Genetic evidence for the involvement of DNA ligase IV in the DNA-PK-dependent pathway of non-homologous end joining in mammalian cells

Cells of vertebrates remove DNA double-strand breaks (DSBs) from their genome predominantly utilizing a fast, DNA-PKcs-dependent form of non-homologous end joining (D-NHEJ). Mutants with inactive DNA-PKcs remove the majority of DNA DSBs utilizing a slow, DNA-PKcs-independent pathway that does not utilize genes of the RAD52 epistasis group, is error-prone and can therefore be classified as a form of NHEJ (termed basic or B-NHEJ). We studied the role of DNA ligase IV in these pathways of NHEJ. Although biochemical studies show physical and functional interactions between the DNA-PKcs/Ku and the DNA ligase IV/Xrcc4 complexes suggesting operation within the same pathway, genetic evidence to support this notion is lacking in mammalian cells. Primary human fibroblasts (180BR) with an inactivating mutation in DNA ligase IV, rejoined DNA DSBs predominantly with slow kinetics similar to those observed in cells deficient in DNA-PKcs, or in wild-type cells treated with wortmannin to inactivate DNA-PK. Treatment of 180BR cells with wortmannin had only a small effect on DNA DSB rejoining and no effect on cell radiosensitivity to killing although it sensitized control cells to 180BR levels. This is consistent with DNA ligase IV functioning as a component of the D-NHEJ, and demonstrates the unperturbed operation of the DNA-PKcs-independent pathway (B-NHEJ) at significantly reduced levels of DNA ligase IV. In vitro, extracts of 180BR cells supported end joining of restriction endonuclease-digested plasmid to the same degree as extracts of control cells when tested at 10 mM Mg2+. At 0.5 mM Mg2+, where only DNA ligase IV is expected to retain activity, low levels of end joining (∼10% of 10 mM) were seen in the control but there was no detectable activity in 180BR cells. Antibodies raised against DNA ligase IV did not measurably inhibit end joining at 10 mM Mg2+ in either cell line. Thus, in contrast to the situation in vivo, end joining in vitro is dominated by pathways with properties similar to B-NHEJ that do not display a strong dependence on DNA ligase IV, with D-NHEJ retaining only a limited contribution. The implications of these observations to studies of NHEJ in vivo and in vitro are discussed.

Homologous DNA recombination in vertebrate cells

The RAD52 epistasis group genes are involved in homologous DNA recombination, and their primary structures are conserved from yeast to humans. Although biochemical studies have suggested that the fundamental mechanism of homologous DNA recombination is conserved from yeast to mammals, recent studies of vertebrate cells deficient in genes of the RAD52 epistasis group reveal that the role of each protein is not necessarily the same as that of the corresponding yeast gene product. This review addresses the roles and mechanisms of homologous recombination-mediated repair with a special emphasis on differences between yeast and vertebrate cells.

Gestational drive and the green-bearded placenta.

A "green beard" refers to a gene, or group of genes, that is able to recognize itself in other individuals and direct benefits to these individuals. Green-beard effects have been dismissed as implausible by authors who have implicitly assumed sophisticated mechanisms of perception and complex behavioral responses. However, many simple mechanisms for genes to "recognize" themselves exist at the maternal-fetal interface of viviparous organisms. Homophilic cell adhesion molecules, for example, are able to interact with copies of themselves on other cells. Thus, the necessary components of a green-beard effect -- feature, recognition, and response -- can be different aspects of the phenotype of a single gene. Other green-beard effects could involve coalitions of genes at closely linked loci. In fact, any form of epistasis between a locus expressed in a mother and a closely linked locus expressed in the fetus has the property of "self-recognition." Green-beard effects have many formal similarities to systems of meiotic drive and, like them, can be a source of intragenomic conflict.

Complex formation in yeast double-strand break repair: participation of Rad51, Rad52, Rad55, and Rad57 proteins.

The repair of DNA double-strand breaks in Saccharomyces cerevisiae requires genes of the RAD52 epistasis group, of which RAD55 and RAD57 are members. Here, we show that the x-ray sensitivity of rad55 and rad57 mutant strains is suppressible by overexpression of RAD51 or RAD52. Virtually complete suppression is provided by the simultaneous overexpression of RAD51 and RAD52. This suppression occurs at 23 degrees C, where these mutants are more sensitive to x-rays, as well as at 30 degrees C and 36 degrees C. In addition, a recombination defect of rad55 and rad57 mutants is similarly suppressed. Direct in vivo interactions between the Rad51 and Rad55 proteins, and between Rad55 and Rad57, have also been identified by using the two-hybrid system. These results indicate that these four proteins constitute part of a complex, a "recombinosome," to effect the recombinational repair of double-strand breaks.

Modularisation épistatique des loci à trait quantitatif associés à la pression artérielle et identification de gènes candidats pour l’hypertension

Problématique: L’hypertension artérielle essentielle, facteur de risque majeur dans le développement des maladies cardiovasculaires, est un trait multigénique complexe dont les connaissances sur le déterminisme génétique nécessitent d’être approfondies. De nombreux loci à trait quantitatif (QTLs); soit des gènes responsables de faire varier la pression artérielle (PA), ont été identifiés chez l’humain et le modèle animal. Cependant, le mystère plane encore sur la façon dont ces gènes fonctionnent ensemble pour réguler la PA. Hypothèse et objectif: Plutôt qu’une addition de QTLs ayant chacun une action infinitésimale sur la PA, une interaction épistatique entre les gènes serait responsable du phénotype hypertendu. Ainsi, l’étude de cette épistasie entre les gènes impliqués, directement ou indirectement, dans l’homéostasie de la PA nous permettrait d’explorer de nouvelles voies de régulation moléculaire en cause dans cette maladie. Méthodes: Via la réalisation de souches congéniques de rats, où un segment chromosomique provenant d’une souche receveuse hypertendue (Dahl Salt Sensitive, SS/Jr) est remplacé par son homologue provenant d’une souche donneuse normotendue (Lewis, LEW), des QTLs peuvent être mis en évidence. Dans ce contexte, la combinaison de QTLs via la création de doubles ou multiples congéniques constitue la première démonstration fonctionnelle des interactions intergéniques. Résultats: Vingt-sept combinaisons au total nous ont menés à l’appréciation d’une modularisation des QTLs. Ces derniers ont été catégorisés selon deux principaux modules épistatiques (EMs) où les QTLs appartenant à un même EM sont épistatiques entre eux et participent à une même voie régulatrice. Les EMs/cascades agissent alors en parallèle pour réguler la PA. Grâce à l’existence de QTLs ayant des effets opposés sur la PA, nous avons pu établir l’ordre hiérarchique entre trois paires de QTLs. Cependant, lorsque cette suite régulatrice ne peut être déterminée, d’autres approches sont nécessaires. Nos travaux nous ont mené à l’identification d’un QTL situé sur le chromosome 16 du rat (C16QTL), appartenant au EM1 et qui révélerait une nouvelle voie de l’homéostasie de la PA. Le gène retinoblastoma-associated protein 140 (Rap140)/family with sequence similarity 208 member A (Fam208a), présentant une mutation non synonyme entre SS/Jr et LEW est le gène candidat le plus plausible pour représenter C16QTL. Celui-ci code pour un facteur de transcription et semblerait influencer l’expression de Solute carrier family 7 (cationic amino acid transporter, y+ system) member 12 (Slc7a12), spécifiquement et significativement sous exprimé dans les reins de la souche congénique portant C16QTL par rapport à la souche SS/Jr. Rap140/Fam208a agirait comme un inhibiteur de la transcription de Slc7a12 menant à une diminution de la pression chez Lewis. Conclusions: L’architecture complexe de la régulation de la PA se dévoile mettant en scène de nouveaux acteurs, pour la plupart inconnus pour leur implication dans la PA. L’étude de la nouvelle voie de signalisation Rap140/Fam208a - Slc7a12 nous permettra d’approfondir nos connaissances quant à l’homéostasie de la pression artérielle et de l’hypertension chez SS/Jr. À long terme, de nouveaux traitements anti-hypertenseurs, ciblant plus d’une voie de régulation à la fois, pourraient voir le jour.

High throughput determination of sequence-function relationships in protein and RNA

Thesis (Ph.D.)--University of Washington, 2016-06