CHARACTERIZATION AND COMPARISON OF SUBHUMAN PRIMATE AND HUMAN TYPE C RETROVIRAL GENE PRODUCTS

The viral specific precursor polyproteins of simian sarcoma/simian associated virus (SiSV/SiAV), baboon endogenous viruses (BaEV), and three human isolate retroviruses, have been analyzed by radioimmunoprecipitation and tryptic peptide mapping. Cells infected with the BaEV isolates are characterized by identical precursor polyproteins: gPr80-85('env), Pr70-71('gag), and Pr65-67('gag). By tryptic digest mapping, m7-BaEV and 455K-BaEV were shown to be highly related. By comparison, mapping studies showed that BILN-BaEV was less highly related to m7-BaEV than was 455K-BaEV. Chase-incubated cells infected with BaEV also contained a stable, p28-related polyprotein termed P72('gag). This polyprotein appeared to arise by posttranslational modification of Pr70-71('gag). Tryptic digest mapping of BaEV and HL23V precursor polyproteins suggested that the BaEV-like component of HL23V was more closely related to m7-BaEV than to 455K-BaEV or BILN-BaEV.^ The intracellular precursor polyproteins of SiSV(SiAV) and gibbon ape leukemia virus (GaLV) were compared to the intracellular proteins of the human retrovirus isolates, HL23V, HEL12V, and A1476V. Cells infected with SiSV(SiAV) were characterized by polyproteins Pr200('gag-pol), gPr80('env), Pr80('gag), pr60('gag), and Pr40('gag). We have found that the human isolates are identical to true SiAV with regard to the size and structure of their precursor polyproteins. Both gPr80('env) and Pr60('gag) of SiAV were identical by tryptic peptide mapping to the respective proteins from the three human retroviral isolates examined. We have also shown that these viruses differ significantly from each of the GaLV isolates studied. Since SiAV differs substantially from any known GaLV isolate, we feel that it is unlikely that SiAV is a subtype of GaLV which exists today in the gibbon gene pool. The experimental evidence suggests that SiAV may be an exogenous human retrovirus that was transmitted originally into the human gene pool in the distant past by cross-species infection with GaLV(,SF) or with the GaLV(,SF) progenitor virus. It is, therefore, quite possible that SiAV expression in the pet woolly monkey arose from a recent infection of that monkey with SiAV from humans.^

Cytochrome P450 4F isoforms in different species: Identification, gene regulation and putative roles in inflammation

The cytochrome P450 4F subfamily comprises a group of enzymes that metabolize derivatives of arachidonic acid such as prostaglandins, lipoxins leukotrienes and hydroxyeicosatetraenoic acids, which are important mediators involved in the inflammatory response. Therefore, we speculate that CYP4Fs might be able to modulate the extent of the inflammation by controlling of the tissue levels of these inflammatory mediators, especially, leukotriene B4. One way to provide support for this hypothesis is to test whether the expression of CYP4Fs changes under inflammatory conditions, since these changes are required to adjust the levels of inflammatory mediators. ^ A lipopolysacchride (LPS) induced rat inflammation model was used to analyze the expressions of rat CYP4F4 and CYP4F5 in liver and kidney. LPS administration did not change the constitutive expression level of CYP4F4 and CYP4F5. In liver, the expressions of CYP4F4 and CYP4F5 decreased to 50–60% of the untreated level. The same effect of LPS on CYP4F4 and CYP4F5 expression can be mimicked in hepatocyte primary cultures treated with LPS, indicating a direct of effect of LPS on hepatocytes. LPS treatment also decreased the activity of liver microsomes towards chlorpromazine, however, antibody inhibition study revealed that liver CYP4Fs are not the only players in metabolizing chlorpromazine. To study further the underlying mechanism, CYP4F5 gene was isolated, characterized, and the promoter region was defined. ^ Accumulating evidence showed that peroxisome proliferator-activated receptors (PPARs) play an active role in inflammation. To investigate the possible role of PPARα in regulating CYP4F expression by inflammation or by clofibrate treatment, the expressions of two new mouse 4F isoforms were analyzed in PPARα knockout mice upon LPS or clofibrate challenge. A novel induction of CYP4F15 by LPS and clofibrate was observed in kidney, and this effect is totally dependent on the presence of PPARα. Renal CYP4F16 expression was not affected by LPS or clofibrate in both (+/+) and (−/−) mice. In contrast, hepatic expressions of CYP4F15 and CYP4F16 were reduced significantly in (+/+) mice, but much less in (−/−) mice, suggesting that PPARα is partially responsible for this down-regulation. Clofibrate treatment reduced the expression of CYP4F16 in liver, but has no effect on CYP4F15 and PPARα does not have a role in hepatic CYP4F expression regulated by clofibrate. In general, CYP4Fs are regulated in an isoform-, tissue- and species-specific manner. ^ A human CYP4F isoform, CYP4F11, was isolated. The genomic structure was also solved by using database mining and bioinformatics tools. Localization of CYP4F11 to chromosome 19, 16 kb upstream of CYP4F2, suggests that human CYP4F genes may form a cluster on chromosome 19. This novel human 4F is highly expressed in liver, as well as in kidney, heart and skeletal muscle. Further study of the activity and gene regulation on CYP4F11 will provide us more insights into the physiological functions of CYP4F subfamily. ^

