Genetic mechanisms of Na+, K+-ATPase regulation in hypertension

The spontaneously hypertensive rat (SHR) is a model of essential hypertension. During the early development of hypertension, the SHR demonstrates increased proximal tubule (PT) Na+ reabsorption. I hypothesized that the increased PT Na+ reabsorption exhibited by the young SHR was due to altered sub-cellular distribution of Na+, K +-ATPase compared to the normotensive Wistar Kyoto (WKY). The hypothesis is supported, herein, by observations of greater Na+, K +-ATPase α 1 abundance in PT plasma membrane and lower abundance in late endosomes of 4wk SHR despite no difference in total PT α 1 abundance. There is a greater amount of Ser-18 unphosphorylated α 1 in the 4wk SHR PT. Total PT Na+, K+-ATPase γ abundance is greater in SHR at 4wk and 16wk but γ abundance in plasma membrane is greater only at 4wk. The phosphatase, calcineurin, was chosen for study because it is involved in the stimulation of Na+, K +-ATPase. No difference in calcineurin coding sequence, expression, or activity was observed in SHR. Gene expression arrays were next used to find candidate genes involved in the regulation of Na+, K +-ATPase. The first candidate analyzed was soluble epoxide hydrolase (sEH). The gene encoding sEH (EPHX2) showed lower expression in SHR. There was also a reduction in sEH protein abundance but there was no correlation between protein abundance and blood pressure in F2 progeny. Two EPHX2 alleles were identified, an ancestral allele and a variant allele containing four polymorphisms. sEH activity was greater in animals carrying the variant allele but the inheritance of the variant allele did not correlate with blood pressure. Gene expression arrays also led to the examination of genes involved in redox balance/Na+, K+-ATPase regulation. A pattern of lower expression of genes involved in reactive radical detoxification in SHR was discerned. Six transcription factor binding sites were identified that occurred more often in these genes. Three transcription factors that bind to the HNF1 site were expressed at lower levels in SHR. This points to the HNF1 transcriptional complex as an important trans-acting regulator of a wide range of genes involved in altered redox balance in SHR. ^

Targeted Deletion Of Functionally Validated Enhancers Defines Their Role in Mouse Limb Development

Transcriptional enhancers are genomic DNA sequences that contain clustered transcription factor (TF) binding sites. When combinations of TFs bind to enhancer sequences they act together with basal transcriptional machinery to regulate the timing, location and quantity of gene transcription. Elucidating the genetic mechanisms responsible for differential gene expression, including the role of enhancers, during embryological and postnatal development is essential to an understanding of evolutionary processes and disease etiology. Numerous methods are in use to identify and characterize enhancers. Several high-throughput methods generate large datasets of enhancer sequences with putative roles in embryonic development. However, few enhancers have been deleted from the genome to determine their roles in the development of specific structures, such as the limb. Manipulation of enhancers at their endogenous loci, such as the deletion of such elements, leads to a better understanding of the regulatory interactions, rules and complexities that contribute to faithful and variant gene transcription – the molecular genetic substrate of evolution and disease. To understand the endogenous roles of two distinct enhancers known to be active in the mouse embryo limb bud we deleted them from the mouse genome. I hypothesized that deletion of these enhancers would lead to aberrant limb development. The enhancers were selected because of their association with p300, a protein associated with active transcription, and because the human enhancer sequences drive distinct lacZ expression patterns in limb buds of embryonic day (E) 11.5 transgenic mice. To confirm that the orthologous mouse enhancers, mouse 280 and 1442 (M280 and M1442, respectively), regulate expression in the developing limb we generated stable transgenic lines, and examined lacZ expression. In M280-lacZ mice, expression was detected in E11.5 fore- and hindlimbs in a region that corresponds to digits II-IV. M1442-lacZ mice exhibited lacZ expression in posterior and anterior margins of the fore- and hindlimbs that overlapped with digits I and V and several wrist bones. We generated mice lacking the M280 and M1442 enhancers by gene targeting. Intercrosses between M280 -/+ and M1442 -/+, respectively, generated M280 and M1442 null mice, which are born at expected Mendelian ratios and manifest no gross limb malformations. Quantitative real-time PCR of mutant E11.5 limb buds indicated that significant changes in transcriptional output of enhancer-proximal genes accompanied the deletion of both M280 and M1442. In neonatal null mice we observed that all limb bones are present in their expected positions, an observation also confirmed by histology of E18.5 distal limbs. Fine-scale measurement of E18.5 digit bone lengths found no differences between mutant and control embryos. Furthermore, when the developmental progression of cartilaginous elements was analyzed in M280 and M1442 embryos from E13.5-E15.5, transient development defects were not detected. These results demonstrate that M280 and M1442 are not required for mouse limb development. Though M280 is not required for embryonic limb development it is required for the development and/or maintenance of body size – adult M280 mice are significantly smaller than control littermates. These studies highlight the importance of experiments that manipulate enhancers in situ to understand their contribution to development.

