Gene by smoking interaction in blood pressure

Hypertension (HT) is mediated by the interaction of many genetic and environmental factors. Previous genome-wide linkage analysis studies have found many loci that show linkage to HT or blood pressure (BP) regulation, but the results were generally inconsistent. Gene by environment interaction is among the reasons that potentially explain these inconsistencies between studies. Here we investigate influences of gene by smoking (GxS) interaction on HT and BP in European American (EA), African American (AA) and Mexican American (MA) families from the GENOA study. A variance component-based method was utilized to perform genome-wide linkage analysis of systolic blood pressure (SBP), diastolic blood pressure (DBP), and HT status, as well as bivariate analysis for SBP and DBP for smokers, non-smokers, and combined groups. The most significant results were found for SBP in MA. The strongest signal was for chromosome 17q24 (LOD = 4.2), increased to (LOD = 4.7) in bivariate analysis but there was no evidence of GxS interaction at this locus (p = 0.48). Two signals were identified only in one group: on chromosome 15q26.2 (LOD = 3.37) in non-smokers and chromosome 7q21.11 (LOD = 1.4) in smokers, both of which had strong evidence for GxS interaction (p = 0.00039 and 0.009 respectively). There were also two other signals, one on chromosome 20q12 (LOD = 2.45) in smokers, which became much higher in the combined sample (LOD = 3.53), and one on chromosome 6p22.2 (LOD = 2.06) in non-smokers. Neither peak had very strong evidence for GxS interaction (p = 0.08 and 0.06 respectively). A fine mapping association study was performed using 200 SNPs in 30 genes located under the linkage signals on chromosomes 15 and 17. Under the chromosome 15 peak, the association analysis identified 6 SNPs accounting for a 7 mmHg increase in SBP in MA non-smokers. For the chromosome 17 linkage peak, the association analysis identified 3 SNPs accounting for a 6 mmHg increase in SBP in MA. However, none of these SNPs was significant after correcting for multiple testing, and accounting for them in the linkage analysis produced very small reductions in the linkage signal. ^ The linkage analysis of BP traits considering the smoking status produced very interesting signals for SBP in the MA population. The fine mapping association analysis gave some insight into the contribution of some SNPs to two of the identified signals, but since these SNPs did not remain significant after multiple testing correction and did not explain the linkage peaks, more work is needed to confirm these exploratory results and identify the culprit variations under these linkage peaks. ^

The molecular basis for beta-lactamase gene expression in B. anthracis, B. cereus and B. thuringiensis

The susceptibility of most Bacillus anthracis strains to β-lactam antibiotics is intriguing considering that the B. anthracis genome harbors two β-lactamase genes, bla1 and bla2, and closely-related species, Bacillus cereus and Bacillus thuringiensis, typically produce β-lactamases. This work demonstrates that B. anthracis bla expression is affected by two genes, sigP and rsp, predicted to encode an extracytoplasmic function sigma factor and an antisigma factor, respectively. Deletion of the sigP/rsp locus abolished bla expression in a penicillin-resistant clinical isolate and had no effect on bla expression in a prototypical penicillin-susceptible strain. Complementation with sigP/rsp from the penicillin-resistant strain, but not the penicillin-susceptible strain, conferred β-lactamase activity upon both mutants. These results are attributed to a nucleotide deletion near the 5' end of rsp in the penicillin-resistant strain that is predicted to result in a nonfunctional protein. B. cereus and B. thuringiensis sigP and rsp homologues are required for inducible penicillin resistance in those species. Expression of the B. cereus or B. thuringiensis sigP and rsp genes in a B. anthracis sigP/rsp-null mutant confers resistance to β-lactam antibiotics, suggesting that while B. anthracis contains the genes necessary for sensing β-lactam antibiotics, the B. anthracis sigP/rsp gene products are insufficient for bla induction. ^ Because alternative sigma factors recognize unique promoter sequence, direct targets can be elucidated by comparing transcriptional profiling results with an in silico search using the sigma factor binding sequence. Potential σP -10 and -35 promoter elements were identified upstream from bla1 bla2 and sigP. Results obtained from searching the B. anthracis genome with the conserved sequences were evaluated against transcriptional profiling results comparing B. anthracis 32 and an isogenic sigP/rsp -null strain. Results from these analyses indicate that while the absence of the sigP gene significantly affects the transcript levels of 16 genes, only bla1, bla2 and sigP are directly regulated by σP. The genomes of B. cereus and B. thuringiensis strains were also analyzed for the potential σP binding elements. The sequence was located upstream from the sigP and bla genes, and previously unidentified genes predicted to encode a penicillin-binding protein (PBP) and a D-alanyl-D-alanine carboxypeptidase, indicating that the σ P regulon in these species responds to cell-wall stress caused by β-lactam antibiotics. ^ β-lactam antibiotics prevent attachment of new peptidoglycan to the cell wall by blocking the active site of PBPs. A B. cereus and B. thuringiensis pbp-encoding gene located near bla1 contains a potential σP recognition sequence upstream from the annotated translational start. Deletion of this gene abolished β-lactam resistance in both strains. Mutations in the active site of the PBP were detrimental to β-lactam resistance in B. cereus, but not B. thuringiensis, indicating that the transpeptidase activity is only important in B. cereus. I also found that transcript levels of the PBP-encoding gene are not significantly affected by the presence of β-lactam antibiotic. Based on these data I hypothesize that the gene product acts a sensor of β-lactam antibiotic. ^