Pharmacogenetics of the human glutathione S-transferase P1 gene in the metabolism and therapeutic efficacy of cis-diamminedichloroplatinum

The human GSTP1 gene has been shown, conclusively, to be polymorphic. The three main GSTP1 alleles, GSTP1*A, GSTP1*B, and GSTP1*C, encode proteins which differ in the 3-dimensional structure of their active sites and in their function in phase II metabolism of carcinogens, mutagens, and anticancer agents. Although, it is well established that GSTP1 is over expressed in many human tumors and that the levels of GSTP1 expression correlate directly with tumor resistance to chemotherapy and inversely with patient survival, the significance of the polymorphic GSTP1 gene locus on tumor response to chemotherapy remains unclear. The goal of this project was to define the role and significance of the polymorphic GSTP1 gene locus in GSTP1-based tumor drug resistance and as a determinant of patient response to chemotherapy. The hypothesis to be tested was that the polymorphic GSTP1 gene locus will confer to tumors a differential ability to metabolize cisplatin resulting in a GSTP1 genotype-based sensitivity to cisplatin. The study examined: (a) whether the different GSTP 1 alleles confer different levels of cellular protection against cisplatin-induced cytotoxicity, (b) whether the allelic GSTP1 proteins metabolize cisplatin with different efficiencies, and (c) whether the GSTP1 genotype is a determinant of tumor response to cisplatin therapy. The results demonstrate that the GSTP1 alleles differentially protect tumors against cisplatin-induced apoptosis and clonogenic cell kill in the rank order: GSTP1*C > GSTP1*B > GSTP1*A. The same rank order was observed for the kinetics of GSTP1-catalyzed cisplatin metabolism, both in cell-free and cellular systems, to the rate-limiting monoglutathionyl-platinum metabolite, which was characterized, for the first time, by mass spectral analysis. Finally, this study demonstrates that both GSTP1 genotype and the level of GSTP1 expression significantly contribute to tumor sensitivity to cisplatin treatment. Overall, the results of this project show that the polymorphic GSTP1 gene locus plays a significant role in tumor sensitivity to cisplatin treatment. Furthermore, these studies have contributed to the overall understanding of the significance of the polymorphic GSTP1 gene locus in tumor resistance to cancer chemotherapy and have provided the basis for further investigations into how this can be utilized to optimize and individualize cancer chemotherapy for cancer patients. ^

Studies on gene expression within germinal centers