XENOESTROGEN-SPECIFIC MECHANISMS OF DEVELOPMENTAL REPROGRAMMING CORRELATE WITH GENE EXPRESSION AND TUMOR DEVELOPMENT

Environmental exposures during sensitive windows of development can reprogram normal physiological responses and alter disease susceptibility later in life in a process known as developmental reprogramming. We have shown that neonatal exposure to the xenoestrogen diethylstilbestrol (DES) can developmentally reprogram the reproductive tract in genetically susceptible Eker rats giving rise to complete penetrance of uterine leiomyoma. Based on this, we hypothesized that xenoestrogens, including genistein (GEN) and bisphenol A (BPA), reprogram estrogen-responsive gene expression in the myometrium and promote the development of uterine leiomyoma. We proposed the mechanism that is responsible for the developmental reprogramming of gene expression was through estrogen (E2)/ xenoestrogen inducedrapid ER signaling, which modifies the histone methyltransferase Enhancer of Zeste homolog 2 (EZH2) via activation of the PI3K/AKT pathway. We further hypothesized that there is a xenostrogen-specific effect on this pathway altering patterns of histone modification, DNA methylation and gene expression. In addition to our novel finding that E2/DES-induced phosphorylation of EZH2 by AKT reduces the levels of H3K27me3 in vitro and in vivo, this work demonstrates in vivo that a brief neonatal exposure to GEN, in contrast to BPA, activates the PI3K/AKT pathway to regulate EZH2 and decreases H3K27me3 levels in the neonatal uterus. Given that H3K27me3 is a repressive mark that has been shown to result in DNA methylation and gene silencing we investigated the methylation of developmentally reprogrammed genes. In support of this evidence, we show that neonatal DES exposure in comparison to VEH, leads to hypomethylation of the promoter of a developmentally reprogrammed gene, Gria2, that become hyper-responsive to estrogen in the adult myometrium indicating vi that DES exposure alter gene expression via chromatin remodeling and loss of DNA methylation. In the adult uterus, GEN and BPA exposure developmentally reprogrammed expression of estrogen-responsive genes in a manner opposite of one another, correlating with our previous data. Furthermore, the ability of GEN and BPA to developmental reprogram gene expression correlated with tumor incidence and multiplicity. These data show that xenoestrogens have unique effects on the activation of non-genomic signaling in the developing uterus that promotes epigenetic and genetic alterations, which are predictive of developmental reprogramming and correlate with their ability to modulate hormone-dependent tumor development.

DOUBLE-STRAND BREAK REPAIR PATHWAYS IN DNA STRUCTURE-INDUCED GENETIC INSTABILITY