G870A polymorphism of cyclin D1 gene influences cancer risk - From epidemiological study to mechanism analysis

Alternate splicing of the cyclin D1 gene gives rise to transcript a and b which encode two protein isoforms cyclin D1a and cyclin D1b. Through testing transcript a and transcript b in a series of human samples, we found that cyclin D1 transcript b is ubiquitously expressed as transcript a but in the lower abundance compared to transcript a. Epidemiological studies have reported that the cyclin D1 gene (CCND1) G870A polymorphism influences the risk for a variety of cancer. In this investigation, we examined the cyclin D1b levels in tumor samples with different genotypes and found that higher levels of cyclin D1b are expressed from the A allele than the G allele. Cyclin D1 is known as a cell cycle regulator facilitating the progression of the cell cycle from G1 to S phase in response to the mitogenic signals. It also interacts with several transcription factors and transcriptional coregulators to modulate their activities. It has been reported that cyclin D1a can substitute for estrogen to activate estrogen receptor α (ERα) mediated transcription and can induce the proliferation of estrogen responsive tissues. However the biological role of cyclin D1b in ERα transcriptional regulation has not been previously explored. In this study, we determined that cyclin D1b antagonizes the action of cyclin D1a on ERα mediated transcription. Cell proliferation assays provided the evidence that cyclin D1b negatively regulates estrogen responsive breast cancer cell growth. Taken together, our findings show that the CCND1 G870A polymorphism is correlated with increased levels of cyclin D1b and that cyclin D1b antagonizes the action of cyclin D1a on ERα mediated transcription providing evidence for the mechanism by which the CCND1 G870A polymorphism may be protective in certain types of breast cancer. ^

Genome-wide gene-gene interaction analysis for cardiovascular disease

Numerous studies have been carried out to try to better understand the genetic predisposition for cardiovascular disease. Although it is widely believed that multifactorial diseases such as cardiovascular disease is the result from effects of many genes which working alone or interact with other genes, most genetic studies have been focused on identifying of cardiovascular disease susceptibility genes and usually ignore the effects of gene-gene interactions in the analysis. The current study applies a novel linkage disequilibrium based statistic for testing interactions between two linked loci using data from a genome-wide study of cardiovascular disease. A total of 53,394 single nucleotide polymorphisms (SNPs) are tested for pair-wise interactions, and 8,644 interactions are found to be significant with p-values less than 3.5×10-11. Results indicate that known cardiovascular disease susceptibility genes tend not to have many significantly interactions. One SNP in the CACNG1 (calcium channel, voltage-dependent, gamma subunit 1) gene and one SNP in the IL3RA (interleukin 3 receptor, alpha) gene are found to have the most significant pair-wise interactions. Findings from the current study should be replicated in other independent cohort to eliminate potential false positive results.^

Molecular mechanisms by which p53/LSD1 and ERalpha/Trim24 complexes mediate gene regulation

In this dissertation, I identify two molecular mechanisms by which transcription factors cooperate with their co-regulators to mediate gene regulation. In the first part, I demonstrate that p53 directly recruits LSD1, a histone demethylase, to AFP chromatin to demethylate methylated H3K4 and actively mediate transcription repression. Loss of p53 and LSD1 interaction at chromatin leads to derepression of AFP in hepatic cells. In the second part, I reveal that Trim24 functions as an important co-activator in ERα-mediated gene activation in response to estrogen stimulation. Trim24 is recruited by ligand-bound ERα to chromatin and stabilizes ERα-chromatin interactions by binding to histone H3 via its PHD finger, which preferentially recognizes unmethylated H3K4. ^

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