During T cell dependent immune responses, the acquisition of B cell memory from naïve cells takes place within a highly specialized microenvironment: The germinal centers (GC) of the secondary lymphoid organs. The GC reaction is a tightly regulated process in which the balance between survival and death is mediated by signals transduced through ligation of critical costimulatory molecules such as CD40 and CD154. While most cognate receptor-ligand interactions occur between T-cells and antigen (Ag)-presenting cells (APC) such as B-cells, evidence of homotypic B cell interactions has emerged. Despite the progress in our understanding of the reaction, several questions remain: (1) What determines the concomitant expression of CD40 and its ligand CD154 by GC B-cells? (2) Which molecules are responsible for inducing GC-B cell survival? and (3) how can cognate T-cell help be recruited into the organized structure of GCs? ^ Because the expression of costimulatory and survival molecules is predominant at distinct Ag-dependent maturation stages, we hypothesized the existence of stage specific gene expression responsible for the regulation of the GC reaction. Our studies reveal several novel genes whose expression may be critical for the GC reaction. The discovery of AKNA reveals the mechanism behind homotypic B cell CD40 and CD40 ligand interactions, which can explain the costimulatory signaling to GC B cells in the absence of T cells. Additionally, the identification of the pro-survival molecule PRELI may provide a novel mechanism for the survival of Ag-selected B cells. We propose that PRELI and its phylogenic homologues represent a novel family of proteins responsible for the protection of cells against caspase-independent apoptosis. Furthermore, we show that GC B cells actively participate in the recruitment of T cells through the secretion of CC and CxC chemokines, thus supporting their mutual involvement in cognate interactions. ^

Downstream targets of the Rhox5 homeobox gene

The X-linked mouse Rhox gene cluster contains over 30 homeobox genes that are candidates to regulate multiple steps in male and female gametogenesis. The founding member of the Rhox gene cluster, Rhox5, is an androgen-dependent gene expressed in Sertoli cells that promotes the survival and differentiation of the adjacent male germ cells. To decipher downstream signaling pathways of Rhox5, I used in vivo and in vitro microarray profiling to identify and characterize downstream targets of Rhox5 in the testis. This led to the identification of many Rhox5 -regulated genes, two of which I focused on in more detail. One of them, Unc5c, encodes a pro-apoptotic receptor with tumor suppressor activity that I found is negatively regulated by Rhox5 through a Rhox5-response element in the Unc5c 5' untranslated region (5' UTR). Examination of other mouse Rhox family members revealed that Rhox2 and Rhox3 also have the ability to downregulate Unc5c expression. The human RHOX protein RHOXF2 also had this ability, indicating that Unc5c repression is a conserved Rhox-dependent response. The repression of Unc5c expression by Rhox5 may, in part, mediate Rhox5's pro-survival function in the testis, as I found that Unc5c mutant mice have decreased germ cell apoptosis in the testis. This along with my other data leads me to propose a model in which Rhox5 is a negative regulator upstream of Unc5c in a Sertoli-cell pathway that promotes germ-cell survival. The other Rhox5-regulated gene that I studied in detail is insulin II (Ins2). Several lines of evidence, including electrophoretic mobility shift anaylsis, promoter mutagenesis, and chromatin immuoprecipitation analysis indicated that Ins2 is a direct target of Rhox5. Structure-function analysis identified homeodomain residues and the RHOX5 amino-terminal domain crucial for conferring Ins2 inducibility. Rhox5 regulates not only the Ins2 gene but also genes encoding other secreted proteins regulating metabolism (adiponectin and resistin), the rate-liming enzyme for monosaturated fatty acid biosynthesis (SCD-1), and transcription factors crucial for regulating metabolism (the nuclear hormone receptor PPARγ). I propose that the regulation of some or all of these molecules in Sertoli cells is responsible for the Rhox5-dependent survival of the adjacent germ cells. ^

Regulation of chromatin structure, Smad binding and AFP expression by the Forkhead box transcription factor: Foxa1