Genetic instability in mammalian cells can occur by many different mechanisms. In the absence of exogenous sources of DNA damage, the DNA structure itself has been implicated in genetic instability. When the canonical B-DNA helix is naturally altered to form a non-canonical DNA structure such as a Z-DNA or H-DNA, this can lead to genetic instability in the form of DNA double-strand breaks (DSBs) (1, 2). Our laboratory found that the stability of these non-B DNA structures was different in mammals versus Escherichia coli (E.coli) bacteria (1, 2). One explanation for the difference between these species may be a result of how DSBs are repaired within each species. Non-homologous end-joining (NHEJ) is primed to repair DSBs in mammalian cells, while bacteria that lack NHEJ (such as E.coli), utilize homologous recombination (HR) to repair DSBs. To investigate the role of the error-prone NHEJ repair pathway in DNA structure-induced genetic instability, E.coli cells were modified to express genes to allow for a functional NHEJ system under different HR backgrounds. The Mycobacterium tuberculosis NHEJ sufficient system is composed of Ku and Ligase D (LigD) (3). These inducible NHEJ components were expressed individually and together in E.coli cells, with or without functional HR (RecA/RecB), and the Z-DNA and H-DNA-induced mutations were characterized. The Z-DNA structure gave rise to higher mutation frequencies compared to the controls, regardless of the DSB repair pathway(s) available; however, the type of mutants produced after repair was greatly dictated on the available DSB repair system, indicated by the shift from 2% large-scale deletions in the total mutant population to 24% large-scale deletions when NHEJ was present (4). This suggests that NHEJ has a role in the large deletions induced by Z-DNA-forming sequences. H-DNA structure, however, did not exhibit an increase in mutagenesis in the newly engineered E.coli environment, suggesting the involvement of other factors in regulating H-DNA formation/stability in bacterial cells. Accurate repair by established DNA DSB repair pathways is essential to maintain the stability of eukaryotic and prokaryotic genomes and our results suggest that an error-prone NHEJ pathway was involved in non-B DNA structure-induced mutagenesis in both prokaryotes and eukaryotes.

A Cbfa1-dependent genetic pathway controls bone formation beyond embryonic development.

The molecular mechanisms controlling bone extracellular matrix (ECM) deposition by differentiated osteoblasts in postnatal life, called hereafter bone formation, are unknown. This contrasts with the growing knowledge about the genetic control of osteoblast differentiation during embryonic development. Cbfa1, a transcriptional activator of osteoblast differentiation during embryonic development, is also expressed in differentiated osteoblasts postnatally. The perinatal lethality occurring in Cbfa1-deficient mice has prevented so far the study of its function after birth. To determine if Cbfa1 plays a role during bone formation we generated transgenic mice overexpressing Cbfa1 DNA-binding domain (DeltaCbfa1) in differentiated osteoblasts only postnatally. DeltaCbfa1 has a higher affinity for DNA than Cbfa1 itself, has no transcriptional activity on its own, and can act in a dominant-negative manner in DNA cotransfection assays. DeltaCbfa1-expressing mice have a normal skeleton at birth but develop an osteopenic phenotype thereafter. Dynamic histomorphometric studies show that this phenotype is caused by a major decrease in the bone formation rate in the face of a normal number of osteoblasts thus indicating that once osteoblasts are differentiated Cbfa1 regulates their function. Molecular analyses reveal that the expression of the genes expressed in osteoblasts and encoding bone ECM proteins is nearly abolished in transgenic mice, and ex vivo assays demonstrated that DeltaCbfa1-expressing osteoblasts were less active than wild-type osteoblasts. We also show that Cbfa1 regulates positively the activity of its own promoter, which has the highest affinity Cbfa1-binding sites characterized. This study demonstrates that beyond its differentiation function Cbfa1 is the first transcriptional activator of bone formation identified to date and illustrates that developmentally important genes control physiological processes postnatally.

ATM Signaling to TSC2: Mechanisms and Implications for Cancer Therapy