Transcriptional and translational regulators of gene expression in germ cell development

Germ cell development is a highly coordinated process driven, in part, by regulatory mechanisms that control gene expression. Not only transcription, but also translation, is under regulatory control to direct proper germ cell development. In this dissertation, I have focused on two regulators of germ cell development. One is the homeobox protein RHOX10, which has the potential to be both a transcriptional and translational regulator in mouse male germ cell development. The other is the RNA-binding protein, Hermes, which functions as a translational regulator in Xenopus laevis female germ cell development. ^ Rhox10 is a member of reproductive homeobox gene X-(linked (Rhox) gene cluster, of which expression is developmentally regulated in developing mouse testes. To identify the cell types and developmental stages in which Rhox10 might function, I characterized its temporal and spatial expression pattern in mouse embryonic, neonatal, and adult tissues. Among other things, this analysis revealed that both the level and the subcellular localization of RHOX10 are regulated during germ cell development. To understand the role of Rhox10 in germ cell development, I generated transgenic mice expressing an artificial microRNA (miRNA) targeting Rhox10. While this artificial miRNA robustly downregulated RHOX10 protein expression in vitro, it did not significantly reduce RHOX10 expression in vivo. So I next elected to knockdown RHOX10 levels in spermatogonial stem cells (SSCs), which I found highly express both Rhox10 mRNA and RHOX10 protein. Using a recently developed in vitro culture system for SSCs combined with a short-hairpin RNA (shRNA) approach, I strongly depleted RHOX10 expression in SSCs. These RHOX10-depleted cells exhibited a defect in the ability to form stem cell clusters in vitro. Expression profiling analysis revealed many genes regulated by Rhox10, including many meiotic genes, which could be downstream of Rhox10 in a molecular pathway that controls SSC differentiation. ^ RNA recognition motif (RRM) containing protein, Hermes is localized in germ plasm, where dormant mRNAs are also located, of Xenopus oocytes, which implicates its role in translational regulator. To understand the function of Hermes in oocyte meiosis, I used a morpholino oligonucleotide (MO) based knockdown approach. Microinjection of Hermes MO into fully grown oocytes, which are arrested in meiotic prophase, caused acceleration of oocytes reentry into meiosis (i.e., maturation) upon progesterone induction. Using a candidate approach, I identified at least three targets of Hermes: Ringo/Spy, Xcat2, and Mos. Ringo/Spy and Mos are known to have functions in oocyte maturation, while Ringo/Spy, Xcat2 mRNA are localized in the germ plasm of oocytes, which drives germ cell specification after fertilization. This led me to propose that Hermes functions in both oocyte maturation and germ cell development through its ability to regulate 3 crucial target mRNAs. ^

XIPHOPHORUS ENZYME GENE EXPRESSION: TISSUE SPECIFICITY, GENETIC CONTROL, AND STABILITY IN CULTURED CELLS

Five permanent cell lines were developed from Xiphophorus maculatus, X. helleri, and their hybrids using three tissue sources, including adults and embryos of different stages. To evaluate cell line gene expression for retention of either tissue-of-origin-specific or ontogenetic stage-specific characters, the activity distribution of 44 enzyme loci was determined in 11 X. maculatus tissues, and the developmental genetics of 17 enzyme loci was charted in X. helleri and in helleri x maculatus hybrids using starch gel electrophoresis. In the process, eight new loci were discovered and characterized for Xiphophorus.^ No Xiphophorus cell line showed retention of tissue-of-origin-specific or ontogenetic stage-specific enzyme gene expressional traits. Instead, gene expression was similar among the cell lines. One enzyme, adenosine deaminase (ADA) was an exception. Two adult-origin cell lines expressed ADA, whereas, three cell lines derived independently from embryos did not. ADA('-) expression of Xiphophorus embryo-derived cell lines may represent retention of an embryonic gene expressional trait. In one cell line (T(,3)) derived from 13 pooled interspecific hybrid (F(,2)) embryos, shifts with time were observed at enzyme loci polymorphic between the two species. This suggested shifts in ratios of cells of different genotypes in the population rather than unstable gene expression in one dominant cell type.^ Verification of this hypothesis was attempted by cloning the culture--seeking clones having different genetic signatures. The large number of loci electrophoretically polymorphic between the two species and whose alleles segregated independently into the 13 progeny from which this culture originated almost guaranteed the presence of different genetic signatures (lineages) in T(,3).^ Seven lineages of cells were found within T(,3), each expressing genotypes at some loci not characteristic of the expression of the culture-as-a-whole, supporting the hypothesis tested. Quantitative studies of ADA expression in the whole culture (ADA('-)) and in clones of these seven lineages suggested the predominance in T(,3) of ADA deficient cell lineages, although moderate to high ADA output clones also occurred. Thus, T(,3) has the potential to shift phenotypes from ADA('-) to ADA('+) by simply changing proportions of its constituent cell types, demonstrating that such shifts can occur in any cell culture containing cells of mixed expressional characteristics.^