Transcription factors must be able to access their DNA binding sites to either activate or repress transcription. However, DNA wrapping and compaction into chromatin occludes most binding sites from ready access by proteins. Pioneer transcription factors are capable of binding their DNA elements within a condensed chromatin context and then reducing the level of nucleosome occupancy so that the chromatin structure is more accessible. This altered accessibility increases the probability of other transcription factors binding to their own DNA binding elements. My hypothesis is that Foxa1, a ‘pioneer’ transcription factor, activates alpha-fetoprotein (AFP) expression by binding DNA in a chromatinized environment, reducing the nucleosome occupancy and facilitating binding of additional transcription factors.^ Using retinoic-acid differentiated mouse embryonic stem cells, we illustrate a mechanism for activation of the tumor marker AFP by the pioneer transcription factor Foxa1 and TGF-β downstream effector transcription factors Smad2 and Smad4. In differentiating embryonic stem cells, binding of the Foxa1 forkhead box transcription factor to chromatin reduces nucleosome occupancy and levels of linker histone H1 at the AFP distal promoter. The more accessible DNA is subsequently bound by the Smad2 and Smad4 transcription factors, concurrent with activation of transcription. Chromatin immunoprecipitation analyses combined with siRNA-mediated knockdown indicate that Smad protein binding and the reduction of nucleosome occupancy at the AFP distal promoter is dependent on Foxa1. In addition to facilitating transcription factor binding, Foxa1 is also associated with histone modifications related to active gene expression. Acetylation of lysine 9 on histone H3, a mark that is associated active transcription, is dependent on Foxa1, while methylation of H3K4, also associated with active transcription, is independent of Foxa1. I propose that Foxa1 potentiates a region of chromatin to respond to Smad proteins, leading to active expression of AFP.^ These studies demonstrate one mechanism whereby a transcription factor can alter the accessibility of additional transcription factors to chromatin, by altering nucleosome positions. Specifically, Foxa1 exposes DNA so that Smad4 can bind to its regulatory element and activate transcription of the tumor-marker gene AFP.^

Mechanism of methylation gene reactivation in human cancer cells

Epigenetic silencing of tumor suppressor genes by DNA hypermethylation at promoter regions is a common event in carcinogenesis and tumor progression. Abrogation of methylation and reversal of epigenetic silencing is a very potent way in cancer treatment. However, the reactivation mechanisms are poorly understood. In this study, we first developed a cell line model system named YB5, derived from SW48 cancer cell line, which bears one copy of stably integrated EGFP gene on Chromosome 1p31.1 region. The GFP gene expression is transcriptionally silenced due to the hypermethylated promoter CMV. However, the GFP expression can be restored using demethylating agent 5-aza-2' deoxycytidine (DAC), and detected by FACS and fluorescent microscopy. Using this system, we observed the heterogeneous reactivation induced by DAC treatment. After flow sorting, GFP negative cells exhibited similar level of incomplete demethylation compared to GFP positive cells on repetitive LINE1 element, tumor suppressor genes such as P16, CDH13, and RASSF1a, and CMV promoter as well. However, the local chromatin of CMV-GFP locus altered to an open structure marked by high H3 lysine 9 acetylation and low H3 lysine 27 tri-methylation in GFP positive cells, while the GFP negative cells retained mostly the original repressive marks. Thus, we concluded that DAC induced DNA hypomethylation alone does not directly determine the level of re-expression, and the resetting of the local chromatin structure under hypomethylation environment is required for gene reactivation. Besides, a lentivirus vector-based shRNA screening was performed using the YB5 system. Although it is the rare chance that vector lands in the neighboring region of GFP, we found that the exogenous vector DNA inserted into the upstream region of GFP gene locus led to the promoter demethylation and reactivated the silenced GFP gene. Thus, epigenetic state can be affected by changing of the adjacent nucleic acid sequences. Further, this hypermethylation silenced system was utilized for epigenetic drug screening. We have found that DAC combined with carboplatin would enhance the GFP% yield and increase expression of other tumor suppressor genes than DAC alone, and this synergistic effect may be related to DNA repair process. In summary, these studies reveal that reversing of methylation silencing requires coordinated alterations of DNA methylation, chromatin structure, and local microenvironment. ^

Mouse retinal development: Role of cell adhesion and mechanism of gene regulation