Ataxia telangiectasia mutated (ATM) is a critical component of the cellular response to DNA damage, where it acts as a damage sensor, and signals to a large network of proteins which execute the important tasks involved in responding to the damage, namely inducing cell cycle checkpoints, inducing DNA repair, modulating transcriptional responses, and regulating cell death pathways if the damage cannot be repaired faithfully. We have now discovered that an additional novel component of this ATM-dependent damage response involves induction of autophagy in response to oxidative stress. In contrast to DNA damage-induced ATM activation however, oxidative stress induced ATM, occurs in the cytoplasm, and does not require nuclear-to-cytoplasmic shuttling of ATM. Using several cell culture systems including MCF7 breast carcinoma cells, SKOV3 ovarian cancer cells, and various lineages of mouse embryonic fibroblasts, we showed that once activated by reactive oxygen species (ROS), ATM signals to mTORC1 to induce autophagy via the LKB1-AMPK-TSC2 pathway. Targeting dysregulation of mTORC1 in Atm-deficient mice, which succumb to lymphomagenesis within 3-4 months of age with daily administration of rapamycin, could significantly extend survival and cause regression of tumors, suggesting that pharmacologically targeting this pathway has therapeutic implications in cancer. We also identified a second contrasting pathway for DNA damage-induced mTORC1 repression which does not require AMPK activation, but does require ATM and TSC2. Several potential mechanisms including mTOR localization and p53-mediated pathways were ruled out however we identified that TSC2 may be an additional cytoplasmic direct ATM substrate that is engaged in response to DNA damage specifically. Lastly, a study was performed to examine whether autophagy induced by ovarian cancer therapeutics (focusing on cisplatin, since paclitaxel does not induce autophagy in the SKOV3 cell line model we used) plays a role in resistance to therapy since autophagy can play both pro-survival mechanisms or be a mechanism of cell death. Using a genetic approach to knock-down Atg5 expression with shRNA in SKOV3 ovarian carcinoma cells, we compared the cytotoxicity of cisplatin in vector or Atg5 knock-down cells, and demonstrated that autophagy does not play any significant role in the response to cisplatin in this cell line.

The physical map of {\it Escherichia coli\/} K-12 strain MG1655: Construction, comparison with other strains and implications for mechanisms controlling oligonucleotide frequency

Complete NotI, SfiI, XbaI and BlnI cleavage maps of Escherichia coli K-12 strain MG1655 were constructed. Techniques used included: CHEF pulsed field gel electrophoresis; transposon mutagenesis; fragment hybridization to the ordered $\lambda$ library of Kohara et al.; fragment and cosmid hybridization to Southern blots; correlation of fragments and cleavage sites with EcoMap, a sequence-modified version of the genomic restriction map of Kohara et al.; and correlation of cleavage sites with DNA sequence databases. In all, 105 restriction sites were mapped and correlated with the EcoMap coordinate system.^ NotI, SfiI, XbaI and BlnI restriction patterns of five commonly used E. coli K-12 strains were compared to those of MG1655. The variability between strains, some of which are separated by numerous steps of mutagenic treatment, is readily detectable by pulsed-field gel electrophoresis. A model is presented to account for the difference between the strains on the basis of simple insertions, deletions, and in one case an inversion. Insertions and deletions ranged in size from 1 kb to 86 kb. Several of the larger features have previously been characterized and some of the smaller rearrangements can potentially account for previously reported genetic features of these strains.^ Some aspects of the frequency and distribution of NotI, SfiI, XbaI and BlnI cleavage sites were analyzed using a method based on Markov chain theory. Overlaps of Dam and Dcm methylase sites with XbaI and SfiI cleavage sites were examined. The one XbaI-Dam overlap in the database is in accord with the expected frequency of this overlap. The occurrence of certain types of SfiI-Dcm overlaps are overrepresented. Of the four subtypes of SfiI-Dcm overlap, only one has a partial inhibitory effect on the activity of SfiI. Recognition sites for all four enzymes are rarer than expected based on oligonucleotide frequency data, with this effect being much stronger for XbaI and BlnI than for NotI and SfiI. The latter two enzyme sites are rare mainly due to apparent negative selection against GGCC (both) and CGGCCG (NotI). The former two enzyme sites are rare mainly due to effects of the VSP repair system on certain di-tri- and tetranucleotides, most notably CTAG. Models are proposed to explain several of the anomalies of oligonucleotide distribution in E. coli, and the biological significance of the systems that produce these anomalies is discussed. ^

Accessing novel developmental mechanisms in the mouse by gene trapping