Characterization of cytochrome P450 4F subfamily: Response in traumatic brain injury and gene regulation

CYP4F enzymes metabolize endogenous molecules including arachidonic acid, leukotrienes and prostaglandins. The involvement of these eisosanoids in inflammation has led to the hypothesis that CYP4Fs may modulate inflammatory conditions after traumatic brain injury (TBI). In rat, TBI elicited changes in mRNA expression of CYP4Fs as a function of time in the cerebrum region. These changes in CYP4F mRNA levels inversely correlated with the cerebral leukotriene B4 (LTB4) level following injury at the same time points. TBI also resulted in changes in CYP4F protein expression and localization around the injury site, where CYP4F1 and CYP4F6 immunoreactivity increased in surrounding astrocytes and CYP4F4 immunoreactivity shifted from endothelia of cerebral vessels to astrocytes. The study with rat primary astrocytes indicated that pro-inflammatory cytokines TNFα and IL-1β could affect the transcription of CYP4Fs to a certain degree, whereas the changing pattern in the primary astrocytes appeared to be different from that in the in vivo TBI model.^ In addition, the regulation of CYP4F genes has been an unsolved issue although factors including cytokines and fatty acids appear to affect CYP4Fs expression in multiple models. In this project, HaCaT cells were used as an in vitro cellular model to define signaling pathways involved in the regulation of human CYP4F genes. Retinoic acids inhibited CYP4F11 expression, whereas cytokines TNFα and IL-1β induced transcription of CYP4F11 in HaCaT cells. The induction of CYP4F11 by both cytokines could be blocked by a JNK specific inhibitor, indicating the involvement of the JNK pathway in the up-regulation of CYP4F11. Retinoic acids are known to function in gene regulation through nuclear receptors RARs and RXRs. The RXR agonist LG268 greatly induced transcription of CYP4F11, whereas RAR agonist TTNPB obviously inhibited CYP4F11 transcription, indicating that the down-regulation of CYP4F11 by retinoic acid was mediated by RARs, and that inhibition of CYP4F11 by retinoic acid may also be related to the competition for RXR receptors. Thus, the CYP4F11 gene is regulated by signaling pathways including the RXR pathway and the JNK pathway. In contrast, the regulation mechanism of other CYP4Fs by retinoic acids appears to be different from that of CYP4F11.^

Comparison of multiple comparison methods for identifying differential gene expression in simulated and real papillary thyroid cancer microarray data

The difficulty of detecting differential gene expression in microarray data has existed for many years. Several correction procedures try to avoid the family-wise error rate in multiple comparison process, including the Bonferroni and Sidak single-step p-value adjustments, Holm's step-down correction method, and Benjamini and Hochberg's false discovery rate (FDR) correction procedure. Each multiple comparison technique has its advantages and weaknesses. We studied each multiple comparison method through numerical studies (simulations) and applied the methods to the real exploratory DNA microarray data, which detect of molecular signatures in papillary thyroid cancer (PTC) patients. According to our results of simulation studies, Benjamini and Hochberg step-up FDR controlling procedure is the best process among these multiple comparison methods and we discovered 1277 potential biomarkers among 54675 probe sets after applying the Benjamini and Hochberg's method to PTC microarray data.^

Estrogen-induced gene expression and patient survival in high-grade serous carcinoma of the ovary