Cell differentiation and pattern formation are fundamental processes in animal development that are under intense investigation. The mouse retina is a good model to study these processes because it has seven distinct cell types, and three well-laminated nuclear layers that form during embryonic and postnatal life. β-catenin functions as both the nuclear effector for the canonical Wnt pathway and a cell adhesion molecule, and is required for the development of various organs. To study the function of β-catenin in retinal development, I used a Cre-loxP system to conditionally ablate β-catenin in the developing retina. Deletion of β-catenin led to disrupted laminar structure but did not affect the differentiation of any of the seven cell types. Eliminating β-catenin did not reduce progenitor cell proliferation, although enhanced apoptosis was observed. Further analysis showed that disruption of cell adhesion was the major cause of the observed patterning defects. Overexpression of β-catenin during retinal development also disrupted the normal retinal lamination and caused a transdifferentiation of neurons into pigmented cells. The results indicate that β-catenin functions as a cell adhesion molecule but not as a Wnt pathway component during retinal neurogenesis, and is essential for lamination but not cell differentiation. The results further imply that retinal lamination and cell differentiation are genetically separable processes. ^ Sonic hedgehog (shh) is expressed in retinal ganglion cells under the control of transcription factor Pou4f2 during retinal development. Previous studies identified a phylogenetically conserved region in the first intron of shh containing a Pou4f2 binding site. Transgenic reporter mice in which reporter gene expression was driven by this region showed that this element can direct gene expression specifically in the retina, but expression was not limited to the ganglion cells. From these data I hypothesized that this element is required for shh expression in the retina but is not sufficient for specific ganglion cell expression. To further test this hypothesis, I created a conditional allele by flanking this region with two loxP sites. Lines carrying this allele will be crossed with retinal-specific Cre lines to remove this element in the retina. My hypothesis predicts that alteration in shh expression and subsequent retinal defects will occur in the retinas of these mice. ^

Systems Biology Approaches to Probe Gene Regulation in Bacteria

Mechanisms that allow pathogens to colonize the host are not the product of isolated genes, but instead emerge from the concerted operation of regulatory networks. Therefore, identifying components and the systemic behavior of networks is necessary to a better understanding of gene regulation and pathogenesis. To this end, I have developed systems biology approaches to study transcriptional and post-transcriptional gene regulation in bacteria, with an emphasis in the human pathogen Mycobacterium tuberculosis (Mtb). First, I developed a network response method to identify parts of the Mtb global transcriptional regulatory network utilized by the pathogen to counteract phagosomal stresses and survive within resting macrophages. As a result, the method unveiled transcriptional regulators and associated regulons utilized by Mtb to establish a successful infection of macrophages throughout the first 14 days of infection. Additionally, this network-based analysis identified the production of Fe-S proteins coupled to lipid metabolism through the alkane hydroxylase complex as a possible strategy employed by Mtb to survive in the host. Second, I developed a network inference method to infer the small non-coding RNA (sRNA) regulatory network in Mtb. The method identifies sRNA-mRNA interactions by integrating a priori knowledge of possible binding sites with structure-driven identification of binding sites. The reconstructed network was useful to predict functional roles for the multitude of sRNAs recently discovered in the pathogen, being that several sRNAs were postulated to be involved in virulence-related processes. Finally, I applied a combined experimental and computational approach to study post-transcriptional repression mediated by small non-coding RNAs in bacteria. Specifically, a probabilistic ranking methodology termed rank-conciliation was developed to infer sRNA-mRNA interactions based on multiple types of data. The method was shown to improve target prediction in Escherichia coli, and therefore is useful to prioritize candidate targets for experimental validation.

Structure-function analysis of human Integrator subunit-4