Genetic analysis is a powerful method for analyzing the function of specific genes in development. I sought to identify novel genes in the mouse using a genetic analysis relying on the expression pattern and phenotype of mutated genes. To this end, I have conducted a gene trap screen using the vector $\rm SA\beta geo,$ a promoterless DNA construct that encodes a fusion protein with lacZ and neomycin resistance activities. Productive integration and expression of the $\beta$geo protein in embryonic stem (ES) cells requires integration into an active transcription unit. The endogenous regulatory elements direct reporter gene expression which reflects the expression of the endogenous gene. Of eight mouse lines generated from gene trap ES cell clones, four showed differential regulation of $\beta$geo activity during embryogenesis. These four were analyzed in more detail.^ Three of the lines RNA 1, RNA2 and RNA 3 had similar expression patterns, within subsets of cells in sites of embryonic hematopoiesis. Cloning of the trapped genes revealed that all three integrations had occurred within 45S rRNA precursor transcription units. These results imply that there exists in these cells some mechanism responsible for the efficient production of the $\beta$geo protein from an RNA polymerase I transcript that is not present in most of the cells in the embryo.^ The fourth line, GT-2, showed widespread, dynamic expression. Many of the sites of expression were important classic embryonic induction systems. Cloning of the sequences fused to the $5\sp\prime$ end of the $\beta$geo sequence revealed that the trapped gene contained significant sequence homology with a previously identified human sequence HumORF5. An open reading frame of this sequence is homologous to a group of eukaryotic proteins that are members of the RNA helicase superfamily I.^ Analysis of the gene trap lines suggests that potentially novel developmental mechanisms have been uncovered. In the case of RNA 1, 2 and 3, the differential production of ribosomal RNAs may be required for differentiation or function of the $\beta$geo positive hematopoietic cells. In the GT-2 line, a previously unsuspected temporal and spatial regulation of a putative RNA helicase implies a role for this activity during specific aspects of mouse development. ^

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. ^

DNA incisions stimulate genetic instability of triplet repeat sequences related to human hereditary neurological diseases

The molecular mechanisms responsible for the expansion and deletion of trinucleotide repeat sequences (TRS) are the focus of our studies. Several hereditary neurological diseases including Huntington's disease, myotonic dystrophy, and fragile X syndrome are associated with the instability of TRS. Using the well defined and controllable model system of Escherichia coli, the influences of three types of DNA incisions on genetic instability of CTG•CAG repeats were studied: DNA double-strand breaks (DSB), single-strand nicks, and single-strand gaps. The DNA incisions were generated in pUC19 derivatives by in vitro cleavage with restriction endonucleases. The cleaved DNA was then transformed into E. coli parental and mutant strains. Double-strand breaks induced deletions throughout the TRS region in an orientation dependent manner relative to the origin of replication. The extent of instability was enhanced by the repeat length and sequence (CTG•CAG vs. CGG•CCG). Mutations in recA and recBC increased deletions, mutations in recF stabilized the TRS, whereas mutations in ruvA had no effect. DSB were repaired by intramolecular recombination, versus an intermolecular gene conversion or crossover mechanism. 30 nt gaps formed a distinct 30 nt deletion product, whereas single strand nicks and gaps of 15 nts did not induce expansions or deletions. Formation of this deletion product required the CTG•CAG repeats to be present in the single-stranded region and was stimulated by E. coli DNA ligase, but was not dependent upon the RecFOR pathway. Models are presented to explain the DSB induced instabilities and formation of the 30 nucleotide deletion product. In addition to the in vitro creation of DSBs, several attempts to generate this incision in vivo with the use of EcoR I restriction modification systems were conducted. ^

BIOLOGICAL MECHANISMS AND CLINICAL IMPLICATIONS OF BCR-ABL-INDUCED MITOCHONDRIAL OXIDATIVE STRESS AND CELL SURVIVAL IN CHRONIC MYELOID LEUKEMIA