Background. Assessment of estrogen receptor (ER) expression has inconsistent utility as a prognostic marker in epithelial ovarian carcinoma. In breast and endometrial cancers, the use of estrogen-induced gene panels, rather than ER expression alone, has shown improved prognostic capability. Specifically, over-expression of estrogen-induced genes in these tumors is associated with a better prognosis and signifies estrogen sensitivity that can be exploited with hormone antagonizing agents. It was therefore hypothesized that estrogen-induced gene expression in ovarian carcinoma would successfully predict outcomes and differentiate between tumors of varying estrogen sensitivities. Methods. Two hundred nineteen (219) patients with ovarian cancer who underwent surgery at M. D. Anderson between 2004 and 2007 were identified. Of these, eighty-three (83) patients were selected for inclusion because they had advanced stage, high-grade serous carcinoma of the ovary or peritoneum, had not received neoadjuvant chemotherapy, and had readily available frozen tissue for study. All patients had also received adjuvant treatment with platinum and taxane agents. The expression of seven genes known to be induced by estrogen in the female reproductive tract (EIG121, sFRP1, sFRP4, RALDH2, PR, IGF-1, and ER) was measured using qRT-PCR. Unsupervised cluster analyses of multiple gene permutations were used to categorize patients as high or low estrogen-induced gene expressors. QPCR gene expression results were then compared to ER and PR immunohistochemical (IHC) expression. Cox proportional hazards models were used to evaluate the effects of both individual genes and selected gene clusters on patient survival. Results. Median follow-up time was 38.7 months (range 1-68 months). In a multivariate model, overall survival was predicted by sFRP1 expression (HR 1.10 [1.02-1.19], p=0.01) and EIG121 expression (HR 1.28 [1.10-1.49], p<0.01). A cluster defined by EIG121 and ER was further examined because that combination appeared to reasonably segregate tumors into distinct groups of high and low estrogen-induced gene expressors. Shorter overall survival was associated with high estrogen-induced gene expressors (HR 2.84 [1.11-7.30], p=0.03), even after adjustment for race, age, body mass index, and residual disease at debulking. No difference in IHC ER or PR expression was noted between gene clusters. Conclusion. In sharp contrast to breast and endometrial cancers, high estrogen-induced gene expression predicts shorter overall survival in patients with high-grade serous ovarian carcinoma. An estrogen-induced gene biomarker panel may have utility as prognostic indicator and may be useful to guide management with estrogen antagonists in this population.^

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

Genetic epidemiology of gallbladder disease in Mexican Americans and cholesterol 7$\alpha$-hydroxylase gene sequence variation

Among Mexican Americans, the second largest minority group in the United States, the prevalence of gallbladder disease is markedly elevated. Previous data from both genetic admixture and family studies indicate that there is a genetic component to the occurrence of gallbladder disease in Mexican Americans. However, prior to this thesis no formal genetic analysis of gallbladder disease had been carried out nor had any contributing genes been identified.^ The results of complex segregation analysis in a sample of 232 Mexican American pedigrees documented the existence of a major gene having two alleles with age- and gender-specific effects influencing the occurrence of gallbladder disease. The estimated frequency of the allele increasing susceptibility was 0.39. The lifetime probabilities that an individual will be affected by gallbladder disease were 1.0, 0.54, and 0.00 for females of genotypes "AA", "Aa", and "aa", respectively, and 0.68, 0.30, and 0.00 for males, respectively. This analysis provided the first conclusive evidence for the existence of a common single gene having a large effect on the occurrence of gallbladder disease.^ Human cholesterol 7$\alpha$-hydroxylase is the rate-limiting enzyme in bile acid synthesis. The results of an association study in both a random sample and a matched case/control sample showed that there is a significant association between cholesterol 7$\alpha$-hydroxylase gene variation and the occurrence of gallbladder disease in Mexican Americans males but not in females. These data have implicated a specific gene, 7$\alpha$-hydroxylase, in the etiology of gallbladder disease in this population.^ Finally, I asked whether the inferred major gene from complex segregation analysis is genetically linked to the cholesterol 7$\alpha$-hydroxylase gene. Three pedigrees predicted to be informative for linkage analysis by virtue of supporting the major gene hypothesis and having parents with informative genotypes and multiple offspring were selected for this linkage analysis. In each of these pedigrees, the recombination fractions maximized at 0 with a positive, albeit low, LOD score. The results of this linkage analysis provide preliminary and suggestive evidence that the cholesterol 7$\alpha$-hydroxylase gene and the inferred gallbladder disease susceptibility gene are genetically linked. ^

Control of the master virulence regulatory gene atxA in Bacillus anthracis