Structure-function analysis of human Integrator subunit 4 Anupama Sataluri Advisor: Eric. J. Wagner, Ph.D. Uridine-rich small nuclear RNAs (U snRNA) are RNA Polymerase-II (RNAPII) transcripts that are ubiquitously expressed and are known to be essential for gene expression. snRNAs play a key role in mRNA splicing and in histone mRNA expression. Inaccurate snRNA biosynthesis can lead to diseases related to defective splicing and histone mRNA expression. Although the 3′ end formation mechanism and processing machinery of other RNAPII transcripts such as mRNA has been well studied, the mechanism of snRNA 3′ end processing has remained a mystery until the recent discovery of the machinery that mediates this process. In 2005, a complex of 14 subunits (the Integrator complex) associated with RNA Polymerase-II was discovered. The 14subunits were annotated Integrator 1-14 based on their size. The subunits of this complex together were found to facilitate 3′ end processing of snRNA. Identification of the Integrator complex propelled research in the direction of understanding the events of snRNA 3’end processing. Recent studies from our lab confirmed that Integrator subunit (IntS) 9 and 11 together perform the endonucleolytic cleavage of the nascent snRNA 3′ end to generate mature snRNA. However, the role of other members of the Integrator complex remains elusive. Current research in our lab is focused on deciphering the role of each subunit within the Integrator complex This work specifically focuses on elucidating the role of human Integrator subunit 4 (IntS4) and understanding how it facilitates the overall function of the complex. IntS4 has structural similarity with a protein called “Symplekin”, which is part of the mRNA 3’end processing machinery. Symplekin has been thoroughly researched in recent years and structure-function correlation studies in the context of mRNA 3’end processing have reported a scaffold function for Symplekin due to the presence of HEAT repeat motifs in its N-terminus. Based upon the structural similarity between IntS4 and Symplekin, we hypothesized that Integrator subunit 4 may be behaving as a Symplekin-like scaffold molecule that facilitates the interaction between other members of the Integrator Complex. To answer this question, the two important goals of this study were to: 1) identify the region of IntS4, which is important for snRNA 3′ end processing and 2) determine binding partners of IntS4 which promote its function as a scaffold. IntS4 structurally consists of a highly conserved N-terminus with 8 HEAT repeats, followed by a nonconserved C- terminus. A series of siRNA resistant N and C-terminus deletion constructs as well as specific point mutants within its N-terminal HEAT repeats were generated for human IntS4 and, utilizing a snRNA transcriptional readthrough GFP-reporter assay, we tested their ability to rescue misprocessing. This assay revealed a possible scaffold like property of IntS4. To probe IntS4 for interaction partners, we performed co-immunoprecipitation on nuclear extracts of IntS4 expressing stable cell lines and identified IntS3 and IntS5 among other Integrator subunits to be binding partners which facilitate the scaffold like function of hIntS4. These findings have established a critical role for IntS4 in snRNA 3′ end processing, identified that both its N and C termini are essential for its function, and mapped putative interaction domains with other Integrator subunits.

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.

Structure and biosynthesis of a pyruvate -dependent enzyme---phosphatidylserine decarboxylase from {\it Escherichia coli\/}

Phosphatidylserine decarboxylase of E. coli, a cytoplasmic membrane protein, catalyzes the formation of phosphatidylethanolamine, the principal phospholipid of the organism. The activity of the enzyme is dependent on a covalently bound pyruvate (Satre and Kennedy (1978) J. Biol. Chem. 253, 479-483). This study shows that the enzyme consists of two nonidentical subunits, $\alpha$ (Mr = 7,332) and $\beta$ (Mr = 28,579), with the pyruvate prosthetic group in amide linkage to the amino-terminus of the $\alpha$ subunit. Partial protein sequence and DNA sequence analysis reveal that the two subunits are derived from a proenzyme ($\pi$ subunit, Mr = 35,893) through a post-translational event. During the conversion of the proenzyme to the $\alpha$ and $\beta$ subunits, the peptide bond between Gly253-Ser254 is cleaved, and Ser254 is converted to the pyruvate prosthetic group at the amino-terminus of the $\alpha$ subunit (Li and Dowhan (1988) J. Biol. Chem. 263, 11516-11522).^ The proenzyme cannot be detected in cells carrying either single or multiple copies of the gene (psd), but can be observed in a T7 RNA polymerase/promoter and transcription-translation system. The cleavage of the wild-type proenzyme occurs rapidly with a half-time on the order of 2 min. Changing of the Ser254 to cysteine (S254C) or threonine (S254T) slows the cleavage rate dramatically and results in mutants with a half-time for processing of around 2-4 h. Change of the Ser254 to alanine (S254A) blocks the cleavage of the proenzyme. The reduced processing rate with the mutations of the proenzyme is consistent with less of the functional enzyme being made. Mutants S254C and S254T produce $\sim$15% and $\sim$1%, respectively, of the activity of the wild-type allele, but can still complement a temperature-sensitive mutant of the psd locus. Neither detectable activity nor complementation is observed by mutant S254A. These results are consistent with the hydroxyl-group of the Ser254 playing a critical role in the cleavage of the peptide bond Gly253-Ser254 of the pro-phosphatidylserine decarboxylase, and support the mechanism proposed by Snell and co-workers (Recsei and Snell (1984) Annu. Rev. Biochem. 53, 357-387) for the formation of the prosthetic group of pyruvate-dependent decarboxylases. ^