Chronic myeloid leukemia (CML), a myeloproliferative disorder, represents approximately 15-20% of all adult leukemia. The development of CML is clearly linked to the constitutively active protein-tyrosine kinase BCR-ABL, which is encoded by BCR-ABL fusion gene as the result of chromosome 9/22 translocation (Philadelphia chromosome). Previous studies have demonstrated that oxidative stress-associated genetic, metabolic and biological alterations contribute to CML cell survival and drug refractory. Mitochondria and NAD(P)H oxidase (NOX) are the major sources of BCR-ABL-induced cellular reactive oxygen species (ROS) production. However, it is still unknown how CML cells maintain the altered redox status, while escaping from the persistent oxidative stress-induced cell death. Therefore, elucidation of the mechanisms by which CML cells cope with oxidative stress will provide new insights into CML leukemogenesis. The major goal of this study is to identify the survival factors protecting CML cells against oxidative stress and develop novel therapeutic strategies to overcome drug resistance. Several experimental models were used to test CML cell redox status and cellular sensitivity to oxidative stress, including BCR-ABL inducible cell lines, BCR-ABL stably transformed cell lines and BCR-ABL-expressing CML blast crisis cells with differential BCL-XL/BCL-2 expressions. Additionally, an artificial CML cell model with heterogenic BCL-XL/BCL-2 expression was established to assess the correlation between differential survival factor expression patterns and cell sensitivity to Imatinib and oxidative stress. In this study, BCL-XL and GSH have been identified as the major survival factors responsive to BCR-ABL-promoted cellular oxidative stress and play a dominant role in regulating the threshold of oxidative stress-induced apoptosis. Cell survival factors BCL-XL and BCL-2 differentially protect mitochondria under oxidative stress. BCL-XL is an essential survival factor in preventing excessive ROS-induced cell death while BCL-2 seems to play a relatively minor role. Furthermore, the redox modulating reagent β-phenethyl isothiocyanate (PEITC) has been found to efficiently deplete GSH and induce potent cell killing effects in drug-resistant CML cells. Combination of PEITC with BCL-XL/BCL2 inhibitor ABT737 or suppression of BCL-XL by BCR-ABL inhibitor Gleevec dramatically sensitizes CML cells to apoptosis. These results have suggested that elevation of BCL-XL and cellular GSH are important for the development of CML, and that redox-directed therapy is worthy of further clinical investigations in CML.

GENETIC ANALYSIS OF THE HIPPO PATHWAY IN MOUSE LIVER

Cancer therapy and tumor treatment remain unsolved puzzles. Genetic screening for tumor suppressor genes in Drosophila revealed the Hippo-signaling pathway as a kinase cascade consisting of five core components. Disrupting the pathway by deleting the main component genes breaks the balance of cell proliferation and apoptosis and results in epithelial tissue tumorigenesis. The pathway is therefore believed to be a tumor suppressor pathway. However, a corresponding role in mammals is yet to be determined. Our lab began to investigate the tumor suppression function of the potent mammalian Hippo pathway by putting floxed alleles into the mouse genome flanking the functional-domain-expressing exons in each component (Mst1, Mst2, Sav1, Lats1 and Lats2). These mice were then crossed with different cre-mouse lines to generate conditional knockout mice. Results indicate a ubiquitous tumor suppression function of these components, predominantly in the liver. A further liver specific analysis of the deletion mutation of these components, as well as the Yap/Taz double deletion mutation, reveals essential roles of the Hippo pathway in regulating hepatic quiescence and embryonic liver development. One of the key cellular mechanisms for the Hippo pathway’s involvement in these liver biological events is likely its cell cycle regulation function. Our work will help to develop potential therapeutic approaches for liver cancer.

Exploration of genetic susceptibility to pancreatic cancer: A GWAS approach