Transcription of the Bacillus anthracis structural genes for the anthrax toxin proteins and biosynthetic operon for capsule are positively regulated by AtxA, a transcription regulator with unique properties. Consistent with the role of atxA in virulence factor expression, a B. anthracis atxA-null mutant is avirulent in a murine model for anthrax. In batch culture, multiple signals impact atxA transcript levels, and the timing and steady state level of atxA expression is critical for optimal toxin and capsule synthesis. Despite the apparent complex control of atxA transcription, only one trans-acting protein, the transition state regulator AbrB, has been demonstrated to directly interact with the atxA promoter. The AbrB-binding site has been described, but additional cis-acting control sequences have not been defined. Using transcriptional lacZ fusions, electrophoretic mobility shift assays, and Western blot analysis, the cis-acting elements and trans-acting factors involved in regulation of atxA in B. anthracis strains containing either both virulence plasmids, pXO1 and pXO2, or only one plasmid, pXO1, were studied. This work demonstrates that atxA transcription from the major start site P1 is dependent upon a consensus sequence for the housekeeping sigma factor SigA, and an A+T-rich upstream element (UP-element) for RNA polymerase (RNAP). In addition, the data show that a trans-acting protein(s) other than AbrB negatively impacts atxA transcription when it binds specifically to a 9-bp palindrome within atxA promoter sequences located downstream of P1. Mutation of the palindrome prevents binding of the trans-acting protein(s) and results in a corresponding increase in AtxA and anthrax toxin production in a strain- and culture-dependent manner. The identity of the trans-acting repressor protein(s) remains elusive; however, phenotypes associated with mutation of the repressor binding site have revealed that the trans-acting repressor protein(s) indirectly controls B. anthracis development. Mutation of the repressor binding site results in misregulation and overexpression of AtxA in conditions conducive for development, leading to a marked sporulation defect that is both atxA- and pXO2-61-dependent. pXO2-61 is homologous to the sensor domain of sporulation sensor histidine kinases and is proposed to titrate an activating signal away from the sporulation phosphorelay when overexpressed by AtxA. These results indicate that AtxA is not only a master virulence regulator, but also a modulator of proper B. anthracis development. Also demonstrated in this work is the impact of the developmental regulators AbrB, Spo0A, and SigH on atxA expression and anthrax toxin production in a genetically incomplete (pXO1+, pXO2-) and genetically complete (pXO1+, pXO2+) strain background. AtxA and anthrax toxin production resulting from deletion of the developmental regulators are strain-dependent suggesting that factors on pXO2 are involved in control of atxA. The only developmental deletion mutant that resulted in a prominent and consistent strain-independent increase in AtxA protein levels was an abrB-null mutant. As a result of increased AtxA levels, there is early and increased production of anthrax toxins in an abrB-null mutant. In addition, the abrB-null mutant exhibited an increase in virulence in a murine model for anthrax. In contrast, virulence of the atxA promoter mutant was unaffected in a murine model for anthrax despite the production of 5-fold more AtxA than the abrB-null mutant. These results imply that AtxA is not the only factor impacting pathogenesis in an abrB-null mutant. Overall, this work highlights the complex regulatory network that governs expression of atxA and provides an additional role for AtxA in B. anthracis development.

DETERMINING THE GENOTYPE-PHENOTYPE CONNECTION IN SYNTHETIC INDUCIBLE GENE EXPRESSION SYSTEMS

Introduction Gene expression is an important process whereby the genotype controls an individual cell’s phenotype. However, even genetically identical cells display a variety of phenotypes, which may be attributed to differences in their environment. Yet, even after controlling for these two factors, individual phenotypes still diverge due to noisy gene expression. Synthetic gene expression systems allow investigators to isolate, control, and measure the effects of noise on cell phenotypes. I used mathematical and computational methods to design, study, and predict the behavior of synthetic gene expression systems in S. cerevisiae, which were affected by noise. Methods I created probabilistic biochemical reaction models from known behaviors of the tetR and rtTA genes, gene products, and their gene architectures. I then simplified these models to account for essential behaviors of gene expression systems. Finally, I used these models to predict behaviors of modified gene expression systems, which were experimentally verified. Results Cell growth, which is often ignored when formulating chemical kinetics models, was essential for understanding gene expression behavior. Models incorporating growth effects were used to explain unexpected reductions in gene expression noise, design a set of gene expression systems with “linear” dose-responses, and quantify the speed with which cells explored their fitness landscapes due to noisy gene expression. Conclusions Models incorporating noisy gene expression and cell division were necessary to design, understand, and predict the behaviors of synthetic gene expression systems. The methods and models developed here will allow investigators to more efficiently design new gene expression systems, and infer gene expression properties of TetR based systems.