The effects of DNA cytosine methylation on H -ras promoter structure and function

The effect of DNA cytosine methylation on H-ras promoter activity was assessed using a transient expression system employing the plasmid H-rasCAT (NaeI H-ras promoter linked to the chloramphenicol acetyltransferase (CAT) gene). This 551 bp promoter is 80% GC rich, enriched with 168 CpG dinucleotides, and contains six functional GC box elements which represent major DNA methylation target sites. Prokaryotic methyltransferases HhaI (CGm$\sp5$CG) and HpaII (Cm$\sp5$CGG) alone or in combination with a human placental methyltransferase (HP MTase) were used to introduce methyl groups at different CpG sites within the promoter. To test for functional promoter activity, the methylated plasmids were introduced into CV-1 cells and CAT activity assessed 48 h post-transfection. Methylation at specific HhaI and HpaII sites reduced CAT expression by 70%, whereas more extensive methylation at generalized CpG sites with HP MTase inactivated the promoter $>$95%. The inhibition of H-ras promoter activity was not attributable to methylation-induced differences in DNA uptake or stability in the cell, topological form of the plasmid, or methylation effects in nonpromoter regions. We also observed that DNA cytosine methylation of a 360 bp promoter fragment by HP MTase induced a local change in DNA conformation. Using three independent methodologies (nitrocellulose filter binding assays, gel mobility shifts, and Southwestern blots), we determined that this change in promoter conformation affected the interaction of nuclear proteins with cis-regulatory sequences residing in the promoter region. The results provide evidence to suggest that DNA methylation may regulate gene expression by inducing changes in local promoter conformation which in turn alters the interactions between DNA and protein factors required for transcription. The results provide supportive evidence for the hypothesis of Cedar and Riggs, who postulated that DNA methylation may regulate gene expression by altering the binding affinities of proteins for DNA. ^

Genetic integration of a nuclear -encoded mitochondrial gene with cardiac function

Cardiovascular disease (CVD) is the leading cause of death in the United States. One manifestation of CVD known to increase mortality is an enlarged, or hypertrophic heart. Hypertrophic cardiomyocytes adapt to increased contractile demand at the genetic level with a re-emergence of the fetal gene program and a downregulation of fatty acid oxidation genes with concomitant increased reliance on glucose-based metabolism. To understand the transcriptional regulatory pathways that implement hypertrophic directives we analyzed the upstream promoter region of the muscle specific isoform of the nuclear-encoded mitochondrial gene, carnitine palmitoyltransferase-1β (CPT-1β) in cultured rat neonatal cardiac myocytes. This enzyme catalyzes the rate-limiting step of fatty acid entry into β-oxidation and is downregulated in cardiac hypertrophy and failure, making it an attractive model for the study of hypertrophic gene regulation and metabolic adaptations. We demonstrate that the muscle-enriched transcription factors GATA-4 and SRF synergistically activate CPT-1β; moreover, DNA binding to cognate sites and intact protein structure are required. This mechanism coordinates upregulation of energy generating processes with activation of the energy consuming contractile promoter for cardiac α-actin. We hypothesized that fatty acid or glucose responsive transcription factors may also regulate CPT-1β. Oleate weakly stimulates CPT-1β activity; in contrast, the glucose responsive Upstream Stimulatory Factors (USF) dramatically depresses the CPT-1β reporter. USF regulates CPT-1β through a novel physical interaction with the cofactor PGC-1 and abrogation of MEF2A/PGC-1 synergistic stimulation. In this way, USF can inversely regulate metabolic gene programs and may play a role in the shift of metabolic substrate preference seen in hypertrophy. Failing hearts have elevated expression of the nuclear hormone receptor COUP-TF. We report that COUP-TF significantly suppresses reporter transcription independent of DNA binding and specific interactions with GATA-4, Nkx2.5 or USF. In summary, CPT-1β transcriptional regulation integrates mitochondrial gene expression with two essential cardiac functions: contraction and metabolic substrate oxidation. ^