Pancreatic cancer is the 4th most common cause for cancer death in the United States, accompanied by less than 5% five-year survival rate based on current treatments, particularly because it is usually detected at a late stage. Identifying a high-risk population to launch an effective preventive strategy and intervention to control this highly lethal disease is desperately needed. The genetic etiology of pancreatic cancer has not been well profiled. We hypothesized that unidentified genetic variants by previous genome-wide association study (GWAS) for pancreatic cancer, due to stringent statistical threshold or missing interaction analysis, may be unveiled using alternative approaches. To achieve this aim, we explored genetic susceptibility to pancreatic cancer in terms of marginal associations of pathway and genes, as well as their interactions with risk factors. We conducted pathway- and gene-based analysis using GWAS data from 3141 pancreatic cancer patients and 3367 controls with European ancestry. Using the gene set ridge regression in association studies (GRASS) method, we analyzed 197 pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Using the logistic kernel machine (LKM) test, we analyzed 17906 genes defined by University of California Santa Cruz (UCSC) database. Using the likelihood ratio test (LRT) in a logistic regression model, we analyzed 177 pathways and 17906 genes for interactions with risk factors in 2028 pancreatic cancer patients and 2109 controls with European ancestry. After adjusting for multiple comparisons, six pathways were marginally associated with risk of pancreatic cancer ( P < 0.00025): Fc epsilon RI signaling, maturity onset diabetes of the young, neuroactive ligand-receptor interaction, long-term depression (Ps < 0.0002), and the olfactory transduction and vascular smooth muscle contraction pathways (P = 0.0002; Nine genes were marginally associated with pancreatic cancer risk (P < 2.62 × 10−5), including five reported genes (ABO, HNF1A, CLPTM1L, SHH and MYC), as well as four novel genes (OR13C4, OR 13C3, KCNA6 and HNF4 G); three pathways significantly interacted with risk factors on modifying the risk of pancreatic cancer (P < 2.82 × 10−4): chemokine signaling pathway with obesity ( P < 1.43 × 10−4), calcium signaling pathway (P < 2.27 × 10−4) and MAPK signaling pathway with diabetes (P < 2.77 × 10−4). However, none of the 17906 genes tested for interactions survived the multiple comparisons corrections. In summary, our current GWAS study unveiled unidentified genetic susceptibility to pancreatic cancer using alternative methods. These novel findings provide new perspectives on genetic susceptibility to and molecular mechanisms of pancreatic cancer, once confirmed, will shed promising light on the prevention and treatment of this disease. ^

Molecular Mechanisms of Vascular Disease in Patients with Rare Variants in MYH11

Thoracic aortic aneurysms and dissections (TAAD) are the primary disease affecting the thoracic ascending aorta, with an incidence rate of 10.4/100,000. Although about 20% of patients carry a mutation in a single gene that causes their disease, the remaining 80% of patients may also have genetic factors that increase their risk for developing TAAD. Many of the genes that predispose to TAAD encode proteins involved in smooth muscle cell (SMC) contraction and the disease-causing mutations are predicted to disrupt contractile function. SMCs are the predominant cell type in the ascending aortic wall. Mutations in MYH11, encoding the smooth muscle specific myosin heavy chain, are a rare cause of inherited TAAD. However, rare but recurrent non-synonymous variants in MYH11 are present in the general population but do not cause inherited TAAD. The goal of this study was to assess the potential role of these rare variants in vascular diseases. Two distinct variants were selected: the most commonly seen rare variant, MYH11 R247C, and a duplication of the chromosomal region spanning the MYH11 locus at 16p13.1. Genetic analyses indicated that both of these variants were significantly enriched in patients with TAAD compared with controls. A knock-in mouse model of the Myh11 R247C rare variant was generated, and these mice survive and reproduce normally. They have no structural abnormalities of the aorta or signs of aortic disease, but do have decreased aortic contractility. Myh11R247C/R247C mice also have increased proliferative response to vascular injury in vivo and increased proliferation of SMCs in vitro. Myh11R247C/R247C SMCs have decreased contractile gene and protein expression and are dedifferentiated. In fibroblasts, myosin force generation is required for maturation of focal adhesions, and enhancers of RhoA activity replace enhancers of Rac1 activity as maturation occurs. Consistent with these previous findings, focal adhesions are smaller in Myh11R247C/R247C SMCs, and there is decreased RhoA activation. A RhoA activator (CN03) rescues the dedifferentiated phenotype of Myh11R247C/R247C SMCs. Myh11R247C/R247C mice were bred with an existing murine model of aneurysm formation, the Acta2-/- mouse. Over time, mice carrying the R247C allele in conjunction with heterozygous or homozygous loss of Acta2 had significantly increased aortic diameter, and a more rapid accumulation of pathologic markers. These results suggest that the Myh11 R247C rare variant acts as a modifier gene increasing the risk for and severity of TAAD in mice. In patients with 16p13.1 duplications, aortic MYH11 expression is increased, but there is no corresponding increase in smooth muscle myosin heavy chain protein. Using SMCs that overexpress Myh11, we identified alterations in SMC phenotype leading to excessive protein turnover. All contractile proteins, not just myosin, are affected, and the proteins are turned over by autophagic degradation. Surprisingly, these cells are also more contractile compared with wild-type SMCs. The results described in this dissertation firmly establish that rare variants in MYH11 significantly affect the phenotype of SMCs. Further, the data suggests that these rare variants do increase the risk of TAAD via pathways involving altered SMC phenotype and contraction. Therefore, this study validates that these rare genetic variants alter vascular SMCs and provides model systems to explore the contribution of rare variants to disease.